1
|
Ancona N, Bastola A, Alexov E. PKAD-2: New entries and expansion of functionalities of the database of experimentally measured pKa's of proteins. JOURNAL OF COMPUTATIONAL BIOPHYSICS AND CHEMISTRY 2023; 22:515-524. [PMID: 37520074 PMCID: PMC10373500 DOI: 10.1142/s2737416523500230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/01/2023]
Abstract
Almost all biological reactions are pH dependent and understanding the origin of pH dependence requires knowledge of the pKa's of ionizable groups. Here we report a new edition of PKAD, the PKAD-2, which is a database of experimentally measured pKa's of proteins, both wild type and mutant proteins. The new additions include 117 wild type and 54 mutant pKa values, resulting in total 1742 experimentally measured pKa's. The new edition of PKAD-2 includes 8 new wild type and 12 new mutant proteins, resulting in total of 220 proteins. This new edition incorporates a visual 3D image of the highlighted residue of interest within the corresponding protein or protein complex. Hydrogen bonds were identified, counted, and implemented as a search feature. Other new search features include the number of neighboring residues <4A from the heaviest atom of the side chain of a given amino acid. Here, we present PKAD-2 with the intention to continuously incorporate novel features and current data with the goal to be used as benchmark for computational methods.
Collapse
Affiliation(s)
- Nicolas Ancona
- Department of Biological Sciences, College of Science, Clemson University, 105 Sikes Hall, Address, Clemson, SC 29634, United States of America
| | - Ananta Bastola
- School of Computing, College of Engineering, Computing and Applied Sciences, Clemson University, 105 Sikes Hall, SC 29634, United States of America
| | - Emil Alexov
- Department of Physics, College of Science, Clemson University, 105 Sikes Hall, Address, Clemson, SC 29634, United States of America
| |
Collapse
|
2
|
Hofer F, Fischer AL, Kamenik AS, Waibl F, Fernández-Quintero ML, Liedl KR. pH-dependent structural diversity of profilin allergens determines thermal stability. FRONTIERS IN ALLERGY 2022; 3:1007000. [PMID: 36324331 PMCID: PMC9618696 DOI: 10.3389/falgy.2022.1007000] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Accepted: 09/20/2022] [Indexed: 11/05/2022] Open
Abstract
The family of profilin allergens is a common class of proteins found in plants, viruses and various eukaryotes including mammals. Profilins are characterized by an evolutionary conserved structural fold, which is responsible for their cross-reactive nature of Immunoglobulin E (IgE) antibodies. Despite their high overall structural similarity, they exhibit substantial differences in their biophysical properties, such as thermal and pH stability. To understand the origin of these functional differences of Amb a 8, Art v 4 and Bet v 2, we performed constant pH molecular dynamics simulation in combination with Gaussian accelerated MD simulations. Depending on the respective protonation at different pH levels, we find distinct differences in conformational flexibility, which are consistent with experimentally determined melting temperatures. These variations in flexibility are accompanied by ensemble shifts in the conformational landscape and quantified and localized by residue-wise B-factors and dihedral entropies. These findings strengthen the link between flexibility of profilin allergens and their thermal stability. Thus, our results clearly show the importance of considering protonation dependent conformational ensembles in solution to elucidate biophysical differences between these structurally similar allergens.
Collapse
|
3
|
Hervø-Hansen S, Højgaard C, Johansson KE, Wang Y, Wahni K, Young D, Messens J, Teilum K, Lindorff-Larsen K, Winther JR. Charge Interactions in a Highly Charge-Depleted Protein. J Am Chem Soc 2021; 143:2500-2508. [PMID: 33529004 DOI: 10.1021/jacs.0c10789] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Electrostatic forces are important for protein folding and are favored targets of protein engineering. However, interactions between charged residues are difficult to study because of the complex network of interactions found in most proteins. We have designed a purposely simple system to investigate this problem by systematically introducing individual and pairs of charged and titratable residues in a protein otherwise free of such residues. We used constant pH molecular dynamics simulations, NMR spectroscopy, and thermodynamic double mutant cycles to probe the structure and energetics of the interaction between the charged residues. We found that the partial burial of surface charges contributes to a shift in pKa value, causing an aspartate to titrate in the neutral pH range. Additionally, the interaction between pairs of residues was found to be highly context dependent, with some pairs having no apparent preferential interaction, while other pairs would engage in coupled titration forming a highly stabilized salt bridge. We find good agreement between experiments and simulations and use the simulations to rationalize our observations and to provide a detailed mechanistic understanding of the electrostatic interactions.
Collapse
Affiliation(s)
- Stefan Hervø-Hansen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark.,Division of Theoretical Chemistry, Department of Chemistry, Lund University, SE 221 00 Lund, Sweden
| | - Casper Højgaard
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Kristoffer Enøe Johansson
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Yong Wang
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Khadija Wahni
- VIB-VUB Center for Structural Biology, Vlaams Instituut voor Biotechnologie, Vrije Universiteit Brussel, B-1050 Brussels, Belgium.,Brussels Center for Redox Biology, Vrije Universiteit Brussel, B-1050 Brussels, Belgium.,Structural Biology Brussels, Vrije Universiteit Brussel, B-1050 Brussels, Belgium
| | - David Young
- VIB-VUB Center for Structural Biology, Vlaams Instituut voor Biotechnologie, Vrije Universiteit Brussel, B-1050 Brussels, Belgium.,Brussels Center for Redox Biology, Vrije Universiteit Brussel, B-1050 Brussels, Belgium.,Structural Biology Brussels, Vrije Universiteit Brussel, B-1050 Brussels, Belgium
| | - Joris Messens
- VIB-VUB Center for Structural Biology, Vlaams Instituut voor Biotechnologie, Vrije Universiteit Brussel, B-1050 Brussels, Belgium.,Brussels Center for Redox Biology, Vrije Universiteit Brussel, B-1050 Brussels, Belgium.,Structural Biology Brussels, Vrije Universiteit Brussel, B-1050 Brussels, Belgium
| | - Kaare Teilum
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Jakob Rahr Winther
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| |
Collapse
|
4
|
Hofer F, Kamenik AS, Fernández-Quintero ML, Kraml J, Liedl KR. pH-Induced Local Unfolding of the Phl p 6 Pollen Allergen From cpH-MD. Front Mol Biosci 2021; 7:603644. [PMID: 33511157 PMCID: PMC7835895 DOI: 10.3389/fmolb.2020.603644] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 12/10/2020] [Indexed: 11/13/2022] Open
Abstract
Susceptibility to endosomal degradation is a decisive contribution to a protein's immunogenicity. It is assumed that the processing kinetics of structured proteins are inherently linked to their probability of local unfolding. In this study, we quantify the impact of endosomal acidification on the conformational stability of the major timothy grass pollen allergen Phl p 6. We use state of the art sampling approaches in combination with constant pH MD techniques to profile pH-dependent local unfolding events in atomistic detail. Integrating our findings into the current view on type 1 allergic sensitization, we characterize local protein dynamics in the context of proteolytic degradation at neutral and acidic pH for the wild type protein and point mutants with varying proteolytic stability. We analyze extensive simulation data using Markov state models and retrieve highly reliable thermodynamic and kinetic information at varying pH levels. Thereby we capture the impact of endolysosomal acidification on the structure and dynamics of the Phl p 6 mutants. We find that upon protonation at lower pH values, the conformational flexibilities in key areas of the wild type protein, i.e., T-cell epitopes and early proteolytic cleavage sites, increase significantly. A decrease of the pH even leads to local unfolding in otherwise stable secondary structure elements, which is a prerequisite for proteolytic cleavage. This effect is even more pronounced in the destabilized mutant, while no unfolding was observed for the stabilized mutant. In summary, we report detailed structural models which rationalize the experimentally observed cleavage pattern during endosomal acidification.
Collapse
|
5
|
Hofer F, Kraml J, Kahler U, Kamenik AS, Liedl KR. Catalytic Site p Ka Values of Aspartic, Cysteine, and Serine Proteases: Constant pH MD Simulations. J Chem Inf Model 2020; 60:3030-3042. [PMID: 32348143 PMCID: PMC7312390 DOI: 10.1021/acs.jcim.0c00190] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
![]()
Enzymatic function and activity of
proteases is closely controlled
by the pH value. The protonation states of titratable residues in
the active site react to changes in the pH value, according to their
pKa, and thereby determine the functionality
of the enzyme. Knowledge of the titration behavior of these residues
is crucial for the development of drugs targeting the active site
residues. However, experimental pKa data
are scarce, since the systems’ size and complexity make determination
of these pKa values inherently difficult.
In this study, we use single pH constant pH MD simulations as a fast
and robust tool to estimate the active site pKa values of a set of aspartic, cysteine, and serine proteases.
We capture characteristic pKa shifts of
the active site residues, which dictate the experimentally determined
activity profiles of the respective protease family. We find clear
differences of active site pKa values
within the respective families, which closely match the experimentally
determined pH preferences of the respective proteases. These shifts
are caused by a distinct network of electrostatic interactions characteristic
for each protease family. While we find convincing agreement with
experimental data for serine and aspartic proteases, we observe clear
deficiencies in the description of the titration behavior of cysteines
within the constant pH MD framework and highlight opportunities for
improvement. Consequently, with this work, we provide a concise set
of active site pKa values of aspartic
and serine proteases, which could serve as reference for future theoretical
as well as experimental studies.
Collapse
Affiliation(s)
- Florian Hofer
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Johannes Kraml
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Ursula Kahler
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Anna S Kamenik
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Klaus R Liedl
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| |
Collapse
|
6
|
Kamenik AS, Hofer F, Handle PH, Liedl KR. Dynamics Rationalize Proteolytic Susceptibility of the Major Birch Pollen Allergen Bet v 1. Front Mol Biosci 2020; 7:18. [PMID: 32154264 PMCID: PMC7045072 DOI: 10.3389/fmolb.2020.00018] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2019] [Accepted: 01/31/2020] [Indexed: 12/21/2022] Open
Abstract
Proteolytic susceptibility during endolysosomal degradation is decisive for allergic sensitization. In the major birch pollen allergen Bet v 1 most protease cleavage sites are located within its secondary structure elements, which are inherently inaccessible to proteases. The allergen thus must unfold locally, exposing the cleavage sites to become susceptible to proteolysis. Hence, allergen cleavage rates are presumed to be linked to their fold stability, i.e., unfolding probability. Yet, these locally unfolded structures have neither been captured in experiment nor simulation due to limitations in resolution and sampling time, respectively. Here, we perform classic and enhanced molecular dynamics (MD) simulations to quantify fold dynamics on extended timescales of Bet v 1a and two variants with higher and lower cleavage rates. Already at the nanosecond-timescale we observe a significantly higher flexibility for the destabilized variant compared to Bet v 1a and the proteolytically stabilized mutant. Estimating the thermodynamics and kinetics of local unfolding around an initial cleavage site, we find that the Bet v 1 variant with the highest cleavage rate also shows the highest probability for local unfolding. For the stabilized mutant on the other hand we only find minimal unfolding probability. These results strengthen the link between the conformational dynamics of allergen proteins and their stability during endolysosomal degradation. The presented approach further allows atomistic insights in the conformational ensemble of allergen proteins and provides probability estimates below experimental detection limits.
Collapse
Affiliation(s)
| | | | | | - Klaus R. Liedl
- Center for Molecular Biosciences Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| |
Collapse
|