Accurate proton affinity and gas-phase basicity values for molecules important in biocatalysis.

Benchmark quantum calculations of proton affinities and gas-phase basicities of molecules relevant to biochemical processes, particularly acid/base catalysis, are presented and compared for a variety of multilevel and density functional quantum models. Included are nucleic acid bases in both keto and enol tautomeric forms, ribose in B-form and A-form sugar pucker conformations, amino acid side chains and backbone molecules, and various phosphates and phosphoranes, including thio substitutions. This work presents a high-level thermodynamic characterization of biologically relevant protonation states and provides a benchmark database for development of next-generation semiempirical and approximate density functional quantum models and parametrization of methods to predict pK(a) values and relative solvation energies.

ESEEM studies of peptide nitrogen hyperfine coupling in tyrosyl radicals and model peptides.

Tyrosyl radicals are important in long-range electron transfer in several enzymes, but the protein environmental factors that control midpoint potential and electron transfer rate are not well understood. To develop a more detailed understanding of the effect of protein sequence, we have performed 14N and 15N electron spin echo envelope modulation (ESEEM) measurements on tyrosyl radical, generated either in polycrystalline tyrosinate or in its 15N-labeled isotopomer, by UV photolysis. 14N-ESEEM was also performed on tyrosyl radical generated in tyrosine-containing pentapeptide samples. Simulation of the 14N- and 15N-tyrosyl radical ESEEM measurements yielded no significant isotropic hyperfine splitting to the amine or amide nitrogen; the amplitude of the anisotropic, nitrogen hyperfine coupling (0.21 MHz) was consistent with a dipole-dipole distance of 3.0 A. Density functional theory was used to calculate the isotropic and anisotropic hyperfine couplings to the amino nitrogen in four different tyrosyl radical conformers. Comparison with the simulated data suggested that the lowest energy radical conformer, generated in tyrosine at pH 11, has a 76 degrees Calpha-Cbeta-C1'-C2' ring and a -73 degrees C-Calpha-Cbeta-C1' backbone dihedral angle. In addition, the magnitude, orientation, and asymmetry of the nuclear quadrupole coupling tensor were derived from analysis of the tyrosyl radical 14N-ESEEM. The simulations showed differences in the coupling and orientation of the nuclear quadrupole tensor, when the tyrosinate and pentapeptide samples were compared. These results suggest sequence- or conformation-induced changes in the ionic character of the NH bond in different tyrosine-containing peptides.

Proton-coupled electron transfer in a biomimetic peptide as a model of enzyme regulatory mechanisms.

Proton-coupled electron-transfer reactions are central to enzymatic mechanism in many proteins. In several enzymes, essential electron-transfer reactions involve oxidation and reduction of tyrosine side chains. For these redox-active tyrosines, proton transfer couples with electron transfer, because the phenolic pKA of the tyrosine is altered by changes in the tyrosine redox state. To develop an experimentally tractable peptide system in which the effect of proton and electron coupling can be investigated, we have designed a novel amino acid sequence that contains one tyrosine residue. The tyrosine can be oxidized by ultraviolet photolysis or electrochemical methods and has a potential cross-strand interaction with a histidine residue. NMR spectroscopy shows that the peptide forms a beta-hairpin with several interstrand dipolar contacts between the histidine and tyrosine side chains. The effect of the cross-strand interaction was probed by electron paramagnetic resonance and electrochemistry. The data are consistent with an increase in histidine pKA when the tyrosine is oxidized; the effect of this thermodynamic coupling is to increase tyrosyl radical yield at low pH. The coupling mechanism is attributed to an interstrand pi-cation interaction, which stabilizes the tyrosyl radical. A similar interaction between histidine and tyrosine in enzymes provides a regulatory mechanism for enzymatic electron-transfer reactions.

Normal modes of redox-active tyrosine: conformation dependence and comparison to experiment.

Redox-active tyrosine residues play important roles in long-distance electron reactions in enzymes such as prostaglandin H synthase, ribonucleotide reductase, and photosystem II (PSII). Spectroscopic characterization of tyrosyl radicals in these systems provides a powerful experimental probe into the role of the enzyme in mediation of long-range electron transfer processes. Interpretation of such data, however, relies critically on first establishing a spectroscopic fingerprint of isotopically labeled tyrosinate and tyrosyl radicals in nonenzymatic environments. In this report, FT-IR results obtained from tyrosinate, tyrosyl radical (produced by ultraviolet photolysis of polycrystalline tyrosinate), and their isotopologues at 77 K are presented. Assignment of peaks and isotope shifts is aided by density-functional B3LYP/6-311++G(3df,2p)//B3LYP/6-31++G(d,p) calculations of tyrosine and tyrosyl radical in several different charge and protonation states. In addition, characterization of the potential energy surfaces of tyrosinate and tyrosyl radical as a function of the backbone and ring torsion angles provides detailed insight into the sensitivity of the vibrational frequencies to conformational changes. These results provide a detailed spectroscopic interpretation, which will elucidate the structures of redox-active tyrosine residues in complex protein environments. Specific application of these data is made to enzymatic systems.

QCRNA 1.0: a database of quantum calculations for RNA catalysis.

This work outlines a new on-line database of quantum calculations for RNA catalysis (QCRNA) available via the worldwide web at http://theory.chem.umn.edu/QCRNA. The database contains high-level density functional calculations for a large range of molecules, complexes and chemical mechanisms important to phosphoryl transfer reactions and RNA catalysis. Calculations are performed using a strict, consistent protocol such that a wealth of cross-comparisons can be made to elucidate meaningful trends in biological phosphate reactivity. Currently, around 2000 molecules have been collected in varying charge states in the gas phase and in solution. Solvation was treated with both the PCM and COSMO continuum solvation models. The data can be used to study important trends in reactivity of biological phosphates, or used as benchmark data for the design of new semiempirical quantum models for hybrid quantum mechanical/molecular mechanical simulations.

Multilevel and density functional electronic structure calculations of proton affinities and gas-phase basicities involved in biological phosphoryl transfer.

Five multilevel model chemistries (CBS-QB3, G3B3, G3MP2B3, MCG3/3, and MC-QCISD/3) and seven hybrid density functional methods (PBE0, B1B95, B3LYP, MPW1KCIS, PBE1KCIS, and MPW1B95) have been applied to the calculation of gas-phase basicity and proton affinity values for a series of 17 molecules relevant to the study of biological phosphoryl transfer. In addition, W1 calculations were performed on a subset of molecules. The accuracy of the methods was assessed and the nature of systematic errors was explored, leading to the introduction of a set of effective bond enthalpy and entropy correction terms. The multicoefficient correlation methods (MCG3/3 and MC-QCISD), with inclusion of specific zero-point scale factors, slightly outperform the other multilevel methods tested (CBS-QB3, G3B3, and G3MP2B3), with significantly less computational cost, and in the case of MC-QCISD, slightly less severe scaling. Four density functional methods, PBE1KCIS, MPW1B95, PBE0, and B1B95 perform nearly as well as the multilevel methods. These results provide an important set of benchmarks relevant to biological phosphoryl transfer reactions.

Benchmark calculations of proton affinities and gas-phase basicities of molecules important in the study of biological phosphoryl transfer.

Benchmark calculations of proton affinities and gas-phase basicities of molecules most relevant to biological phosphoryl transfer reactions are presented and compared with available experimental results. The accuracy of proton affinity and gas-phase basicity results obtained from several multi-level model chemistries (CBS-QB3, G3B3, and G3MP2B3) and density-functional quantum models (PBE0, B1B95, and B3LYP) are assessed and compared. From these data, a set of empirical bond enthalpy, entropy, and free energy corrections are introduced that considerably improve the accuracy and predictive capability of the methods. These corrections are applied to the prediction of proton affinity and gas-phase basicity values of important biological phosphates and phosphoranes for which experimental data does not currently exist. Comparison is made with results from semiempirical quantum models that are commonly employed in hybrid quantum mechanical/molecular mechanical simulations. Data suggest that the design of improved semiempirical quantum models with increased accuracy for relative proton affinity values is necessary to obtain quantitative accuracy for phosphoryl transfer reactions in solution, enzymes, and ribozymes.

The contribution of phosphate-phosphate repulsions to the free energy of DNA bending.

DNA bending is important for the packaging of genetic material, regulation of gene expression and interaction of nucleic acids with proteins. Consequently, it is of considerable interest to quantify the energetic factors that must be overcome to induce bending of DNA, such as base stacking and phosphate-phosphate repulsions. In the present work, the electrostatic contribution of phosphate-phosphate repulsions to the free energy of bending DNA is examined for 71 bp linear and bent-form model structures. The bent DNA model was based on the crystallographic structure of a full turn of DNA in a nucleosome core particle. A Green's function approach based on a linear-scaling smooth conductor-like screening model was applied to ascertain the contribution of individual phosphate-phosphate repulsions and overall electrostatic stabilization in aqueous solution. The effect of charge neutralization by site-bound ions was considered using Monte Carlo simulation to characterize the distribution of ion occupations and contribution of phosphate repulsions to the free energy of bending as a function of counterion load. The calculations predict that the phosphate-phosphate repulsions account for approximately 30% of the total free energy required to bend DNA from canonical linear B-form into the conformation found in the nucleosome core particle.

Structure and binding of Mg(II) ions and di-metal bridge complexes with biological phosphates and phosphoranes.

Divalent Mg(2+) ions often serve as cofactors in enzyme or ribozyme-catalyzed phosphoryl transfer reactions. In this work, the interaction of Mg(2+) ions and di-metal bridge complexes with phosphates, phosphoranes, and other biological ligands relevant to RNA catalysis are characterized with density functional methods. The effect of bulk solvent is treated with two continuum solvation methods (PCM and COSMO) for comparison. The relative binding affinity for different biological ligands to Mg(2+) are quantified in different protonation states. The structure and stability of the single-metal and di-metal complexes are characterized, and the changes in phosphate and phosphorane geometry induced by metal ion binding are discussed. Di-metal bridge complexes are a ubiquitous motif and the key factors governing their electrostatic stabilization are outlined. The results presented here provide quantitative characterization of metal ion binding to ligands of importance to RNA catalysis, and lay the groundwork for design of new generation quantum models that can be applied to the full biological enzymatic systems.

The structure and stability of biological metaphosphate, phosphate, and phosphorane compounds in the gas phase and in solution.

Density functional calculations of a series of metaphosphates, acyclic and cyclic phosphates and phosphoranes relevant to RNA catalysis are presented. Solvent effects calculated with three well-established solvation models are analyzed and compared. The structure and stability of the compounds are characterized in terms of thermodynamic quantities for isomerization and ligand substitution reactions, gas-phase proton affinities, and microscopic solution pK(a)() values. The large dataset of compounds allows the estimation of bond energies to determine the relative strengths of axial and equatorial P-O phosphorane single bonds and P-O single and double bonds in metaphosphates and phosphates. The relative apicophilicty of hydroxyl and methoxy ligands in phosphoranes are characterized. The results presented here provide quantitative insight into RNA catalysis and serve as a first step toward the construction of a high-level quantum database for development of new semiempirical Hamiltonian models for biological reactions

Spectroscopic properties of tyrosyl radicals in dipeptides.

Redox-active tyrosine residues play important roles in long-distance electron reactions in enzymes, including prostaglandin H synthase, galactose oxidase, ribonucleotide reductase, and photosystem II. Magnetic resonance and vibrational spectroscopy provide methods with which to study the structures of redox-active amino acids in proteins. In this report, ultraviolet photolysis was used to generate tyrosyl radicals from polycrystalline tyrosinate or dipeptides, and the structure of the radical was investigated with EPR and reaction-induced FT-IR spectroscopy at 77 K. Photolysis at 77 K is expected to generate a neutral tyrosyl radical through oxidation of the aromatic ring. EPR and FT-IR results obtained from (13)C-labeled tyrosine were consistent with that expectation. Surprisingly, labeling of the tyrosyl amino group with (15)N also resulted in isotope-shifted bands in the photolysis spectrum. The force constant of a NH deformation mode increased when the tyrosyl radical was generated. These data suggest an interaction between the pi system of the tyrosyl radical and the amino group. In spectra acquired from the dipeptides, evidence for a sequence-dependent interaction between the tyrosyl radical and the amide bond of the dipeptide was also obtained. We postulate that perturbation of the amino or the amide/imide groups may occur through a spin polarization mechanism, which is indirectly detected as a change in NH force constant. This conclusion is supported by density functional calculations, which suggest a conformationally sensitive delocalization of spin density onto the amino and carboxylate groups of the tyrosyl radical. These experiments provide a step toward a detailed spectral interpretation for protein-based tyrosyl radicals.