Allosteric inactivation of an engineered optogenetic GTPase.

Optogenetics is a technique for establishing direct spatiotemporal control over molecular function within living cells using light. Light application induces conformational changes within targeted proteins that produce changes in function. One of the applications of optogenetic tools is an allosteric control of proteins via light-sensing domain (LOV2), which allows direct and robust control of protein function. Computational studies supported by cellular imaging demonstrated that application of light allosterically inhibited signaling proteins Vav2, ITSN, and Rac1, but the structural and dynamic basis of such control has yet to be elucidated by experiment. Here, using NMR spectroscopy, we discover principles of action of allosteric control of cell division control protein 42 (CDC42), a small GTPase involved in cell signaling. Both LOV2 and Cdc42 employ flexibility in their function to switch between "dark"/"lit" or active/inactive states, respectively. By conjoining Cdc42 and phototropin1 LOV2 domains into the bi-switchable fusion Cdc42Lov, application of light-or alternatively, mutation in LOV2 to mimic light absorption-allosterically inhibits Cdc42 downstream signaling. The flow and patterning of allosteric transduction in this flexible system are well suited to observation by NMR. Close monitoring of the structural and dynamic properties of dark versus "lit" states of Cdc42Lov revealed lit-induced allosteric perturbations that extend to Cdc42's downstream effector binding site. Chemical shift perturbations for lit mimic, I539E, have distinct regions of sensitivity, and both the domains are coupled together, leading to bidirectional interdomain signaling. Insights gained from this optoallosteric design will increase our ability to control response sensitivity in future designs.

Mixed, nonclassical behavior in a classic allosteric protein.

Allostery is a major driver of biological processes requiring coordination. Thus, it is one of the most fundamental and remarkable phenomena in nature, and there is motivation to understand and manipulate it to a multitude of ends. Today, it is often described in terms of two phenomenological models proposed more than a half-century ago involving only T(tense) or R(relaxed) conformations. Here, methyl-based NMR provides extensive detail on a dynamic T to R switch in the classical dimeric allosteric protein, yeast chorismate mutase (CM), that occurs in the absence of substrate, but only with the activator bound. Switching of individual subunits is uncoupled based on direct observation of mixed TR states in the dimer. This unique finding excludes both classic models and solves the paradox of a coexisting hyperbolic binding curve and highly skewed substrate-free T-R equilibrium. Surprisingly, structures of the activator-bound and effector-free forms of CM appear the same by NMR, providing another example of the need to account for dynamic ensembles. The apo enzyme, which has a sigmoidal activity profile, is shown to switch, not to R, but to a related high-energy state. Thus, the conformational repertoire of CM does not just change as a matter of degree depending on the allosteric input, be it effector and/or substrate. Rather, the allosteric model appears to completely change in different contexts, which is only consistent with modern ensemble-based frameworks.

Positive Cooperativity in Substrate Binding by Human Thymidylate Synthase.

Thymidylate synthase (TS) catalyzes the production of the nucleotide dTMP from deoxyuridine monophosphate (dUMP), making the enzyme necessary for DNA replication and consequently a target for cancer therapeutics. TSs are homodimers with active sites separated by ∼30 Å. Reports of half-the-sites activity in TSs from multiple species demonstrate the presence of allosteric communication between the active sites of this enzyme. A simple explanation for the negative allosteric regulation occurring in half-the-sites activity would be that the two substrates bind with negative cooperativity. However, previous work on Escherichia coli TS revealed that dUMP substrate binds without cooperativity. To gain further insight into TS allosteric function, binding cooperativity in human TS is examined here. Isothermal titration calorimetry and two-dimensional lineshape analysis of NMR titration spectra are used to characterize the thermodynamics of dUMP binding, with a focus on quantification of cooperativity between the two substrate binding events. We find that human TS binds dUMP with ∼9-fold entropically driven positive cooperativity (ρITC = 9 ± 1, ρNMR = 7 ± 1), in contrast to the apparent strong negative cooperativity reported previously. Our work further demonstrates the necessity of globally fitting isotherms collected under various conditions, as well as accurate determination of binding competent protein concentration, for calorimetric characterization of homotropic cooperative binding. Notably, an initial curvature of the isotherm is found to be indicative of positively cooperative binding. Two-dimensional lineshape analysis NMR is also found to be an informative tool for quantifying binding cooperativity, particularly in cases in which bound intermediates yield unique resonances.

Inter-Active Site Communication Mediated by the Dimer Interface ß-Sheet in the Half-the-Sites Enzyme, Thymidylate Synthase.

Thymidylate synthase (TS) is a dimeric enzyme conserved in all life forms that exhibits the allosteric feature of half-the-sites activity. Neither the reason for nor the mechanism of this phenomenon is understood. We used a combined nuclear magnetic resonance (NMR) and molecular dynamics approach to study a stable intermediate preceding hydride transfer, which is the rate-limiting and half-the-sites step. In NMR titrations with ligands leading to this intermediate, we measured chemical shifts of the apoenzyme (lig0), the saturated holoenzyme (lig2), and the typically elusive singly bound (lig1) states. Approximately 40 amides showed quartet patterns providing direct NMR evidence of coupling between the active site and probes >30 Å away in the distal subunit. Quartet peak patterns have symmetrical character, indicating reciprocity in communicating the first and second binding events to the distal protomer. Quartets include key catalytic residues and map to the dimer interface ß-sheet, which also represents the shortest path between the two active sites. Simulations corroborate the coupling observed in solution in that there is excellent overlap between quartet residues and main-chain atoms having intersubunit cross-correlated motions. Simulations identify five hot spot residues, three of which lie at the kink in the unique ß-bulge abutting the active sites on either end of the sheet. Interstrand cross-correlated motions become more organized and pronounced as the enzyme progresses from lig0 to lig1 and ultimately lig2. Coupling in the apparently symmetrical complex has implications for half-the-sites reactivity and potentially resolves the paradox of inequivalent TS active sites despite the vast majority of X-ray structures appearing to be symmetrical.

Chemical shift imprint of intersubunit communication in a symmetric homodimer.

Allosteric communication is critical for protein function and cellular homeostasis, and it can be exploited as a strategy for drug design. However, unlike many protein-ligand interactions, the structural basis for the long-range communication that underlies allostery is not well understood. This lack of understanding is most evident in the case of classical allostery, in which a binding event in one protomer is sensed by a second symmetric protomer. A primary reason why study of interdomain signaling is challenging in oligomeric proteins is the difficulty in characterizing intermediate, singly bound species. Here, we use an NMR approach to isolate and characterize a singly ligated state ("lig1") of a homodimeric enzyme that is otherwise obscured by rapid exchange with apo and saturated forms. Mixed labeled dimers were prepared that simultaneously permit full population of the lig1 state and isotopic labeling of either protomer. Direct visualization of peaks from lig1 yielded site-specific ligand-state multiplets that provide a convenient format for assessing mechanisms of intersubunit communication from a variety of NMR measurements. We demonstrate this approach on thymidylate synthase from Escherichia coli, a homodimeric enzyme known to be half-the-sites reactive. Resolving the dUMP1 state shows that active site communication occurs not upon the first dUMP binding, but upon the second. Surprisingly, for many sites, dUMP1 peaks are found beyond the limits set by apo and dUMP2 peaks, indicating that binding the first dUMP pushes the enzyme ensemble to further conformational extremes than the apo or saturated forms. The approach used here should be generally applicable to homodimers.

Protein Mass Effects on Formate Dehydrogenase.

Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.

Widespread Perturbation of Function, Structure, and Dynamics by a Conservative Single-Atom Substitution in Thymidylate Synthase.

Thymidylate synthase (TSase) is responsible for synthesizing the sole de novo source of dTMP in all organisms. TSase is a drug target, and as such, it has been well studied in terms of both structure and reaction mechanism. Cysteine 146 in Escherichia coli TSase is universally conserved because it serves as the nucleophile in the enzyme mechanism. Here we use the C146S mutation to probe the role of the sulfur atom in early events in the catalytic cycle beyond serving as the nucleophile. Surprisingly, the single-atom substitution severely decreases substrate binding affinity, and the unfavorable ΔΔG°bind is comprised of roughly equal enthalpic and entropic components at 25 °C. Chemical shifts in the free and dUMP-bound states show the mutation causes perturbations throughout TSase, including regions important for complex stability, in agreement with a less favorable enthalpy change. We measured the nuclear magnetic resonance methyl symmetry axis order parameter (S2axis), a proxy for conformational entropy, for TSase at all vertices of the dUMP binding/C146S mutation thermodynamic cycle and found that the calculated TΔΔS°conf is similar in sign and magnitude to the calorimetric TΔΔS°. Further, we ascribed minor resonances in wild-type-dUMP spectra to a state with a covalent bond between Sγ of C146 and C6 of dUMP and find S2axis values are unaffected by covalent bond formation, indicating this reaction step is neutral with respect to ΔS°conf. Lastly, the C146S mutation allowed us to measure cofactor analog binding by isothermal titration calorimetry without the confounding heat signature of covalent bond formation. Raltitrexed binds free and singly bound TSase with similar affinities, yet the two binding events have different enthalpy changes, providing further evidence of communication between the two active sites.

The Effect of Protein Mass Modulation on Human Dihydrofolate Reductase.

Dihydrofolate reductase (DHFR) from Escherichia coli has long served as a model enzyme with which to elucidate possible links between protein dynamics and the catalyzed reaction. Such physical properties of its human counterpart have not been rigorously studied so far, but recent computer-based simulations suggest that these two DHFRs differ significantly in how closely coupled the protein dynamics and the catalyzed C-H → C hydride transfer step are. To test this prediction, two contemporary probes for studying the effect of protein dynamics on catalysis were combined here: temperature dependence of intrinsic kinetic isotope effects (KIEs), which are sensitive to the physical nature of the chemical step, and protein mass modulation, which slows down fast dynamics (femto- to picosecond time scale) throughout the protein. The intrinsic H/T KIEs of human DHFR, like those of E. coli DHFR, are shown to be temperature-independent in the range from 5 to 45 °C, indicating fast sampling of donor and acceptor distances (DADs) at the reaction's transition state (or tunneling ready state, TRS). Mass modulation of these enzymes through isotopic labeling with (13)C, (15)N, and (2)H at nonexchangeable hydrogens yields an 11% heavier enzyme. The additional mass has no effect on the intrinsic KIEs of the human enzyme. This finding indicates that the mass modulation of the human DHFR affects neither DAD distribution nor the DAD's conformational sampling dynamics. Furthermore, reduction in the enzymatic turnover number and the dissociation rate constant for the product indicate that the isotopic substitution affects kinetic steps that are not the catalyzed C-H → C hydride transfer. The findings are discussed in terms of fast dynamics and their role in catalysis, the comparison of calculations and experiments, and the interpretation of isotopically modulated heavy enzymes in general.

Bacterial Thymidylate Synthase Binds Two Molecules of Substrate and Cofactor without Cooperativity.

Thymidylate synthase (TSase) is a clinically important enzyme because it catalyzes synthesis of the sole de novo source of deoxy-thymidylate. Without this enzyme, cells die a "thymineless death" since they are starved of a crucial DNA synthesis precursor. As a drug target, TSase is well studied in terms of its structure and reaction mechanism. An interesting mechanistic feature of dimeric TSase is that it is "half-the-sites reactive", which is a form of negative cooperativity. Yet, the basis for this is not well-understood. Some experiments point to cooperativity at the binding steps of the reaction cycle as being responsible for the phenomenon, but the literature contains conflicting reports. Here we use ITC and NMR to resolve these inconsistencies. This first detailed thermodynamic dissection of multisite binding of dUMP to E. coli TSase shows the nucleotide binds to the free and singly bound forms of the enzyme with nearly equal affinity over a broad range of temperatures and in multiple buffers. While small but significant differences in ΔC°P for the two binding events show that the active sites are not formally equivalent, there is little-to-no allostery at the level of ΔG°bind. In addition NMR titration data reveal that there is minor intersubunit cooperativity in formation of a ternary complex with the mechanism based inhibitor, 5F-dUMP, and cofactor. Taken together, the data show that functional communication between subunits is minimal for both binding steps of the reaction coordinate.

Evidence for dynamics in proteins as a mechanism for ligand dissociation.

Signal transduction, regulatory processes and pharmaceutical responses are highly dependent upon ligand residence times. Gaining insight into how physical factors influence residence times (1/k(off)) should enhance our ability to manipulate biological interactions. We report experiments that yield structural insight into k(off) involving a series of eight 2,4-diaminopyrimidine inhibitors of dihydrofolate reductase whose binding affinities vary by six orders of magnitude. NMR relaxation-dispersion experiments revealed a common set of residues near the binding site that undergo a concerted millisecond-timescale switching event to a previously unidentified conformation. The rate of switching from ground to excited conformations correlates exponentially with the binding affinity K(i) and k(off), suggesting that protein dynamics serves as a mechanical initiator of ligand dissociation within this series and potentially for other macromolecule-ligand systems. Although the forward rate of conformational exchange, k(conf,forward), is faster than k(off), the use of the ligand series allowed for connections to be drawn between kinetic events on different timescales.

Mg2+ binds to the surface of thymidylate synthase and affects hydride transfer at the interior active site.

Thymidylate synthase (TSase) produces the sole intracellular de novo source of thymidine (i.e., the DNA base T) and thus is a common target for antibiotic and anticancer drugs. Mg(2+) has been reported to affect TSase activity, but the mechanism of this interaction has not been investigated. Here we show that Mg(2+) binds to the surface of Escherichia coli TSase and affects the kinetics of hydride transfer at the interior active site (16 Å away). Examination of the crystal structures identifies a Mg(2+) near the glutamyl moiety of the folate cofactor, providing the first structural evidence for Mg(2+) binding to TSase. The kinetics and NMR relaxation experiments suggest that the weak binding of Mg(2+) to the protein surface stabilizes the closed conformation of the ternary enzyme complex and reduces the entropy of activation on the hydride transfer step. Mg(2+) accelerates the hydride transfer by ~7-fold but does not affect the magnitude or temperature dependence of the intrinsic kinetic isotope effect. These results suggest that Mg(2+) facilitates the protein motions that bring the hydride donor and acceptor together, but it does not change the tunneling ready state of the hydride transfer. These findings highlight how variations in cellular Mg(2+) concentration can modulate enzyme activity through long-range interactions in the protein, rather than binding at the active site. The interaction of Mg(2+) with the glutamyl tail of the folate cofactor and nonconserved residues of bacterial TSase may assist in designing antifolates with polyglutamyl substitutes as species-specific antibiotic drugs.

Phosphorylation of a PDZ domain extension modulates binding affinity and interdomain interactions in postsynaptic density-95 (PSD-95) protein, a membrane-associated guanylate kinase (MAGUK).

Postsynaptic density-95 is a multidomain scaffolding protein that recruits glutamate receptors to postsynaptic sites and facilitates signal processing and connection to the cytoskeleton. It is the leading member of the membrane-associated guanylate kinase family of proteins, which are defined by the PSD-95/Discs large/ZO-1 (PDZ)-Src homology 3 (SH3)-guanylate kinase domain sequence. We used NMR to show that phosphorylation of conserved tyrosine 397, which occurs in vivo and is located in an atypical helical extension (α3), initiates a rapid equilibrium of docked and undocked conformations. Undocking reduced ligand binding affinity allosterically and weakened the interaction of PDZ3 with SH3 even though these domains are separated by a ~25-residue linker. Additional phosphorylation at two linker sites further disrupted the interaction, implicating α3 and the linker in tuning interdomain communication. These experiments revealed a novel mode of regulation by a detachable PDZ element and offer a first glimpse at the dynamic interaction of PDZ and SH3-guanylate kinase domains in membrane-associated guanylate kinases.

Hidden dynamic allostery in a PDZ domain.

Structure-function relationships in proteins are predicated on the spatial proximity of noncovalently interacting groups of atoms. Thus, structural elements located away from a protein's active site are typically presumed to serve a stabilizing or scaffolding role for the larger structure. Here we report a functional role for a distal structural element in a PDZ domain, even though it is not required to maintain PDZ structure. The third PDZ domain from PSD-95/SAP90 (PDZ3) has an unusual additional third alpha helix (alpha3) that packs in contiguous fashion against the globular domain. Although alpha3 lies outside the active site and does not make direct contact with C-terminal peptide ligand, removal of alpha3 reduces ligand affinity by 21-fold. Further investigation revealed that the difference in binding free energies between the full-length and truncated constructs is predominantly entropic in nature and that without alpha3, picosecond-nanosecond side-chain dynamics are enhanced throughout the domain, as determined by (2)H methyl NMR relaxation. Thus, the distal modulation of binding function appears to occur via a delocalized conformational entropy mechanism. Without removal of alpha3 and characterization of side-chain dynamics, this dynamic allostery would have gone unnoticed. Moreover, what appeared at first to be an artificial modification of PDZ3 has been corroborated by experimentally verified phosphorylation of alpha3, revealing a tangible biological mechanism for this novel regulatory scheme. This hidden dynamic allostery raises the possibility of as-yet unidentified or untapped allosteric regulation in this PDZ domain and is a very clear example of function arising from dynamics rather than from structure.

Dynamic allostery in substrate binding by human thymidylate synthase.

Human thymidylate synthase (hTS) is essential for DNA replication and therefore a therapeutic target for cancer. Effective targeting requires knowledge of the mechanism(s) of regulation of this 72 kDa homodimeric enzyme. Here, we investigate the mechanism of binding cooperativity of the nucleotide substrate. We have employed exquisitely sensitive methyl-based CPMG and CEST NMR experiments enabling us to identify residues undergoing bifurcated linear 3-state exchange, including concerted switching between active and inactive conformations in the apo enzyme. The inactive state is populated to only ~1.3%, indicating that conformational selection contributes negligibly to the cooperativity. Instead, methyl rotation axis order parameters, determined by 2H transverse relaxation rates, suggest that rigidification of the enzyme upon substrate binding is responsible for the entropically-driven cooperativity. Lack of the rigidification in product binding and substrate binding to an N-terminally truncated enzyme, both non-cooperative, support this idea. In addition, the lack of this rigidification in the N-terminal truncation indicates that interactions between the flexible N-terminus and the rest of the protein, which are perturbed by substrate binding, play a significant role in the cooperativity-a novel mechanism of dynamic allostery. Together, these findings yield a rare depth of insight into the substrate binding cooperativity of an essential enzyme.

Exploring the role of structure and dynamics in the function of chymotrypsin inhibitor 2.

Increasing awareness of the possible role of internal dynamics in protein function has led to the development of new methods for experimentally characterizing protein dynamics across multiple time scales, especially using NMR spectroscopy. A few analyses of the conformational dynamics of proteins ranging from nonallosteric single domains to multidomain allosteric enzymes are now available; however, demonstrating a connection between dynamics and function remains difficult on account of the comparative lack of studies examining both changes in dynamics and changes in function in response to the same perturbations. In previous work, we characterized changes in structure and dynamics on the ps­ns time scale resulting from hydrophobic core mutations in chymotrypsin inhibitor 2 and found that there are moderate, persistent global changes in dynamics in the absence of gross structural changes (Whitley et al., Biochemistry 2008;47:8566­8576). Here, we assay those and additional mutants for inhibitory ability toward the serine proteases elastase and chymotrypsin to determine the effects of mutation on function. Results indicate that core mutation has only a subtle effect on CI2 function. Using chemical shifts, we also studied the effect of complex formation on CI2 structure and found that perturbations are greatest at the complex interface but also propagate toward CI2's hydrophobic core. The structure­dynamics­function data set completed here suggests that dynamics plays a limited role in the function of this small model system, although we do observe a correlation between nanosecond-scale reactive loop motions and inhibitory ability for mutations at one key position in the hydrophobic core.

Direct detection of structurally resolved dynamics in a multiconformation receptor-ligand complex.

Structure-based drug design relies on static protein structures despite significant evidence for the need to include protein dynamics as a serious consideration. In practice, dynamic motions are neglected because they are not understood well enough to model, a situation resulting from a lack of explicit experimental examples of dynamic receptor-ligand complexes. Here, we report high-resolution details of pronounced ~1 ms time scale motions of a receptor-small molecule complex using a combination of NMR and X-ray crystallography. Large conformational dynamics in Escherichia coli dihydrofolate reductase are driven by internal switching motions of the drug-like, nanomolar-affinity inhibitor. Carr-Purcell-Meiboom-Gill relaxation dispersion experiments and NOEs revealed the crystal structure to contain critical elements of the high energy protein-ligand conformation. The availability of accurate, structurally resolved dynamics in a protein-ligand complex should serve as a valuable benchmark for modeling dynamics in other receptor-ligand complexes and prediction of binding affinities.

Dynamic dysfunction in dihydrofolate reductase results from antifolate drug binding: modulation of dynamics within a structural state.

The arduous task of rationally designing small-molecule enzyme inhibitors is complicated by the inherent flexibility of the protein scaffold. To gain insight into the changes in dynamics associated with small-molecule-based inhibition, we have characterized, using NMR spectroscopy, Escherichia coli dihydrofolate reductase in complex with two drugs: methotrexate and trimethoprim. The complexes allowed the intrinsic dynamic effects of drug binding to be revealed within the context of the "closed" structural ensemble. Binding of both drugs results in an identical decoupling of global motion on the micro- to millisecond timescale. Consistent with a change in overall dynamic character, the drugs' perturbations to pico- to nanosecond backbone and side-chain methyl dynamics are also highly similar. These data show that the inhibitors simultaneously modulate slow concerted switching and fast motions at distal regions of dihydrofolate reductase, providing a dynamic link between the substrate binding site and distal loop residues known to affect catalysis.

Backbone and ILVM methyl resonance assignments of human thymidylate synthase in apo and substrate bound forms.

Human thymidylate synthase (hTS) is a 72 kDa homodimeric enzyme responsible for the conversion of deoxyuridine monophosphate (dUMP) to deoxythymidine monophosphate (dTMP), making it the sole source of de novo dTMP in human cells. As a result, hTS is an attractive anti-cancer therapeutic target. Additionally, hTS is known to possess a number of interesting biophysical features, including adoption of active and inactive conformations, positively cooperative substrate binding, half-the-sites activity, and interacting with its own mRNA. The physical mechanisms underlying these properties, and how they may be leveraged to guide therapeutic development, are yet to be fully explored. Here, as a preface to detailed NMR characterization, we present backbone amide and ILVM methyl resonance assignments for hTS in apo and dUMP bound forms. In addition, we present backbone amide resonance assignments for hTS bound to a substrate analog and the native cofactor.

Visualizing an Allosteric Intermediate Using CuAAC Stabilization of an NMR Mixed Labeled Dimer.

Homodimers are the most abundant type of enzyme in cells, and as such, they represent the most elemental system for studying the phenomenon of allostery. In these systems, in which the allosteric features are manifest by the effect of the first binding event on a similar event at the second site, the most informative state is the asymmetric singly bound (lig1) form, yet it tends to be thermodynamically elusive. Here we obtain milligram quantities of lig1 of the allosteric homodimer, chorismate mutase, in the form of a mixed isotopically labeled dimer stabilized by Cu(I)-catalyzed azide-alkyne cycloaddition (CuAAC) between the subunits. Below, we outline several critical steps required to generate high yields of both types of unnatural amino acid-containing proteins and overcome multiple pitfalls intrinsic to CuAAC to obtain high yields of a highly purified, fully intact, active mixed labeled dimer, which provides the first glimpse of the lig1 intermediate. These data not only will make possible NMR-based investigations of allostery envisioned by us but also should facilitate other structural applications in which specific linkage of proteins is helpful.