CR1/CR2 interactions modulate the functions of the cell surface epidermal growth factor receptor

Recent crystallographic data on the isolated extracellular domain of the epidermal growth factor receptor (EGFR) have suggested a model for its activation by ligand. We have tested this model in the context of the full-length EGFR displayed at the cell surface, by introducing mutations in two regions (CR1 and CR2) of the extracellular domain thought to be critical for regulation of receptor activation. Mutations in the CR1 and CR2 domains have opposing effects on ligand binding affinity, receptor dimerization, tyrosine kinase activation, and signaling competence. Tyr(246) is a critical residue in the CR1 loop, which is implicated in the positioning and stabilization of the receptor dimer interface after ligand binding; mutations of Tyr(246) impair or abolish receptor function. Mutations in CR2, which weaken the interaction that restricts the receptor to the tethered (inactive) state, enhance responsiveness to EGF by increasing affinity for the ligand. However, weakening of the CR1/CR2 interaction does not result in spontaneous activation of the receptors' kinase. We have used an antibody (mAb 806), which recognizes a transition state of the EGF receptor between the negatively constrained, tethered state and the fully active back-to-back dimer conformation, to follow conformational changes in the wild-type and mutant EGF receptors after ligand binding. Our results suggest that EGFR on the cell surface can be untethered, but this form is inactive; thus, untethering of the receptor is not sufficient for activation, and ligand binding is essential for the correct positioning of the two receptor subunits to achieve kinase activation.     

Evidence for extended YFP-EGFR dimers in the absence of ligand on the surface of living cells

The epidermal growth factor receptor (EGFR) is a member of the erbB tyrosine kinase family of receptors. Structural studies have revealed two distinct conformations of the ectodomain of the EGFR: a compact, tethered, conformation and an untethered extended conformation. In the context of a monomer-dimer transition model, ligand binding is thought to untether the monomeric receptor leading to exposure of a dimerization arm which then facilitates receptor dimerization, kinase activation and signaling. For receptors directed orthogonal to the local plane of the membrane surface, this would lead to a large change in the distance of the receptor N-terminus from the membrane surface. To investigate this experimentally, we produced stable BaF/3 cell lines expressing a biochemically functional yellow fluorescent protein (YFP)-EGFR chimera and determined the vertical separation of the N-terminal YFP tag from the membrane using fluorescence resonance energy transfer (FRET) techniques. Homo-FRET/rFLIM was employed to determine the presence of unliganded dimers and to measure the average distance between the N-terminal tags in those dimers. The results suggest that EGF-induced activation occurs within or between pre-formed and extended dimers with very little change in the extension of the N-terminii from the membrane surface. These results provide constraints on possible models for EGFR activation.     

Cell-free expression and purification of the fragments of the receptor tyrosine kinases of the EGFR family,containing the transmembrane domain with the juxtamembrane region,for structural studies

The EGFR/HER receptor family of an epidermal growth factor represents an important class of the receptor tyrosine kinases playing the key role in the control of cell growth and differentiation in mammalian cells, as well as in the development of a number of pathological processes, including oncogenesis. Binding of a ligand to the extracellular domains initiates switching of the EGFR/HER receptor between the alternative dimeric states that causes the allosteric activation of kinase domains in cell cytoplasm. The transmembrane (TM) domain and adjacent flexible regions alternatively interacting with either membrane surface or kinase domains are directly involved in the complex conformational transition in EGFR/HERs. Here we report on a highly efficient system of the cell free production of the EGFR/HER TM domains with functionally important juxtamembrane (JM) regions for the investigation of the molecular basis of biochemical signal transduction across the cell membrane. To increase the efficiency of synthesis of the EGFR/HER TM-JM fragments of the receptors, we used two N-terminal expression tags, which significantly increased the protein yield. In the case of the TM-JM fragments of EGFR (residues 638–692) and HER2 (residues 644–700), the method allowed us to obtain milligram quantities of the ¹³C,¹⁵N-labeled protein for structural and biophysical investigations in the membrane-mimicking environments using high-resolution heteronuclear NMR spectroscopy.     

mAb806 binding to epidermal growth factor receptor: a computational study

The epidermal growth factor receptor (EGFR) is an important target in the treatment of cancer. A very potent antibody, mAb806, has been developed against overexpressed EGFR and was found to be particularly active in brain tumors. Structural studies reveal that it binds to an epitope on the extracellular region of the EGFR. However, this epitope is cryptic/buried in crystal structures of the active (untethered) and inactive (tethered) EGFR, and it is unclear as to how the antibody interacts with this region. To explore this interaction, we combined molecular docking, steered molecular dynamics, and equilibrium molecular dynamics simulations. Our computational models reveal that the antibody induces local unfolding around the epitope to form the antibody–EGFR complex. In addition, regions in the vicinity of the epitope also modulate the interaction, which are in accordance with several other known antibody–antigen interactions, and offers new possibilities for the design of antibodies with increased potency and specificity for this receptor. Proteins 2015; 83:153–168. © 2014 Wiley Periodicals, Inc.     

Molecular Determinants of Epidermal Growth Factor Binding: A Molecular Dynamics Study

The epidermal growth factor receptor (EGFR) is a member of the receptor tyrosine kinase family that plays a role in multiple cellular processes. Activation of EGFR requires binding of a ligand on the extracellular domain to promote conformational changes leading to dimerization and transphosphorylation of intracellular kinase domains. Seven ligands are known to bind EGFR with affinities ranging from sub-nanomolar to near micromolar dissociation constants. In the case of EGFR, distinct conformational states assumed upon binding a ligand is thought to be a determining factor in activation of a downstream signaling network. Previous biochemical studies suggest the existence of both low affinity and high affinity EGFR ligands. While these studies have identified functional effects of ligand binding, high-resolution structural data are lacking. To gain a better understanding of the molecular basis of EGFR binding affinities, we docked each EGFR ligand to the putative active state extracellular domain dimer and 25.0 ns molecular dynamics simulations were performed. MM-PBSA/GBSA are efficient computational approaches to approximate free energies of protein-protein interactions and decompose the free energy at the amino acid level. We applied these methods to the last 6.0 ns of each ligand-receptor simulation. MM-PBSA calculations were able to successfully rank all seven of the EGFR ligands based on the two affinity classes: EGF>HB-EGF>TGF-α>BTC>EPR>EPG>AR. Results from energy decomposition identified several interactions that are common among binding ligands. These findings reveal that while several residues are conserved among the EGFR ligand family, no single set of residues determines the affinity class. Instead we found heterogeneous sets of interactions that were driven primarily by electrostatic and Van der Waals forces. These results not only illustrate the complexity of EGFR dynamics but also pave the way for structure-based design of therapeutics targeting EGF ligands or the receptor itself.     

Crystal structure of a truncated epidermal growth factor receptor extracellular domain bound to transforming growth factor alpha

We report the crystal structure, at 2.5 A resolution, of a truncated human EGFR ectodomain bound to TGFalpha. TGFalpha interacts with both L1 and L2 domains of EGFR, making many main chain contacts with L1 and interacting with L2 via key conserved residues. The results indicate how EGFR family members can bind a family of highly variable ligands. In the 2:2 TGFalpha:sEGFR501 complex, each ligand interacts with only one receptor molecule. There are two types of dimers in the asymmetric unit: a head-to-head dimer involving contacts between the L1 and L2 domains and a back-to-back dimer dominated by interactions between the CR1 domains of each receptor. Based on sequence conservation, buried surface area, and mutagenesis experiments, the back-to-back dimer is favored to be biologically relevant.     

The dependence of EGFR oligomerization on environment and structure: A camera-based N&B study

Number and brightness (N&B) analysis is a fluorescence spectroscopy technique to quantify oligomerization of the mobile fraction of proteins. Accurate results, however, rely on a good knowledge of nonfluorescent states of the fluorescent labels, especially of fluorescent proteins, which are widely used in biology. Fluorescent proteins have been characterized for confocal, but not camera-based, N&B, which allows, in principle, faster measurements over larger areas. Here, we calibrate camera-based N&B implemented on a total internal reflection fluorescence microscope for various fluorescent proteins by determining their propensity to be fluorescent. We then apply camera-based N&B in live CHO-K1 cells to determine the oligomerization state of the epidermal growth factor receptor (EGFR), a transmembrane receptor tyrosine kinase that is a crucial regulator of cell proliferation and survival with implications in many cancers. EGFR oligomerization in resting cells and its regulation by the plasma membrane microenvironment are still under debate. Therefore, we investigate the effects of extrinsic factors, including membrane organization, cytoskeletal structure, and ligand stimulation, and intrinsic factors, including mutations in various EGFR domains, on the receptor’s oligomerization. Our results demonstrate that EGFR oligomerization increases with removal of cholesterol or sphingolipids or the disruption of GM3-EGFR interactions, indicating raft association. However, oligomerization is not significantly influenced by the cytoskeleton. Mutations in either I706/V948 residues or E685/E687/E690 residues in the kinase and juxtamembrane domains, respectively, lead to a decrease in oligomerization, indicating their necessity for EGFR dimerization. Finally, EGFR phosphorylation is oligomerization dependent, involving the extracellular domain (550–580 residues). Coupled with biochemical investigations, camera-based N&B indicates that EGFR oligomerization and phosphorylation are the outcomes of several molecular interactions involving the lipid content and structure of the cell membrane and multiple residues in the kinase, juxtamembrane, and extracellular domains.     

Elucidation of interactions regulating conformational stability and dynamics of SARS-CoV-2 S-protein

The ongoing COVID-19 pandemic caused by the new coronavirus, SARS-CoV-2, calls for urgent developments of vaccines and antiviral drugs. The spike protein of SARS-CoV-2 (S-protein), which consists of trimeric polypeptide chains with glycosylated residues on the surface, triggers the virus entry into a host cell. Extensive structural and functional studies on this protein have rapidly advanced our understanding of the S-protein structure at atomic resolutions, although most of these structural studies overlook the effect of glycans attached to the S-protein on the conformational stability and functional motions between the inactive down and active up forms. Here, we performed all-atom molecular dynamics simulations of both down and up forms of a fully glycosylated S-protein in solution as well as targeted molecular dynamics simulations between them to elucidate key interdomain interactions for stabilizing each form and inducing the large-scale conformational transitions. The residue-level interaction analysis of the simulation trajectories detects distinct amino acid residues and N-glycans as determinants on conformational stability of each form. During the conformational transitions between them, interdomain interactions mediated by glycosylated residues are switched to play key roles on the stabilization of another form. Electrostatic interactions, as well as hydrogen bonds between the three receptor binding domains, work as driving forces to initiate the conformational transitions toward the active form. This study sheds light on the mechanisms underlying conformational stability and functional motions of the S-protein, which are relevant for vaccine and antiviral drug developments.     

Modeling Signal Propagation Mechanisms and Ligand-Based Conformational Dynamics of the Hsp90 Molecular Chaperone Full-Length Dimer

Hsp90 is a molecular chaperone essential for protein folding and activation in normal homeostasis and stress response. ATP binding and hydrolysis facilitate Hsp90 conformational changes required for client activation. Hsp90 plays an important role in disease states, particularly in cancer, where chaperoning of the mutated and overexpressed oncoproteins is important for function. Recent studies have illuminated mechanisms related to the chaperone function. However, an atomic resolution view of Hsp90 conformational dynamics, determined by the presence of different binding partners, is critical to define communication pathways between remote residues in different domains intimately affecting the chaperone cycle. Here, we present a computational analysis of signal propagation and long-range communication pathways in Hsp90. We carried out molecular dynamics simulations of the full-length Hsp90 dimer, combined with essential dynamics, correlation analysis, and a signal propagation model. All-atom MD simulations with timescales of 70 ns have been performed for complexes with the natural substrates ATP and ADP and for the unliganded dimer. We elucidate the mechanisms of signal propagation and determine “hot spots” involved in interdomain communication pathways from the nucleotide-binding site to the C-terminal domain interface. A comprehensive computational analysis of the Hsp90 communication pathways and dynamics at atomic resolution has revealed the role of the nucleotide in effecting conformational changes, elucidating the mechanisms of signal propagation. Functionally important residues and secondary structure elements emerge as effective mediators of communication between the nucleotide-binding site and the C-terminal interface. Furthermore, we show that specific interdomain signal propagation pathways may be activated as a function of the ligand. Our results support a “conformational selection model” of the Hsp90 mechanism, whereby the protein may exist in a dynamic equilibrium between different conformational states available on the energy landscape and binding of a specific partner can bias the equilibrium toward functionally relevant complexes.     

Identification of a determinant of epidermal growth factor receptor ligand-binding specificity using a truncated, high-affinity form of the ectodomain

Murine and human epidermal growth factor receptors (EGFRs) bind human EGF (hEGF), mouse EGF (mEGF), and human transforming growth factor alpha (hTGF-alpha) with high affinity despite the significant differences in the amino acid sequences of the ligands and the receptors. In contrast, the chicken EGFR can discriminate between mEGF (and hEGF) and hTGF-alpha and binds the EGFs with approximately 100-fold lower affinity. The regions responsible for this poor binding are known to be Arg(45) in hEGF and the L2 domain in the chicken EGFR. In this study we have produced a truncated form of the hEGFR ectodomain comprising residues 1-501 (sEGFR501), which, unlike the full-length hEGFR ectodomain (residues 1-621, sEGFR621), binds hEGF and hTGF-alpha with high affinity (K(D) = 13-21 and 35-40 nM, respectively). sEGFR501 was a competitive inhibitor of EGF-stimulated mitogenesis, being almost 10-fold more effective than the full-length EGFR ectodomain and three times more potent than the neutralizing anti-EGFR monoclonal antibody Mab528. Analytical ultracentrifugation showed that the primary EGF binding sites on sEGFR501 were saturated at an equimolar ratio of ligand and receptor, leading to the formation of a 2:2 EGF:sEGFR501 dimer complex. We have used sEGFR501 to generate three mutants with single position substitutions at Glu(367), Gly(441), or Glu(472) to Lys, the residue found in the corresponding positions in the chicken EGFR. All three mutants bound hTGF-alpha and were recognized by Mab528. However, mutant Gly(441)Lys showed markedly reduced binding to hEGF, implicating Gly(441), in the L2 domain, as part of the binding site that recognizes Arg(45) of hEGF.     

Ligand-induced changes in membrane-bound acetylcholine receptor observed by ethidium fluorescence. 3. Stopped-flow studies with histrionicotoxin.

Rapid kinetic studies of histrionicotoxin interactions with membrane-bound acetylcholine-receptor showed a conformational change in the receptor-histironicotoxin complex as reflected by a decrease in fluorescence intensity of the extrinsic probe ethidium. The simplest kinetic mechanism consistent with the observed data is one in which a rapid preequiliibrium exists between receptor and toxin (K = 3.33 micrometers), followed by a slow conformational change (k1 congruent to 2 X 10(-2) s-1 and k-1 congruent to 1.5 X 10(-3) s-1). The overall equilibrium constant (Kov) determined from a fit of the amplitude dependence on toxin concentration had a value of 0.25 micrometer. The data preclude kinetic mechanisms where histrionicotoxin acts as an effector, shifting equilibria between preexisting, discrete, and slowly interconverting receptor forms.     

Escherichia coli cAMP receptor protein: evidence for three protein conformational states with different promoter binding affinities

Cyclic AMP receptor protein (CRP) from Escherichia coli is assumed to exist in two states, namely, those represented by the free protein and that of the ligand-protein complex. To establish a quantitative structure-function relation between cAMP binding and the cAMP-induced conformational changes in the receptor, protein conformational change was quantitated as a function of cAMP concentration up to 10 mM. The protein conformation was monitored by four different methods at pH 7.8 and 23 degrees C, namely, rate of proteolytic digestion by subtilisin, rate of chemical modification of Cys-178, tryptophan fluorescence, and fluorescence of the extrinsic fluorescence probe 8-anilino-1-naphthalenesulfonic acid (ANS). Each of these techniques reveals a biphasic dependence of protein conformation on cAMP concentration. At low cAMP concentrations ranging from 0 to 200 microM, the rates of proteolytic digestion and that of Cys-178 modification increase, whereas the fluorescence intensity of the ANS-protein complex is quenched, and there is no change in the fluorescence intensity of the tryptophan residues in the protein. At higher cAMP concentrations, the rates of proteolytic and chemical modification of the protein decrease, while the fluorescence intensity of the ANS-protein complex is further quenched but there is an increase in the intensity of tryptophan fluorescence. These results show unequivocally that there are at least three conformational states of the protein. The association constants for the formation of CRP-cAMP and CRP-(cAMP)2 complexes derived from conformational studies are in good agreement with those determined by equilibrium dialysis, nonequilibrium dialysis, and ultrafiltration. Therefore, the simplest explanation would be that the protein exhibits three conformational states, free CRP and two cAMP-dependent states, which correspond to the CRP-cAMP and CRP-(cAMP)2 complexes. The binding properties of CRP-cAMP and CRP-(cAMP)2 to the lac promoter were studied by using the gel retardation technique. At a high concentration of cAMP which favors the formation of the CRP-(cAMP)2 complex, binding of the protein to DNA is decreased. This, together with conformational data, strongly suggests that only the CRP-cAMP complex is active in specific DNA binding whereas CRP and CRP-(cAMP)2 are not.     

Crystal structures of the maltodextrin/maltose-binding protein complexed with reduced oligosaccharides: flexibility of tertiary structure and ligand binding

The structure of the maltodextrin or maltose-binding protein, an initial receptor for bacterial ABC-type active transport and chemotaxis, consists of two globular domains that are separated by a groove wherein the ligand is bound and enclosed by an inter-domain rotation. Here, we report the determination of the crystal structures of the protein complexed with reduced maltooligosaccharides (maltotriitol and maltotetraitol) in both the "closed" and "open" forms. Although these modified sugars bind to the receptor, they are not transported by the wild-type transporter. In the closed structures, the reduced sugars are buried in the groove and bound by both domains, one domain mainly by hydrogen-bonding interactions and the other domain primarily by non-polar interactions with aromatic side-chains. In the open structures, which abrogate both cellular activities of active transport and chemotaxis because of the large separation between the two domains, the sugars are bound almost exclusively to the domain rich in aromatic residues. The binding site for the open chain glucitol residue extends to a subsite that is distinct from those for the glucose residues that were uncovered in prior structural studies of the binding of active linear maltooligosaccharides. Occupation of this subsite may also account for the inability of the reduced oligosaccharides to be transported. The structures reported here, combined with those previously determined for several other complexes with active oligosaccharides in the closed form and with cyclodextrin in the open form, revealed at least four distinct modes of ligand binding but with only one being functionally active. This versatility reflects the flexibility of the protein, from very large motions of interdomain rotation to more localized side-chain conformational changes, and adaptation by the oligosaccharides as well.     

Activation of Zap-70 tyrosine kinase due to a structural rearrangement induced by tyrosine phosphorylation and/or ITAM binding

The protein tyrosine kinase ZAP-70 is implicated in the early steps of the T-cell antigen receptor (TCR) signaling. Binding of ZAP-70 to the phosphorylated immunoreceptor tyrosine-based activation motifs (ITAMs) of the TCR zeta chain through its two src-homology 2 (SH2) domains results in its activation coupled to phosphorylation on multiple tyrosine residues, mediated by Src kinases including Lck as well as by autophosphorylation. The mechanism of ZAP-70 activation following receptor binding is still not completely understood. Here we investigated the effect of intramolecular interactions and autophosphorylation by following the kinetics of recombinant ZAP-70 activation in a spectrophotometric substrate phosphorylation assay. Under these conditions, we observed a lag phase of several minutes before full ZAP-70 activation, which was not observed using a truncated form lacking the first 254 residues, suggesting that it might be due to an intramolecular interaction involving the interdomain A and SH2 region. Accordingly, the lag phase could be reproduced by testing the truncated form in the presence of recombinant SH2 domains and was abolished by the addition of diphosphorylated ITAM peptide. Preincubation with ATP or phosphorylation by Lck also abolished the lag phase and resulted in a more active enzyme. The same results were obtained using a ZAP-70 mutant lacking the interdomain B tyrosines. These findings are consistent with a mechanism in which ZAP-70 phosphorylation/autophosphorylation on tyrosine(s) other than 292, 315, and 319, as well as engagement of the SH2 domains by the phosphorylated TCR, can induce a conformational change leading to accelerated enzyme kinetics and higher catalytic efficiency.     

Determination of the Individual Roles of the Linker Residues in the Interdomain Motions of Calmodulin Using NMR Chemical Shifts

Many protein molecules are formed by two or more domains whose structures and dynamics are closely related to their biological functions. It is thus important to develop methods to determine the structural properties of these multidomain proteins. Here, we characterize the interdomain motions in the calcium-bound state of calmodulin (Ca^2 +-CaM) using NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations. We find that the conformational fluctuations of the interdomain linker, which are largely responsible for the overall interdomain motions of CaM, can be well described by exploiting the information provided by chemical shifts. We thus identify 10 residues in the interdomain linker region that change their conformations upon substrate binding. Five of these residues (Met76, Lys77, Thr79, Asp80 and Ser81) are highly flexible and cover the range of conformations observed in the substrate-bound state, while the remaining five (Arg74, Lys75, Asp78, Glu82 and Glu83) are much more rigid and do not populate conformations typical of the substrate-bound form. The ensemble of conformations representing the Ca^2 +-CaM state obtained in this study is in good agreement with residual dipolar coupling, paramagnetic resonance enhancement, small-angle X-ray scattering and fluorescence resonance energy transfer measurements, which were not used as restraints in the calculations. These results provide initial evidence that chemical shifts can be used to characterize the conformational fluctuations of multidomain proteins.     

Conformation of a peptide ligand bound to its G-protein coupled receptor

Many peptide hormones elicit a wide array of physiological effects by binding to G-protein coupled receptors. We have determined the conformation of pituitary adenylate cyclase activating polypeptide, PACAP(1--21)NH(2), bound to a PACAP-specific receptor by NMR spectroscopy. Residues 3--7 form a unique beta-coil structure that is preceded by an N-terminal extended tail. This beta-coil creates a patch of hydrophobic residues that is important for receptor binding. In contrast, the C-terminal region (residues 8--21) forms an alpha-helix, similar to that in the micelle-bound PACAP. Thus, the conformational difference between PACAP in the receptor-bound and the micelle-bound states is limited to the N-terminal seven residues. This observation is consistent with the two-step ligand transportation model in which PACAP first binds to the membrane nonspecifically and then diffuses two-dimensionally in search of its receptor; a conformational change at the N-terminal region then allows specific interactions between the ligand and the receptor.     

Role of the Src homology 2 domains and interdomain regions in ZAP-70 phosphorylation and enzymatic activity.

The protein tyrosine kinase ZAP-70, which mediates T-cell antigen receptor (TCR) signalling, contains three distinct functional modules, two tandemly arranged SH2 domains, a kinase domain and a linker region (interdomain B) that connects them. ZAP-70 enzymatic activation is strictly dependent on the binding, via its SH2 domains, to the triggered TCR and on tyrosine phosphorylation. Here we utilized recombinant ZAP-70 and carried out a mutational analysis to understand the structural requirements for its activation. We show that deletion of both SH2 domains corresponding to the first 254 residues moderately increases ZAP-70 enzymatic activity on an exogenous substrate in vitro, results in increased tyrosine phosphorylation and produces subtle conformational changes, as judged by altered SDS/PAGE migration. Mutation of Tyr292, 315 and 319 to Phe in the interdomain B region, which constitute the major phosphorylation sites both in vitro and in vivo, did not affect ZAP-70 enzymatic activity. Moreover, deletion analysis of the interdomain B region established residues 320-619 as a minimal region endowed with full kinase activity. We propose that binding of ZAP-70 to the TCR promotes, through conformational changes, its extensive phosphorylation on tyrosine. However, Tyr292, 315 and 319 do not affect ZAP-70 enzymatic activity and may influence ZAP-70 signalling only indirectly by mediating its association with intracellular transducers.     

Structural model of the mAb 806-EGFR complex using computational docking followed by computational and experimental mutagenesis

In this work, we combined computational protein-protein docking with computational and experimental mutagenesis to predict the structure of the complex formed by monoclonal antibody 806 (mAb 806) and the epidermal growth factor receptor (EGFR). We docked mAb 806, an antitumor antibody, to its epitope of EGFR residues 287-302. Potential mAb 806-EGFR orientations were generated, and computational mutagenesis was used to filter them according to their agreement with experimental mutagenesis data. Further computational mutagenesis suggested additional mutations, which were tested to arrive at a final structure that was most consistent with experimental mutagenesis data. We propose that this is the EGFR-mAb 806 structure, in which mAb 806 binds to an untethered form of the receptor, consistent with published experimental results. The steric hindrance created by the antibody near the EGFR dimer interface interferes with receptor dimerization, and we postulate this as the structural origin for the antitumor effect of mAb 806.     

Stability of hairpin ribozyme tertiary structure is governed by the interdomain junction.

The equilibrium distributions of hairpin ribozyme conformational isomers have been examined by time-resolved fluorescence resonance energy transfer. Ribozymes partition between active (docked) and inactive (extended) conformers, characterized by unique interdomain distance distributions, which define differences in folding free energy. The active tertiary structure is stabilized both by specific interactions between the catalytic and the substrate-binding domains and by the structure of the intervening helical junction. Under physiological conditions, the docking equilibrium of the natural four-way junction dramatically favors the active conformer, while those of a three-way and the two-way junction used in gene therapy applications favor the inactive conformer.