Computational investigation of proton transfer,pKa shifts and pH‐optimum of protein–DNA and protein–RNA complexes

Protein–nucleic acid interactions play a crucial role in many biological processes. This work investigates the changes of pKa values and protonation states of ionizable groups (including nucleic acid bases) that may occur at protein–nucleic acid binding. Taking advantage of the recently developed pKa calculation tool DelphiPka, we utilize the large protein–nucleic acid interaction database (NPIDB database) to model pKa shifts caused by binding. It has been found that the protein's interfacial basic residues experience favorable electrostatic interactions while the protein acidic residues undergo proton uptake to reduce the energy cost upon the binding. This is in contrast with observations made for protein–protein complexes. In terms of DNA/RNA, both base groups and phosphate groups of nucleotides are found to participate in binding. Some DNA/RNA bases undergo pKa shifts at complex formation, with the binding process tending to suppress charged states of nucleic acid bases. In addition, a weak correlation is found between the pH‐optimum of protein–DNA/RNA binding free energy and the pH‐optimum of protein folding free energy. Overall, the pH‐dependence of protein–nucleic acid binding is not predicted to be as significant as that of protein–protein association. Proteins 2017; 85:282–295. © 2016 Wiley Periodicals, Inc.     

Analysis of difference spectra of protonated DNA: determination of degree of protonation of nitrogen bases and the fractions of disordered nucleotide pairs.

The titration curves of nitrogen bases and fractions of disordered nucleotide pairs are obtained during DNA protonation. It is shown that purine bases are the first sites of the DNA double helix protonation. The cytosine protonation is due to proton-induced conformational transition within GC pairs with the sequence proton transfer from (N-7) of guanine to (N-3) of cytosine. Within DNA with unwound regions the bases are protonated in the following order: cytosine, adenine, guanine. It is shown that GC pairs are the primary centres in which the unwinding of protonated DNAs occurs.     

Protonation free energy levels in complex molecular systems

All proteins, nucleic acids, and other biomolecules contain residues capable of exchanging protons with their environment. These proton transfer phenomena lead to pH sensitivity of many molecular processes underlying biological phenomena. In the course of biological evolution, Nature has invented some mechanisms to use pH gradients to regulate biomolecular processes inside cells or in interstitial fluids. Therefore, an ability to model protonation equilibria in molecular systems accurately would be of enormous value for our understanding of biological processes and for possible rational influence on them, like in developing pH dependent drugs to treat particular diseases. This work presents a derivation, by thermodynamic and statistical mechanical methods, of an expression for the free energy of a complex molecular system at arbitrary ionization state of its titratable residues. This constitutes one of the elements of modeling protonation equilibria. Starting from a consideration of a simple acid-base equilibrium of a model compound with a single tritratable group, we arrive at an expression which is of general validity for complex systems. The only approximation used in this derivation is the postulating that the interaction energy between any pair of titratable sites does not depend on the protonation states of all the remaining ionizable groups.     

Acid-induced exchange of the imino proton in G.C pairs.

Acid-induced catalysis of imino proton exchange in G.C pairs of DNA duplexes is surprisingly fast, being nearly as fast as for the isolated nucleoside, despite base-pair dissociation constants in the range of 10(-5) at neutral or basic pH. It is also observed in terminal G.C pairs of duplexes and in base pairs of drug-DNA complexes. We have measured imino proton exchange in deoxyguanosine and in the duplex (ATATAGATCTATAT) as a function of pH. We show that acid-induced exchange can be assigned to proton transfer from N7-protonated guanosine to cytidine in the open state of the pair. This is faster than transfer from neutral guanosine (the process of intrinsic catalysis previously characterized at neutral ph) due to the lower imino proton pK of the protonated form, 7.2 instead of 9.4. Other interpretations are excluded by a study of exchange catalysis by formiate and cytidine as exchange catalysts. The cross-over pH between the regimes of pH-independent and acid-induced exchange rates is more basic in the case of base pairs than in the mononucleoside, suggestive of an increase by one to two decades in the dissociation constant of the base pair upon N7 protonation of G. Acid-induced catalysis is much weaker in A.T base pairs, as expected in view of the low pK for protonation of thymidine.     

Models of interaction between nucleic acids and proteins. Hydrogen bonding of arginine with nucleic acid bases, phosphate groups and carboxylic acids.

Complex formation between the side chain of arginine and nucleic acid bases has been investigated by proton magnetic resonance in dimethylsulfoxide. Simultaneous formation of two hydrogen bonds leads to a selectivity of arginine interaction towards cytosine and guanine. A comparison is made of the interaction of arginine side chain with nucleic acid bases, phosphate and carboxylate anions. It is shown that interaction between carboxylate and arginine is stronger than between phosphate and arginine. These results are discussed with respect to the selective recognition of nucleic acid bases by arginine side chains and by the arginyl-glutamyl ion pair which could form in proteins interacting with nucleic acids.     

The protonation state of the Glu-71/Asp-80 residues in the KcsA potassium channel: a first-principles QM/MM molecular dynamics study

Although a few x-ray structures of the KcsA K(+) channel have been crystallized several issues concerning the mechanisms of the ionic permeation and the protonation state of the selectivity filter ionizable side chains are still open. Using a first-principles quantum mechanical/molecular mechanical simulation approach, we have investigated the protonation state of Glu-71 and Asp-80, two important residues located in the vicinity of the selectivity filter. Results from the dynamics show that a proton is shared between the two residues, with a slight preference for Glu-71. The proton is found to exchange on the picosecond timescale, an interesting phenomenon that cannot be observed in classical molecular dynamics. Simulations of different ionic loading states of the filter show that the probability for the proton transfer is correlated with the filter occupancy. In addition, the Glu-71/Asp-80 pair is able to modulate the potential energy profile experienced by a K(+) ion as it translates along the pore axis. These theoretical predictions, along with recent experimental results, suggest that changes of the filter structure could be associated with a shift in the Glu-Asp protonation state, which in turn would influence the ion translocation.     

Monte Carlo simulation of hydration of the nucleic acid fragments

Systems containing a base or a base pair and 25 water molecules, as well as a helical stack and 30 water molecules per base pair, have been simulated. Changes in the base hydration shell structure, after the bases have been included into the pair and then into the base pair stack, are discussed. Hydration shells of several configurations of the base pair stacks are discussed. Probabilities of formation of the hydrogen-bonded bridges of 1, 2 and 3 water molecules between hydrophilic centres have been estimated. The hydration shell structure was shown to depend on the nature of the base pair and on the stack configuration, while dependence of the global hydration shell characteristics on the stack configuration has been proved to be rather slight. The most typical structural elements of hydration shells, in the glycosidic (minor in B-like conformation) and non-glycosidic (major) grooves, for different configurations of AU and GC stacks, have been found and discussed. The number of hydrogen bonds between water molecules and bases per water molecule was shown to change upon transformation of the stack from A to B configuration. This result is discussed in connection with the reasons for B to A conformational transition and the concept of "water economy". Hydration shell patterns of NH2-groups of AU and GC helical stacks differ significantly.     

Effect of intramolecular interaction on the electron excitation state of nucleic acid components

A theoretical and experimental investigation of absorption and luminescence features of crystals and aggregates of nucleic acid bases were carried out. The long wavelength low intensity bands in UV-absorption nd excitation spectra, bathochromic shift of fluorescence spectra, the change of correlation between the intensity of fluorescence and phosphorescence spectra were obtained. The interpretation of these experimental results was proposed on the basis of pair interaction calculations (exciton-resonance and charge-resonance) in different conformations of cytosine dimers. The energy transfer after excitation at lambda 280 and 312 nm was investigated in nucleic acid base aggregates.     

Catalytic roles for proton transfer and protonation in ribozymes

Utilization of proton transfer in catalysis, which is well known in the mechanisms of protein enzymes, has been described only relatively recently for RNA enzymes. In this article, we present a current understanding of proton transfer by nucleic acids. Rate enhancement and specificity conferred by general acid-base catalysis are discussed. We also present possibilities for electrostatic catalysis from general acids and bases as well as cationic base pairs. The microenvironments of a large RNA provide the possibility of histidine-like pK(a)s for proton transfer, as well as lysine- and arginine-like pK(a)s for electrostatic catalysis. Discussion on proton transfer focuses on the hepatitis delta virus (HDV) and hairpin ribozymes, with select examples drawn from the protein literature. Discussion on electrostatic catalysis also draws on these two ribozymes, and a postulate for electrostatic catalysis by a cationic base pair in the mechanism of peptidyl transfer in the ribosome is presented. We also provide a perspective on possibilities for phosphoryl transfer mechanisms involving phosphorane intermediates and unusual tautomeric forms of the bases. Lastly, a distinction is made between ground state and "transition state" pK(a)s. We favor a model in which changes in pH lead to changes in the distribution of reactive and nonreactive ionizations of the ribozyme molecules in the ground state, and therefore suggest that "pK(a) changes in the transition state" do not provide an acceptable explanation for observed pH-rate profiles.     

Computational delineation of the catalytic step of a high‐fidelity DNA polymerase

The Bacillus fragment, belonging to a class of high‐fidelity polymerases, demonstrates high processivity (adding ～115 bases per DNA binding event) and exceptional accuracy (1 error in 10⁶ nucleotide incorporations) during DNA replication. We present analysis of structural rearrangements and energetics just before and during the chemical step (phosphodiester bond formation) using a combination of classical molecular dynamics, mixed quantum mechanics molecular mechanics simulations, and free energy computations. We find that the reaction is associative, proceeding via the two‐metal‐ion mechanism, and requiring the proton on the terminal primer O3′ to transfer to the pyrophosphate tail of the incoming nucleotide before the formation of the pentacovalent transition state. Different protonation states for key active site residues direct the system to alternative pathways of catalysis and we estimate a free energy barrier of ～12 kcal/mol for the chemical step. We propose that the protonation of a highly conserved catalytic aspartic acid residue is essential for the high processivity demonstrated by the enzyme and suggest that global motions could be part of the reaction free energy landscape.     

Conformational characteristics of the hexanucleoside pentaphosphate AUAUAU: a 2D NMR study at 500 MHz.

All 36 ribose proton resonances and most of the base proton resonances of the hexanucleoside pentaphosphate AUAUAU have been assigned unequivocally using 2D J-resolved spectroscopy, spin echo correlated spectroscopy (SECSY) and 2D NOE spectroscopy (NOESY). The NMR parameters of AUAUAU are compared with those of smaller fragments that contain methylated adenine bases: m62AU, m62AUm62A, m62AUm62AU and m62AUm62AUm62A. Previous studies on this series of compounds have shown that in all these cases purine-pyrimidine-purine sequences prefer to adopt a mixture of states which have as common feature that the interior pyrimidine residues are bulged out, whereas the purine residues stack upon each other. Chemical shift data, proton-proton coupling constants, as well as the observation of imino-proton resonances for AUAUAU show unambiguously that upon lowering the temperature the high-temperature "bulged out" situation reverts to a normal A-RNA-like double helix.     

Protonation changes upon ligand binding to trypsin and thrombin: structural interpretation based on pK(a) calculations and ITC experiments

The protonation states of a protein and a ligand can be altered upon complex formation. Such changes can be detected experimentally by isothermal titration calorimetry (ITC). For a series of ligands binding to the serine proteases trypsin and thrombin, we previously performed an extensive ITC and crystallographic study and were able to identify protonation changes for four complexes. However, since ITC measures only the overall proton exchange, it does not provide structural insights into the functional groups involved in the proton transfer. Using Poisson-Boltzmann calculations based on our recently developed PEOE_PB charges, we compute pK(a) values for all complexes of our former study in order to reveal the residues with altered protonation states. The results indicate that His57, a member of the catalytic triad, is responsible for the most relevant pK(a) shifts leading to the experimentally detected protonation changes. This finding is in contrast to our previous assumption that the observed protonation changes occur at the carboxylic group of the ligands. The newly detected proton acceptor is used for a revised factorization of the ITC data, which is necessary whenever the protonation inventory changes upon complexation. The pK(a) values of complexes showing no protonation change in the ITC experiment are reliably predicted in most cases, whereas predictions of strongly coupled systems remain problematic.     

Fluorescence and binding properties of phenazine derivatives in complexes with polynucleotides of various base compositions and secondary structures

The interactions of two phenazine derivatives, one with a neutral chromophore (glycoside) and the other with a cationic one (quaternary salt), with various synthetic single- and double-stranded polynucleotides and natural DNA were studied by fluorescence techniques, conducting measurements of steady-state fluorescence intensity and polarization degree as well as fluorescence lifetime. These dyes show fluorescence quenching upon intercalation into the GC sequences of the double-stranded nucleic acids and an increase in fluorescence emission and lifetime upon incorporation into the AT and AU sequences. GC base pairs in continuous deoxynucleotide sequences were found to be preferred as binding sites for both phenazines, in contrast to AT base pairs. On the contrary, the continuous ribonucleotide GC sequence binds the phenazines more weakly than does the AU sequence. With regard to the interaction of the phenazines with single-stranded polynucleotides, a stacking interaction of the dye chromophores with the nucleic bases was observed. In that case the guanine residue quenches the cationic phenazine fluorescence, while the stacking interaction with the other bases results in an increase in the fluorescence quantum yield. Unlike the cationic dye, the fluorescence of the neutral phenazine was quenched by both purine bases.     

Monte Carlo Simulation of Hydration of the Nucleic Acid Fragments

Abstract

Systems containing a base or a base pair and 25 water molecules, as well as a helical stack and 30 water molecules per base pair, have been simulated. Changes in the base hydration shell structure, after the bases have been included into the pair and then into the base pair stack are discussed. Hydration shells of several configurations of the base pair stacks are discussed. Probabilities of formation of the hydrogen-bonded bridges of 1, 2 and 3 water molecules between hydrophilic centres have been estimated. The hydration shell structure was shown to depend on the nature of the base pair and on the stack configuration, while dependence of the global hydration shell characteristics on the stack configuration has been proved to be rather slight. The most typical structural elements of hydration shells, in the glycosidic (minor in B-like conformation) and non-glycosidic (major) grooves, for different configurations of AU and GC stacks, have been found and discussed. The number of hydrogen bonds between water molecules and bases per water molecule was shown to change upon transformation of the stack from A to B configuration. This result is discussed in connection with the reasons for B to A conformational transition and the concept of “water economy”. Hydration shell patterns of NH₂-groups of AU and GC helical stacks differ significantly.     

Theoretical insights into the protonation states of active site cysteine and citrullination mechanism of Porphyromonas gingivalis peptidylarginine deiminase

Porphyromonas gingivalis peptidylarginine deiminase (PPAD) catalyzes the citrullination of peptidylarginine, which plays a critical role in the rheumatoid arthritis (RA) and gene regulation. For a better understanding of citrullination mechanism of PPAD, it is required to establish the protonation states of active site cysteine, which is still a controversial issue for the members of guanidino‐group‐modifying enzyme superfamily. In this work, we first explored the transformation between the two states: State N (both C351 and H236 are neutral) and State I (both residues exist as a thiolate–imidazolium ion pair), and then investigated the citrullination reaction of peptidylarginine, using a combined QM/MM approach. State N is calculated to be more stable than State I by 8.46 kcal/mol, and State N can transform to State I via two steps of substrate‐assisted proton transfer. Citrullination of the peptidylarginine contains deamination and hydrolysis. Starting from State N, the deamination reaction corresponds to an energy barrier of 18.82 kcal/mol. The deprotonated C351 initiates the nucleophilic attack to the substrate, which is the key step for deamination reaction. The hydrolysis reaction contains two chemical steps. Both the deprotonated D238 and H236 can act as the bases to activate the hydrolytic water, which correspond to similar energy barriers (～17 kcal/mol). On the basis of our calculations, C351, D238, and H236 constitute a catalytic triad, and their protonation states are critical for both the deamination and hydrolysis processes. In view of the sequence similarity, these findings may be shared with human PAD1–PAD4 and other guanidino‐group‐modifying enzymes. Proteins 2017; 85:1518–1528. © 2017 Wiley Periodicals, Inc.     

Electronic and molecular structure of M-DNA fragments

To study M-DNA molecular structure (such DNA with transition metal ions placed between the nucleic bases is able to conduct the electric current) and its conductivity mechanisms, we carried out ab initio quantum-mechanical calculations of electronic and spatial structures, thermodynamic characteristics of adenine-thymine (АТ) and guanine-cytosine (GC) base pair complexes with Zn²⁺ and Ni²⁺. To take into account the influence of the alkaline environment, calculations for these complexes were also carried out with hydroxyl and two water molecules. Computations were performed at MP2 level of theory using 6–31+G* basis set. Analogous calculations were carried out for (AC)(TG) stacking dimer of nucleic acid base pairs with two Zn²⁺. The calculation of the interaction energy in complexes has shown the preference of locating the metal ion (instead of the imino proton) between bases in M-DNA. The electronic transition energy calculation has revealed the reduction of the first singlet transition energy in АТ and GC complexes with Ni²⁺ from 4.5 eV to 0.4 - 0.6 eV. Ni²⁺ orbitals take part in the formation of HOMO and LUMO on the complexes investigated. It was shown that charges of metal ions incorporated into complexes with nucleic bases and in dimer decrease significantly.     

Computational and experimental studies on the catalytic mechanism of biliverdin-IXbeta reductase

BVR-B (biliverdin-IXbeta reductase) also known as FR (flavin reductase) is a promiscuous enzyme catalysing the pyridine-nucleotide-dependent reduction of a variety of flavins, biliverdins, PQQ (pyrroloquinoline quinone) and ferric ion. Mechanistically it is a good model for BVR-A (biliverdin-IXalpha reductase), a potential pharmacological target for neonatal jaundice and also a potential target for adjunct therapy to maintain protective levels of biliverdin-IXalpha during organ transplantation. In a commentary on the structure of BVR-B it was noted that one outstanding issue remained: whether the mechanism was a concerted hydride transfer followed by protonation of a pyrrolic anion or protonation of the pyrrole followed by hydride transfer. In the present study we have attempted to address this question using QM/MM (quantum mechanics/molecular mechanics) calculations. QM/MM potential energy surfaces show that the lowest energy pathway proceeds with a positively charged pyrrole intermediate via two transition states. These initial calculations were performed with His(153) as the source of the proton. However site-directed mutagenesis studies with both the H153A and the H153N mutant reveal that His(153) is not required for catalytic activity. We have repeated the calculation with a solvent hydroxonium donor and obtain a similar energy landscape indicating that protonation of the pyrrole is the most likely first step followed by hydride transfer and that the required proton may come from bulk solvent. The implications of the present study for the design of inhibitors of BVR-A are discussed.     

Water-hydroxide exchange reactions at the catalytic site of heme-copper oxidases

Membrane-bound heme-copper oxidases catalyze the reduction of O(2) to water. Part of the free energy associated with this process is used to pump protons across the membrane. The O(2) reduction reaction results in formation of high-pK(a) protonatable groups at the catalytic site. The free energy associated with protonation of these groups is used for proton pumping. One of these protonatable groups is OH(-), coordinated to the heme and Cu(B) at the catalytic site. Here we present results from EPR experiments on the Rhodobacter sphaeroides cytochrome c oxidase, which show that at high pH (9) approximately 50% of oxidized heme a(3) is hydroxide-ligated, while at low pH (6.5), no hydroxide is bound to heme a(3). The kinetics of hydroxide binding to heme a(3) were investigated after dissociation of CO from heme a(3) in the enzyme in which the heme a(3)-Cu(B) center was reduced while the remaining redox sites were oxidized. The dissociation of CO results in a decrease of the midpoint potential of heme a(3), which results in electron transfer (tau approximately equal 3 micros) from heme a(3) to heme a in approximately 100% of the enzyme population. At pH >7.5, the electron transfer is followed by proton release from a H(2)O molecule to the bulk solution (tau approximately equal 2 ms at pH 9). This reaction is also associated with absorbance changes of heme a(3), which on the basis of the results from the EPR experiments are attributed to formation of hydroxide-ligated heme a(3). The OH(-) bound to heme a(3) under equilibrium conditions at high pH is also formed transiently after O(2) reduction at low pH. It is proposed that the free energy associated with electron transfer to the binuclear center and protonation of this OH(-) upon reduction of the recently oxidized enzyme provides the driving force for the pumping of one proton.     

Unraveling the mechanism of proton translocation in the extracellular half‐channel of bacteriorhodopsin

Bacteriorhodopsin, a light activated protein that creates a proton gradient in halobacteria, has long served as a simple model of proton pumps. Within bacteriorhodopsin, several key sites undergo protonation changes during the photocycle, moving protons from the higher pH cytoplasm to the lower pH extracellular side. The mechanism underlying the long‐range proton translocation between the central (the retinal Schiff base SB216, D85, and D212) and exit clusters (E194 and E204) remains elusive. To obtain a dynamic view of the key factors controlling proton translocation, a systematic study using molecular dynamics simulation was performed for eight bacteriorhodopsin models varying in retinal isomer and protonation states of the SB216, D85, D212, and E204. The side‐chain orientation of R82 is determined primarily by the protonation states of the residues in the EC. The side‐chain reorientation of R82 modulates the hydrogen‐bond network and consequently possible pathways of proton transfer. Quantum mechanical intrinsic reaction coordinate calculations of proton‐transfer in the methyl guanidinium‐hydronium‐hydroxide model system show that proton transfer via a guanidinium group requires an initial geometry permitting proton donation and acceptance by the same amine. In all the bacteriorhodopsin models, R82 can form proton wires with both the CC and the EC connected by the same amine. Alternatively, rare proton wires for proton transfer from the CC to the EC without involving R82 were found in an O′ state where the proton on D85 is transferred to D212. Proteins 2016; 84:639–654. © 2016 Wiley Periodicals, Inc.     

Buffer catalysis of amino proton exchange in compounds of adenosine, cytidine and their endocyclic N-methylated derivatives

The use of buffer catalysts having a wide range of pK (dissociation) values (4-12) provides the first estimates of two generally useful empirical parameters of amino proton exchange in compounds of adenine and cytosine. These are a nucleobase amino group dissociation constant (pKD) and the 'encounter frequency' for proton transfer (kD), which can be used to predict amino proton exchange rates. Values of amino pKD fall in the range 8.6-9.4 for the unsubstituted nucleobases and their endocyclic N-methylated derivatives. Similar values of kD are obtained for all nucleobases (1 X 10(8) M-1 s-1). These constants were obtained from a statistical fit of second-order catalytic rate constants for amino proton exchange, measured by amino 1H-NMR lineshape at varying field frequencies (100, 300 and 360 MHz). These results confirm the requirement for buffer conjugate base formation and nucleobase protonation, but point to a different mechanism of exchange at low pH; most probably direct amino protonation for adenine, but not for cytosine compounds. Anionic buffer conjugate bases (phosphate and acetate) show a greater catalytic effect than neutral (nitrogen) bases, especially with cytosine compounds. The use of high concentrations of sodium perchlorate to sharpen amino 1H resonances of 1-methyladenosine is examined, with respect to chemical and rotational exchange and NMR line broadening.