Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.     

Structural basis for the aldolase and epimerase activities of Staphylococcus aureus dihydroneopterin aldolase

Dihydroneopterin aldolase (DHNA) catalyzes the conversion of 7,8-dihydroneopterin (DHNP) to 6-hydroxymethyl-7,8-dihydropterin (HP) and the epimerization of DHNP to 7,8-dihydromonopterin (DHMP). Although crystal structures of the enzyme from several microorganisms have been reported, no structural information is available about the critical interactions between DHNA and the trihydroxypropyl moiety of the substrate, which undergoes bond cleavage and formation. Here, we present the structures of Staphylococcus aureus DHNA (SaDHNA) in complex with neopterin (NP, an analog of DHNP) and with monapterin (MP, an analog of DHMP), filling the gap in the structural analysis of the enzyme. In combination with previously reported SaDHNA structures in its ligand-free form (PDB entry 1DHN) and in complex with HP (PDB entry 2DHN), four snapshots for the catalytic center assembly along the reaction pathway can be derived, advancing our knowledge about the molecular mechanism of SaDHNA-catalyzed reactions. An additional step appears to be necessary for the epimerization of DHMP to DHNP. Three active site residues (E22, K100, and Y54) function coordinately during catalysis: together, they organize the catalytic center assembly, and individually, each plays a central role at different stages of the catalytic cycle.     

Crystal structures of Trypanosoma brucei and Staphylococcus aureus mevalonate diphosphate decarboxylase inform on the determinants of specificity and reactivity

Mevalonate diphosphate decarboxylase (MDD) catalyzes the ATP-dependent decarboxylation of mevalonate 5-diphosphate (MDP) to form isopentenyl pyrophosphate, a ubiquitous precursor for isoprenoid biosynthesis. MDD is a poorly understood component of this important metabolic pathway. Complementation of a temperature-sensitive yeast mutant by the putative mdd genes of Trypanosoma brucei and Staphylococcus aureus provides proof-of-function. Crystal structures of MDD from T. brucei (TbMDD, at 1.8 A resolution) and S. aureus (SaMDD, in two distinct crystal forms, each diffracting to 2.3 A resolution) have been determined. Gel-filtration chromatography and analytical ultracentrifugation experiments indicate that TbMDD is predominantly monomeric in solution while SaMDD is dimeric. The new crystal structures and comparison with that of the yeast Saccharomyces cerevisiae enzyme (ScMDD) reveal the structural basis for this variance in quaternary structure. The presence of an ordered sulfate in the structure of TbMDD reveals for the first time details of a ligand binding in the MDD active site and, in conjunction with well-ordered water molecules, comparisons with the related enzyme mevalonate kinase, structural and biochemical data derived on ScMDD and SaMDD, allows us to model a ternary complex with MDP and ATP. This model facilitates discussion of the molecular determinants of substrate recognition and contributions made by specific residues to the enzyme mechanism.     

Structural basis for bile salt inhibition of pancreatic phospholipase A2

Bile salt interactions with phospholipid monolayers of fat emulsions are known to regulate the actions of gastrointestinal lipolytic enzymes in order to control the uptake of dietary fat. Specifically, on the lipid/aqueous interface of fat emulsions, the anionic portions of amphipathic bile salts have been thought to interact with and activate the enzyme group-IB phospholipase A2 (PLA2) derived from the pancreas. To explore this regulatory process, we have determined the crystal structures of the complexes of pancreatic PLA2 with the naturally occurring bile salts: cholate, glycocholate, taurocholate, glycochenodeoxycholate, and taurochenodeoxycholate. The five PLA2-bile salt complexes each result in a partly occluded active site, and the resulting ligand binding displays specific hydrogen bonding interactions and extensive hydrophobic packing. The amphipathic bile salts are bound to PLA2 with their polar hydroxyl and sulfate/carboxy groups oriented away from the enzyme's hydrophobic core. The impaired catalytic and interface binding functions implied by these structures provide a basis for the previous numerous observations of a biphasic dependence of the rate of PLA2 catalyzed hydrolysis of zwitterionic glycerophospholipids in the presence of bile salts. The rising or activation phase is consistent with enhanced binding and activation of the bound PLA2 by the bile salt induced anionic charge in a zwitterionic interface. The falling or inhibitory phase can be explained by the formation of a catalytically inert stoichiometric complex between PLA2 and any bile salts in which it forms a stable complex. The model provides new insight into the regulatory role that specific PLA2-bile salt interactions are likely to play in fat metabolism.     

Shape variation in protein binding pockets and their ligands

A common assumption about the shape of protein binding pockets is that they are related to the shape of the small ligand molecules that can bind there. But to what extent is that assumption true? Here we use a recently developed shape matching method to compare the shapes of protein binding pockets to the shapes of their ligands. We find that pockets binding the same ligand show greater variation in their shapes than can be accounted for by the conformational variability of the ligand. This suggests that geometrical complementarity in general is not sufficient to drive molecular recognition. Nevertheless, we show when considering only shape and size that a significant proportion of the recognition power of a binding pocket for its ligand resides in its shape. Additionally, we observe a "buffer zone" or a region of free space between the ligand and protein, which results in binding pockets being on average three times larger than the ligand that they bind.     

Crystal structures of the luciferase and green fluorescent protein from Renilla reniformis

Due to its ability to emit light, the luciferase from Renilla reniformis (RLuc) is widely employed in molecular biology as a reporter gene in cell culture experiments and small animal imaging. To accomplish this bioluminescence, the 37-kDa enzyme catalyzes the degradation of its substrate coelenterazine in the presence of molecular oxygen, resulting in the product coelenteramide, carbon dioxide, and the desired photon of light. We successfully crystallized a stabilized variant of this important protein (RLuc8) and herein present the first structures for any coelenterazine-using luciferase. These structures are based on high-resolution data measured to 1.4 Å and demonstrate a classic α/β-hydrolase fold. We also present data of a coelenteramide-bound luciferase and reason that this structure represents a secondary conformational form following shift of the product out of the primary active site. During the course of this work, the structure of the luciferase's accessory green fluorescent protein (RrGFP) was also determined and shown to be highly similar to that of Aequorea victoria GFP.     

Structural analysis and dynamics of retinal chromophore in dark and meta I states of rhodopsin from 2H NMR of aligned membranes

Rhodopsin is a prototype for G protein-coupled receptors (GPCRs) that are implicated in many biological responses in humans. A site-directed (2)H NMR approach was used for structural analysis of retinal within its binding cavity in the dark and pre-activated meta I states. Retinal was labeled with (2)H at the C5, C9, or C13 methyl groups by total synthesis, and was used to regenerate the opsin apoprotein. Solid-state (2)H NMR spectra were acquired for aligned membranes in the low-temperature lipid gel phase versus the tilt angle to the magnetic field. Data reduction assumed a static uniaxial distribution, and gave the retinylidene methyl bond orientations plus the alignment disorder (mosaic spread). The dark-state (2)H NMR structure of 11-cis-retinal shows torsional twisting of the polyene chain and the beta-ionone ring. The ligand undergoes restricted motion, as evinced by order parameters of approximately 0.9 for the spinning C-C(2)H(3) groups, with off-axial fluctuations of approximately 15 degrees . Retinal is accommodated within the rhodopsin binding pocket with a negative pre-twist about the C11=C12 double bond that explains its rapid photochemistry and the trajectory of 11-cis to trans isomerization. In the cryo-trapped meta I state, the (2)H NMR structure shows a reduction of the polyene strain, while torsional twisting of the beta-ionone ring is maintained. Distortion of the retinal conformation is interpreted through substituent control of receptor activation. Steric hindrance between trans retinal and Trp265 can trigger formation of the subsequent activated meta II state. Our results are pertinent to quantum and molecular mechanics simulations of ligands bound to GPCRs, and illustrate how (2)H NMR can be applied to study their biological mechanisms of action.     

Crystal structure of pyruvate dehydrogenase phosphatase 1 and its functional implications

Pyruvate dehydrogenase phosphatase 1 (PDP1) catalyzes dephosphorylation of pyruvate dehydrogenase (E1) in the mammalian pyruvate dehydrogenase complex (PDC), whose activity is regulated by the phosphorylation-dephosphorylation cycle by the corresponding protein kinases (PDHKs) and phosphatases. The activity of PDP1 is greatly enhanced through Ca2+ -dependent binding of the catalytic subunit (PDP1c) to the L2 (inner lipoyl) domain of dihydrolipoyl acetyltransferase (E2), which is also integrated in PDC. Here, we report the crystal structure of the rat PDP1c at 1.8 A resolution. The structure reveals that PDP1 belongs to the PPM family of protein serine/threonine phosphatases, which, in spite of a low level of sequence identity, share the structural core consisting of the central beta-sandwich flanked on both sides by loops and alpha-helices. Consistent with the previous studies, two well-fixed magnesium ions are coordinated by five active site residues and five water molecules in the PDP1c catalytic center. Structural analysis indicates that, while the central portion of the PDP1c molecule is highly conserved among the members of the PPM protein family, a number of structural insertions and deletions located at the periphery of PDP1c likely define its functional specificity towards the PDC. One notable feature of PDP1c is a long insertion (residues 98-151) forming a unique hydrophobic pocket on the surface that likely accommodates the lipoyl moiety of the E2 domain in a fashion similar to that of PDHKs. The cavity, however, appears more open than in PDHK, suggesting that its closure may be required to achieve tight, specific binding of the lipoic acid. We propose a mechanism in which the closure of the lipoic acid binding site is triggered by the formation of the intermolecular (PDP1c/L2) Ca2+ binding site in a manner reminiscent of the Ca2+ -induced closure of the regulatory domain of troponin C.     

Crystal structure of tryptophanyl-tRNA synthetase complexed with adenosine-5' tetraphosphate: evidence for distributed use of catalytic binding energy in amino acid activation by class I aminoacyl-tRNA synthetases

Tryptophanyl-tRNA synthetase (TrpRS) is a functionally dimeric ligase, which specifically couples hydrolysis of ATP to AMP and pyrophosphate to the formation of an ester bond between tryptophan and the cognate tRNA. TrpRS from Bacillus stearothermophilus binds the ATP analogue, adenosine-5' tetraphosphate (AQP) competitively with ATP during pyrophosphate exchange. Estimates of binding affinity from this competitive inhibition and from isothermal titration calorimetry show that AQP binds 200 times more tightly than ATP both under conditions of induced-fit, where binding is coupled to an unfavorable conformational change, and under exchange conditions, where there is no conformational change. These binding data provide an indirect experimental measurement of +3.0 kcal/mol for the conformational free energy change associated with induced-fit assembly of the active site. Thermodynamic parameters derived from the calorimetry reveal very modest enthalpic changes, consistent with binding driven largely by a favorable entropy change. The 2.5 A structure of the TrpRS:AQP complex, determined de novo by X-ray crystallography, resembles that of the previously described, pre-transition state TrpRS:ATP complexes. The anticodon-binding domain untwists relative to the Rossmann-fold domain by 20% of the way toward the orientation observed for the Products complex. An unexpected tetraphosphate conformation allows the gamma and deltad phosphate groups to occupy positions equivalent to those occupied by the beta and gamma phosphates of ATP. The beta-phosphate effects a 1.11 A extension that relocates the alpha-phosphate toward the tryptophan carboxylate while the PPi mimic moves deeper into the KMSKS loop. This configuration improves interactions between enzyme and nucleotide significantly and uniformly in the adenosine and PPi binding subsites. A new hydrogen bond forms between S194 from the class I KMSKS signature sequence and the PPi mimic. These complementary thermodynamic and structural data are all consistent with the conclusion that the tetraphosphate mimics a transition-state in which the KMSKS loop develops increasingly tight bonds to the PPi leaving group, weakening linkage to the Palpha as it is relocated by an energetically favorable domain movement. Consistent with extensive mutational data on Tyrosyl-tRNA synthetase, this aspect of the mechanism develops high transition-state affinity for the adenosine and pyrophosphate moieties, which move significantly, relative to one another, during the catalytic step.