Use of intrinsic binding energy for catalysis by an RNA enzyme

The contribution of several individual ribozyme⋅substrate base pairs to binding and catalysis has been investigated using hammerhead ribozyme substrates that were truncated at their 3′ or 5′ ends. The base pairs at positions 1.1–2.1 and 15.2–16.2, which flank the conserved core, each contribute 104-fold in the chemical step, without affecting substrate binding. In contrast, base pairs distal to the core contribute to substrate binding but have no effect on the chemical step. These results suggest a “fraying model” in which each ribozyme⋅substrate helix can exist in either an unpaired (“open”) state or a helical (“closed”) state, with the closed state required for catalysis. The base pairs directly adjacent to the conserved core contribute to catalysis by allowing the closed state to form. Once the number of base pairs is sufficient to ensure that the closed helical state predominates, additional residues provide stabilization of the helix, and therefore increase binding, but have no further effect on the chemical step. Remarkably, the >5 kcal/mol free energy contribution to catalysis from each of the internal base pairs is considerably greater than the free energy expected for formation of a base pair. It is suggested that this unusually large energetic contribution arises because free energy that is typically lost in constraining residues within a base pair is expressed in the transition state, where it is used for positioning. This extends the concept of “intrinsic binding energy” from protein to RNA enzymes, suggesting that intrinsic binding energy is a fundamental feature of biological catalysis.

Use of binding energy by an RNA enzyme for catalysis by positioning and substrate destabilization.

A fundamental catalytic principle for protein enzymes in the use of binding interactions away from the site of chemical transformation for catalysis. We have compared the binding and reactivity of a series of oligonucleotide substrates and products of the Tetrahymena ribozyme, which catalyzes a site-specific phosphodiester cleavage reaction: CCCUCUpA+G<-->CCCUCU-OH+GpA. The results suggest that this RNA enzyme, like protein enzymes, can utilize binding interactions to achieve substantial catalysis via entropic fixation and substrate destabilization. The stronger binding of the all-ribose oligonucleotide product compared to an analog with a terminal 3' deoxyribose residue gives an effective concentration of 2200 M for the 3' hydroxyl group, a value approaching those obtained with protein enzymes and suggesting the presence of a structurally well defined active site capable of precise positioning. The stabilization from tertiary binding interactions is 40-fold less for the oligonucleotide substrate than the oligonucleotide product, despite the presence of the reactive phosphoryl group in the substrate. This destabilization is accounted for by a model in which tertiary interactions away from the site of bond cleavage position the electron-deficient 3' bridging phosphoryl oxygen of the oligonucleotide substrate next to an electropositive Mg ion. As the phosphodiester bond breaks and this 3' oxygen atom develops a negative charge in the transition state, the weak interaction of the substrate with Mg2+ becomes strong. These strategies of "substrate destabilization" and "transition state stabilization" provide estimated rate enhancements of approximately 280- and approximately 60-fold, respectively. Analogous substrate destabilization by a metal ion or hydrogen bond donor may be used more generally by RNA and protein enzymes catalyzing reactions of phosphate esters.

The change in hydrogen bond strength accompanying charge rearrangement: Implications for enzymatic catalysis

The equilibrium for formation of the intramolecular hydrogen bond (KHB) in a series of substituted salicylate monoanions was investigated as a function of ΔpKa, the difference between the pKa values of the hydrogen bond donor and acceptor, in both water and dimethyl sulfoxide. The dependence of log KHB upon ΔpKa is linear in both solvents, but is steeper in dimethyl sulfoxide (slope = 0.73) than in water (slope = 0.05). Thus, hydrogen bond strength can undergo substantially larger increases in nonaqueous media than aqueous solutions as the charge density on the donor or acceptor atom increases. These results support a general mechanism for enzymatic catalysis, in which hydrogen bonding to a substrate is strengthened as charge rearranges in going from the ground state to the transition state; the strengthening of the hydrogen bond would be greater in a nonaqueous enzymatic active site than in water, thus providing a rate enhancement for an enzymatic reaction relative to the solution reaction. We suggest that binding energy of an enzyme is used to fix the substrate in the low-dielectric active site, where the strengthening of the hydrogen bond in the course of a reaction is increased.

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding.

We propose a general mean field model of ligand-protein interactions to determine the thermodynamic equilibrium of a system at finite temperature. The method is employed in structural assessments of two human immuno-deficiency virus type 1 protease complexes where the gross effects of protein flexibility are incorporated by utilizing a data base of crystal structures. Analysis of the energy spectra for these complexes has revealed that structural and thermo-dynamic aspects of molecular recognition can be rationalized on the basis of the extent of frustration in the binding energy landscape. In particular, the relationship between receptor-specific binding of these ligands to human immunodeficiency virus type 1 protease and a minimal frustration principle is analyzed.

Promotion of purine nucleotide binding to thymidylate synthase by a potent folate analogue inhibitor, 1843U89.

A folate analogue, 1843U89 (U89), with potential as a chemotherapeutic agent due to its potent and specific inhibition of thymidylate synthase (TS; EC 2.1.1.45), greatly enhances not only the binding of 5-fluoro-2'-deoxyuridine 5'-monophosphate (FdUMP) and dUMP to Escherichia coli TS but also that of dGMP, GMP, dIMP, and IMP. Guanine nucleotide binding was first detected by CD analysis, which revealed a unique spectrum for the TS-dGMP-U89 ternary complex. The quantitative binding of dGMP relative to GMP, FdUMP, and dUMP was determined in the presence and absence of U89 by ultrafiltration analysis, which revealed that although the binding of GMP and dGMP could not be detected in the absence of U89 both were bound in its presence. The Kd for dGMP was about the same as that for dUMP and FdUMP, with binding of the latter two nucleotides being increased by two orders of magnitude by U89. An explanation for the binding of dGMP was provided by x-ray diffraction studies that revealed an extensive stacking interaction between the guanine of dGMP and the benzoquinazoline ring of U89 and hydrogen bonds similar to those involved in dUMP binding. In addition, binding energy was provided through a water molecule that formed hydrogen bonds to both N7 of dGMP and the hydroxyl of Tyr-94. Accommodation of the larger dGMP molecule was accomplished through a distortion of the active site and a shift of the deoxyribose moiety to a new position. These rearrangements also enabled the binding of GMP to occur by creating a pocket for the ribose 2' hydroxyl group, overcoming the normal TS discrimination against nucleotides containing the 2' hydroxyl.

Construction of a high-affinity receptor site for dihydropyridine agonists and antagonists by single amino acid substitutions in a non-l-type Ca2+ channel

The activity of l-type Ca2+ channels is increased by dihydropyridine (DHP) agonists and inhibited by DHP antagonists, which are widely used in the therapy of cardiovascular disease. These drugs bind to the pore-forming α1 subunits of l-type Ca2+ channels. To define the minimal requirements for DHP binding and action, we constructed a high-affinity DHP receptor site by substituting a total of nine amino acid residues from DHP-sensitive l-type α1 subunits into the S5 and S6 transmembrane segments of domain III and the S6 transmembrane segment of domain IV of the DHP-insensitive P/Q-type α1A subunit. The resulting chimeric α1A/DHPS subunit bound DHP antagonists with high affinity in radioligand binding assays and was inhibited by DHP antagonists with high affinity in voltage clamp experiments. Substitution of these nine amino acid residues yielded 86% of the binding energy of the l-type α1C subunit and 92% of the binding energy of the l-type α1S subunit for the high-affinity DHP antagonist PN200–110. The activity of chimeric Ca2+ channels containing α1A/DHPS was increased 3.5 ± 0.7-fold by the DHP agonist (−)Bay K8644. The effect of this agonist was stereoselective as in l-type Ca2+ channels since (+) Bay K8644 inhibited the activity of α1A/DHPS. The results show conclusively that DHP agonists and antagonists bind to a single receptor site at which they have opposite effects on Ca2+ channel activity. This site contains essential components from both domains III and IV, consistent with a domain interface model for binding and allosteric modulation of Ca2+ channel activity by DHPs.

Uracil-DNA glycosylase–DNA substrate and product structures: Conformational strain promotes catalytic efficiency by coupled stereoelectronic effects

Enzymatic transformations of macromolecular substrates such as DNA repair enzyme/DNA transformations are commonly interpreted primarily by active-site functional-group chemistry that ignores their extensive interfaces. Yet human uracil–DNA glycosylase (UDG), an archetypical enzyme that initiates DNA base-excision repair, efficiently excises the damaged base uracil resulting from cytosine deamination even when active-site functional groups are deleted by mutagenesis. The 1.8-Å resolution substrate analogue and 2.0-Å resolution cleaved product cocrystal structures of UDG bound to double-stranded DNA suggest enzyme–DNA substrate-binding energy from the macromolecular interface is funneled into catalytic power at the active site. The architecturally stabilized closing of UDG enforces distortions of the uracil and deoxyribose in the flipped-out nucleotide substrate that are relieved by glycosylic bond cleavage in the product complex. This experimentally defined substrate stereochemistry implies the enzyme alters the orientation of three orthogonal electron orbitals to favor electron transpositions for glycosylic bond cleavage. By revealing the coupling of this anomeric effect to a delocalization of the glycosylic bond electrons into the uracil aromatic system, this structurally implicated mechanism resolves apparent paradoxes concerning the transpositions of electrons among orthogonal orbitals and the retention of catalytic efficiency despite mutational removal of active-site functional groups. These UDG/DNA structures and their implied dissociative excision chemistry suggest biology favors a chemistry for base-excision repair initiation that optimizes pathway coordination by product binding to avoid the release of cytotoxic and mutagenic intermediates. Similar excision chemistry may apply to other biological reaction pathways requiring the coordination of complex multistep chemical transformations.

A host–guest system to study structure–function relationships of membrane fusion peptides

We designed a host–guest fusion peptide system, which is completely soluble in water and has a high affinity for biological and lipid model membranes. The guest sequences are those of the fusion peptides of influenza hemagglutinin, which are solubilized by a highly charged unstructured C-terminal host sequence. These peptides partition to the surface of negatively charged liposomes or erythrocytes and elicit membrane fusion or hemolysis. They undergo a conformational change from random coil to an obliquely inserted (≈33° from the surface) α-helix on binding to model membranes. Partition coefficients for membrane insertion were measured for influenza fusion peptides of increasing lengths (n = 8, 13, 16, and 20). The hydrophobic contribution to the free energy of binding of the 20-residue fusion peptide at pH 5.0 is −7.6 kcal/mol (1 cal = 4.18 J). This energy is sufficient to stabilize a “stalk” intermediate if a typical number of fusion peptides assemble at the site of membrane fusion. The fusion activity of the fusion peptides increases with each increment in length, and this increase strictly correlates with the hydrophobic binding energy and the angle of insertion.

Energy transfer to a proton-transfer fluorescence probe: Tryptophan to a flavonol in human serum albumin

A protein fluorescence probe system, coupling excited-state intermolecular Förster energy transfer and intramolecular proton transfer (PT), is presented. As an energy donor for this system, we used tryptophan, which transfers its excitation energy to 3-hydroxyflavone (3-HF) as a flavonol prototype, an acceptor exhibiting excited-state intramolecular PT. We demonstrate such a coupling in human serum albumin–3-HF complexes, excited via the single intrinsic tryptophan (Trp-214). Besides the PT tautomer fluorescence (λmax = 526 nm), these protein–probe complexes exhibit a 3-HF anion emission (λmax = 500 nm). Analysis of spectroscopic data leads to the conclusion that two binding sites are involved in the human serum albumin–3-HF interaction. The 3-HF molecule bound in the higher affinity binding site, located in the IIIA subdomain, has the association constant (k1) of 7.2 × 105 M−1 and predominantly exists as an anion. The lower affinity site (k2 = 2.5 × 105 M−1), situated in the IIA subdomain, is occupied by the neutral form of 3-HF (normal tautomer). Since Trp-214 is situated in the immediate vicinity of the 3-HF normal tautomer bound in the IIA subdomain, the intermolecular energy transfer for this donor/acceptor pair has a 100% efficiency and is followed by the PT tautomer fluorescence. Intermolecular energy transfer from the Trp-214 to the 3-HF anion bound in the IIIA subdomain is less efficient and has the rate of 1.61 × 108 s−1, thus giving for the donor/acceptor distance a value of 25.5 Å.

Preassembly of interleukin 2 (IL-2) receptor subunits on resting Kit 225 K6 T cells and their modulation by IL-2, IL-7, and IL-15: A fluorescence resonance energy transfer study

Assembly and mutual proximities of α, β, and γc subunits of the interleukin 2 receptors (IL-2R) in plasma membranes of Kit 225 K6 T lymphoma cells were investigated by fluorescence resonance energy transfer (FRET) using fluorescein isothiocyanate- and Cy3-conjugated monoclonal antibodies (mAbs) that were directed against the IL-2Rα, IL-2Rβ, and γc subunits of IL-2R. The cell-surface distribution of subunits was analyzed at the nanometer scale (2–10 nm) by FRET on a cell-by-cell basis. The cells were probed in resting phase and after coculture with saturating concentrations of IL-2, IL-7, and IL-15. FRET data from donor- and acceptor-labeled IL-2Rβ-α, γ-α, and γ-β pairs demonstrated close proximity of all subunits to each other in the plasma membrane of resting T cells. These mutual proximities do not appear to represent mAb-induced microaggregation, because FRET measurements with Fab fragments of the mAbs gave similar results. The relative proximities were meaningfully modulated by binding of IL-2, IL-7, and IL-15. Based on FRET analysis the topology of the three subunits at the surface of resting cells can be best described by a “triangular model” in the absence of added interleukins. IL-2 strengthens the bridges between the subunits, making the triangle more compact. IL-7 and IL-15 act in the opposite direction by opening the triangle possibly because they associate their private specific α receptors with the β and/or γc subunits of the IL-2R complex. These data suggest that IL-2R subunits are already colocalized in resting T cells and do not require cytokine-induced redistribution. This colocalization is significantly modulated by binding of relevant interleukins in a cytokine-specific manner.

Identification of the proton pathway in bacterial reaction centers: Inhibition of proton transfer by binding of Zn2+ or Cd2+

The reaction center (RC) from Rhodobacter sphaeroides converts light into chemical energy through the light induced two-electron, two-proton reduction of a bound quinone molecule QB (the secondary quinone acceptor). A unique pathway for proton transfer to the QB site had so far not been determined. To study the molecular basis for proton transfer, we investigated the effects of exogenous metal ion binding on the kinetics of the proton-assisted electron transfer kAB(2) (QA−•QB−• + H+ → QA(QBH)−, where QA is the primary quinone acceptor). Zn2+ and Cd2+ bound stoichiometrically to the RC (KD ≤ 0.5 μM) and reduced the observed value of kAB(2) 10-fold and 20-fold (pH 8.0), respectively. The bound metal changed the mechanism of the kAB(2) reaction. In native RCs, kAB(2) was previously shown to be rate-limited by electron transfer based on the dependence of kAB(2) on the driving force for electron transfer. Upon addition of Zn2+ or Cd2+, kAB(2) became approximately independent of the electron driving force, implying that the rate of proton transfer was reduced (≥ 102-fold) and has become the rate-limiting step. The lack of an effect of the metal binding on the charge recombination reaction D+•QAQB−• → DQAQB suggests that the binding site is located far (>10 Å) from QB. This hypothesis is confirmed by preliminary x-ray structure analysis. The large change in the rate of proton transfer caused by the stoichiometric binding of the metal ion shows that there is one dominant site of proton entry into the RC from which proton transfer to QB−• occurs.

Probing the interaction between two single molecules: fluorescence resonance energy transfer between a single donor and a single acceptor.

We extend the sensitivity of fluorescence resonance energy transfer (FRET) to the single molecule level by measuring energy transfer between a single donor fluorophore and a single acceptor fluorophore. Near-field scanning optical microscopy (NSOM) is used to obtain simultaneous dual color images and emission spectra from donor and acceptor fluorophores linked by a short DNA molecule. Photodestruction dynamics of the donor or acceptor are used to determine the presence and efficiency of energy transfer. The classical equations used to measure energy transfer on ensembles of fluorophores are modified for single-molecule measurements. In contrast to ensemble measurements, dynamic events on a molecular scale are observable in single pair FRET measurements because they are not canceled out by random averaging. Monitoring conformational changes, such as rotations and distance changes on a nanometer scale, within single biological macromolecules, may be possible with single pair FRET.

Single-pair fluorescence resonance energy transfer on freely diffusing molecules: Observation of Förster distance dependence and subpopulations

Photon bursts from single diffusing donor-acceptor labeled macromolecules were used to measure intramolecular distances and identify subpopulations of freely diffusing macromolecules in a heterogeneous ensemble. By using DNA as a rigid spacer, a series of constructs with varying intramolecular donor-acceptor spacings were used to measure the mean and distribution width of fluorescence resonance energy transfer (FRET) efficiencies as a function of distance. The mean single-pair FRET efficiencies qualitatively follow the distance dependence predicted by Förster theory. Possible contributions to the widths of the FRET efficiency distributions are discussed, and potential applications in the study of biopolymer conformational dynamics are suggested. The ability to measure intramolecular (and intermolecular) distances for single molecules implies the ability to distinguish and monitor subpopulations of molecules in a mixture with different distances or conformational states. This is demonstrated by monitoring substrate and product subpopulations before and after a restriction endonuclease cleavage reaction. Distance measurements at single-molecule resolution also should facilitate the study of complex reactions such as biopolymer folding. To this end, the denaturation of a DNA hairpin was examined by using single-pair FRET.