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.

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.

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.

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.

Antigen binding forces of individually addressed single-chain Fv antibody molecules

Antibody single-chain Fv fragment (scFv) molecules that are specific for fluorescein have been engineered with a C-terminal cysteine for a directed immobilization on a flat gold surface. Individual scFv molecules can be identified by atomic force microscopy. For selected molecules the antigen binding forces are then determined by using a tip modified with covalently immobilized antigen. An scFv mutant of 12% lower free energy for ligand binding exhibits a statistically significant 20% lower binding force. This strategy of covalent immobilization and measuring well separated single molecules allows the characterization of ligand binding forces in molecular repertoires at the single molecule level and will provide a deeper insight into biorecognition processes.

Structural changes in hemoglobin during adsorption to solid surfaces: Effects of pH, ionic strength, and ligand binding

We have studied the adsorption of two structurally similar forms of hemoglobin (met-Hb and HbCO) to a hydrophobic self-assembled methyl-terminated thiol monolayer on a gold surface, by using a Quartz Crystal Microbalance (QCM) technique. This technique allows time-resolved simultaneous measurements of changes in frequency (f) (c.f. mass) and energy dissipation (D) (c.f. rigidity/viscoelastic properties) of the QCM during the adsorption process, which makes it possible to investigate the viscoelastic properties of the different protein layers during the adsorption process. Below the isoelectric points of both met-Hb and HbCO, the ΔD vs. Δf graphs displayed two phases with significantly different slopes, which indicates two states of the adsorbed proteins with different visco-elastic properties. The slope of the first phase was smaller than that of the second phase, which indicates that the first phase was associated with binding of a more rigidly attached, presumably denatured protein layer, whereas the second phase was associated with formation of a second layer of more loosely bound proteins. This second layer desorbed, e.g., upon reduction of Fe3+ of adsorbed met-Hb and subsequent binding of carbon monoxide (CO) forming HbCO. Thus, the results suggest that the adsorbed proteins in the second layer were in a native-like state. This information could only be obtained from simultaneous, time-resolved measurements of changes in both D and f, demonstrating that the QCM technique provides unique information about the mechanisms of protein adsorption to solid surfaces.

The Aer protein and the serine chemoreceptor Tsr independently sense intracellular energy levels and transduce oxygen, redox, and energy signals for Escherichia coli behavior

We identified a protein, Aer, as a signal transducer that senses intracellular energy levels rather than the external environment and that transduces signals for aerotaxis (taxis to oxygen) and other energy-dependent behavioral responses in Escherichia coli. Domains in Aer are similar to the signaling domain in chemotaxis receptors and the putative oxygen-sensing domain of some transcriptional activators. A putative FAD-binding site in the N-terminal domain of Aer shares a consensus sequence with the NifL, Bat, and Wc-1 signal-transducing proteins that regulate gene expression in response to redox changes, oxygen, and blue light, respectively. A double mutant deficient in aer and tsr, which codes for the serine chemoreceptor, was negative for aerotaxis, redox taxis, and glycerol taxis, each of which requires the proton motive force and/or electron transport system for signaling. We propose that Aer and Tsr sense the proton motive force or cellular redox state and thereby integrate diverse signals that guide E. coli to environments where maximal energy is available for growth.

“Switch I” mutant forms of the bacterial enhancer-binding protein NtrC that perturb the response to DNA

NtrC (nitrogen regulatory protein C) is a bacterial enhancer-binding protein of 469 residues that activates transcription by σ54-holoenzyme. A region of its transcriptional activation (central) domain that is highly conserved among homologous activators of σ54-holoenzyme—residues 206–220—is essential for interaction with this RNA polymerase: it is required for contact with the polymerase and/or for coupling the energy from ATP hydrolysis to a change in the conformation of the polymerase that allows it to form transcriptionally productive open complexes. Several mutant NtrC proteins with amino acid substitutions in this region, including NtrCA216V and NtrCG219K, have normal ATPase activity but fail in transcriptional activation. We now report that other mutant forms carrying amino acid substitutions at these same positions, NtrCA216C and NtrCG219C, are capable of activating transcription when they are not bound to a DNA template (non-DNA-binding derivatives with an altered helix–turn–helix DNA-binding motif at the C terminus of the protein) but are unable to do so when they are bound to a DNA template, whether or not it carries a specific enhancer. Enhancer DNA remains a positive allosteric effector of ATP hydrolysis, as it is for wild-type NtrC but, surprisingly, appears to have become a negative allosteric effector for some aspect of interaction with σ54-holoenzyme. The conserved region in which these amino acid substitutions occur (206–220) is equivalent to the Switch I region of a large group of purine nucleotide-binding proteins. Interesting analogies can be drawn between the Switch I region of NtrC and that of p21ras.

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 Å.