Characterization of hydrogen bonding in the complex of adenosine deaminase with a transition state analogue: a Raman spectroscopic study

Computational Studies of Catalytic Loop Dynamics in Yersinia Protein Tyrosine Phosphatase Using Pathway Optimization Methods

Yersinia Protein Tyrosine Phosphatase (YopH) is the most efficient enzyme among all known PTPases and relies on its catalytic loop movements for substrate binding and catalysis. Fluorescence, NMR, and UV resonance Raman (UVRR) techniques have been used to study the thermodynamic and dynamic properties of the loop motions. In this study, a computational approach based on the pathway refinement methods nudged elastic band (NEB) and harmonic Fourier beads (HFB) has been developed to provide structural interpretations for the experimentally observed kinetic processes. In this approach, the minimum potential energy pathways for the loop open/closure conformational changes were determined by NEB using a one-dimensional global coordinate. Two dimensional data analyses of the NEB results were performed as an efficient method to qualitatively evaluate the energetics of transitions along several specific physical coordinates. The free energy barriers for these transitions were then determined more precisely using the HFB method. Kinetic parameters were estimated from the energy barriers using transition state theory and compared against experimentally determined kinetic parameters. When the calculated energy barriers are calibrated by a simple "scaling factor", as have been done in our previous vibrational frequency calculations to explain the ligand frequency shift upon its binding to protein, it is possible to make structural interpretations of several observed enzyme dynamic rates. For example, the nanosecond kinetics observed by fluorescence anisotropy may be assigned to the translational motion of the catalytic loop and microsecond kinetics observed in fluorescence T-jump can be assigned to the loop backbone dihedral angle flipping. Furthermore, we can predict that a Trp354 conformational conversion associated with the loop movements would occur on the tens of nanoseconds time scale, to be verified by future UVRR T-jump studies.

Exquisite Modulation of the Active Site of Methanocaldococcus jannaschii Adenylosuccinate Synthetase in Forward Reaction Complexes

In enzymes that conduct complex reactions involving several substrates and chemical transformations, the active site must reorganize at each step to complement the transition state of that chemical step. Adenylosuccinate synthetase (ADSS) utilizes a molecule each of guanosine 5'-monophosphate (GTP) and aspartate to convert inosine 5'-monophosphate (IMP) into succinyl adenosine 5'-monophosphate (sAMP) through several kinetic intermediates. Here we followed catalysis by ADSS through high-resolution vibrational spectral fingerprints of each substrate and intermediate involved in the forward reaction. Vibrational spectra show differential ligand distortion at each step of catalysis, and band positions of substrates are influenced by binding of cosubstrates. We found that the bound IMP is distorted toward its N1-deprotonated form even in the absence of any other ligands. Several specific interactions between GTP and active-site amino acid residues result in large Raman shifts and contribute substantially to intrinsic binding energy. When both IMP and GTP are simultaneously bound to ADSS, IMP is converted into an intermediate 6-phosphoryl inosine 5'-monophosphate (6-pIMP). The 6-pIMP·ADSS complex was found to be stable upon binding of the third ligand, hadacidin (HDA), an analogue of l-aspartate. We find that in the absence of HDA, 6-pIMP is quickly released from ADSS, is unstable in solution, and converts back into IMP. HDA allosterically stabilizes ADSS through local conformational rearrangements. We captured this complex and determined the spectra and structure of 6-pIMP in its enzyme-bound state. These results provide important insights into the exquisite tuning of active-site interactions with changing substrate at each kinetic step of catalysis.

Solution structure of ligands involved in purine salvage pathway

Analogues of intermediates involved in the purine salvage pathway can be exploited as potential drug molecules against enzymes of protozoan parasites. To develop such analogues we need knowledge of the solution structures, predominant tautomer at physiological pH and protonation-state of the corresponding natural ligand. In this regard, we have employed ultraviolet resonance Raman spectroscopy (UVRR) in combination with density functional theory (DFT) to study the solution structures of two relatively unexplored intermediates, 6-phosphoryl IMP (6-pIMP) and succinyl adenosine-5'-monophosphate (sAMP), of purine salvage pathway. These molecules are intermediates in a two step enzymatic process that converts inosine-5'-monpophosphate (IMP) to adenosine-5'-monophosphate (AMP). Experimental data on the molecular structure of these ligands is lacking. We report UVRR spectra of these two ligands, obtained at an excitation wavelength of 260 nm. Using isotope induced shifts and DFT calculations we assigned observed spectra to computed normal modes. We find that sAMP exists as neutral species at physiological pH and the predominant tautomer in solution bears proton at N10 position of purine ring. Though transient in solution, 6-pIMP is captured in the enzyme-bound form. This work provides the structural information of these ligands in solution state at physiological pH. We further compare these structures with the structures of AMP and IMP. Despite the presence of similar purine rings in AMP and sAMP, their UVRR spectra are found to be very different. Similarly, though the purine ring in 6-pIMP resembles that of IMP, UVRR spectra of the two molecules are distinct. These differences in the vibrational spectra provide direct information on the effects of exocyclic groups on the skeletal structures of these molecules. Our results identify key bands in the vibrational spectra of these ligands which may serve as markers of hydrogen bonding interactions upon binding to the active-sites of enzymes.

8-Azapurines as isosteric purine fluorescent probes for nucleic acid and enzymatic research

The 8-azapurines, and their 7-deaza and 9-deaza congeners, represent a unique class of isosteric (isomorphic) analogues of the natural purines, frequently capable of substituting for the latter in many biochemical processes. Particularly interesting is their propensity to exhibit pH-dependent room-temperature fluorescence in aqueous medium, and in non-polar media. We herein review the physico-chemical properties of this class of compounds, with particular emphasis on the fluorescence emission properties of their neutral and/or ionic species, which has led to their widespread use as fluorescent probes in enzymology, including enzymes involved in purine metabolism, agonists/antagonists of adenosine receptors, mechanisms of catalytic RNAs, RNA editing, etc. They are also exceptionally useful fluorescent probes for analytical and clinical applications in crude cell homogenates.

Enzyme active site interactions by Raman/FTIR, NMR, and ab initio calculations

Characterization of enzyme active site structure and interactions at high resolution is important for the understanding of the enzyme catalysis. Vibrational frequency and NMR chemical shift measurements of enzyme-bound ligands are often used for such purpose when X-ray structures are not available or when higher resolution active site structures are desired. This review is focused on how ab initio calculations may be integrated with vibrational and NMR chemical shift measurements to quantitatively determine high-resolution ligand structures (up to 0.001 Å for bond length and 0.01 Å for hydrogen bonding distance) and how interaction energies between bound ligand and its surroundings at the active site may be determined. Quantitative characterization of substrate ionic states, bond polarizations, tautomeric forms, conformational changes and its interactions with surroundings in enzyme complexes that mimic ground state or transition state can provide snapshots for visualizing the substrate structural evolution along enzyme-catalyzed reaction pathway. Our results have shown that the integration of spectroscopic studies with theoretical computation greatly enhances our ability to interpret experimental data and significantly increases the reliability of the theoretical analysis.

Molecular recognition of adenosine deaminase: 15N NMR studies

The elucidation of the molecular recognition of adenosine deaminase (ADA), the interpretation of the catalytic mechanism, and the design of novel inhibitors are based mostly on data obtained for the crystalline state of the enzyme. To obtain evidence for molecular recognition of the physiologically relevant soluble enzyme, we studied its interactions with the in situ formed inhibitor, 6-OH-purine riboside (HDPR), by 1D-15N- and 2D-(1H-15N)- NMR using the labeled primary inhibitor [15N4]-PR. We synthesized both [15N4]-PR and an [15N4]-HDPR model, from relatively inexpensive 15N sources. The [15N4]-HDPR model was used to simulate H-bonding and possible Zn2+-coordination of HDPR with ADA. We also explored possible ionic interactions between PR and ADA by 15N-NMR monitored pH-titrations of [15N4]-PR. Finally, we investigated the [15N4]-PR-ADA 1:1 complex by 2D-(1H-15N) NMR. We found that HDPR recognition determinants in ADA do not include any ionic-interactions. HDPR N1 H is an H-bond acceptor, and not an H-bond donor. Despite the proximity of N7 to the Zn2+-ion, no coordination occurs; instead, N7 is an H-bond acceptor. We found an overall agreement between the crystallographic data for the crystallized ADA:HDPR complex and the 15N-NMR signals for the corresponding soluble complex. This finding justifies the use of ADA's crystallographic data for the design of novel inhibitors.

Raman spectroscopy of protein and nucleic acid assemblies

The Raman spectrum of a protein or nucleic acid consists of numerous discrete bands representing molecular normal modes of vibration and serves as a sensitive and selective fingerprint of three-dimensional structure, intermolecular interactions, and dynamics. Recent improvements in instrumentation, coupled with innovative approaches in experimental design, dramatically increase the power and scope of the method, particularly for investigations of large supramolecular assemblies. Applications are considered that involve the use of (a) time-resolved Raman spectroscopy to elucidate assembly pathways in icosahedral viruses, (b) polarized Raman microspectroscopy to determine detailed structural parameters in filamentous viruses, (c) ultraviolet-resonance Raman spectroscopy to probe selective DNA and protein residues in nucleoprotein complexes, and (d) difference Raman methods to understand mechanisms of protein/DNA recognition in gene regulatory and chromosomal complexes.

Computational studies on nonenzymatic and enzymatic pyridoxal phosphate catalyzed decarboxylations of 2-aminoisobutyrate

A computational study of nonenzymatic and enzymatic pyridoxal phosphate-catalyzed decarboxylation of 2-aminoisobutyrate (AIB) is presented. Four prototropic isomers of a model aldimine between AIB and 5'-deoxypyridoxal, with acetate interacting with the pyridine nitrogen, were employed in calculations of both gas phase and water model (PM3 and PM3-SM3) decarboxylation reaction paths. Calculations employing the transition state structures obtained for the four isomers allow the demonstration of stereoelectronic effects in transition state stabilization as well as a separation of the contributions of the Schiff base and pyridine ring moieties to this stabilization. The unprotonated Schiff base contribution (approximately 16 kcal/mol) is larger than that of the pyridine ring even when it is protonated (approximately 10 kcal/mol), providing an explanation of the catalytic power of pyruvoyl-dependent amino acid decarboxylases. An active site model of dialkylglycine decarboxylase was constructed and validated, and enzymatic decarboxylation reaction paths were calculated. The reaction coordinate is shown to be complex, with proton transfer from Lys272 to the coenzyme C4' likely simultaneous with C alpha--CO(2)(-) bond cleavage. The proposed concerted decarboxylation/proton-transfer mechanism provides a simple explanation for the observed specificity of this enzyme toward oxidative decarboxylation.

A vibrational structure of 7,8-dihydrobiopterin bound to dihydroneopterin aldolase

Dihydroneopterin aldolase (DHNA) catalyzes the conversion of 7, 8-dihydroneopterin to 6-hydroxymethyl-7,8-dihydropterin and glycolaldehyde. An inhibitor of the enzyme, 7,8-dihydrobiopterin, free in solution and bound in its complex with the enzyme has been studied by Raman difference spectroscopy. By using isotopically labeled 7,8-dihydrobiopterin and normal mode analyses based on ab initio quantum mechanic methods, we have positively identified some of the Raman bands in the enzyme-bound inhibitor, particularly the important N5=C6 stretch mode. The spectrum of the enzyme-bound inhibitor shows that the pK(a) of N5 is not significantly increased in the complex. This result suggests that N5 of 7,8-dihydroneopterin is not protonated before the bond cleavage of 7,8-dihydroneopterin during the DHNA-catalyzed reaction as has been suggested. Our results also show that the N5=C6 stretch mode of 7, 8-dihydrobiopterin shifts 19 cm(-)(1) upon binding to DHNA. Various possibilities on how the enzyme can bring about such large frequency change of the N5=C6 stretch mode are discussed.