1
|
Radhakrishnan K, Singh SP. Explicit characterization of counterion dynamics around a flexible polyelectrolyte. Phys Rev E 2022; 105:044501. [PMID: 35590562 DOI: 10.1103/physreve.105.044501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Accepted: 03/04/2022] [Indexed: 06/15/2023]
Abstract
The article presents a comprehensive study of counterion dynamics around a generic linear polyelectrolyte chain with the help of coarse-grained computer simulations. The ion-chain coupling is discussed in the form of binding time, mean-square displacement (MSD) relative to the chain, local ion transport coefficient, and spatiotemporal correlations in the effective charge. We have shown that a counterion exhibits subdiffusive behavior 〈δR^{2}〉∼t^{δ}, δ≈0.9 w.r.t. chain's center of mass. The MSD of ions perpendicularly outward from the chain segment exhibits a smaller subdiffusive exponent compared to the one along the chain backbone. The effective diffusivity of ion is the lowest in chain's close proximity, extending up to the length-scale of radius of gyration R_{g}. Beyond R_{g} at larger distances, they attain diffusivity of free ion with a smooth cross-over from the adsorbed regime to the free ion regime. We have shown that the effective diffusivity drastically decreases for the multivalent ions, while the crossover length scale remains the same. Conversely, with increasing salt concentration the coupling-length scale reduces, while the diffusivity remains unaltered. The effective diffusivity of adsorbed-ion reveals an exponential reduction with electrostatic interaction strength. We further corroborate this from the binding time of ions on the chain, which also grows exponentially with the coupling strength of the ion-polymer duo. Moreover, the binding time of ions exhibits a weak dependence with salt concentration for the monovalent salt, while for multivalent salts the binding time decreases dramatically with concentration. Our work also elucidates fluctuations in the effective charge per site, where it exhibits strong negative correlations at short length-scales.
Collapse
Affiliation(s)
- Keerthi Radhakrishnan
- Department of Physics, Indian Institute of Science Education and Research, Bhopal 462 066, Madhya Pradesh, India
| | - Sunil P Singh
- Department of Physics, Indian Institute of Science Education and Research, Bhopal 462 066, Madhya Pradesh, India
| |
Collapse
|
2
|
da Rocha L, Baptista AM, Campos SRR. Approach to Study pH-Dependent Protein Association Using Constant-pH Molecular Dynamics: Application to the Dimerization of β-Lactoglobulin. J Chem Theory Comput 2022; 18:1982-2001. [PMID: 35171602 PMCID: PMC9775224 DOI: 10.1021/acs.jctc.1c01187] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Protein-protein association is often mediated by electrostatic interactions and modulated by pH. However, experimental and computational studies have often overlooked the effect of association on the protonation state of the protein. In this work, we present a methodological approach based on constant-pH molecular dynamics (MD), which aims to provide a detailed description of a pH-dependent protein-protein association, and apply it to the dimerization of β-lactoglobulin (BLG). A selection of analyses is performed using the data generated by constant-pH MD simulations of monomeric and dimeric forms of bovine BLG, in the pH range 3-8. First, we estimate free energies of dimerization using a computationally inexpensive approach based on the Wyman-Tanford linkage theory, calculated in a new way through the use of thermodynamically based splines. The individual free energy contribution of each titratable site is also calculated, allowing for identification of relevant residues. Second, the correlations between the proton occupancies of pairs of sites are calculated (using the Pearson coefficient), and extensive networks of correlated sites are observed at acidic pH values, sometimes involving distant pairs. In general, strongly correlated sites are also slow proton exchangers and contribute significantly to the pH-dependency of the dimerization free energy. Third, we use ionic density as a fingerprint of protein charge distribution and observe electrostatic complementarity between the monomer faces that form the dimer interface, more markedly at the isoionic point (where maximum dimerization occurs) than at other pH values, which might contribute to guide the association. Finally, the pH-dependent dimerization modes are inspected using PCA, among other analyses, and two states are identified: a relaxed state at pH 4-8 (with the typical alignment of the crystallographic structure) and a compact state at pH 3-4 (with a tighter association and rotated alignment). This work shows that an approach based on constant-pH MD simulations can produce rich detailed pictures of pH-dependent protein associations, as illustrated for BLG dimerization.
Collapse
|
3
|
Antosiewicz JM, Długosz M. Constant-pH Brownian Dynamics Simulations of a Protein near a Charged Surface. ACS OMEGA 2020; 5:30282-30298. [PMID: 33251463 PMCID: PMC7689933 DOI: 10.1021/acsomega.0c04817] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2020] [Accepted: 11/04/2020] [Indexed: 05/04/2023]
Abstract
We have developed a rigid-body Brownian dynamics algorithm that allows for simulations of a globular protein suspended in an ionic solution confined by a charged planar boundary, with an explicit treatment of pH-dependent protein protonation equilibria and their couplings to the electrostatic potential of the plane. Electrostatic interactions are described within a framework of the continuum Poisson-Boltzmann model, whereas protein-plane hydrodynamic interactions are evaluated based on analytical expressions for the position- and orientation-dependent near-wall friction tensor of a spheroid. The algorithm was applied to simulate near-surface diffusion of lysozyme in solutions having pH in the range 4-10 and ionic strengths of 10 and 150 mM. As a reference, we performed Brownian dynamics simulations in which the protein is assigned a fixed, most probable protonation state, appropriate for given solution conditions and unaffected by the presence of the charged plane, and Brownian dynamics simulations in which the protein probes possible protonation states with the pH-dependent probability, but these variations are not coupled to the electric field generated by the boundary. We show that electrostatic interactions with the negatively charged plane substantially modify probabilities of different protonation states of lysozyme and shift protonation equilibria of both acidic and basic amino acid side chains toward higher pH values. Consequently, equilibrium energy distributions, equilibrium position-orientation distributions, and functions that characterize rotational dynamics, which for a protein with multiple ionization sites, such as lysozyme, in the presence of a charged obstacle are pH-dependent, are significantly affected by the approach taken to incorporate the solution pH into simulations.
Collapse
|
4
|
Kumawat A, Chakrabarty S. Protonation-Induced Dynamic Allostery in PDZ Domain: Evidence of Perturbation-Independent Universal Response Network. J Phys Chem Lett 2020; 11:9026-9031. [PMID: 33043672 DOI: 10.1021/acs.jpclett.0c02885] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Dynamic allostery is a relatively new paradigm where certain external perturbations may lead to modulation of conformational dynamics at a distant part of a protein without significant changes in the overall structure. While most well-characterized examples of dynamic allostery involve binding with other entities like small molecules, peptides, or nucleic acids, in this work we demonstrate that chemical modifications like protonation may lead to significant dynamical allosteric response in a PDZ domain protein. Tuning the protonation states of two histidine residues (H317 and H372), we identify the allosteric pathways responsible for the dynamic response. Interestingly, the same set of residues that constitute the allosteric response network upon ligand binding seem to be responsible for protonation-induced dynamic allostery. Thus, we propose the existence of an inherent universal response network in signaling proteins, where the same set of residues can respond to varying types of external perturbations in terms of rearrangement of hydrogen-bonded network and redistribution of electrostatic interaction energies.
Collapse
Affiliation(s)
- Amit Kumawat
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
- Physical and Materials Chemistry Division, CSIR-National Chemical Laboratory, Pune 411008, India
| | - Suman Chakrabarty
- Department of Chemical, Biological & Macromolecular Sciences, S. N. Bose National Centre for Basic Sciences, Kolkata 700106, India
| |
Collapse
|
5
|
Bakhshandeh A, Frydel D, Levin Y. Charge regulation of colloidal particles in aqueous solutions. Phys Chem Chem Phys 2020; 22:24712-24728. [PMID: 33104140 DOI: 10.1039/d0cp03633a] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We study the charge regulation of colloidal particles inside aqueous electrolyte solutions. To stabilize a colloidal suspension against precipitation, colloidal particles are synthesized with either acidic or basic groups on their surface. On contact with water, these surface groups undergo proton transfer reactions, resulting in colloidal surface charge. The charge is determined by the condition of local chemical equilibrium between hydronium ions inside the solution and at the colloidal surface. We use a model of Baxter sticky spheres to explicitly calculate the equilibrium dissociation constants and to construct a theory which is able to quantitatively predict the effective charge of colloidal particles with either acidic or basic surface groups. The predictions of the theory for the model are found to be in excellent agreement with the results of Monte Carlo simulations. This theory is further extended to treat colloidal particles with a mixture of both acidic and basic surface groups.
Collapse
Affiliation(s)
- Amin Bakhshandeh
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil.
| | | | | |
Collapse
|
6
|
Kawade R, Kuroda D, Tsumoto K. How the protonation state of a phosphorylated amino acid governs molecular recognition: insights from classical molecular dynamics simulations. FEBS Lett 2019; 594:903-912. [DOI: 10.1002/1873-3468.13674] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2019] [Revised: 10/07/2019] [Accepted: 11/06/2019] [Indexed: 11/05/2022]
Affiliation(s)
- Raiji Kawade
- Department of Bioengineering School of Engineering The University of Tokyo Japan
| | - Daisuke Kuroda
- Department of Bioengineering School of Engineering The University of Tokyo Japan
- Medical Device Development and Regulation Research Center School of Engineering The University of Tokyo Japan
| | - Kouhei Tsumoto
- Department of Bioengineering School of Engineering The University of Tokyo Japan
- Medical Device Development and Regulation Research Center School of Engineering The University of Tokyo Japan
- Laboratory of Medical Proteomics The Institute of Medical Science The University of Tokyo Japan
| |
Collapse
|
7
|
Blanco PM, Madurga S, Narambuena CF, Mas F, Garcés JL. Role of Charge Regulation and Fluctuations in the Conformational and Mechanical Properties of Weak Flexible Polyelectrolytes. Polymers (Basel) 2019; 11:polym11121962. [PMID: 31795443 PMCID: PMC6960815 DOI: 10.3390/polym11121962] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2019] [Revised: 11/21/2019] [Accepted: 11/25/2019] [Indexed: 01/29/2023] Open
Abstract
This work addresses the role of charge regulation (CR) and the associated fluctuations in the conformational and mechanical properties of weak polyelectrolytes. Due to CR, changes in the pH-value modifies the average macromolecular charge and conformational equilibria. A second effect is that, for a given average charge per site, fluctuations can alter the intensity of the interactions by means of correlation between binding sites. We investigate both effects by means of Monte Carlo simulations at constant pH-value, so that the charge is a fluctuating quantity. Once the average charge per site is available, we turn off the fluctuations by assigning the same average charge to every site. A constant charge MC simulation is then performed. We make use of a model which accounts for the main fundamental aspects of a linear flexible polyelectrolyte that is, proton binding, angle internal rotation, bond stretching and bending. Steric excluded volume and differentiated treatment for short-range and long-range interactions are also included. This model can be regarded as a kind of “minimal” in the sense that it contains a minimum number of parameters but still preserving the atomistic detail. It is shown that, if fluctuations are activated, gauche state bond probabilities increase and the persistence length decreases, so that the polymer becomes more folded. Macromolecular stretching is also analyzed in presence of CR (the charge depends on the applied force) and without CR (the charge is fixed to the value at zero force). The analysis of the low force scaling behavior concludes that Pincus exponent becomes pH-dependent. Both, with and without CR, a transition from 1/2 at high pH-values (phantom chain) to 3/5 at low pH-values (Pincus regime) is observed. Finally, the intermediate force stretching regime is investigated. It is found that CR induces a moderate influence in the force-extension curves and persistence length (which in this force regime becomes force-dependent). It is thus concluded that the effect of CR on the stretching curves is mainly due to the changes in the average charge at zero force. It is also found that, for the cases studied, the effect of steric excluded volume is almost irrelevant compared to electrostatic interactions.
Collapse
Affiliation(s)
- Pablo M. Blanco
- Physical Chemistry Unit, Materials Science and Physical Chemistry Department & Research Institute of Theoretical and Computational Chemistry (IQTCUB) of Barcelona University (UB), 08028 Barcelona, Catalonia, Spain;
- Correspondence: (P.M.B.); (F.M.)
| | - Sergio Madurga
- Physical Chemistry Unit, Materials Science and Physical Chemistry Department & Research Institute of Theoretical and Computational Chemistry (IQTCUB) of Barcelona University (UB), 08028 Barcelona, Catalonia, Spain;
| | - Claudio F. Narambuena
- Facultad Regional San Rafael, Universidad Tecnológica Nacional & Instituto de Física Aplicada (INFAP), Universidad Nacional de San Luis-CONICET, 5600 San Rafael, Argentina;
| | - Francesc Mas
- Physical Chemistry Unit, Materials Science and Physical Chemistry Department & Research Institute of Theoretical and Computational Chemistry (IQTCUB) of Barcelona University (UB), 08028 Barcelona, Catalonia, Spain;
- Correspondence: (P.M.B.); (F.M.)
| | - Josep L. Garcés
- Chemistry Department, Technical School of Agricultural Engineering & AGROTECNIO of Lleida University (UdL), 25003 Lleida, Catalonia, Spain;
| |
Collapse
|
8
|
Zhou HX, Pang X. Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation. Chem Rev 2018; 118:1691-1741. [PMID: 29319301 DOI: 10.1021/acs.chemrev.7b00305] [Citation(s) in RCA: 454] [Impact Index Per Article: 75.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.
Collapse
Affiliation(s)
- Huan-Xiang Zhou
- Department of Chemistry and Department of Physics, University of Illinois at Chicago , Chicago, Illinois 60607, United States.,Department of Physics and Institute of Molecular Biophysics, Florida State University , Tallahassee, Florida 32306, United States
| | - Xiaodong Pang
- Department of Physics and Institute of Molecular Biophysics, Florida State University , Tallahassee, Florida 32306, United States
| |
Collapse
|
9
|
Wahle CW, Martini KM, Hollenbeck DM, Langner A, Ross DS, Hamilton JF, Thurston GM. Model for screened, charge-regulated electrostatics of an eye lens protein: Bovine gammaB-crystallin. Phys Rev E 2017; 96:032415. [PMID: 29346981 PMCID: PMC5830141 DOI: 10.1103/physreve.96.032415] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2016] [Indexed: 06/07/2023]
Abstract
We model screened, site-specific charge regulation of the eye lens protein bovine gammaB-crystallin (γB) and study the probability distributions of its proton occupancy patterns. Using a simplified dielectric model, we solve the linearized Poisson-Boltzmann equation to calculate a 54×54 work-of-charging matrix, each entry being the modeled voltage at a given titratable site, due to an elementary charge at another site. The matrix quantifies interactions within patches of sites, including γB charge pairs. We model intrinsic pK values that would occur hypothetically in the absence of other charges, with use of experimental data on the dependence of pK values on aqueous solution conditions, the dielectric model, and literature values. We use Monte Carlo simulations to calculate a model grand-canonical partition function that incorporates both the work-of-charging and the intrinsic pK values for isolated γB molecules and we calculate the probabilities of leading proton occupancy configurations, for 4<pH<8 and Debye screening lengths from 6 to 20 Å. We select the interior dielectric value to model γB titration data. At pH 7.1 and Debye length 6.0 Å, on a given γB molecule the predicted top occupancy pattern is present nearly 20% of the time, and 90% of the time one or another of the first 100 patterns will be present. Many of these occupancy patterns differ in net charge sign as well as in surface voltage profile. We illustrate how charge pattern probabilities deviate from the multinomial distribution that would result from use of effective pK values alone and estimate the extents to which γB charge pattern distributions broaden at lower pH and narrow as ionic strength is lowered. These results suggest that for accurate modeling of orientation-dependent γB-γB interactions, consideration of numerous pairs of proton occupancy patterns will be needed.
Collapse
Affiliation(s)
- Christopher W. Wahle
- School of Mathematical Sciences, Rochester Institute of Technology, Rochester, New York 14623, USA
| | - K. Michael Martini
- School of Physics and Astronomy, Rochester Institute of Technology, Rochester, New York 14623, USA
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana-Champaign, Illinois 61801, USA
| | - Dawn M. Hollenbeck
- School of Physics and Astronomy, Rochester Institute of Technology, Rochester, New York 14623, USA
| | - Andreas Langner
- School of Chemistry and Materials Science, Rochester Institute of Technology, Rochester, New York 14623, USA
| | - David S. Ross
- School of Mathematical Sciences, Rochester Institute of Technology, Rochester, New York 14623, USA
| | - John F. Hamilton
- School of Mathematical Sciences, Rochester Institute of Technology, Rochester, New York 14623, USA
| | - George M. Thurston
- School of Physics and Astronomy, Rochester Institute of Technology, Rochester, New York 14623, USA
| |
Collapse
|
10
|
Kim MO, McCammon JA. Computation of pH-dependent binding free energies. Biopolymers 2016; 105:43-9. [PMID: 26202905 PMCID: PMC4623928 DOI: 10.1002/bip.22702] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2015] [Accepted: 07/20/2015] [Indexed: 01/21/2023]
Abstract
Protein-ligand binding accompanies changes in the surrounding electrostatic environments of the two binding partners and may lead to changes in protonation upon binding. In cases where the complex formation results in a net transfer of protons, the binding process is pH-dependent. However, conventional free energy computations or molecular docking protocols typically employ fixed protonation states for the titratable groups in both binding partners set a priori, which are identical for the free and bound states. In this review, we draw attention to these important yet largely ignored binding-induced protonation changes in protein-ligand association by outlining physical origins and prevalence of the protonation changes upon binding. Following a summary of various theoretical methods for pKa prediction, we discuss the theoretical framework to examine the pH dependence of protein-ligand binding processes.
Collapse
Affiliation(s)
- M. Olivia Kim
- Department of Pharmacology, University of California San Diego, La Jolla, CA 92093, USA
| | - J. Andrew McCammon
- Department of Pharmacology, University of California San Diego, La Jolla, CA 92093, USA
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
- Howard Hughes Medical Institute, University of California San Diego, La Jolla, CA 92093, USA
- National Biomedical Computation Resource, University of California San Diego, La Jolla, CA 92093, USA
| |
Collapse
|
11
|
Conformational Dynamics and Binding Free Energies of Inhibitors of BACE-1: From the Perspective of Protonation Equilibria. PLoS Comput Biol 2015; 11:e1004341. [PMID: 26506513 PMCID: PMC4623973 DOI: 10.1371/journal.pcbi.1004341] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2015] [Accepted: 05/17/2015] [Indexed: 11/19/2022] Open
Abstract
BACE-1 is the β-secretase responsible for the initial amyloidogenesis in Alzheimer’s disease, catalyzing hydrolytic cleavage of substrate in a pH-sensitive manner. The catalytic mechanism of BACE-1 requires water-mediated proton transfer from aspartyl dyad to the substrate, as well as structural flexibility in the flap region. Thus, the coupling of protonation and conformational equilibria is essential to a full in silico characterization of BACE-1. In this work, we perform constant pH replica exchange molecular dynamics simulations on both apo BACE-1 and five BACE-1-inhibitor complexes to examine the effect of pH on dynamics and inhibitor binding properties of BACE-1. In our simulations, we find that solution pH controls the conformational flexibility of apo BACE-1, whereas bound inhibitors largely limit the motions of the holo enzyme at all levels of pH. The microscopic pKa values of titratable residues in BACE-1 including its aspartyl dyad are computed and compared between apo and inhibitor-bound states. Changes in protonation between the apo and holo forms suggest a thermodynamic linkage between binding of inhibitors and protons localized at the dyad. Utilizing our recently developed computational protocol applying the binding polynomial formalism to the constant pH molecular dynamics (CpHMD) framework, we are able to obtain the pH-dependent binding free energy profiles for various BACE-1-inhibitor complexes. Our results highlight the importance of correctly addressing the binding-induced protonation changes in protein-ligand systems where binding accompanies a net proton transfer. This work comprises the first application of our CpHMD-based free energy computational method to protein-ligand complexes and illustrates the value of CpHMD as an all-purpose tool for obtaining pH-dependent dynamics and binding free energies of biological systems. Formation of insoluble amyloid plaques in the vascular and hippocampal areas of the brain characterizes Alzheimer’s disease, a devastating neurodegenerative disorder causing dementia. Site-specific hydrolytic catalysis of β-secretase, or BACE-1, is responsible for production of oligomerative amyloid β-peptide. As the catalytic activity of BACE-1 is pH-dependent and its structural dynamics are intrinsic to the catalysis, we examine the dependence of dynamics of BACE-1 on solution pH and its implications on the catalytic mechanism of BACE-1. Also, we highlight the importance of accurate description of protonation states of the titratable groups in computer-aided drug discovery targeting BACE-1. We hope the understanding of pH dependence of the dynamics and inhibitor binding properties of BACE-1 will aid the structure-based inhibitor design efforts against Alzheimer’s disease.
Collapse
|
12
|
Lee J, Miller BT, Brooks BR. Computational scheme for pH-dependent binding free energy calculation with explicit solvent. Protein Sci 2015; 25:231-43. [PMID: 26189656 PMCID: PMC4815317 DOI: 10.1002/pro.2755] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2015] [Accepted: 07/16/2015] [Indexed: 01/25/2023]
Abstract
We present a computational scheme to compute the pH‐dependence of binding free energy with explicit solvent. Despite the importance of pH, the effect of pH has been generally neglected in binding free energy calculations because of a lack of accurate methods to model it. To address this limitation, we use a constant‐pH methodology to obtain a true ensemble of multiple protonation states of a titratable system at a given pH and analyze the ensemble using the Bennett acceptance ratio (BAR) method. The constant pH method is based on the combination of enveloping distribution sampling (EDS) with the Hamiltonian replica exchange method (HREM), which yields an accurate semi‐grand canonical ensemble of a titratable system. By considering the free energy change of constraining multiple protonation states to a single state or releasing a single protonation state to multiple states, the pH dependent binding free energy profile can be obtained. We perform benchmark simulations of a host‐guest system: cucurbit[7]uril (CB[7]) and benzimidazole (BZ). BZ experiences a large pKa shift upon complex formation. The pH‐dependent binding free energy profiles of the benchmark system are obtained with three different long‐range interaction calculation schemes: a cutoff, the particle mesh Ewald (PME), and the isotropic periodic sum (IPS) method. Our scheme captures the pH‐dependent behavior of binding free energy successfully. Absolute binding free energy values obtained with the PME and IPS methods are consistent, while cutoff method results are off by 2 kcal mol−1. We also discuss the characteristics of three long‐range interaction calculation methods for constant‐pH simulations.
Collapse
Affiliation(s)
- Juyong Lee
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute (NHLBI), National Institutes of Health (NIH), Bethesda, Maryland, 20892
| | - Benjamin T Miller
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute (NHLBI), National Institutes of Health (NIH), Bethesda, Maryland, 20892
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute (NHLBI), National Institutes of Health (NIH), Bethesda, Maryland, 20892
| |
Collapse
|
13
|
Kim MO, Blachly PG, Kaus JW, McCammon JA. Protocols utilizing constant pH molecular dynamics to compute pH-dependent binding free energies. J Phys Chem B 2014; 119:861-72. [PMID: 25134690 PMCID: PMC4306499 DOI: 10.1021/jp505777n] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
![]()
In protein–ligand binding,
the electrostatic environments
of the two binding partners may vary significantly in bound and unbound
states, which may lead to protonation changes upon binding. In cases
where ligand binding results in a net uptake or release of protons,
the free energy of binding is pH-dependent. Nevertheless, conventional
free energy calculations and molecular docking protocols typically
do not rigorously account for changes in protonation that may occur
upon ligand binding. To address these shortcomings, we present a simple
methodology based on Wyman’s binding polynomial formalism to
account for the pH dependence of binding free energies and demonstrate
its use on cucurbit[7]uril (CB[7]) host–guest systems. Using
constant pH molecular dynamics and a reference binding free energy
that is taken either from experiment or from thermodynamic integration
computations, the pH-dependent binding free energy is determined.
This computational protocol accurately captures the large pKa shifts observed experimentally upon CB[7]:guest
association and reproduces experimental binding free energies at different
levels of pH. We show that incorrect assignment of fixed protonation
states in free energy computations can give errors of >2 kcal/mol
in these host–guest systems. Use of the methods presented here
avoids such errors, thus suggesting their utility in computing proton-linked
binding free energies for protein–ligand complexes.
Collapse
Affiliation(s)
- M Olivia Kim
- Department of Chemistry and Biochemistry, University of California San Diego , La Jolla, California 92093, United States
| | | | | | | |
Collapse
|
14
|
Abstract
AbstractProteins and other biomolecules contain acidic and basic titratable groups that give rise to intricate charge distributions and control electrostatic interactions. ‘Charge regulation’ concerns how the proton equilibria of these sites are perturbed when approached by alien molecular matter such as other proteins, surfaces and membranes, DNA, polyelectrolytes etc. Importantly, this perturbation generates a charge response that leads to attractive intermolecular interactions that can be conveniently described by a single molecular property – the charge capacitance. The capacitance quantifies molecular charge fluctuations, i.e. it is the variance of the mean charge and is an intrinsic property on par with the net charge and the dipole moment. It directly enters the free energy expression for intermolecular interactions and can be obtained experimentally from the derivative of the titration curve or theoretically from simulations. In this review, we focus on the capacitance concept as a predictive parameter for charge regulation and demonstrate how it can be used to estimate the interaction of a protein with other proteins, polyelectrolytes, membranes as well as with ligands.
Collapse
|
15
|
Abstract
Formation of protein-ligand complexes causes various changes in both the receptor and the ligand. This review focuses on changes in pK and protonation states of ionizable groups that accompany protein-ligand binding. Physical origins of these effects are outlined, followed by a brief overview of the computational methods to predict them and the associated corrections to receptor-ligand binding affinities. Statistical prevalence, magnitude and spatial distribution of the pK and protonation state changes in protein-ligand binding are discussed in detail, based on both experimental and theoretical studies. While there is no doubt that these changes occur, they do not occur all the time; the estimated prevalence varies, both between individual complexes and by method. The changes occur not only in the immediate vicinity of the interface but also sometimes far away. When receptor-ligand binding is associated with protonation state change at particular pH, the binding becomes pH dependent: we review the interplay between sub-cellular characteristic pH and optimum pH of receptor-ligand binding. It is pointed out that there is a tendency for protonation state changes upon binding to be minimal at physiologically relevant pH for each complex (no net proton uptake/release), suggesting that native receptor-ligand interactions have evolved to reduce the energy cost associated with ionization changes. As a result, previously reported statistical prevalence of these changes - typically computed at the same pH for all complexes - may be higher than what may be expected at optimum pH specific to each complex. We also discuss whether proper account of protonation state changes appears to improve practical docking and scoring outcomes relevant to structure-based drug design. An overview of some of the existing challenges in the field is provided in conclusion.
Collapse
Affiliation(s)
- Alexey V Onufriev
- Department of Computer Science and Physics, 2050 Torgersen Hall, Virginia Tech, Blacksburg, VA 24061, USA.
| | | |
Collapse
|
16
|
Schönichen A, Webb BA, Jacobson MP, Barber DL. Considering protonation as a posttranslational modification regulating protein structure and function. Annu Rev Biophys 2013; 42:289-314. [PMID: 23451893 DOI: 10.1146/annurev-biophys-050511-102349] [Citation(s) in RCA: 116] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Posttranslational modification is an evolutionarily conserved mechanism for regulating protein activity, binding affinity, and stability. Compared with established posttranslational modifications such as phosphorylation or ubiquitination, posttranslational modification by protons within physiological pH ranges is a less recognized mechanism for regulating protein function. By changing the charge of amino acid side chains, posttranslational modification by protons can drive dynamic changes in protein conformation and function. Addition and removal of a proton is rapid and reversible and, in contrast to most other posttranslational modifications, does not require an enzyme. Signaling specificity is achieved by only a minority of sites in proteins titrating within the physiological pH range. Here, we examine the structural mechanisms and functional consequences of proton posttranslational modification of pH-sensing proteins regulating different cellular processes.
Collapse
Affiliation(s)
- André Schönichen
- Department of Cell and Tissue Biology, University of California, San Francisco, USA
| | | | | | | |
Collapse
|
17
|
Roos G, Foloppe N, Messens J. Understanding the pK(a) of redox cysteines: the key role of hydrogen bonding. Antioxid Redox Signal 2013; 18:94-127. [PMID: 22746677 DOI: 10.1089/ars.2012.4521] [Citation(s) in RCA: 168] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Many cellular functions involve cysteine chemistry via thiol-disulfide exchange pathways. The nucleophilic cysteines of the enzymes involved are activated as thiolate. A thiolate is much more reactive than a neutral thiol. Therefore, determining and understanding the pK(a)s of functional cysteines are important aspects of biochemistry and molecular biology with direct implications for redox signaling. Here, we describe the experimental and theoretical methods to determine cysteine pK(a) values, and we examine the factors that control these pK(a)s. Drawing largely on experience gained with the thioredoxin superfamily, we examine the roles of solvation, charge-charge, helix macrodipole, and hydrogen bonding interactions as pK(a)-modulating factors. The contributions of these factors in influencing cysteine pK(a)s and the associated chemistry, including the relevance for the reaction kinetics and thermodynamics, are discussed. This analysis highlights the critical role of direct hydrogen bonding to the cysteine sulfur as a key factor modulating the equilibrium between thiol S-H and thiolate S(-). This role is easily understood intuitively and provides a framework for biochemical functional insights.
Collapse
Affiliation(s)
- Goedele Roos
- General Chemistry, Vrije University Brussel, Brussels, Belgium
| | | | | |
Collapse
|
18
|
Coussens NP, Schuck P, Zhao H. Strategies for assessing proton linkage to bimolecular interactions by global analysis of isothermal titration calorimetry data. THE JOURNAL OF CHEMICAL THERMODYNAMICS 2012; 52:95-107. [PMID: 22773848 PMCID: PMC3388511 DOI: 10.1016/j.jct.2012.02.008] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Isothermal titration calorimetry (ITC) is a traditional and powerful method for studying the linkage of ligand binding to proton uptake or release. The theoretical framework has been developed for more than two decades and numerous applications have appeared. In the current work, we explored strategic aspects of experimental design. To this end, we simulated families of ITC data sets that embed different strategies with regard to the number of experiments, range of experimental pH, buffer ionization enthalpy, and temperature. We then re-analyzed the families of data sets in the context of global analysis, employing a proton linkage binding model implemented in the global data analysis platform SEDPHAT, and examined the information content of all data sets by a detailed statistical error analysis of the parameter estimates. In particular, we studied the impact of different assumptions about the knowledge of the exact concentrations of the components, which in practice presents an experimental limitation for many systems. For example, the uncertainty in concentration may reflect imperfectly known extinction coefficients and stock concentrations or may account for different extents of partial inactivation when working with proteins at different pH values. Our results show that the global analysis can yield reliable estimates of the thermodynamic parameters for intrinsic binding and protonation, and that in the context of the global analysis the exact molecular component concentrations may not be required. Additionally, a comparison of data from different experimental strategies illustrates the benefit of conducting experiments at a range of temperatures.
Collapse
Affiliation(s)
- Nathan P. Coussens
- Laboratory of Cellular and Molecular Biology, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, U.S.A
| | - Peter Schuck
- Dynamics of Macromolecular Assembly, Laboratory of Cellular Imaging and Macromolecular Biophysics, National Institute of Biomedical Imaging and Bioengineering, National Institutes of Health, Bethesda, MD 20892, U.S.A
| | - Huaying Zhao
- Dynamics of Macromolecular Assembly, Laboratory of Cellular Imaging and Macromolecular Biophysics, National Institute of Biomedical Imaging and Bioengineering, National Institutes of Health, Bethesda, MD 20892, U.S.A
| |
Collapse
|
19
|
Yuan S, Le Roy K, Venken T, Lammens W, Van den Ende W, De Maeyer M. pKa modulation of the acid/base catalyst within GH32 and GH68: a role in substrate/inhibitor specificity? PLoS One 2012; 7:e37453. [PMID: 22662155 PMCID: PMC3360783 DOI: 10.1371/journal.pone.0037453] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2011] [Accepted: 04/21/2012] [Indexed: 11/18/2022] Open
Abstract
Glycoside hydrolases of families 32 (GH32) and 68 (GH68) belong to clan GH-J, containing hydrolytic enzymes (sucrose/fructans as donor substrates) and fructosyltransferases (sucrose/fructans as donor and acceptor substrates). In GH32 members, some of the sugar substrates can also function as inhibitors, this regulatory aspect further adding to the complexity in enzyme functionalities within this family. Although 3D structural information becomes increasingly available within this clan and huge progress has been made on structure-function relationships, it is not clear why some sugars bind as inhibitors without being catalyzed. Conserved aspartate and glutamate residues are well known to act as nucleophile and acid/bases within this clan. Based on the available 3D structures of enzymes and enzyme-ligand complexes as well as docking simulations, we calculated the pKa of the acid-base before and after substrate binding. The obtained results strongly suggest that most GH-J members show an acid-base catalyst that is not sufficiently protonated before ligand entrance, while the acid-base can be fully protonated when a substrate, but not an inhibitor, enters the catalytic pocket. This provides a new mechanistic insight aiming at understanding the complex substrate and inhibitor specificities observed within the GH-J clan. Moreover, besides the effect of substrate entrance on its own, we strongly suggest that a highly conserved arginine residue (in the RDP motif) rather than the previously proposed Tyr motif (not conserved) provides the proton to increase the pKa of the acid-base catalyst.
Collapse
Affiliation(s)
- Shuguang Yuan
- Laboratory of Molecular Plant Physiology, Institute of Botany and Microbiology, KU Leuven, Heverlee, Belgium
- Laboratory for Biomolecular Modelling, Department of Chemistry, Division of Biochemistry, Molecular and Structural Biology, KU Leuven, Heverlee, Belgium
| | - Katrien Le Roy
- Laboratory of Molecular Plant Physiology, Institute of Botany and Microbiology, KU Leuven, Heverlee, Belgium
| | - Tom Venken
- Laboratory for Biomolecular Modelling, Department of Chemistry, Division of Biochemistry, Molecular and Structural Biology, KU Leuven, Heverlee, Belgium
| | - Willem Lammens
- Laboratory of Molecular Plant Physiology, Institute of Botany and Microbiology, KU Leuven, Heverlee, Belgium
| | - Wim Van den Ende
- Laboratory of Molecular Plant Physiology, Institute of Botany and Microbiology, KU Leuven, Heverlee, Belgium
- * E-mail: (WV); (MDM)
| | - Marc De Maeyer
- Laboratory for Biomolecular Modelling, Department of Chemistry, Division of Biochemistry, Molecular and Structural Biology, KU Leuven, Heverlee, Belgium
- * E-mail: (WV); (MDM)
| |
Collapse
|
20
|
Vreven T, Hwang H, Pierce BG, Weng Z. Prediction of protein-protein binding free energies. Protein Sci 2012; 21:396-404. [PMID: 22238219 DOI: 10.1002/pro.2027] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2011] [Revised: 12/23/2011] [Accepted: 01/04/2012] [Indexed: 11/09/2022]
Abstract
We present an energy function for predicting binding free energies of protein-protein complexes, using the three-dimensional structures of the complex and unbound proteins as input. Our function is a linear combination of nine terms and achieves a correlation coefficient of 0.63 with experimental measurements when tested on a benchmark of 144 complexes using leave-one-out cross validation. Although we systematically tested both atomic and residue-based scoring functions, the selected function is dominated by residue-based terms. Our function is stable for subsets of the benchmark stratified by experimental pH and extent of conformational change upon complex formation, with correlation coefficients ranging from 0.61 to 0.66.
Collapse
Affiliation(s)
- Thom Vreven
- Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | | | | | | |
Collapse
|
21
|
Meyer T, Kieseritzky G, Knapp EW. Electrostatic pK
a
computations in proteins: Role of internal cavities. Proteins 2011; 79:3320-32. [DOI: 10.1002/prot.23092] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2011] [Revised: 05/03/2011] [Accepted: 05/09/2011] [Indexed: 11/09/2022]
|
22
|
Abstract
Accurate computational methods for predicting electrostatic energies are of major importance for our understanding of protein energetics in general for computer-aided drug design as well as for the design of novel biocatalysts and protein therapeutics. Electrostatic energies are of particular importance in such applications as virtual screening, drug design and protein-protein docking due to the high charge density of protein ligands and small-molecule drugs, and the frequent protonation state changes observed when drugs bind to their protein targets. Therefore, the development of a reliable and fast algorithm for the evaluation of electrostatic free energies, as an important contributor to the overall protein energy function, has been the focus for many scientists over the past three decades. In this review we describe the current state-of-the-art in modeling electrostatic effects in proteins and protein-ligand complexes. We focus mainly on the merits and drawbacks of the continuum methodology, and speculate on future directions in refining algorithms for calculating electrostatic energies in proteins using experimental data.
Collapse
|
23
|
Park MS, Gao C, Stern HA. Estimating binding affinities by docking/scoring methods using variable protonation states. Proteins 2011; 79:304-14. [PMID: 21058298 DOI: 10.1002/prot.22883] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
To investigate the effects of multiple protonation states on protein-ligand recognition, we generated alternative protonation states for selected titratable groups of ligands and receptors. The selection of states was based on the predicted pK(a) of the unbound receptor and ligand and the proximity of titratable groups of the receptor to the binding site. Various ligand tautomer states were also considered. An independent docking calculation was run for each state. Several protocols were examined: using an ensemble of all generated states of ligand and receptor, using only the most probable state of the unbound ligand/receptor, and using only the state giving the most favorable docking score. The accuracies of these approaches were compared, using a set of 176 protein-ligand complexes (15 receptors) for which crystal structures and measured binding affinities are available. The best agreement with experiment was obtained when ligand poses from experimental crystal structures were used. For 9 of 15 receptors, using an ensemble of all generated protonation states of the ligand and receptor gave the best correlation between calculated and measured affinities.
Collapse
Affiliation(s)
- Min-Sun Park
- Department of Biochemistry and Biophysics, University of Rochester, Rochester, New York 14642, USA
| | | | | |
Collapse
|
24
|
Abstract
The role of electrostatics in protein-protein interactions and binding is reviewed in this paper. A brief outline of the computational modeling, in the framework of continuum electrostatics, is presented and the basic electrostatic effects occurring upon the formation of the complex are discussed. The effect of the salt concentration and pH of the water phase on protein-protein binding free energy is demonstrated which indicates that the increase of the salt concentration tends to weaken the binding, an observation that is attributed to the optimization of the charge-charge interactions across the interface. It is pointed out that the pH-optimum (pH of optimal binding affinity) varies among the protein-protein complexes, and perhaps is a result of their adaptation to particular subcellular compartments. The similarities and differences between hetero- and homo-complexes are outlined and discussed with respect to the binding mode and charge complementarity.
Collapse
Affiliation(s)
- Zhe Zhang
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson,SC 29634, USA
| | | | | |
Collapse
|
25
|
Bjarnadottir U, Nielsen JE. Calculating pKa values in the cAMP-dependent protein kinase: the effect of conformational change and ligand binding. Protein Sci 2011; 19:2485-97. [PMID: 20954248 DOI: 10.1002/pro.530] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The conformational change observed upon ligand binding and phosphorylation for the cAMP-dependent protein kinase (protein kinase A-PKA) is of high importance for the regulation of its activity. We calculate pKa values and net charges for 18 3D structures of PKA in various conformations and liganded states to examine the role of electrostatics in ligand binding and activation. We find that the conformational change of PKA takes place without any significant net proton uptake/release at all pH values, thus indicating that PKA has evolved to reduce any pH-dependent barriers to the conformational motion. We furthermore find that the binding of ligands induces large changes in the net charge of PKA at most pH values, but significantly, we find that the net charge difference at physiological pH is close to zero, thus indicating that the active-site pKa values have been preorganized for substrate binding. We are unable to unequivocally resolve the identity of the groups responsible for determining the pH-activity profile of PKA but speculate that the titration of Lys 168 or the titration of ATP itself could be responsible for the loss of activity at high pH values. Finally, we examine the effect of point mutations on the pKa values of the PKA catalytic residues and find these to be relatively insensitive to both noncharge-altering and charge-altering mutations.
Collapse
Affiliation(s)
- Una Bjarnadottir
- School of Biomolecular and Biomedical Science, Centre for Synthesis and Chemical Biology, UCD Conway institute, University College Dublin, Belfield, Dublin 4, Ireland
| | | |
Collapse
|
26
|
Antosiewicz JM, Shugar D. Poisson–Boltzmann continuum-solvation models: applications to pH-dependent properties of biomolecules. MOLECULAR BIOSYSTEMS 2011; 7:2923-49. [DOI: 10.1039/c1mb05170a] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
27
|
Mitra RC, Zhang Z, Alexov E. In silico modeling of pH-optimum of protein-protein binding. Proteins 2010; 79:925-36. [PMID: 21287623 DOI: 10.1002/prot.22931] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2010] [Revised: 10/12/2010] [Accepted: 10/29/2010] [Indexed: 01/05/2023]
Abstract
Protein-protein association is a pH-dependent process and thus the binding affinity depends on the local pH. In vivo the association occurs in a particular cellular compartment, where the individual monomers are supposed to meet and form a complex. Since the monomers and the complex exist in the same micro environment, it is plausible that they coevolved toward its properties, in particular, toward the characteristic subcellular pH. Here we show that the pH at which the monomers are most stable (pH-optimum) or the pH at which stability is almost pH-independent (pH-flat) of monomers are correlated with the pH-optimum of maximal affinity (pH-optimum of binding) or pH interval at which affinity is almost pH-independent (pH-flat of binding) of the complexes made of the corresponding monomers. The analysis of interfacial properties of protein complexes demonstrates that pH-dependent properties can be roughly estimated using the interface charge alone. In addition, we introduce a parameter beta, proportional to the square root of the absolute product of the net charges of monomers, and show that protein complexes characterized with small or very large beta tend to have neutral pH-optimum. Further more, protein complexes made of monomers carrying the same polarity net charge at neutral pH have either very low or very high pH-optimum of binding. These findings are used to propose empirical rule for predicting pH-optimum of binding provided that the amino acid compositions of the corresponding monomers are available.
Collapse
Affiliation(s)
- Rooplekha C Mitra
- Physics Department, Computational Biophysics and Bioinformatics, Clemson University, Clemson, South Carolina 29634, USA
| | | | | |
Collapse
|
28
|
Fiorucci S, Zacharias M. Prediction of protein-protein interaction sites using electrostatic desolvation profiles. Biophys J 2010; 98:1921-30. [PMID: 20441756 DOI: 10.1016/j.bpj.2009.12.4332] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2009] [Revised: 11/30/2009] [Accepted: 12/23/2009] [Indexed: 01/15/2023] Open
Abstract
Protein-protein complex formation involves removal of water from the interface region. Surface regions with a small free energy penalty for water removal or desolvation may correspond to preferred interaction sites. A method to calculate the electrostatic free energy of placing a neutral low-dielectric probe at various protein surface positions has been designed and applied to characterize putative interaction sites. Based on solutions of the finite-difference Poisson equation, this method also includes long-range electrostatic contributions and the protein solvent boundary shape in contrast to accessible-surface-area-based solvation energies. Calculations on a large set of proteins indicate that in many cases (>90%), the known binding site overlaps with one of the six regions of lowest electrostatic desolvation penalty (overlap with the lowest desolvation region for 48% of proteins). Since the onset of electrostatic desolvation occurs even before direct protein-protein contact formation, it may help guide proteins toward the binding region in the final stage of complex formation. It is interesting that the probe desolvation properties associated with residue types were found to depend to some degree on whether the residue was outside of or part of a binding site. The probe desolvation penalty was on average smaller if the residue was part of a binding site compared to other surface locations. Applications to several antigen-antibody complexes demonstrated that the approach might be useful not only to predict protein interaction sites in general but to map potential antigenic epitopes on protein surfaces.
Collapse
Affiliation(s)
- Sébastien Fiorucci
- School of Engineering and Science, Jacobs University Bremen, Bremen, Germany.
| | | |
Collapse
|
29
|
Palpant NJ, Houang EM, Delport W, Hastings KEM, Onufriev AV, Sham YY, Metzger JM. Pathogenic peptide deviations support a model of adaptive evolution of chordate cardiac performance by troponin mutations. Physiol Genomics 2010; 42:287-99. [PMID: 20423961 DOI: 10.1152/physiolgenomics.00033.2010] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
In cardiac muscle, the troponin (cTn) complex is a key regulator of myofilament calcium sensitivity because it serves as a molecular switch required for translating myocyte calcium fluxes into sarcomeric contraction and relaxation. Studies of several species suggest that ectotherm chordates have myofilaments with heightened calcium responsiveness. However, genetic polymorphisms in cTn that cause increased myofilament sensitivity to activating calcium in mammals result in cardiac disease including arrhythmias, diastolic dysfunction, and increased susceptibility to sudden cardiac death. We hypothesized that specific residue modifications in the regulatory arm of troponin I (TnI) were critical in mediating the observed decrease in myofilament calcium sensitivity within the mammalian taxa. We performed large-scale phylogenetic analysis, atomic resolution molecular dynamics simulations and modeling, and computational alanine scanning. This study provides evidence that a His to Ala substitution within mammalian cardiac TnI (cTnI) reduced the thermodynamic potential at the interface between cTnI and cardiac TnC (cTnC) in the calcium-saturated state by disrupting a strong intermolecular electrostatic interaction. This key residue modification reduced myofilament calcium sensitivity by making cTnI molecularly untethered from cTnC. To meet the requirements for refined mammalian adult cardiac performance, we propose that compensatory evolutionary pressures favored mutations that enhanced the relaxation properties of cTn by decreasing its sensitivity to activating calcium.
Collapse
Affiliation(s)
- Nathan J Palpant
- Department of Integrative Biology and Physiology, University of Minnesota Academic Health Center, 321 Church Street SE, Minneapolis, MN 55455, USA
| | | | | | | | | | | | | |
Collapse
|
30
|
Aguilar B, Anandakrishnan R, Ruscio JZ, Onufriev AV. Statistics and physical origins of pK and ionization state changes upon protein-ligand binding. Biophys J 2010; 98:872-80. [PMID: 20197041 PMCID: PMC2830434 DOI: 10.1016/j.bpj.2009.11.016] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2009] [Revised: 10/21/2009] [Accepted: 11/04/2009] [Indexed: 10/19/2022] Open
Abstract
This work investigates statistical prevalence and overall physical origins of changes in charge states of receptor proteins upon ligand binding. These changes are explored as a function of the ligand type (small molecule, protein, and nucleic acid), and distance from the binding region. Standard continuum solvent methodology is used to compute, on an equal footing, pK changes upon ligand binding for a total of 5899 ionizable residues in 20 protein-protein, 20 protein-small molecule, and 20 protein-nucleic acid high-resolution complexes. The size of the data set combined with an extensive error and sensitivity analysis allows us to make statistically justified and conservative conclusions: in 60% of all protein-small molecule, 90% of all protein-protein, and 85% of all protein-nucleic acid complexes there exists at least one ionizable residue that changes its charge state upon ligand binding at physiological conditions (pH = 6.5). Considering the most biologically relevant pH range of 4-8, the number of ionizable residues that experience substantial pK changes (DeltapK > 1.0) due to ligand binding is appreciable: on average, 6% of all ionizable residues in protein-small molecule complexes, 9% in protein-protein, and 12% in protein-nucleic acid complexes experience a substantial pK change upon ligand binding. These changes are safely above the statistical false-positive noise level. Most of the changes occur in the immediate binding interface region, where approximately one out of five ionizable residues experiences substantial pK change regardless of the ligand type. However, the physical origins of the change differ between the types: in protein-nucleic acid complexes, the pK values of interface residues are predominantly affected by electrostatic effects, whereas in protein-protein and protein-small molecule complexes, structural changes due to the induced-fit effect play an equally important role. In protein-protein and protein-nucleic acid complexes, there is a statistically significant number of substantial pK perturbations, mostly due to the induced-fit structural changes, in regions far from the binding interface.
Collapse
Affiliation(s)
- Boris Aguilar
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia
| | | | - Jory Z. Ruscio
- Department of Bioengineering, University of California, Berkeley, California
| | - Alexey V. Onufriev
- Departments of Computer Science and Physics, Virginia Tech, Blacksburg, Virginia
| |
Collapse
|
31
|
Długosz M, Antosiewicz JM, Trylska J. pH-dependent association of proteins. The test case of monoclonal antibody HyHEL-5 and its antigen hen egg white lysozyme. J Phys Chem B 2010; 113:15662-9. [PMID: 19883097 DOI: 10.1021/jp906829z] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
We describe a method for determining diffusion-controlled rate constants for protein-protein association that explicitly includes the solution pH. The method combines the transient-complex theory for computing electrostatically enhanced association rates with an approach based on a rigorous thermodynamic cycle and partition functions for energy levels characterizing protonation states of associating proteins and their complexes. To test our method, we determine the pH-dependent kinetics of association of the HyHEL-5 antibody with its antigen hen egg white lysozyme. It was shown experimentally that their association rate constant depends on pH, increasing linearly in the pH range 6-8 and saturating or even exhibiting a flat maximum in the pH range 8-10. The presented methodology leads to a qualitative agreement with the experimental data. Our approach allows one to study diffusion-controlled protein-protein association under different pH conditions by taking into account the ensembles of protonation states rather than just the most probable protonation state of each protein.
Collapse
Affiliation(s)
- Maciej Długosz
- Interdisciplinary Centre for Mathematical and Computational Modelling, University of Warsaw, Zwirki i Wigury 93, Warsaw 02-089, Poland.
| | | | | |
Collapse
|
32
|
Rapp CS, Schonbrun C, Jacobson MP, Kalyanaraman C, Huang N. Automated site preparation in physics-based rescoring of receptor ligand complexes. Proteins 2009; 77:52-61. [PMID: 19382204 PMCID: PMC2744578 DOI: 10.1002/prot.22415] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Hydrogen atoms are not typically observable in X-ray crystal structures, but inferring their locations is often important in structure-based drug design. In addition, protonation states of the protein can change in response to ligand binding, as can the orientations of OH groups, a subtle form of "induced fit." We implement and evaluate an automated procedure for optimizing polar hydrogens in protein-binding sites in complex with ligands. Specifically, we apply the previously described Independent Cluster Decomposition Algorithm (ICDA) algorithm (Li et al., Proteins 2007;66:824-837), which assigns the ionization states of titratable residues, the amide orientations of Asn/Gln side chains, the imidazole ring orientation in His, and the orientations of OH/SH groups, in a unified algorithm. We test the utility of this method for identifying nativelike ligand poses using 247 protein-ligand complexes from an established database of docked decoys. Pose selection is performed with a physics-based scoring function based on a molecular mechanics energy function and a Generalized Born implicit solvent model. The use of the ICDA receptor preparation protocol, implemented with no knowledge of the native ligand pose, increases the accuracy of pose selection significantly, with the average RMSD over all complexes decreasing from 2.7 to 1.5 A when applying ICDA. Large improvements are seen for specific classes of binding sites with titratable groups, such as aspartyl proteases.
Collapse
Affiliation(s)
- Chaya S Rapp
- Department of Chemistry, Stern College for Women, Yeshiva University, New York, New York 10016, USA.
| | | | | | | | | |
Collapse
|
33
|
Talley K, Ng C, Shoppell M, Kundrotas P, Alexov E. On the electrostatic component of protein-protein binding free energy. PMC BIOPHYSICS 2008; 1:2. [PMID: 19351424 PMCID: PMC2666630 DOI: 10.1186/1757-5036-1-2] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/25/2008] [Accepted: 11/05/2008] [Indexed: 01/02/2023]
Abstract
Calculations of electrostatic properties of protein-protein complexes are usually done within framework of a model with a certain set of parameters. In this paper we present a comprehensive statistical analysis of the sensitivity of the electrostatic component of binding free energy (DeltaDeltaGel) with respect with different force fields (Charmm, Amber, and OPLS), different values of the internal dielectric constant, and different presentations of molecular surface (different values of the probe radius). The study was done using the largest so far set of entries comprising 260 hetero and 2148 homo protein-protein complexes extracted from a previously developed database of protein complexes (ProtCom). To test the sensitivity of the energy calculations with respect to the structural details, all structures were energy minimized with corresponding force field, and the energies were recalculated. The results indicate that the absolute value of the electrostatic component of the binding free energy (DeltaDeltaGel) is very sensitive to the force field parameters, the minimization procedure, the values of the internal dielectric constant, and the probe radius. Nevertheless our results indicate that certain trends in DeltaDeltaGel behavior are much less sensitive to the calculation parameters. For instance, the fraction of the homo-complexes, for which the electrostatics was found to oppose binding, is 80% regardless of the force fields and parameters used. For the hetero-complexes, however, the percentage of the cases for which electrostatics opposed binding varied from 43% to 85%, depending on the protocol and parameters employed. A significant correlation was found between the effects caused by raising the internal dielectric constant and decreasing the probe radius. Correlations were also found among the results obtained with different force fields. However, despite of the correlations found, the absolute DeltaDeltaGel calculated with different force field parameters could differ more than tens of kcal/mol in some cases. Set of rules of obtaining confident predictions of absolute DeltaDeltaGel and DeltaDeltaGel sign are provided in the conclusion section.PACS codes: 87.15.A-, 87.15. km.
Collapse
Affiliation(s)
- Kemper Talley
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson, SC 29634, USA
| | - Carmen Ng
- James Byrnes High School, Duncan, SC 29334, USA
| | - Michael Shoppell
- South Carolina Governor School for Science and Mathematics, Hartsville, SC 29550, USA
| | - Petras Kundrotas
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson, SC 29634, USA
- Center for Bioinformatics, The University of Kansas, Lawrence, KS 66047, USA
| | - Emil Alexov
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson, SC 29634, USA
| |
Collapse
|
34
|
Piłat Z, Antosiewicz JM. Multiple Protonation Equilibria in Electrostatics of Protein−Protein Binding. J Phys Chem B 2008; 112:15074-85. [DOI: 10.1021/jp8029659] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zofia Piłat
- Division of Biophysics, Institute of Experimental Physics, Department of Physics, University of Warsaw, Zwirki i Wigury 93 St., Warsaw 02-089, Poland
| | - Jan M. Antosiewicz
- Division of Biophysics, Institute of Experimental Physics, Department of Physics, University of Warsaw, Zwirki i Wigury 93 St., Warsaw 02-089, Poland
| |
Collapse
|
35
|
Bas DC, Rogers DM, Jensen JH. Very fast prediction and rationalization of pKa values for protein-ligand complexes. Proteins 2008; 73:765-83. [PMID: 18498103 DOI: 10.1002/prot.22102] [Citation(s) in RCA: 882] [Impact Index Per Article: 55.1] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Delphine C Bas
- Equipe de Chimie et Biochimie Théoriques, UMR 7565 - CNRS, Université Henri Poincaré, Nancy I, Boulevard des Aiguillettes BP 239, 54506 Vandoeuvre-lès-Nancy Cedex, France
| | | | | |
Collapse
|