FRET-FCS detection of intralobe dynamics in calmodulin

Detection and quantification of Na,K-ATPase dimers in the plasma membrane of living cells by FRET-FCS

The sodium potassium pump, Na,K-ATPase (NKA), is an integral plasma membrane protein, expressed in all eukaryotic cells. It is responsible for maintaining the transmembrane Na+ gradient and is the major determinant of the membrane potential. Self-interaction and oligomerization of NKA in cell membranes has been proposed and discussed but is still an open question. Here, we have used a combination of FRET and Fluorescence Correlation Spectroscopy, FRET-FCS, to analyze NKA in the plasma membrane of living cells. Click chemistry was used to conjugate the fluorescent labels Alexa 488 and Alexa 647 to non-canonical amino acids introduced in the NKA α1 and β1 subunits. We demonstrate that FRET-FCS can detect an order of magnitude lower concentration of green-red labeled protein pairs in a single-labeled red and green background than what is possible with cross-correlation (FCCS). We show that a significant fraction of NKA is expressed as a dimer in the plasma membrane. We also introduce a method to estimate not only the number of single and double labeled NKA, but the number of unlabeled, endogenous NKA and estimate the density of endogenous NKA at the plasma membrane to 1400 ± 800 enzymes/μm2.

Energetic and structural insights behind calcium induced conformational transition in calmodulin

Calmodulin (CaM) is a key signaling protein that triggers several cellular and physiological processes inside the cell. Upon binding with calcium ion, CaM undergoes large scale conformational transition from a closed state to an open state that facilitates its interaction with various target protein and regulates their activity. This work explores the origin of the energetic and structural variation of the wild type and mutated CaM and explores the molecular origin for the structural differences between them. We first calculated the sequential calcium binding energy to CaM using the PDLD/S-LRA/β approach. This study  shows a very good correlation with experimental calcium binding energies. Next we calculated the calcium binding energies to the wild type CaM and several mutated CaM systems which were reported experimentally. On the structural aspect, it has been reported experimentally that certain mutation (Q41L-K75I) in calcium bound CaM leads to complete conformational transition from an open to a closed state. By using equilibrium molecular dynamics simulation, free energy calculation and contact frequency map analysis, we have shown that the formation of a cluster of long-range hydrophobic contacts, initiated by the Q41L-K75I CaM variant is the driving force behind its closing motion. This study unravels the energetics and structural aspects behind calcium ion induced conformational changes in wild type CaM and its variant.

Simple and Efficient Detection Scheme of Two-Color Fluorescence Correlation Spectroscopy for Protein Dynamics Investigation from Nanoseconds to Milliseconds

Nanosecond resolved fluorescence correlation spectroscopy (ns-FCS) based on two-color fluorescence detection is a powerful strategy for investigating the fast dynamics of biological macromolecules labeled with donor and acceptor fluorophores. The standard methods of ns-FCS use two single-photon avalanche diodes (SPADs) for the detection of single-color signals (four SPADs for two-color signals) to eliminate the afterpulse artifacts of SPAD at the expense of the efficiency of utilizing photon data in the calculation of correlograms. Herein, we demonstrated that hybrid photodetectors (HPDs) enable the recording of fluorescence photons in ns-FCS based on the minimal system using two HPDs for the detection of two-color signals. However, HPD exhibited afterpulses at a yield with respect to the rate of photodetection (<10-4) much lower than that of SPADs (∼10-2), which could still hamper correlation measurements. We demonstrated that the simple subtraction procedure could eliminate afterpulse artifacts. While the quantum efficiency of photodetection for HPDs is lower than that for high-performance SPADs, the developed system can be practically used for two-color ns-FCS in a time domain longer than a few nanoseconds. The fast chain dynamics of the B domain of protein A in the unfolded state was observed using the new method.

DNA facilitates heterodimerization between human transcription factors FoxP1 and FoxP2 by increasing their conformational flexibility

Transcription factors regulate gene expression by binding to DNA. They have disordered regions and specific DNA-binding domains. Binding to DNA causes structural changes, including folding and interactions with other molecules. The FoxP subfamily of transcription factors in humans is unique because they can form heterotypic interactions without DNA. However, it is unclear how they form heterodimers and how DNA binding affects their function. We used computational and experimental methods to study the structural changes in FoxP1's DNA-binding domain when it forms a heterodimer with FoxP2. We found that FoxP1 has complex and diverse conformational dynamics, transitioning between compact and extended states. Surprisingly, DNA binding increases the flexibility of FoxP1, contrary to the typical folding-upon-binding mechanism. In addition, we observed a 3-fold increase in the rate of heterodimerization after FoxP1 binds to DNA. These findings emphasize the importance of structural flexibility in promoting heterodimerization to form transcriptional complexes.

Role of water-bridged interactions in metal ion coupled protein allostery

Allosteric communication between distant parts of proteins controls many cellular functions, in which metal ions are widely utilized as effectors to trigger the allosteric cascade. Due to the involvement of strong coordination interactions, the energy landscape dictating the metal ion binding is intrinsically rugged. How metal ions achieve fast binding by overcoming the landscape ruggedness and thereby efficiently mediate protein allostery is elusive. By performing molecular dynamics simulations for the Ca²⁺ binding mediated allostery of the calmodulin (CaM) domains, each containing two Ca²⁺ binding helix-loop-helix motifs (EF-hands), we revealed the key role of water-bridged interactions in Ca²⁺ binding and protein allostery. The bridging water molecules between Ca²⁺ and binding residue reduces the ruggedness of ligand exchange landscape by acting as a lubricant, facilitating the Ca²⁺ coupled protein allostery. Calcium-induced rotation of the helices in the EF-hands, with the hydrophobic core serving as the pivot, leads to exposure of hydrophobic sites for target binding. Intriguingly, despite being structurally similar, the response of the two symmetrically arranged EF-hands upon Ca²⁺ binding is asymmetric. Breakage of symmetry is needed for efficient allosteric communication between the EF-hands. The key roles that water molecules play in driving allosteric transitions are likely to be general in other metal ion mediated protein allostery.

Natural proteins often utilize allostery in executing a variety of functions. Metal ions are typical cofactors to trigger the allosteric cascade. In this work, using the Ca²⁺ sensor protein calmodulin as the model system, we revealed crucial roles of water-bridged interactions in the metal ion coupled protein allostery. The coordination of the Ca²⁺ to the binding site involves an intermediate in which the water molecule bridges the Ca²⁺ and the liganding residue. The bridging water reduces the free energy barrier height of ligand exchange, therefore facilitating the ligand exchange and allosteric coupling by acting as a lubricant. We also showed that the response of the two symmetrically arranged EF-hand motifs of CaM domains upon Ca²⁺ binding is asymmetric, which is directly attributed to the differing dehydration process of the Ca²⁺ ions and is needed for efficient allosteric communication.

Response of Terahertz Protein Vibrations to Ligand Binding: Calmodulin-Peptide Complexes as a Case Study

Terahertz vibrations are sensitive reporters of the structure and interactions of proteins. Ligand binding alters the nature and distribution of these collective vibrations. The ligand-induced changes in the terahertz protein vibrations contribute to the binding entropy and to the overall thermodynamic stability of the resultant protein-ligand complexes. Here, we have examined the response of the low-frequency (below 6 terahertz) collective vibrations of the calcium-loaded calmodulin (CaM) to binding to five different ligands, both in the presence and absence of water, using normal-mode analysis and molecular dynamics simulations. A comparison of the vibrational spectra of hydrated and dry systems reveals that protein-solvent interactions stiffen the terahertz protein vibrations and that these solvent-coupled collective vibrations contribute significantly to the hydration-sensitive variation in the vibrational entropy of CaM. In the absence of water, the low-frequency vibrations of CaM are stiffened by ligand binding. On the contrary, the number and the cumulative vibrational entropy of low-frequency vibrational modes (ω < 200 cm^-1) of the hydrated CaM are increased noticeably after binding to the peptides, indicating binding-induced softening of collective vibrations of the protein. Although the calculated and experimental binding affinities of the chosen complexes correlated reasonably well, no systematic correlation was observed between the protein vibrational entropy and the binding affinity. The results underscored the importance of the interplay of protein-ligand and solvent interactions in modulating the low-frequency vibrations of proteins.

Molecular Brightness Approach for FRET Analysis of Donor-Linker-Acceptor Constructs at the Single Molecule Level: A Concept

In this report, we have developed a simple approach using single-detector fluorescence autocorrelation spectroscopy (FCS) to investigate the Förster resonance energy transfer (FRET) of genetically encoded, freely diffusing crTC2.1 (mTurquoise2.1-linker-mCitrine) at the single molecule level. We hypothesize that the molecular brightness of the freely diffusing donor (mTurquoise2.1) in the presence of the acceptor (mCitrine) is lower than that of the donor alone due to FRET. To test this hypothesis, the fluorescence fluctuation signal and number of molecules of freely diffusing construct were measured using FCS to calculate the molecular brightness of the donor, excited at 405 nm and detected at 475/50 nm, in the presence and absence of the acceptor. Our results indicate that the molecular brightness of cleaved crTC2.1 in a buffer is larger than that of the intact counterpart under 405-nm excitation. The energy transfer efficiency at the single molecule level is larger and more spread in values as compared with the ensemble-averaging time-resolved fluorescence measurements. In contrast, the molecular brightness of the intact crTC2.1, under 488 nm excitation of the acceptor (531/40 nm detection), is the same or slightly larger than that of the cleaved counterpart. These FCS-FRET measurements on freely diffusing donor-acceptor pairs are independent of the precise time constants associated with autocorrelation curves due to the presence of potential photophysical processes. Ultimately, when used in living cells, the proposed approach would only require a low expression level of these genetically encoded constructs, helping to limit potential interference with the cell machinery.

Clusterin inhibits Aβ42 aggregation through a "strawberry model" as detected by FRET-FCS

Extracellular plaque deposits of β-amyloid peptide (Aβ) are one of the main pathological features of Alzheimer's disease (AD). The aggregation of Aβ₄₂ species, especially Aβ₄₂ oligomers, is still an active research field in AD pathogenesis. Secretory clusterin protein (sCLU), an extracellular chaperone, plays an important role in AD pathogenesis. Although sCLU interacts directly with Aβ₄₂ in vitro and in vivo, the mechanism is not clear. In this paper, His-tagged sCLU (sCLU-His) was cloned, expressed and purified, and we applied florescence resonance energy transfer-fluorescence correlation spectroscopy (FRET-FCS) to investigate the direct interaction of sCLU-His and Aβ₄₂ at the single-molecule fluorescence level in vitro. Here, we chose four different fluorescently labeled Aβ₄₂ oligomers to form two different groups of aggregation models, easy or difficult to aggregate. The results showed that sCLU-His could form complexes with both aggregation models, and sCLU-His inhibited the aggregation of Aβ_42/RB  ~ Aβ_42/Atto647 (easy to aggregate model). The complexes were produced as the Aβ_42/Label adhered to the sCLU-His, which is similar to a "strawberry model," as strawberry seeds are dotted on the outer surface of strawberries. This work provided additional insight into the interaction mechanism of sCLU and Aβ₄₂ .

Intrinsically Disordered Regions of the DNA-Binding Domain of Human FoxP1 Facilitate Domain Swapping

Forkhead box P (FoxP) proteins are unique transcription factors that spatiotemporally regulate gene expression by tethering two chromosome loci together via functional domain-swapped dimers formed through their DNA-binding domains. Further, the differential kinetics on this dimerization mechanism underlie an intricate gene regulation network at physiological conditions. Nonetheless, poor understanding of the structural dynamics and steps of the association process impedes to link the functional domain swapping to human-associated diseases. Here, we have characterized the DNA-binding domain of human FoxP1 by integrating single-molecule Förster resonance energy transfer and hydrogen-deuterium exchange mass spectrometry data with molecular dynamics simulations. Our results confirm the formation of a previously postulated domain-swapped (DS) FoxP1 dimer in solution and reveal the presence of highly populated, heterogeneous, and locally disordered dimeric intermediates along the dimer dissociation pathway. The unique features of FoxP1 provide a glimpse of how intrinsically disordered regions can facilitate domain swapping oligomerization and other tightly regulated association mechanisms relevant in biological processes.

Stimulation-induced changes in diffusion and structure of calmodulin and calmodulin-dependent protein kinase II proteins in neurons

Calcium/calmodulin-dependent protein kinase II (CaMKII) and calmodulin (CaM) play essential roles in synaptic plasticity, which is an elementary process of learning and memory. In this study, fluorescence correlation spectroscopy (FCS) revealed diffusion properties of CaM, CaMKIIα and CaMKIIβ proteins in human embryonic kidney 293 (HEK293) cells and hippocampal neurons. A simultaneous multiple-point FCS recording system was developed on a random-access two-photon microscope, which facilitated efficient analysis of molecular dynamics in neuronal compartments. The diffusion of CaM in neurons was slower than that in HEK293 cells at rest, while the diffusion in stimulated neurons was accelerated and indistinguishable from that in HEK293 cells. This implied that activity-dependent binding partners of CaM exist in neurons, which slow down the diffusion at rest. Diffusion properties of CaMKIIα and β proteins implied that major populations of these proteins exist as holoenzymatic forms. Upon stimulation of neurons, the diffusion of CaMKIIα and β proteins became faster with reduced particle brightness, indicating drastic structural changes of the proteins such as dismissal from holoenzyme structure and further fragmentation.

The use of fluorescence correlation spectroscopy to characterize the molecular mobility of fluorescently labelled G protein-coupled receptors

The membranes of living cells have been shown to be highly organized into distinct microdomains, which has spatial and temporal consequences for the interaction of membrane bound receptors and their signalling partners as complexes. Fluorescence correlation spectroscopy (FCS) is a technique with single cell sensitivity that sheds light on the molecular dynamics of fluorescently labelled receptors, ligands or signalling complexes within small plasma membrane regions of living cells. This review provides an overview of the use of FCS to probe the real time quantification of the diffusion and concentration of G protein-coupled receptors (GPCRs), primarily to gain insights into ligand–receptor interactions and the molecular composition of signalling complexes. In addition we document the use of photon counting histogram (PCH) analysis to investigate how changes in molecular brightness (ε) can be a sensitive indicator of changes in molecular mass of fluorescently labelled moieties.

Comparing allosteric transitions in the domains of calmodulin through coarse-grained simulations

Calmodulin (CaM) is a ubiquitous Ca(2+)-binding protein consisting of two structurally similar domains with distinct stabilities, binding affinities, and flexibilities. We present coarse grained simulations that suggest that the mechanism for the domain's allosteric transitions between the open and closed conformations depends on subtle differences in the folded state topology of the two domains. Throughout a wide temperature range, the simulated transition mechanism of the N-terminal domain (nCaM) follows a two-state transition mechanism while domain opening in the C-terminal domain (cCaM) involves unfolding and refolding of the tertiary structure. The appearance of the unfolded intermediate occurs at a higher temperature in nCaM than it does in cCaM consistent with nCaM's higher thermal stability. Under approximate physiological conditions, the simulated unfolded state population of cCaM accounts for 10% of the population with nearly all of the sampled transitions (approximately 95%) unfolding and refolding during the conformational change. Transient unfolding significantly slows the domain opening and closing rates of cCaM, which can potentially influence its Ca(2+)-binding mechanism.

Enhanced conformational sampling technique provides an energy landscape view of large-scale protein conformational transitions

A novel in silico approach (NMA–ITS) is introduced to rapidly and effectively sample the configuration space and give quantitative data for exploring the conformational changes of proteins.

Highly Sensitive FRET-FCS Detects Amyloid β-Peptide Oligomers in Solution at Physiological Concentrations

Oligomers formed by the amyloid β-peptide (Aβ) are pathogens in Alzheimer's disease. Increased knowledge on the oligomerization process is crucial for understanding the disease and for finding treatments. Ideally, Aβ oligomerization should be studied in solution and at physiologically relevant concentrations, but most popular techniques of today are not capable of such analyses. We demonstrate here that the combination of Förster Resonance Energy Transfer and Fluorescence Correlation Spectroscopy (FRET-FCS) has a unique ability to detect small subpopulations of FRET-active molecules and oligomers. FRET-FCS could readily detect a FRET-active oligonucleotide present at levels as low as 0.5% compared to FRET-inactive dye molecules. In contrast, three established fluorescence fluctuation techniques (FCS, FCCS, and PCH) required fractions between 7 and 11%. When applied to the analysis of Aβ, FRET-FCS detected oligomers consisting of less than 10 Aβ molecules, which coexisted with the monomers at fractions as low as 2 ± 2%. Thus, we demonstrate for the first time direct detection of small fractions of Aβ oligomers in solution at physiological concentrations. This ability of FRET-FCS could be an indispensable tool for studying biological oligomerization processes, in general, and for finding therapeutically useful oligomerization inhibitors.

Microsecond protein dynamics observed at the single-molecule level

How polypeptide chains acquire specific conformations to realize unique biological functions is a central problem of protein science. Single-molecule spectroscopy, combined with fluorescence resonance energy transfer, is utilized to study the conformational heterogeneity and the state-to-state transition dynamics of proteins on the submillisecond to second timescales. However, observation of the dynamics on the microsecond timescale is still very challenging. This timescale is important because the elementary processes of protein dynamics take place and direct comparison between experiment and simulation is possible. Here we report a new single-molecule technique to reveal the microsecond structural dynamics of proteins through correlation of the fluorescence lifetime. This method, two-dimensional fluorescence lifetime correlation spectroscopy, is applied to clarify the conformational dynamics of cytochrome c. Three conformational ensembles and the microsecond transitions in each ensemble are indicated from the correlation signal, demonstrating the importance of quantifying microsecond dynamics of proteins on the folding free energy landscape.

Single molecule spectroscopy can visualise dynamic changes in protein conformation on the submillisecond timescale. Here, Otosu et al. apply two-dimensional fluorescence lifetime correlation spectroscopy to visualise dynamics between seven conformers of cytochrome c on the microsecond timescale.

The effect of intrachain electrostatic repulsion on conformational disorder and dynamics of the Sic1 protein

The yeast cyclin-dependent kinase inhibitor Sic1 is a disordered protein that, upon multisite phosphorylation, forms a dynamic complex with the Cdc4 subunit of an SCF ubiquitin ligase. To understand the multisite phosphorylation dependence of the Sic1:Cdc4 interaction, which ultimately leads to a sharp cell cycle transition, the conformational properties of the disordered Sic1 N-terminal targeting region were studied using single-molecule fluorescence spectroscopy. Multiple conformational populations with different sensitivities to charge screening were identified by performing experiments in nondenaturing salts and ionic denaturants. Both the end-to-end distance and the hydrodynamic radius decrease monotonically with increasing the salt concentration, and a rollover of the chain dimensions in high denaturant conditions is observed. The data were fit to the polyelectrolyte binding-screening model, yielding parameters such as the excluded volume of the uncharged chain and the binding constant to denaturant. An overall scaling factor of ∼1.2 was needed for fitting the data, which implies that Sic1 cannot be approximated by a random Gaussian chain. Fluorescence correlation spectroscopy reveals Sic1 structure fluctuations occurring on both fast (10-100 ns) and slow (∼10 ms) time scales, with the fast phase absent in low salt solutions. The results of this study provide direct evidence that long-range intrachain electrostatic repulsions are a significant factor for the conformational landscape of Sic1, and support the role of electrostatics in determining the overall shape and hydrodynamic properties of intrinsically disordered proteins.

Sub-diffusion decays in fluorescence correlation spectroscopy: dye photophysics or protein dynamics?

Transitions between bright and dark fluorescent states of several rhodamine dyes were investigated by fluorescence correlation spectroscopy. We resolved two sub-diffusion exponential decays for free rhodamines in aqueous solutions, of which the slower component scales linearly with the viscosity of the solution. Correlation data for proteins and DNA labeled with tetramethylrhodamine were fitted with three to four exponential decays describing flickering dynamics on a time scale between 0.5 and 100 μs. We investigated the nature of these processes by performing experiments under different experimental conditions and for different samples. On the basis of how their population and lifetime change with viscosity, the oxygen content of the solution, the laser irradiance, and the detection geometry, we assigned these states, in the order of increasing lifetimes, to a triplet state, a hybrid between twisted-intramolecular-charge-transfer state and a ground state lactonic state, a lactonic state, and a photoionized state, respectively. Our data suggests that none of the observed sub-diffusion correlation decays can be directly assigned to the intramolecular dynamics of the labeled biomolecules. However, we found evidence that the intrinsic conformational dynamics of the biomolecule appears in the correlation curves as a modulation of the photophysics of the dye label. This shows the importance of accurate control measurements and appropriate modeling of the dye photophysics in fluorescence correlation studies, and it cautions against direct assignments of dark-state relaxation times to folding kinetics in proteins and nucleic acids.

Analyzing Förster resonance energy transfer with fluctuation algorithms

Fluorescence correlation spectroscopy (FCS) in combination with Förster resonance energy transfer (FRET) has been developed to a powerful statistical tool, which allows for the analysis of FRET fluctuations in the huge time of nanoseconds to seconds. FRET-FCS utilizes the strong distance dependence of the FRET efficiency on the donor (D)-acceptor (A) distance so that it developed to a perfect method for studying structural fluctuation in biomolecules involved in conformational flexibility, structural dynamics, complex formation, folding, and catalysis. Structural fluctuations thereby result in anticorrelated donor and acceptor signals, which are analyzed by FRET-FCS in order to characterize underlying structural dynamics. Simulated and experimental examples are discussed. First, we review experimental implementations of FRET-FCS and present theory for a two-state interconverting system. Additionally, we consider a very common case of FRET dynamics in the presence of donor-only labeled species. We demonstrate that the mean relaxation time for the structural dynamics can be easily obtained in most of cases, whereas extracting meaningful information from correlation amplitudes can be challenging. We present a strategy to avoid a fit with an underdetermined model function by restraining the D and A brightnesses of the at least one involved state, so that both FRET efficiencies and both rate constants (i.e., the equilibrium constant) can be determined. For samples containing several fluorescent species, the use of pulsed polarized excitation with multiparameter fluorescence detection allows for filtered FCS (fFCS), where species-specific correlation functions can be obtained, which can be directly interpreted. The species selection is achieved by filtering using fluorescence decays of individual species. Analytical functions for species auto- and cross-correlation functions are given. Moreover, fFCS is less affected by photophysical artifacts and often offers higher contrast, which effectively increases its time resolution and significantly enhances its capability to resolve multistate kinetics. fFCS can also differentiate between species even when their brightnesses are the same and thus opens up new possibilities to characterize complex dynamics. Alternative fluctuation algorithms to study FRET dynamics are also briefly reviewed.

Direct observation of T4 lysozyme hinge-bending motion by fluorescence correlation spectroscopy

Bacteriophage T4 Lysozyme (T4L) catalyzes the hydrolysis of the peptidoglycan layer of the bacterial cell wall late in the infection cycle. It has long been postulated that equilibrium dynamics enable substrate access to the active site located at the interface between the N- and C-terminal domains. Crystal structures of WT-T4L and point mutants captured a range of conformations that differ by the hinge-bending angle between the two domains. Evidence of equilibrium between open and closed conformations in solution was gleaned from distance measurements between the two domains but the nature of the equilibrium and the timescale of the underlying motion have not been investigated. Here, we used fluorescence fluctuation spectroscopy to directly detect T4L equilibrium conformational fluctuations in solution. For this purpose, Tetramethylrhodamine probes were introduced at pairs of cysteines in regions of the molecule that undergo relative displacement upon transition from open to closed conformations. Correlation analysis of Tetramethylrhodamine intensity fluctuations reveals hinge-bending motion that changes the relative distance and orientation of the N- and C-terminal domains with ≅ 15 μs relaxation time. That this motion involves interconversion between open and closed conformations was further confirmed by the dampening of its amplitude upon covalent substrate trapping. In contrast to the prevalent two-state model of T4L equilibrium, molecular brightness and number of particles obtained from cumulant analysis suggest that T4L populates multiple intermediate states, consistent with the wide range of hinge-bending angles trapped in the crystal structure of T4L mutants.

Conformational flexibility and the mechanisms of allosteric transitions in topologically similar proteins

Conformational flexibility plays a central role in allosteric transition of proteins. In this paper, we extend the analysis of our previous study [S. Tripathi and J. J. Portman, Proc. Natl. Acad. Sci. U.S.A. 106, 2104 (2009)] to investigate how relatively minor structural changes of the meta-stable states can significantly influence the conformational flexibility and allosteric transition mechanism. We use the allosteric transitions of the domains of calmodulin as an example system to highlight the relationship between the transition mechanism and the inter-residue contacts present in the meta-stable states. In particular, we focus on the origin of transient local unfolding (cracking), a mechanism that can lower free energy barriers of allosteric transitions, in terms of the inter-residue contacts of the meta-stable states and the pattern of local strain that develops during the transition. We find that the magnitude of the local strain in the protein is not the sole factor determining whether a region will ultimately crack during the transition. These results emphasize that the residue interactions found exclusively in one of the two meta-stable states is the key in understanding the mechanism of allosteric conformational change.