1
|
Zanetti-Polzi L, Daidone I, Amadei A. Fully Atomistic Multiscale Approach for p Ka Prediction. J Phys Chem B 2020; 124:4712-4722. [PMID: 32427481 DOI: 10.1021/acs.jpcb.0c01752] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The ionization state of titratable amino acids strongly affects proteins structure and functioning in a large number of biological processes. It is therefore essential to be able to characterize the pKa of ionizable groups inside proteins and to understand its microscopic determinants in order to gain insights into many functional properties of proteins. A big effort has been devoted to the development of theoretical approaches for the prediction of deprotonation free energies, yet the accurate theoretical/computational calculation of pKa values is recognized as a current challenge. A methodology based on a hybrid quantum/classical approach is here proposed for the computation of deprotonation free energies. The method is applied to calculate the pKa of formic acid, methylammonium, and methanethiol, providing results in good agreement with the corresponding experimental estimates. The pKa is also calculated for aspartic acid and lysine as single residues in solution and for three aspartic/glutamic acids inside a well-characterized protein: hen egg white lysozyme. While for small molecules the method is able to deal with multiple protonation states of all titratable groups, this becomes computationally very expensive for proteins. The calculated pKa values for the single amino acids (except for the zwitterionic aspartic acid) and inside the protein display a systematic shift with respect to the experimental values that suggests that the fine balance between hydrophobic and polar interactions might be not accurately reproduced by the usual classical force-fields, thus affecting the computation of deprotonation free energies. The calculated pKa shifts inside the protein are in good agreement with the corresponding experimental ones (within 1 pKa unit), well reproducing the pKa changes due to the protein environment even in the case of large pKa shifts.
Collapse
Affiliation(s)
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, Via Vetoio, I-67010 L'Aquila, Italy
| | - Andrea Amadei
- Department of Chemical and Technological Sciences, University of Rome "Tor Vergata", Via della Ricerca Scientifica, I-00185 Rome, Italy
| |
Collapse
|
2
|
|
3
|
Vorobjev YN. An effective molecular blocker of ion channel of M2 protein as anti-influenza a drug. J Biomol Struct Dyn 2020; 39:2352-2363. [PMID: 32212957 DOI: 10.1080/07391102.2020.1747550] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Design of a drug compound that can effectively bind to the M2 ion channel and block the diffusion of hydrogen ions (H+) through and inhibit influenza A virus replication is an important task. Known anti-influenza drugs amantadine and rimantadine have a weak effect on influenza A virus. A new class of positively charged, +2 e.u., molecules is proposed here to block diffusion of H+ ion through the M2 channel. Several drug candidates, derivatives of a lead compound (diazabicyclooctane), were proposed and investigated. Molecular dynamics of thermal fluctuations of M2 protein structure and ionization-conformation coupling of all the ionizable residues were simulated at physiological pH. The influence of the most probable mutations of key drug-binding amino acid residues in the M2 ion channel were investigated too. It is shown that the suggested new blocker has high binding affinity for the M2 channel. There are two in-channel binding sites of high affinity, the first one has H-bonds with two of four serine residues Ser-31A (B) or Ser-31C(D), and the second one has H-bonds with two of four histidine residues His-37A (B), or His-37C(D). The main advantage of the new drug molecule is the positive charge, +2 e.u., which creates a positive electrostatic potential barrier (in addition to a steric one) for a transfer of H+ ion through M2 channel and may serve as an effective anti-influenza A virus drug.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Yury N Vorobjev
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences, Novosibirsk, Russia
| |
Collapse
|
4
|
Han CT, Song J, Chan T, Pruett C, Han S. Electrostatic Environment of Proteorhodopsin Affects the pKa of Its Buried Primary Proton Acceptor. Biophys J 2020; 118:1838-1849. [PMID: 32197061 DOI: 10.1016/j.bpj.2020.02.027] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Revised: 01/27/2020] [Accepted: 02/27/2020] [Indexed: 01/18/2023] Open
Abstract
The protonation state of embedded charged residues in transmembrane proteins (TMPs) can control the onset of protein function. It is understood that interactions between an embedded charged residue and other charged or polar residues in the moiety would influence its pKa, but how the surrounding environment in which the TMP resides affects the pKa of these residues is unclear. Proteorhodopsin (PR), a light-responsive proton pump from marine bacteria, was used as a model to examine externally accessible factors that tune the pKa of its embedded charged residue, specifically its primary proton acceptor D97. The pKa of D97 was compared between PR reconstituted in liposomes with different net headgroup charges and equilibrated in buffer with different ion concentrations. For PR reconstituted in net positively charged compared to net negatively charged liposomes in low-salt buffer solutions, a drop of the apparent pKa from 7.6 to 5.6 was observed, whereas intrinsic pKa modeled with surface pH calculated from Gouy-Chapman predictions found an opposite trend for the pKa change, suggesting that surface pH does not account for the main changes observed in the apparent pKa. This difference in the pKa of D97 observed from PR reconstituted in oppositely charged liposome environments disappeared when the NaCl concentration was increased to 150 mM. We suggest that protein-intrinsic structural properties must play a role in adjusting the local microenvironment around D97 to affect its pKa, as corroborated with observations of changes in protein side-chain and hydration dynamics around the E-F loop of PR. Understanding the effect of externally controllable factors in tuning the pKa of TMP-embedded charged residues is important for bioengineering and biomedical applications relying on TMP systems, in which the onset of functions can be controlled by the protonation state of embedded residues.
Collapse
Affiliation(s)
- Chung-Ta Han
- Department of Chemical Engineering, University of California, Santa Barbara, California
| | - Jichao Song
- Department of Chemical Engineering, University of California, Santa Barbara, California
| | - Tristan Chan
- Department of Chemistry, University of California, Santa Barbara, California
| | - Christine Pruett
- Department of Chemical Engineering, University of California, Santa Barbara, California
| | - Songi Han
- Department of Chemical Engineering, University of California, Santa Barbara, California; Department of Chemistry, University of California, Santa Barbara, California.
| |
Collapse
|
5
|
Alekseeva IV, Bakman AS, Vorobjev YN, Fedorova OS, Kuznetsov NA. Role of Ionizing Amino Acid Residues in the Process of DNA Binding by Human AP Endonuclease 1 and in Its Catalysis. J Phys Chem B 2019; 123:9546-9556. [PMID: 31633353 DOI: 10.1021/acs.jpcb.9b07150] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In the repair of the damage to bases, human apurinic/apyrimidinic (AP) endonuclease 1 (APE1) is a key participant via the DNA base excision repair pathway. APE1 cleaves AP sites in DNA, which are potentially cytotoxic and highly mutagenic if left unrepaired. According to existing structural data, this enzyme's active site contains many polar amino acid residues, which form extensive contacts with a DNA substrate. A few alternative catalytic mechanisms of the phosphodiester bond hydrolysis by APE1 have been reported. Here, the kinetics of conformational changes of the enzyme and of DNA substrate molecules were studied during the recognition and cleavage of the abasic site in the pH range from 5.5 to 9.0 using stopped-flow fluorescence techniques. The activity of APE1 increased with an increase in pH because of acceleration of the rates of catalytic complex formation and of the catalytic reaction. Molecular dynamics simulation uncovered a significant increase in the pKa of His-309 located in the active site of the enzyme. This finding revealed that the observed enhancement of enzymatic activity with pH could be associated with deprotonation of not only Tyr-171 but also His-309. The obtained data allowed us to hypothesize that the ionized state of these residues could be a molecular switch between the alternative catalytic mechanisms, which involve different functionalities of these residues throughout the reaction.
Collapse
Affiliation(s)
- Irina V Alekseeva
- Institute of Chemical Biology and Fundamental Medicine , Novosibirsk 630090 , Russia
| | - Artemiy S Bakman
- Institute of Chemical Biology and Fundamental Medicine , Novosibirsk 630090 , Russia
| | - Yury N Vorobjev
- Institute of Chemical Biology and Fundamental Medicine , Novosibirsk 630090 , Russia
| | - Olga S Fedorova
- Institute of Chemical Biology and Fundamental Medicine , Novosibirsk 630090 , Russia.,Department of Natural Sciences , Novosibirsk State University , Novosibirsk 630090 , Russia
| | - Nikita A Kuznetsov
- Institute of Chemical Biology and Fundamental Medicine , Novosibirsk 630090 , Russia.,Department of Natural Sciences , Novosibirsk State University , Novosibirsk 630090 , Russia
| |
Collapse
|
6
|
Martin OA, Vorobjev Y, Scheraga HA, Vila JA. Outline of an experimental design aimed to detect a protein A mirror image in solution. PEERJ PHYSICAL CHEMISTRY 2019; 1. [PMID: 34079958 DOI: 10.7717/peerj-pchem.2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
There is abundant theoretical evidence indicating that a mirror image of Protein A may occur during the protein folding process. However, as to whether such mirror image exists in solution is an unsolved issue. Here we provide outline of an experimental design aimed to detect the mirror image of Protein A in solution. The proposal is based on computational simulations indicating that the use of a mutant of protein A, namely Q10H, could be used to detect the mirror image conformation in solution. Our results indicate that the native conformation of the protein A should have a pKa, for the Q10H mutant, at ≈6.2, while the mirror-image conformation should have a pKa close to ≈7.3. Naturally, if all the population is in the native state for the Q10H mutant, the pKa should be ≈6.2, while, if all are in the mirror-image state, it would be ≈7.3, and, if it is a mixture, the pKa should be largerthan 6.2, presumably in proportion to the mirror population. In addition, evidence is provided indicating the tautomeric distribution of H10 must also change between the native and mirror conformations. Although this may not be completely relevant for the purpose of determining whether the protein A mirror image exists in solution, it could provide valuable information to validate the pKa findings. We hope this proposal will foster experimental work on this problem either by direct application of our proposed experimental design or serving as inspiration and motivation for other experiments.
Collapse
Affiliation(s)
- Osvaldo A Martin
- Instituto de Matemática Aplicada San Luis, UNSL-CONICET, San Luis, Argentina
| | - Yury Vorobjev
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Science, Novosibirsk, Russia
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, United States
| | - Jorge A Vila
- Instituto de Matemática Aplicada San Luis, UNSL-CONICET, San Luis, Argentina.,Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, United States
| |
Collapse
|
7
|
Sakipov SN, Flores-Canales JC, Kurnikova MG. A Hierarchical Approach to Predict Conformation-Dependent Histidine Protonation States in Stable and Flexible Proteins. J Phys Chem B 2019; 123:5024-5034. [PMID: 31095377 DOI: 10.1021/acs.jpcb.9b00656] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Solution acidity measured by pH is an important environmental factor that affects protein structure. It influences the protonation state of protein residues, which in turn may be coupled to protein conformational changes, unfolding, and ligand binding. It remains difficult to compute and measure the p Ka of individual residues, as well as to relate them to pH-dependent protein transitions. This paper presents a hierarchical approach to compute the p Ka of individual protonatable residues, specifically histidines, coupled with underlying structural changes of a protein. A fast and efficient free energy perturbation (FEP) algorithm has also been developed utilizing a fast implementation of standard molecular dynamics (MD) algorithms. Specifically, a CUDA version of the AMBER MD engine is used in this paper. Eight histidine p Ka's are computed in a diverse set of pH stable proteins to demonstrate the proposed approach's utility and assess the predictive quality of the AMBER FF99SB force field. A reference molecule is carefully selected and tested for convergence. A hierarchical approach is used to model p Ka's of the six histidine residues of the diphtheria toxin translocation domain (DTT), which exhibits a diverse ensemble of individual conformations and pH-dependent unfolding. The hierarchical approach consists of first sampling equilibrium conformational ensembles of a protein with protonated and neutral histidine residues via long equilibrium MD simulations (Flores-Canales, J. C.; et al. bioRxiv, 2019, 572040). A clustering method is then used to identify sampled protein conformations, and p Ka's of histidines in each protein conformation are computed. Finally, an ensemble averaging formalism is developed to compute weighted average histidine p Ka's. These can be compared with an apparent experimentally measured p Ka of the DTT protein and thus allows us to propose a mechanism of pH-dependent unfolding of the DTT protein.
Collapse
Affiliation(s)
- Serzhan N Sakipov
- Chemistry Department , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Jose C Flores-Canales
- Chemistry Department , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Maria G Kurnikova
- Chemistry Department , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| |
Collapse
|
8
|
Vorobjev YN, Scheraga HA, Vila JA. A comprehensive analysis of the computed tautomer fractions of the imidazole ring of histidines in Loligo vulgaris. J Biomol Struct Dyn 2018; 36:3094-3105. [PMID: 28884632 PMCID: PMC6102090 DOI: 10.1080/07391102.2017.1377636] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2017] [Accepted: 09/04/2017] [Indexed: 10/18/2022]
Abstract
A recently introduced electrostatic-based method to determine the pKa values of ionizable residues and fractions of ionized and tautomeric forms of histidine (His) and acid residues in proteins, at a given fixed pH, is applied here to the analysis of a His-rich protein, namely Loligo vulgaris (pdb id 1E1A), a 314-residue all-β protein. The average tautomeric fractions for the imidazole ring of each of the six histidines in the sequence were computed using an approach that includes, but is not limited to, molecular dynamic simulations coupled with calculations of the ionization states for all 94 ionizable residues of protein 1E1A in water at pH 6.5 and 300 K. The electrostatic-calculated tautomeric fractions of the imidazole ring of His were compared with predictions obtained from an existent NMR-based methodology. Our results indicate that: (i) the averaged electrostatic-based tautomeric predictions for the imidazole ring of all histidines of Loligo vulgaris are dominated by the Nε2-H rather than the Nδ1-H form, although such preferences from the NMR-based methodology are not so well defined; (ii) the computed average absolute difference between the electrostatic- and the NMR-based tautomeric predictions among all six histidines vary among 0% to 17%; (iii) for the His showing the largest fraction of the neutral form (81%), the absolute difference between the NMR- and electrostatic-based computed tautomeric predictions is only 3%; and (iv) the tautomeric predictions for the imidazole ring of His computed with the NMR-based methodology are stable within a certain, well-defined, range of variations of a tautomer-related parameter.
Collapse
Affiliation(s)
- Yury N Vorobjev
- a Institute of Chemical Biology and Fundamental Medicine of the Siberian Branch of the Russian Academy of Science , Lavrentiev Avenue 8, Novosibirsk 630090 , Russia
- b Department of Chemistry , Novosibirsk State University , Novosibirsk 630090 , Russia
- c Baker Laboratory of Chemistry and Chemical Biology , Cornell University , Ithaca , NY 14853-1301 , USA
| | - Harold A Scheraga
- c Baker Laboratory of Chemistry and Chemical Biology , Cornell University , Ithaca , NY 14853-1301 , USA
| | - Jorge A Vila
- d IMASL-CONICET , Universidad Nacional de San Luis , Ejército de Los Andes 950, 5700 San Luis , Argentina
| |
Collapse
|