Dataset Information

PKa determination of histidine residues in ?-conotoxin MII peptides by 1H NMR and constant pH molecular dynamics simulation.

ABSTRACT: ?-Conotoxin MII (?-CTxMII) is a potent and selective peptide antagonist of neuronal nicotinic acetylcholine receptors (nAChR's). Studies have shown that His9 and His12 are significant determinants of toxin binding affinity for nAChR, while Glu11 may dictate differential toxin affinity between nAChR isoforms. The protonation state of these histidine residues and therefore the charge on the ?-CTx may contribute to the observed differences in binding affinity and selectivity. In this study, we assess the pH dependence of the protonation state of His9 and His12 by (1)H NMR spectroscopy and constant pH molecular dynamics (CpHMD) in ?-CTxMII, ?-CTxMII[E11A], and the triple mutant, ?-CTxMII[N5R:E11A:H12K]. The E11A mutation does not significantly perturb the pKa of His9 or His12, while N5R:E11A:H12K results in a significant decrease in the pKa value of His9. The pKa values predicted by CpHMD simulations are in good agreement with (1)H NMR spectroscopy, with a mean absolute deviation from experiment of 0.3 pKa units. These results support the use of CpHMD as an efficient and inexpensive predictive tool to determine pKa values and structural features of small peptides critical to their function.

SUBMITTER: McDougal OM

PROVIDER: S-EPMC3632324 | biostudies-literature | 2013 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

pKa determination of histidine residues in α-conotoxin MII peptides by 1H NMR and constant pH molecular dynamics simulation.

McDougal Owen M OM Granum David M DM Swartz Mark M Rohleder Conrad C Maupin C Mark CM

The journal of physical chemistry. B 20130225 9

α-Conotoxin MII (α-CTxMII) is a potent and selective peptide antagonist of neuronal nicotinic acetylcholine receptors (nAChR's). Studies have shown that His9 and His12 are significant determinants of toxin binding affinity for nAChR, while Glu11 may dictate differential toxin affinity between nAChR isoforms. The protonation state of these histidine residues and therefore the charge on the α-CTx may contribute to the observed differences in binding affinity and selectivity. In this study, we asse ...[more]

PMID: 23336579

Similar Datasets

Project description:We present a GPU implementation of the continuous constant pH molecular dynamics (CpHMD) based on the most recent generalized Born implicit-solvent model in the pmemd engine of the Amber molecular dynamics package. To test the accuracy of the tool for rapid pKa predictions, a series of 2 ns single-pH simulations were performed for over 120 titratable residues in 10 benchmark proteins that were previously used to test the various continuous CpHMD methods. The calculated pKa's showed a root-mean-square deviation of 0.80 and correlation coefficient of 0.83 with respect to experiment. Also, 90% of the pKa's were converged with estimated errors below 0.1 pH units. Surprisingly, this level of accuracy is similar to our previous replica-exchange simulations with 2 ns per replica and an exchange attempt frequency of 2 ps-1 (Huang, Harris, and Shen J. Chem. Inf. Model. 2018 , 58 , 1372 - 1383 ). Interestingly, for the linked titration sites in two enzymes, although residue-specific protonation state sampling in the single-pH simulations was not converged within 2 ns, the protonation fraction of the linked residues appeared to be largely converged, and the experimental macroscopic pKa values were reproduced to within 1 pH unit. Comparison with replica-exchange simulations with different exchange attempt frequencies showed that the splitting between the two macroscopic pKa's is underestimated with frequent exchange attempts such as 2 ps-1, while single-pH simulations overestimate the splitting. The same trend is seen for the single-pH vs replica-exchange simulations of a hydrogen-bonded aspartyl dyad in a much larger protein. A 2 ns single-pH simulation of a 400-residue protein takes about 1 h on a single NVIDIA GeForce RTX 2080 graphics card, which is over 1000 times faster than a CpHMD run on a single CPU core of a high-performance computing cluster node. Thus, we envision that GPU-accelerated continuous CpHMD may be used in routine pKa predictions for a variety of applications, from assisting MD simulations with protonation state assignment to offering pH-dependent corrections of binding free energies and identifying reactive hot spots for covalent drug design.

Project description:The serine proteases constitute a group of endopeptidases whose members owe their catalytic activity to the presence of a catalytic triad of amino acids consisting of a serine, a histidine and an aspartate. The pK(a) values for this histidine have been determined for several cases in which there is a negative charge installed at the serine to mimic the oxyanionic intermediate and related transition state for the catalytic pathway. Instances from this laboratory include (1) replacement of the serine by a cysteine in subtilisin to create a thiolate; (2) formation of monoisopropylphosphoryl-Ser 195 monoanionic phosphodiesters (in trypsin and chymotrypsin, Ser 221 in subtilisins); and (3) tetrahedral boronates formed with peptide boronic acids. The nuclear magnetic resonance (NMR) signals pertinent to this histidine, or signals indirectly reflecting the state of ionization of this histidine, have been used effectively to monitor changes in the active center ionization state. In every case studied, there is elevation of the pK(a) at the histidine when the negative charge is installed at the serine position. Herein is reported the first NMR measurement of the active center His 63 pK(a) in thiolsubtilisin Carlsberg; it is elevated by 3 units compared with the parent enzyme. Using a numerical solution (finite difference) of the Poisson-Boltzmann equation, a protein dielectric constant of 4 provides a good estimate of the experimentally observed pK(a) elevations. Very significantly, a very low protein dielectric constant (epsilon(p) = 3-5) is required in all of the comparisons, and for all three enzymes used (chymotrypsin, trypsin, and subtilisin). Finally, we discuss why the electrostatic perturbation sensed at His of the active center is more amplified by a negative charge on the Ser side than the same charge on the Asp side. A plausible explanation is that the positive charge on the imidazolium ring of the His is localized, with the N(delta 1) carrying a smaller fraction, the N(epsilon 2) carrying the bulk of the positive charge.

Dataset Information

PKa determination of histidine residues in ?-conotoxin MII peptides by 1H NMR and constant pH molecular dynamics simulation.

Publications

pKa determination of histidine residues in α-conotoxin MII peptides by 1H NMR and constant pH molecular dynamics simulation.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets