Rational design of a calcium-binding protein

Effects of selection stringency on the outcomes of directed evolution

Directed evolution makes mutant lineages compete in climbing complicated sequence-function landscapes. Given this underlying complexity it is unclear how selection stringency, a ubiquitous parameter of directed evolution, impacts the outcome. Here we approach this question in terms of the fitnesses of the candidate variants at each round and the heterogeneity of their distributions of fitness effects. We show that even if the fittest mutant is most likely to yield the fittest mutants in the next round of selection, diversification can improve outcomes by sampling a larger variety of fitness effects. We find that heterogeneity in fitness effects between variants, larger population sizes, and evolution over a greater number of rounds all encourage diversification.

Protein MRI Contrast Agents as an Effective Approach for Precision Molecular Imaging

ABSTRACT

Cancer and other acute and chronic diseases are results of perturbations of common molecular determinants in key biological and signaling processes. Imaging is critical for characterizing dynamic changes in tumors and metastases, the tumor microenvironment, tumor-stroma interactions, and drug targets, at multiscale levels. Magnetic resonance imaging (MRI) has emerged to be a primary imaging modality for both clinical and preclinical applications due to its advantages over other modalities, including sensitivity to soft tissues, nondepth limitations, and the use of nonionizing radiation. However, extending the application of MRI to achieve both qualitative and quantitative precise molecular imaging with the capability to quantify molecular biomarkers for early detection, staging, and monitoring therapeutic treatment requires the capacity to overcome several major challenges including the trade-off between metal-binding affinity and relaxivity, which is an issue frequently associated with small chelator contrast agents. In this review, we will introduce the criteria of ideal contrast agents for precision molecular imaging and discuss the relaxivity of current contrast agents with defined first shell coordination water molecules. We will then report our advances in creating a new class of protein-targeted MRI contrast agents (ProCAs) with contributions to relaxivity largely derived from the secondary sphere and correlation time. We will summarize our rationale, design strategy, and approaches to the development and optimization of our pioneering ProCAs with desired high relaxivity, metal stability, and molecular biomarker-targeting capability, for precision MRI. From first generation (ProCA1) to third generation (ProCA32), we have achieved dual high r1 and r2 values that are 6- to 10-fold higher than clinically approved contrast agents at magnetic fields of 1.5 T, and their relaxivity values at high field are also significantly higher, which enables high resolution during small animal imaging. Further engineering of multiple targeting moieties enables ProCA32 agents that have strong biomarker-binding affinity and specificity for an array of key molecular biomarkers associated with various chronic diseases, while maintaining relaxation and exceptional metal-binding and selectivity, serum stability, and resistance to transmetallation, which are critical in mitigating risks associated with metal toxicity. Our leading product ProCA32.collagen has enabled the first early detection of liver metastasis from multiple cancers at early stages by mapping the tumor environment and early stage of fibrosis from liver and lung in vivo, with strong translational potential to extend to precision MRI for preclinical and clinical applications for precision diagnosis and treatment.

Design of Diverse Asymmetric Pockets in De Novo Homo-oligomeric Proteins

A challenge for design of protein-small-molecule recognition is that incorporation of cavities with size, shape, and composition suitable for specific recognition can considerably destabilize protein monomers. This challenge can be overcome through binding pockets formed at homo-oligomeric interfaces between folded monomers. Interfaces surrounding the central homo-oligomer symmetry axes necessarily have the same symmetry and so may not be well suited to binding asymmetric molecules. To enable general recognition of arbitrary asymmetric substrates and small molecules, we developed an approach to designing asymmetric interfaces at off-axis sites on homo-oligomers, analogous to those found in native homo-oligomeric proteins such as glutamine synthetase. We symmetrically dock curved helical repeat proteins such that they form pockets at the asymmetric interface of the oligomer with sizes ranging from several angstroms, appropriate for binding a single ion, to up to more than 20 Å across. Of the 133 proteins tested, 84 had soluble expression in E. coli, 47 had correct oligomeric states in solution, 35 had small-angle X-ray scattering (SAXS) data largely consistent with design models, and 8 had negative-stain electron microscopy (nsEM) 2D class averages showing the structures coming together as designed. Both an X-ray crystal structure and a cryogenic electron microscopy (cryoEM) structure are close to the computational design models. The nature of these proteins as homo-oligomers allows them to be readily built into higher-order structures such as nanocages, and the asymmetric pockets of these structures open rich possibilities for small-molecule binder design free from the constraints associated with monomer destabilization.

Design of Calcium-Binding Proteins to Sense Calcium

Calcium controls numerous biological processes by interacting with different classes of calcium binding proteins (CaBP’s), with different affinities, metal selectivities, kinetics, and calcium dependent conformational changes. Due to the diverse coordination chemistry of calcium, and complexity associated with protein folding and binding cooperativity, the rational design of CaBP’s was anticipated to present multiple challenges. In this paper we will first discuss applications of statistical analysis of calcium binding sites in proteins and subsequent development of algorithms to predict and identify calcium binding proteins. Next, we report efforts to identify key determinants for calcium binding affinity, cooperativity and calcium dependent conformational changes using grafting and protein design. Finally, we report recent advances in designing protein calcium sensors to capture calcium dynamics in various cellular environments.

A new approach to biomining: Bioengineering surfaces for metal recovery from aqueous solutions

Electronics waste production has been fueled by economic growth and the demand for faster, more efficient consumer electronics. The glass and metals in end-of-life electronics components can be reused or recycled; however, conventional extraction methods rely on energy-intensive processes that are inefficient when applied to recycling e-waste that contains mixed materials and small amounts of metals. To make e-waste recycling economically viable and competitive with obtaining raw materials, recovery methods that lower the cost of metal reclamation and minimize environmental impact need to be developed. Microbial surface adsorption can aid in metal recovery with lower costs and energy requirements than traditional metal-extraction approaches. We introduce a novel method for metal recovery by utilizing metal-binding peptides to functionalize fungal mycelia and enhance metal recovery from aqueous solutions such as those found in bioremediation or biomining processes. Using copper-binding as a proof-of-concept, we compared binding parameters between natural motifs and those derived in silico, and found comparable binding affinity and specificity for Cu. We then combined metal-binding peptides with chitin-binding domains to functionalize a mycelium-based filter to enhance metal recovery from a Cu-rich solution. This finding suggests that engineered peptides could be used to functionalize biological surfaces to recover metals of economic interest and allow for metal recovery from metal-rich effluent with a low environmental footprint, at ambient temperatures, and under circumneutral pH.

Redox-Responsive Protein Design: Design of a Small Protein Motif Dependent on Glutathionylation

Cysteine S-glutathionylation is a protein post-translational modification that promotes cellular responses to changes in oxidative conditions. The design of protein motifs that directly depend on defined changes to protein side chains provides new methods for probing diverse protein post-translational modifications. A canonical, 12-residue EF-hand motif was redesigned to be responsive to cysteine glutathionylation. The key design principle was the replacement of the metal-binding Glu12 carboxylate of an EF-hand with a motif capable of metal binding via a free carboxylate in the glutathione-conjugated peptide. In the optimized peptide (DKDADGWCG), metal binding and terbium luminescence were dependent on glutathionylation, with weaker metal binding in the presence of reduced cysteine but increased metal affinity and a 3.5-fold increase in terbium luminescence at 544 nm when cysteine was glutathionylated. Nuclear magnetic resonance spectroscopy indicated that the structure at all residues of the glutathionylated peptide changed in the presence of metal, with chemical shift changes consistent with the adoption of an EF-hand-like structure in the metal-bound glutathionylated peptide. This small protein motif consists of canonical amino acids and is thus genetically encodable, for its potential use as a localized tag to probe protein glutathionylation.

Many-body effect determines the selectivity for Ca2+ and Mg2+ in proteins

Calcium ion is a versatile messenger in many cell-signaling processes. To achieve their functions, calcium-binding proteins selectively bind Ca2+ against a background of competing ions such as Mg2+ The high specificity of calcium-binding proteins has been intriguing since Mg2+ has a higher charge density than Ca2+ and is expected to bind more tightly to the carboxylate groups in calcium-binding pockets. Here, we showed that the specificity for Ca2+ is dictated by the many-body polarization effect, which is an energetic cost arising from the dense packing of multiple residues around the metal ion. Since polarization has stronger distance dependence compared with permanent electrostatics, the cost associated with the smaller Mg2+ is much higher than that with Ca2+ and outweighs the electrostatic attraction favorable for Mg2+ With the AMOEBA (atomic multipole optimized energetics for biomolecular simulation) polarizable force field, our simulations captured the relative binding free energy between Ca2+ and Mg2+ for proteins with various types of binding pockets and explained the nonmonotonic size dependence of the binding free energy in EF-hand proteins. Without electronic polarization, the smaller ions are always favored over larger ions and the relative binding free energy is roughly proportional to the net charge of the pocket. The many-body effect depends on both the number and the arrangement of charged residues. Fine-tuning of the ion selectivity could be achieved by combining the many-body effect and geometric constraint.

Fluorescent analysis of bioactive molecules in single cells based on microfluidic chips

Single-cell analysis of bioactive molecules is an essential strategy for a better understanding of cell biology, exploring cell heterogeneity, and improvement of the ability to detect early diseases. In single-cell analysis, highly efficient single-cell manipulation techniques and high-sensitive detection schemes are in urgent need. The rapid development of fluorescent analysis techniques combined with microfluidic chips have offered a widely applicable solution. Thus, in this review, we mainly focus on the application of fluorescence methods in components analysis on microchips at a single-cell level. By targeting different types of biological molecules in cells such as nucleic acids, proteins, and active small molecules, we specially introduce and comment on their corresponding fluorescent probes, fluorescence labelling and sensing strategies, and different fluorescence detection instruments used in single-cell analysis on a microfluidic chip. We hope that through this review, readers will have a better understanding of single-cell fluorescence analysis, especially for single-cell component fluorescence analysis based on microfluidic chips.

Benchmarking a computational design method for the incorporation of metal ion-binding sites at symmetric protein interfaces

The design of novel metal-ion binding sites along symmetric axes in protein oligomers could provide new avenues for metalloenzyme design, construction of protein-based nanomaterials and novel ion transport systems. Here, we describe a computational design method, symmetric protein recursive ion-cofactor sampling (SyPRIS), for locating constellations of backbone positions within oligomeric protein structures that are capable of supporting desired symmetrically coordinated metal ion(s) chelated by sidechains (chelant model). Using SyPRIS on a curated benchmark set of protein structures with symmetric metal binding sites, we found high recovery of native metal coordinating rotamers: in 65 of the 67 (97.0%) cases, native rotamers featured in the best scoring model while in the remaining cases native rotamers were found within the top three scoring models. In a second test, chelant models were crossmatched against protein structures with identical cyclic symmetry. In addition to recovering all native placements, 10.4% (8939/86013) of the non-native placements, had acceptable geometric compatibility scores. Discrimination between native and non-native metal site placements was further enhanced upon constrained energy minimization using the Rosetta energy function. Upon sequence design of the surrounding first-shell residues, we found further stabilization of native placements and a small but significant (1.7%) number of non-native placement-based sites with favorable Rosetta energies, indicating their designability in existing protein interfaces. The generality of the SyPRIS approach allows design of novel symmetric metal sites including with non-natural amino acid sidechains, and should enable the predictive incorporation of a variety of metal-containing cofactors at symmetric protein interfaces.

Calcium Dynamics Mediated by the Endoplasmic/Sarcoplasmic Reticulum and Related Diseases

The flow of intracellular calcium (Ca2+) is critical for the activation and regulation of important biological events that are required in living organisms. As the major Ca2+ repositories inside the cell, the endoplasmic reticulum (ER) and the sarcoplasmic reticulum (SR) of muscle cells are central in maintaining and amplifying the intracellular Ca2+ signal. The morphology of these organelles, along with the distribution of key calcium-binding proteins (CaBPs), regulatory proteins, pumps, and receptors fundamentally impact the local and global differences in Ca2+ release kinetics. In this review, we will discuss the structural and morphological differences between the ER and SR and how they influence localized Ca2+ release, related diseases, and the need for targeted genetically encoded calcium indicators (GECIs) to study these events.

High Resolution Prediction of Calcium-Binding Sites in 3D Protein Structures Using FEATURE

Metal-binding proteins are ubiquitous in biological systems ranging from enzymes to cell surface receptors. Among the various biologically active metal ions, calcium plays a large role in regulating cellular and physiological changes. With the increasing number of high-quality crystal structures of proteins associated with their metal ion ligands, many groups have built models to identify Ca(2+) sites in proteins, utilizing information such as structure, geometry, or homology to do the inference. We present a FEATURE-based approach in building such a model and show that our model is able to discriminate between nonsites and calcium-binding sites with a very high precision of more than 98%. We demonstrate the high specificity of our model by applying it to test sets constructed from other ions. We also introduce an algorithm to convert high scoring regions into specific site predictions and demonstrate the usage by scanning a test set of 91 calcium-binding protein structures (190 calcium sites). The algorithm has a recall of more than 93% on the test set with predictions found within 3 Å of the actual sites.

Fast kinetics of calcium signaling and sensor design

Fast calcium signaling is regulated by numerous calcium channels exhibiting high spatiotemporal profiles which are currently measured by fluorescent calcium sensors. There is still a strong need to improve the kinetics of genetically encoded calcium indicators (sensors) to capture calcium dynamics in the millisecond time frame. In this review, we summarize several major fast calcium signaling pathways and discuss the recent developments and application of genetically encoded calcium indicators to detect these pathways. A new class of genetically encoded calcium indicators designed with site-directed mutagenesis on the surface of beta-barrel fluorescent proteins to form a pentagonal bipyramidal-like calcium binding domain dramatically accelerates calcium binding kinetics. Furthermore, novel genetically encoded calcium indicators with significantly increased fluorescent lifetime change are advantageous in deep-field imaging with high light-scattering and notable morphology change.

Design of ProCAs (protein-based Gd(3+) MRI contrast agents) with high dose efficiency and capability for molecular imaging of cancer biomarkers

Magnetic resonance imaging (MRI) is the leading imaging technique for disease diagnostics, providing high resolution, three-dimensional images noninvasively. MRI contrast agents are designed to improve the contrast and sensitivity of MRI. However, current clinically used MRI contrast agents have relaxivities far below the theoretical upper limit, which largely prevent advancing molecular imaging of biomarkers with desired sensitivity and specificity. This review describes current progress in the development of a new class of protein-based MRI contrast agents (ProCAs) with high relaxivity using protein design to optimize the parameters that govern relaxivity. Further, engineering with targeting moiety allows these contrast agents to be applicable for molecular imaging of prostate cancer biomarkers by MRI. The developed protein-based contrast agents also exhibit additional in vitro and in vivo advantages for molecular imaging of disease biomarkers, such as high metal-binding stability and selectivity, reduced toxicity, proper blood circulation time, and higher permeability in tumor tissue in addition to improved relaxivities.

De novo design of peptide scaffolds as novel preorganized ligands for metal-ion coordination

This chapter describes how de novo designed peptides can be used as novel preorganized ligands for metal ion coordination. The focus is on the design of peptides which are programmed to spontaneously self-assemble into α-helical coiled coils in aqueous solution, and how metal ion binding sites can be engineered onto and into these structures. In addition to describing the various design principles, some key examples are covered illustrating the success of this approach, including a more detailed example in the case study.

Molecular imaging of EGFR/HER2 cancer biomarkers by protein MRI contrast agents

Epidermal growth factor receptor (EGFR) and HER2 are major prognosis biomarkers and drug targets overexpressed in various types of cancer cells. There is a pressing need to develop MRI contrast agents capable of enhancing the contrast between normal tissues and tumors with high relaxivity, capable of targeting tumors, and with high intratumoral distribution and minimal toxicity. In this review, we first discuss EGFR signaling and its role in tumor progression as a major drug target. We then report our progress in the development of protein contrast agents with significant improvement of both r1 and r2 relaxivities, pharmacokinetics, in vivo retention time, and in vivo dose efficiency. Finally, we report our effort in the development of EGFR-targeted protein contrast agents with the capability to cross the endothelial boundary and with good tissue distribution across the entire tumor mass. The noninvasive capability of MRI to visualize spatially and temporally the intratumoral distribution as well as quantify the levels of EGFR and HER2 would greatly improve our ability to track changes of the biomarkers during tumor progression, monitor treatment efficacy, aid in patient selection, and further develop novel targeted therapies for clinical application.

Structural analysis of autoinhibition in the Ras-specific exchange factor RasGRP1

RasGRP1 and SOS are Ras-specific nucleotide exchange factors that have distinct roles in lymphocyte development. RasGRP1 is important in some cancers and autoimmune diseases but, in contrast to SOS, its regulatory mechanisms are poorly understood. Activating signals lead to the membrane recruitment of RasGRP1 and Ras engagement, but it is unclear how interactions between RasGRP1 and Ras are suppressed in the absence of such signals. We present a crystal structure of a fragment of RasGRP1 in which the Ras-binding site is blocked by an interdomain linker and the membrane-interaction surface of RasGRP1 is hidden within a dimerization interface that may be stabilized by the C-terminal oligomerization domain. NMR data demonstrate that calcium binding to the regulatory module generates substantial conformational changes that are incompatible with the inactive assembly. These features allow RasGRP1 to be maintained in an inactive state that is poised for activation by calcium and membrane-localization signals.

DOI:http://dx.doi.org/10.7554/eLife.00813.001

Individual cells within the human body must grow, divide or specialize to perform the tasks required of them. The fates of these cells are often directed by proteins in the Ras family, which detect signals from elsewhere in the body and orchestrate responses within each cell. The activities of these proteins must be tightly controlled, because cancers and developmental diseases can result if Ras proteins are not properly regulated.

Binding to the small molecule GTP activates Ras and causes conformational changes that allow it to interact with other proteins in various signaling pathways in the cell. GTP is loaded into Ras by proteins called nucleotide exchange factors, which can replace ‘used’ nucleotides with ‘fresh’ ones to activate Ras.

These nucleotide exchange factors are also tightly regulated. For example, the genes for many exchange factors are only switched on after particular signals are received, which can restrict their presence to defined times and locations (e.g., cells or tissues). Also, when activating signals are absent, nucleotide exchange factors commonly reside in the cytoplasm, whereas the Ras proteins remain bound to lipid membranes inside the cell.

RasGRP1 is a nucleotide exchange factor that controls the development of immune cells, and leukemia and lupus can result if it is not regulated correctly. However, many questions about RasGRP1 remain unanswered, including how it is able to remain inactive, and how it is activated by various different signals.

Iwig et al. have now revealed the mechanisms through which RasGRP1 suppresses Ras signaling in immune cells by solving the structures of two fragments of RasGRP1 and then using a combination of structural, biochemical and cell-based methods to explore how it is activated. These analyses revealed that inactive RasGRP1 adopts a conformation in which one of its regulatory elements blocks access to the Ras binding site. Surprisingly, RasGRP1 can form dimers; this hides the portions of the protein that associate with the membrane and thereby keeps RasGRP1 away from Ras. Iwig et al. also found that two signals, calcium ions and a lipid called diacylglycerol, overcome these inhibitory mechanisms by changing the conformation of RasGRP1 and recruiting it to the membrane.

These studies provide a framework for understanding how disease-associated mutations in RasGRP1 bypass the regulatory mechanisms that insure proper immune cell development, and will be critical for developing therapeutic agents that inhibit RasGRP1 activity.

DOI:http://dx.doi.org/10.7554/eLife.00813.002

Computational enzyme design

Recent developments in computational chemistry and biology have come together in the "inside-out" approach to enzyme engineering. Proteins have been designed to catalyze reactions not previously accelerated in nature. Some of these proteins fold and act as catalysts, but the success rate is still low. The achievements and limitations of the current technology are highlighted and contrasted to other protein engineering techniques. On its own, computational "inside-out" design can lead to the production of catalytically active and selective proteins, but their kinetic performances fall short of natural enzymes. When combined with directed evolution, molecular dynamics simulations, and crowd-sourced structure-prediction approaches, however, computational designs can be significantly improved in terms of binding, turnover, and thermal stability.

Design of a novel class of protein-based magnetic resonance imaging contrast agents for the molecular imaging of cancer biomarkers

Magnetic resonance imaging (MRI) of disease biomarkers, especially cancer biomarkers, could potentially improve our understanding of the disease and drug activity during preclinical and clinical drug treatment and patient stratification. MRI contrast agents with high relaxivity and targeting capability to tumor biomarkers are highly required. Extensive work has been done to develop MRI contrast agents. However, only a few limited literatures report that protein residues can function as ligands to bind Gd(3+) with high binding affinity, selectivity, and relaxivity. In this paper, we focus on reporting our current progress on designing a novel class of protein-based Gd(3+) MRI contrast agents (ProCAs) equipped with several desirable capabilities for in vivo application of MRI of tumor biomarkers. We will first discuss our strategy for improving the relaxivity by a novel protein-based design. We then discuss the effect of increased relaxivity of ProCAs on improving the detection limits for MRI contrast agent, especially for in vivo application. We will further report our efforts to improve in vivo imaging capability and our achievement in molecular imaging of cancer biomarkers with potential preclinical and clinical applications.

The calmodulin regulator protein, PEP-19, sensitizes ATP-induced Ca2+ release

PEP-19 is a small, intrinsically disordered protein that binds to the C-domain of calmodulin (CaM) via an IQ motif and tunes its Ca(2+) binding properties via an acidic sequence. We show here that the acidic sequence of PEP-19 has intrinsic Ca(2+) binding activity, which may modulate Ca(2+) binding to CaM by stabilizing an initial Ca(2+)-CaM complex or by electrostatically steering Ca(2+) to and from CaM. Because PEP-19 is expressed in cells that exhibit highly active Ca(2+) dynamics, we tested the hypothesis that it influences ligand-dependent Ca(2+) release. We show that PEP-19 increases the sensitivity of HeLa cells to ATP-induced Ca(2+) release to greatly increase the percentage of cells responding to sub-saturating doses of ATP and increases the frequency of Ca(2+) oscillations. Mutations in the acidic sequence of PEP-19 that inhibit or prevent it from modulating Ca(2+) binding to CaM greatly inhibit its effect on ATP-induced Ca(2+) release. Thus, this cellular effect of PEP-19 does not depend simply on binding to CaM via the IQ motif but requires its acidic metal binding domain. Tuning the activities of Ca(2+) mobilization pathways places PEP-19 at the top of CaM signaling cascades, with great potential to exert broad effects on downstream CaM targets, thus expanding the biological significance of this small regulator of CaM signaling.

Engineering of an environmentally responsive beta roll peptide for use as a calcium-dependent cross-linking domain for peptide hydrogel formation

We have created a set of rationally designed peptides that form calcium-dependent hydrogels based on the beta roll peptide domain. In the absence of calcium, the beta roll domain is intrinsically disordered. Upon the addition of calcium, the peptide forms a beta helix secondary structure. We have designed two variations of our beta roll domain. First, we have mutated one face of the beta roll domain to contain leucine residues so that the calcium-dependent structural formation leads to dimerization through hydrophobic interactions. Second, an α-helical leucine zipper domain is appended to the engineered beta roll domain as an additional means of forming intermolecular cross-links. This full peptide construct forms a hydrogel only in calcium-rich environments. The resulting structural and mechanical properties of the supramolecular assemblies are compared with the wild-type domain using several biophysical techniques including circular dichroism, FRET, bis-ANS binding and microrheology. The calcium responsiveness and rheological properties of the leucine beta roll containing construct confirm the potential of this allosterically regulated scaffold to serve as a cross-linking domain for stimulus-responsive biomaterials development.

Identification of regions responsible for the open conformation of S100A10 using chimaeric S100A11-S100A10 proteins

S100A11 is a dimeric EF-hand calcium-binding protein. Calcium binding to S100A11 results in a large conformational change that uncovers a broad hydrophobic surface used to interact with phospholipid-binding proteins (annexins A1 and A2) and facilitate membrane vesiculation events. In contrast with other S100 proteins, S100A10 is unable to bind calcium due to deletion and substitution of calcium-ligating residues. Despite this, calcium-free S100A10 assumes an 'open' conformation that is very similar to S100A11 in its calcium-bound state. To understand how S100A10 is able to adopt an open conformation in the absence of calcium, seven chimaeric proteins were constructed where regions from calcium-binding sites I and II, and helices II-IV in S100A11 were replaced with the corresponding regions of S100A10. The chimaeric proteins having substitutions in calcium-binding site II displayed increased hydrophobic surface exposure as assessed by bis-ANS (4,4'-dianilino-1,1'-binaphthyl-5,5'disulfonic acid, dipotassium salt) fluorescence and phenyl-Sepharose binding in the absence of calcium. This response is similar to that observed for Ca2+-S100A11 and calcium-free S100A10. Further, this substitution resulted in calcium-insensitive binding to annexin A2 for one chimaeric protein. The results indicate that residues within site II are important in stabilizing the open conformation of S100A10 and presentation of its target binding site. In contrast, S100A11 chimaeric proteins with helical substitutions displayed poorer hydrophobic surface exposure and, consequently, unobservable annexin A2 binding. The present study represents a first attempt to systematically understand the molecular basis for the calcium-insensitive open conformation of S100A10.

Design of lanthanide fingers: compact lanthanide-binding metalloproteins

Lanthanides have interesting chemical properties; these include luminescent, magnetic, and catalytic functions. Toward the development of proteins incorporating novel functions, we have designed a new lanthanide-binding motif, lanthanide fingers. These were designed based on the Zif268 zinc finger, which exhibits a beta beta alpha structural motif. Lanthanide fingers utilize an Asp(2)Glu(2) metal-coordination environment to bind lanthanides through a tetracarboxylate peptide ligand. The iterative design of a general lanthanide-binding peptide incorporated the following key elements: 1) residues with high alpha-helix and beta-sheet propensities in the respective secondary structures; 2) an optimized big box alpha-helix N-cap; 3) a Schellman alpha-helix C-cap motif; and 4) an optional D-Pro-Ser type II' beta-turn in the beta-hairpin. The peptides were characterized for lanthanide binding by circular dichroism (CD), NMR, and fluorescence spectroscopy. In all instances, stabilization of the peptide secondary structures resulted in an increase in metal affinity. The optimized protein design was a 25-residue peptide that was a general lanthanide-binding motif; this binds all lanthanides examined in a competitive aqueous environment, with a dissociation constant of 9.3 microM for binding Er(3+). CD spectra of the peptide-lanthanide complexes are similar to those of zinc fingers and other beta beta alpha proteins. Metal binding involves residues from the N-terminal beta-hairpin and the C terminal alpha-helical segments of the peptide. NMR data indicated that metal binding induced a global change in the peptide structure. The D-Pro-Ser type II' beta-turn motif could be replaced by Thr-Ile to generate genetically encodable lanthanide fingers. Replacement of the central Phe with Trp generated genetically encodable lanthanide fingers that exhibited terbium luminescence greater than that of an EF-hand peptide.

Remeasuring HEWL pK(a) values by NMR spectroscopy: methods, analysis, accuracy, and implications for theoretical pK(a) calculations

Site-specific pK(a) values measured by NMR spectroscopy provide essential information on protein electrostatics, the pH-dependence of protein structure, dynamics and function, and constitute an important benchmark for protein pK(a) calculation algorithms. Titration curves can be measured by tracking the NMR chemical shifts of several reporter nuclei versus sample pH. However, careful analysis of these curves is needed to extract residue-specific pK(a) values since pH-dependent chemical shift changes can arise from many sources, including through-bond inductive effects, through-space electric field effects, and conformational changes. We have re-measured titration curves for all carboxylates and His 15 in Hen Egg White Lysozyme (HEWL) by recording the pH-dependent chemical shifts of all backbone amide nitrogens and protons, Asp/Glu side chain protons and carboxyl carbons, and imidazole protonated carbons and protons in this protein. We extracted pK(a) values from the resulting titration curves using standard fitting methods, and compared these values to each other, and with those measured previously by ¹H NMR (Bartik et al., Biophys J 1994;66:1180–1184). This analysis gives insights into the true accuracy associated with experimentally measured pK(a) values. We find that apparent pK(a) values frequently differ by 0.5–1.0 units depending upon the nuclei monitored, and that larger differences occasionally can be observed. The variation in measured pK(a) values, which reflects the difficulty in fitting and assigning pH-dependent chemical shifts to specific ionization equilibria, has significant implications for the experimental procedures used for measuring protein pK(a) values, for the benchmarking of protein pK(a) calculation algorithms, and for the understanding of protein electrostatics in general.

Capturing, sharing and analysing biophysical data from protein engineering and protein characterization studies

Large amounts of data are being generated annually on the connection between the sequence, structure and function of proteins using site-directed mutagenesis, protein design and directed evolution techniques. These data provide the fundamental building blocks for our understanding of protein function, molecular biology and living organisms in general. However, much experimental data are never deposited in databases and is thus 'lost' in journal publications or in PhD theses. At the same time theoretical scientists are in need of large amounts of experimental data for benchmarking and calibrating novel predictive algorithms, and theoretical progress is therefore often hampered by the lack of suitable data to validate or disprove a theoretical assumption. We present PEAT (Protein Engineering Analysis Tool), an application that integrates data deposition, storage and analysis for researchers carrying out protein engineering projects or biophysical characterization of proteins. PEAT contains modules for DNA sequence manipulation, primer design, fitting of biophysical characterization data (enzyme kinetics, circular dichroism spectroscopy, NMR titration data, etc.), and facilitates sharing of experimental data and analyses for a typical university-based research group. PEAT is freely available to academic researchers at http://enzyme.ucd.ie/PEAT.

Calcium-induced folding of a beta roll motif requires C-terminal entropic stabilization

Beta roll motifs are associated with several proteins secreted by the type 1 secretion system (T1SS). Located just upstream of the C-terminal T1SS secretion signal, they are believed to act as calcium-induced switches that prevent folding before secretion. Bordetella pertussis adenylate cyclase (CyaA) toxin has five blocks of beta roll motifs (or repeats-in-toxin motifs) separated by linkers. The block V motif on its own has been reported to be non-responsive to calcium. Only when the N- and C-terminal linkers, or flanking groups, were fused did the motif bind calcium and fold. In an effort to understand the requirements for beta roll folding, we have truncated the N- and C-terminal flanks at several locations to determine the minimal essential sequences. Calcium-responsive beta roll folding occurred even in the absence of the natural N-terminal flank. The natural C-terminal flank could not be truncated without decreased calcium affinity and only partially truncated before losing calcium-responsiveness. Globular protein fusion at the C-terminus likewise enabled calcium-induced folding but fusions solely at the N-terminus failed. This demonstrates that calcium-induced folding is an inherent property of the beta roll motif rather than the flanking groups. Given the disparate nature of the observed functional flanking groups, C-terminal fusions appear to confer calcium-responsiveness to the beta roll motif via a non-specific mechanism, suggesting that entropic stabilization of the unstructured C-terminus can enable beta roll folding. Increased calcium affinity was observed when the natural C-terminal flank was used to enable calcium-induced folding, pointing to its cooperative participation in beta roll formation. This work indicates that a general principle of C-terminal entropic stabilization can enable stimulus-responsive repeat protein folding, while the C-terminal flank has a specific role in tuning calcium-responsive beta roll formation. These observations are in stark contrast to what has been reported for other repeat proteins.

Integration of Diverse Research Methods to Analyze and Engineer Ca-Binding Proteins: From Prediction to Production

In recent years, increasingly sophisticated computational and bioinformatics tools have evolved for the analyses of protein structure, function, ligand interactions, modeling and energetics. This includes the development of algorithms to recursively evaluate side-chain rotamer permutations, identify regions in a 3D structure that meet some set of search parameters, calculate and minimize energy values, and provide high-resolution visual tools for theoretical modeling. Here we discuss the interdependency between different areas of bioinformatics, the evolution of different algorithm design approaches, and finally the transition from theoretical models to real-world design and application as they relate to Ca(2+)-binding proteins. Within this context, it has become evident that significant pre-experimental design and calculations can be modeled through computational methods, thus eliminating potentially unproductive research and increasing our confidence in the correlation between real and theoretical models. Moving from prediction to production, it is anticipated that bioinformatics tools will play an increasingly significant role in research and development, improving our ability to both understand the physiological roles of Ca(2+) and other metals and to extend that knowledge to the design of function-specific synthetic proteins capable of fulfilling different roles in medical diagnostics and therapeutics.

Calciomics: prediction and analysis of EF-hand calcium binding proteins by protein engineering

Ca²⁺ plays a pivotal role in the physiology and biochemistry of prokaryotic and mammalian organisms. Viruses also utilize the universal Ca²⁺ signal to create a specific cellular environment to achieve coexistence with the host, and to propagate. In this paper we first describe our development of a grafting approach to understand site-specific Ca²⁺ binding properties of EF-hand proteins with a helix-loop-helix Ca²⁺ binding motif, then summarize our prediction and identification of EF-hand Ca²⁺ binding sites on a genome-wide scale in bacteria and virus, and next report the application of the grafting approach to probe the metal binding capability of predicted EF-hand motifs within the streptococcal hemoprotein receptor (Shr) of Streptococcus pyrogenes and the nonstructural protein 1 (nsP1) of Sindbis virus. When methods such as the grafting approach are developed in conjunction with prediction algorithms we are better able to probe continuous Ca²⁺-binding sites that have been previously underrepresented due to the limitation of conventional methodology.

Docking of calcium ions in proteins with flexible side chains and deformable backbones

A method of docking Ca(2+) ions in proteins with flexible side chains and deformable backbones is proposed. The energy was calculated with the AMBER force field, implicit solvent, and solvent exposure-dependent and distance-dependent dielectric function. Starting structures were generated with Ca(2+) coordinates and side-chain torsions sampled in 1000 A(3) cubes centered at the experimental Ca(2+) positions. The energy was Monte Carlo-minimized. The method was tested on fourteen Ca(2+)-binding sites. For twelve Ca(2+)-binding sites the root mean square (RMS) deviation of the apparent global minimum from the experimental structure was below 1.3 and 1.7 A for Ca(2+) ions and side-chain heavy atoms, respectively. Energies of multiple local minima correlate with the RMS deviations from the X-ray structures. Two Ca(2+)-binding sites at the surface of proteinase K were not predicted, because of underestimation of Ca(2+) hydration energy by the implicit-solvent method.

Design of functional metalloproteins

Metalloproteins catalyse some of the most complex and important processes in nature, such as photosynthesis and water oxidation. An ultimate test of our knowledge of how metalloproteins work is to design new metalloproteins. Doing so not only can reveal hidden structural features that may be missing from studies of native metalloproteins and their variants, but also can result in new metalloenzymes for biotechnological and pharmaceutical applications. Although it is much more challenging to design metalloproteins than non-metalloproteins, much progress has been made in this area, particularly in functional design, owing to recent advances in areas such as computational and structural biology.

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.

Rational design of a conformation-switchable Ca2+- and Tb3+-binding protein without the use of multiple coupled metal-binding sites

Ca2+, as a messenger of signal transduction, regulates numerous target molecules via Ca2+-induced conformational changes. Investigation into the determinants for Ca2+-induced conformational change is often impeded by cooperativity between multiple metal-binding sites or protein oligomerization in naturally occurring proteins. To dissect the relative contributions of key determinants for Ca2+-dependent conformational changes, we report the design of a single-site Ca2+-binding protein (CD2.trigger) created by altering charged residues at an electrostatically sensitive location on the surface of the host protein rat Cluster of Differentiation 2 (CD2).CD2.trigger binds to Tb3+ and Ca2+ with dissociation constants of 0.3 +/- 0.1 and 90 +/- 25 microM, respectively. This protein is largely unfolded in the absence of metal ions at physiological pH, but Tb3+ or Ca2+ binding results in folding of the native-like conformation. Neutralization of the charged coordination residues, either by mutation or protonation, similarly induces folding of the protein. The control of a major conformational change by a single Ca2+ ion, achieved on a protein designed without reliance on sequence similarity to known Ca2+-dependent proteins and coupled metal-binding sites, represents an important step in the design of trigger proteins.

Structural differences between Pb2+- and Ca2+-binding sites in proteins: implications with respect to toxicity

Pb(2+) is known to displace physiologically-relevant metal ions in proteins. To investigate potential relationships between Pb(2+)/protein complexes and toxicity, data from the protein data bank were analyzed to compare structural properties of Pb(2+)- and Ca(2+)-binding sites. Results of this analysis reveal that the majority of Pb(2+) sites (77.1%) involve 2-5 binding ligands, compared with 6+/-2 for non-EF-Hand and 7+/-1 for EF-Hand Ca(2+)-binding sites. The mean net negative charge by site (1.7) fell between values noted for non-EF-Hand (1+/-1) and EF-Hand (3+/-1). Oxygen is the dominant ligand for both Pb(2+) and Ca(2+), but Pb(2+) binds predominantly with sidechain Glu (38.4%), which is less prevalent in both non-EF-Hand (10.4%) and EF-Hand (26.6%) Ca(2+)-binding sites. A comparison of binding geometries where Pb(2+) has replaced Ca(2+) in calmodulin (CaM) and Zn(2+) in 5-aminolaevulinic acid dehydratase (ALAD) revealed protein structural changes that appear to be unrelated to ionic displacement. Structural changes observed with CaM may be related to opportunistic binding of Pb(2+) in regions of high electrostatic charge, whereas ALAD may bind multiple Pb(2+) ions in the active site. These results suggest that Pb(2+) adapts to structurally-diverse binding geometries and that opportunistic binding may play an active role in molecular metal toxicity.

Statistical analysis of structural characteristics of protein Ca2+-binding sites

To better understand the biological significance of Ca(2+), we report a comprehensive statistical analysis of calcium-binding proteins from the Protein Data Bank to identify structural parameters associated with EF-hand and non-EF-hand Ca(2+)-binding sites. Comparatively, non-EF-hand sites utilize lower coordination numbers (6 +/- 2 vs. 7 +/- 1), fewer protein ligands (4 +/- 2 vs. 6 +/- 1), and more water ligands (2 +/- 2 vs. 1 +/- 0) than EF-hand sites. The orders of ligand preference for non-EF-hand and EF-hand sites, respectively, were H(2)O (33.1%) > side-chain Asp (24.5%) > main-chain carbonyl (23.9%) > side-chain Glu (10.4%), and side-chain Asp (29.7%) > side-chain Glu (26.6%) > main-chain carbonyl (21.4%) > H(2)O (13.3%). Less formal negative charge was observed in the non-EF-hand than in the EF-hand binding sites (1 +/- 1 vs. 3 +/- 1). Additionally, over 20% of non-EF-hand sites had formal charge values of zero due to increased utilization of water and carbonyl oxygen ligands. Moreover, the EF-hand sites presented a narrower range of ligand distances and bond angles than non-EF-hand sites, possibly owing to the highly conserved helix-loop-helix motif. Significant differences between ligand types (carbonyl, side chain, bidentate) demonstrated that angles associated with each type must be classified separately, and the EF-hand side-chain Ca-O-C angles exhibited an unusual bimodal quality consistent with an Asp distribution that differed from the Gaussian model observed for non-EF-hand proteins. The results of this survey more accurately describe differences between EF-hand and non-EF-hand proteins and provide new parameters for the prediction and design of different classes of Ca(2+)-binding proteins.

Rational design of protein-based MRI contrast agents

We describe the rational design of a novel class of magnetic resonance imaging (MRI) contrast agents with engineered proteins (CAi.CD2, i = 1, 2,..., 9) chelated with gadolinium. The design of protein-based contrast agents involves creating high-coordination Gd(3+) binding sites in a stable host protein using amino acid residues and water molecules as metal coordinating ligands. Designed proteins show strong selectivity for Gd(3+) over physiological metal ions such as Ca(2+), Zn(2+), and Mg(2+). These agents exhibit a 20-fold increase in longitudinal and transverse relaxation rate values over the conventional small-molecule contrast agents, e.g., Gd-DTPA (diethylene triamine pentaacetic acid), used clinically. Furthermore, they exhibit much stronger contrast enhancement and much longer blood retention time than Gd-DTPA in mice. With good biocompatibility and potential functionalities, these protein contrast agents may be used as molecular imaging probes to target disease markers, thereby extending applications of MRI.

Rational design of a novel calcium-binding site adjacent to the ligand-binding site on CD2 increases its CD48 affinity

Electrostatic interactions are important for molecular recognition processes including Ca2+-binding and cell adhesion. To understand these processes, we have successfully introduced a novel Ca2+-binding site into the non-Ca2+-dependent cell adhesion protein CD2 using our criteria that are specifically tailored to the structural and functional properties of the protein environment and charged adhesion surface. This designed site with ligand residues exclusively from the beta-sheets selectively binds to Ca2+ and Ln3+ over other mono- and divalent cations. While Ca2+ and Ln3+ binding specifically alters the local environment of the designed Ca2+-binding site, the designed protein undergoes a significantly smaller conformation change compared with those observed in naturally occurring Ca2+-binding sites that are composed of at least part of the flexible loop and helical regions. In addition, the CD2-CD48-binding affinity increased approximately threefold after protein engineering, suggesting that the cell adhesion of CD2 can be modulated by altering the local electrostatic environment. The study provides site-specific information for regulating cell adhesion within CD2 and gives insight into the structural factors required for Ca2+-modulated biological processes.

Identification and dissection of Ca(2+)-binding sites in the extracellular domain of Ca(2+)-sensing receptor

Ca(2+)-sensing receptors (CaSRs) represent a class of receptors that respond to changes in the extracellular Ca(2+) concentration ([Ca(2+)](o)) and activate multiple signaling pathways. A major barrier to advancing our understanding of the role of Ca(2+) in regulating CaSRs is the lack of adequate information about their Ca(2+)-binding locations, which is largely hindered by the lack of a solved three-dimensional structure and rapid off rates due to low Ca(2+)-binding affinities. In this paper, we have reported the identification of three potential Ca(2+)-binding sites in a modeled CaSR structure using computational algorithms based on the geometric description and surface electrostatic potentials. Mutation of the predicted ligand residues in the full-length CaSR caused abnormal responses to [Ca(2+)](o), similar to those observed with naturally occurring activating or inactivating mutations of the CaR, supporting the essential role of these predicted Ca(2+)-binding sites in the sensing capability of the CaSR. In addition, to probe the intrinsic Ca(2+)-binding properties of the predicted sequences, we engineered two predicted continuous Ca(2+)-binding sequences individually into a scaffold protein provided by a non-Ca(2+)-binding protein, CD2. We report herein the estimation of the metal-binding affinities of these predicted sites in the CaSR by monitoring aromatic-sensitized Tb(3+) fluorescence energy transfer. Removing the predicted Ca(2+)-binding ligands resulted in the loss of or significantly weakened cation binding. The potential Ca(2+)-binding residues were shown to be involved in Ca(2+)/Ln(3+) binding by high resolution NMR and site-directed mutagenesis, further validating our prediction of Ca(2+)-binding sites within the extracellular domain of the CaSR.

Predicting calcium-binding sites in proteins - a graph theory and geometry approach

Identifying calcium-binding sites in proteins is one of the first steps towards predicting and understanding the role of calcium in biological systems for protein structure and function studies. Due to the complexity and irregularity of calcium-binding sites, a fast and accurate method for predicting and identifying calcium-binding protein is needed. Here we report our development of a new fast algorithm (GG) to detect calcium-binding sites. The GG algorithm uses a graph theory algorithm to find oxygen clusters of the protein and a geometric algorithm to identify the center of these clusters. A cluster of four or more oxygen atoms has a high potential for calcium binding. High performance with about 90% site sensitivity and 80% site selectivity has been obtained for three datasets containing a total of 123 proteins. The results suggest that a sphere of a certain size with four or more oxygen atoms on the surface and without other atoms inside is necessary and sufficient for quickly identifying the majority of the calcium-binding sites with high accuracy. Our finding opens a new avenue to visualize and analyze calcium-binding sites in proteins facilitating the prediction of functions from structural genomic information.

Prediction of water and metal binding sites and their affinities by using the Fold-X force field

The empirical force field Fold-X was developed previously to allow rapid free energy calculations in proteins. Here, we present an enhanced version of the force field allowing prediction of the position of structural water molecules and metal ions, together called single atom ligands. Fold-X picks up 76% of water molecules found to interact with two or more polar atoms of proteins in high-resolution crystal structures and predicts their position to within 0.8 A on average. The prediction of metal ion-binding sites have success rates between 90% and 97% depending on the metal, with an overall standard deviation on the position of binding of 0.3-0.6 A. The following metals were included in the force field: Mg2+, Ca2+, Zn2+, Mn2+, and Cu2+. As a result, the current version of Fold-X can accurately decorate a protein structure with biologically important ions and water molecules. Additionally, the free energy of binding of Ca2+ and Zn2+ (i.e., the natural logarithm of the dissociation constant) and its dependence on ionic strength correlate reasonably well with the experimental data available in the literature, allowing one to discriminate between high- and low-affinity binding sites. Importantly, the accuracy of the energy prediction presented here is sufficient to efficiently discriminate between Mg2+, Ca2+, and Zn2+ binding.

Altering protein specificity: techniques and applications

Protein engineering constitutes a powerful tool for generating novel proteins that serve as catalysts to induce selective chemical and biological transformations that would not otherwise be possible. Protocols that are commonly employed for altering the substrate specificity and selectivity profiles by mutating known enzymes include rational and random methods as well as techniques that entail evolution, selection and screening. Proteins identified by these techniques play important roles in a variety of industrial and medicinal applications and in the study of protein structure-function relationships. Herein we present a critical overview of methods for creating new functional proteins having altered specificity profiles and some practical case studies in which these techniques have been applied to solving problems in synthetic and medicinal chemistry and to elucidating enzyme function and biological pathways.

A united residue force-field for calcium-protein interactions

United-residue potentials are derived for interactions of the calcium cation with polypeptide chains in energy-based prediction of protein structure with a united-residue (UNRES) force-field. Specific potentials were derived for the interaction of the calcium cation with the Asp, Glu, Asn, and Gln side chains and the peptide group. The analytical expressions for the interaction energies for each of these amino acids were obtained by averaging the electrostatic interaction energy, expressed by a multipole series over the dihedral angles not considered in the united-residue model, that is, the side-chain dihedral angles chi and the dihedral angles lambda for the rotation of peptide groups about the C(alpha)...C(alpha) virtual-bond axes. For the side-chains that do not interact favorably with calcium, simple excluded-volume potentials were introduced. The parameters of the potentials were obtained from ab initio quantum mechanical calculations of model systems at the Restricted Hartree-Fock (RHF) level with the 6-31G(d,p) basis set. The energy surfaces of pairs consisting of Ca(2+)-acetate, Ca(2+)-propionate, Ca(2+)-acetamide, Ca(2+)-propionamide, and Ca(2+)-N-methylacetamide systems (modeling the Ca(2+)-Asp(-), Ca(2+)-Glu(-), Ca(2+)-Asn, Ca(2+)-Gln, and Ca(2+)-peptide group interactions) at different distances and orientations were calculated. For each pair, the restricted free energy (RFE) surfaces were calculated by numerical integration over the degrees of freedom lost when switching from the all-atom model to the united-residue model. Finally, the analytical expressions for each pair were fitted to the RFE surfaces. This force-field was able to distinguish the EF-hand motif from all potential binding sites in the crystal structures of bovine alpha-lactalbumin, whiting parvalbumin, calbindin D9K, and apo-calbindin D9K.

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

The control of cellular physiology and gene expression in response to extracellular signals is a basic property of living systems. We have constructed a synthetic bacterial signal transduction pathway in which gene expression is controlled by extracellular Zn2+. In this system a computationally designed Zn2+-binding periplasmic receptor senses the extracellular solute and triggers a two-component signal transduction pathway via a chimeric transmembrane protein, resulting in transcriptional up-regulation of a beta-galactosidase reporter gene. The Zn2+-binding site in the designed receptor is based on a four-coordinate, tetrahedral primary coordination sphere consisting of histidines and glutamates. In addition, mutations were introduced in a secondary coordination sphere to satisfy the residual hydrogen-bonding potential of the histidines coordinated to the metal. The importance of the secondary shell interactions is demonstrated by their effect on metal affinity and selectivity, as well as protein stability. Three designed protein sequences, comprising two distinct metal-binding positions, were all shown to bind Zn2+ and to function in the cell-based assay, indicating the generality of the design methodology. These experiments demonstrate that biological systems can be manipulated with computationally designed proteins that have drastically altered ligand-binding specificities, thereby extending the repertoire of genetic control by extracellular signals.