skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Amino acid deprotonation rates from classical force fields
Acid ionization constants (pK a ’s) of titratable amino acid side chains have received a large amount of experimental and theoretical attention. In many situations, however, the rates of protonation and deprotonation, k on and k off , may also be important, for example, in understanding the mechanism of action of proton channels or membrane proteins that couple proton transport to other processes. Protonation and deprotonation involve the making and breaking of covalent bonds, which cannot be studied by classical force fields. However, environment effects on the rates should be captured by such methods. Here, we present an approach for estimating deprotonation rates based on Warshel’s extension of Marcus’s theory of electron transfer, with input from molecular simulations. The missing bond dissociation energy is represented by a constant term determined by fitting the pK a value in solution. The statistics of the energy gap between protonated and deprotonated states is used to compute free energy curves of the two states and, thus, free energy barriers, from which the rate can be deduced. The method is applied to Glu, Asp, and His in bulk solution and select membrane proteins: the M2 proton channel, bacteriorhodopsin, and cytochrome c oxidase.  more » « less
Award ID(s):
1855942
PAR ID:
10416411
Author(s) / Creator(s):
;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
157
Issue:
8
ISSN:
0021-9606
Page Range / eLocation ID:
085101
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Lactose permease is a paradigm for the major facilitator superfamily, the largest family of ion-coupled membrane transport proteins known at present. LacY carries out the coupled stoichiometric symport of a galactoside with an H+, using the free energy released from downhill translocation of H+to drive accumulation of galactosides against a concentration gradient. In neutrophilicEscherichia coli, internal pH is kept at ∼7.6 over the physiological range, but the apparent pK (pKapp) for galactoside binding is 10.5. Surface-enhanced infrared absorption spectroscopy (SEIRAS) demonstrates that the high pKais due to Glu325 (helix X), which must be protonated for LacY to bind galactoside effectively. Deprotonation is also obligatory for turnover, however. Here, we utilize SEIRAS to study the effect of mutating residues in the immediate vicinity of Glu325 on its pKa. The results are consistent with the idea that Arg302 (helix IX) is important for deprotonation of Glu325. 
    more » « less
  2. null (Ed.)
    We present a single molecular polyoxometalate cluster (K 41 [(P 2 W 12 Nb 6 O 62 ) 6 {Mn 3 (OH) 3 (H 2 O) 6 } 4 {Mn 3 Na(H 2 O) 16 }]·26H 2 O) with controllable release of a large number of protons (∼40 per molecule) in its aqueous solution upon addition of a base. The deprotonation/protonation process is reversible with the clusters remaining intact. This molecule can also absorb up to 11 protons per cluster when an acid, HCl, was added to its original aqueous solution. To the best of our knowledge, such large proton absorption/release capacity along with excellent stability is unprecedented. 
    more » « less
  3. EmrE is anEscherichia colimultidrug efflux pump and member of the small multidrug resistance (SMR) family that transports drugs as a homodimer by harnessing energy from the proton motive force. SMR family transporters contain a conserved glutamate residue in transmembrane 1 (Glu14 in EmrE) that is required for binding protons and drugs. Yet the mechanism underlying proton-coupled transport by the two glutamate residues in the dimer remains unresolved. Here, we used NMR spectroscopy to determine acid dissociation constants (pKa) for wild-type EmrE and heterodimers containing one or two Glu14 residues in the dimer. For wild-type EmrE, we measured chemical shifts of the carboxyl side chain of Glu14 using solid-state NMR in lipid bilayers and obtained unambiguous evidence on the existence of asymmetric protonation states. Subsequent measurements of pKavalues for heterodimers with a single Glu14 residue showed no significant differences from heterodimers with two Glu14 residues, supporting a model where the two Glu14 residues have independent pKavalues and are not electrostatically coupled. These insights support a transport pathway with well-defined protonation states in each monomer of the dimer, including a preferred cytoplasmic-facing state where Glu14 is deprotonated in monomer A and protonated in monomer B under pH conditions in the cytoplasm ofE. coli. Our findings also lead to a model, hop-free exchange, which proposes how exchangers with conformation-dependent pKavalues reduce proton leakage. This model is relevant to the SMR family and transporters comprised of inverted repeat domains. 
    more » « less
  4. null (Ed.)
    Rotary vacuolar adenosine triphosphatases (V-ATPases) drive transmembrane proton transport through a V o proton channel subcomplex. Despite recent high-resolution structures of several rotary ATPases, the dynamic mechanism of proton pumping remains elusive. Here, we determined a 2.7-Å cryo–electron microscopy (cryo-EM) structure of yeast V o proton channel in nanodisc that reveals the location of ordered water molecules along the proton path, details of specific protein-lipid interactions, and the architecture of the membrane scaffold protein. Moreover, we uncover a state of V o that shows the c -ring rotated by ~14°. Molecular dynamics simulations demonstrate that the two rotary states are in thermal equilibrium and depict how the protonation state of essential glutamic acid residues couples water-mediated proton transfer with c -ring rotation. Our cryo-EM models and simulations also rationalize a mechanism for inhibition of passive proton transport as observed for free V o that is generated as a result of V-ATPase regulation by reversible disassembly in vivo. 
    more » « less
  5. Interactions in enzymes between catalytic and neighboring amino acids and how these interactions facilitate catalysis are examined. In examples from both natural and designed enzymes, it is shown that increases in catalytic rates may be achieved through elongation of the buffer range of the catalytic residues; such perturbations in the protonation equilibria are, in turn, achieved through enhanced coupling of the protonation equilibria of the active ionizable residues with those of other ionizable residues. The strongest coupling between protonation states for a pair of residues that deprotonate to form an anion (or a pair that accept a proton to form a cation) is achieved when the difference in the intrinsic pKas of the two residues is approximately within 1 pH unit. Thus, catalytic aspartates and glutamates are often coupled to nearby acidic residues. For an anion-forming residue coupled to a cation-forming residue, the elongated buffer range is achieved when the intrinsic pKa of the anion-forming residue is higher than the intrinsic pKa of the (conjugate acid of the) cation-forming residue. Therefore, the high pKa, anion-forming residues tyrosine and cysteine make good coupling partners for catalytic lysine residues. For the anion–cation pairs, the optimum difference in intrinsic pKas is a function of the energy of interaction between the residues. For the energy of interaction ε expressed in units of (ln 10)RT, the optimum difference in intrinsic pKas is within ∼1 pH unit of ε. 
    more » « less