Article Information

• PDF (217 kb)
• Full Text with Large Figures
• Supporting Material
• Cited by in Scopus (4)

PubMed

• Articles by David L. Bostick
• Articles by Charles L. Brooks

Related Articles

• Ammonium Recruitment and Ammonia Transport by E. coli Ammoni...

  Ammonium Recruitment and Ammonia Transport by E. coli Ammonia Channel AmtB Biophysical Journal, Volume 91, Issue 12, 15 December 2006, Pages 4401-4412 Thomas P. Nygaard, Carme Rovira, Günther H. Peters and Morten Ø. Jensen Abstract To investigate substrate recruitment and transport across the Ammonia transporter B (AmtB) protein, we performed molecular dynamics simulations of the AmtB trimer. We have identified residues important in recruitment of ammonium and intraluminal binding sites selective of ammonium, which provide a means of cation selectivity. Our results indicate that A162 guides translocation of an extraluminal ammonium into the pore lumen. We propose a mechanism for transporting the intraluminally recruited proton back to periplasm. Our mechanism conforms to net transport of ammonia and can explain why ammonia conduction is lost upon mutation of the conserved residue D160. We unify previous suggestions of D160 having either a structural or an ammonium binding function. Finally, our simulations show that the channel lumen is hydrated from the cytoplasmic side via the formation of single file water, while the F107/F215 stack at the inner-most part of the periplasmic vestibule constitutes a hydrophobic filter preventing AmtB from conducting water. Abstract | Full Text | PDF (1357 kb)

• Detailed Mechanism for AmtB Conducting NH4/NH3: Molecular Dy...

  Detailed Mechanism for AmtB Conducting NH4/NH3: Molecular Dynamics Simulations Biophysical Journal, Volume 92, Issue 3, 1 February 2007, Pages 877-885 Huaiyu Yang, Yechun Xu, Weiliang Zhu, Kaixian Chen and Hualiang Jiang Abstract The mechanism by which the ammonium transporter, AmtB, conducts into the cytoplasm was investigated using conventional molecular dynamics (MD) simulations. These simulations revealed that the neutral molecule, NH, passes automatically through the channel upon its arrival at the Am2 site and that the function of the channel as a one-way valve for passage of NH is not determined by the cytoplasmic exit gate but, rather, by the periplasmic entrance gate of the channel. The NH, produced by deprotonation of at the entrance gate, is spontaneously conveyed to the central region of the channel via a hydrogenbond network comprising His-168, His-318, Tyr-32, and the NH molecule. Finally, the NH molecule exits the channel through the exit gate formed by Phe-31, Ile-266, Val-314, and His-318. In addition, Ser-263 is shown to play a critical role in the final stages, acting as a pivoting arm to shunt the NH molecule from the cytoplasmic exit gate of the channel out into the cytoplasm. This finding is further supported by another simulation which shows that NH fails to be translocated through the channel formed by the Ser-263–Ala mutation. Thus, this study casts new insights on the mechanism of AmtB-mediated passage of NH across cellular membranes. Abstract | Full Text | PDF (1147 kb)

• A Stable Water Chain in the Hydrophobic Pore of the AmtB Amm...

  A Stable Water Chain in the Hydrophobic Pore of the AmtB Ammonium Transporter Biophysical Journal, Volume 92, Issue 9, 1 May 2007, Pages L82-L84 Guillaume Lamoureux, Michael L. Klein and Simon Bernèche Abstract The accessibility of water molecules to the pore of the AmtB ammonium transporter is studied using molecular dynamics simulations. Free energy calculations show that the so-called hydrophobic pore can stabilize a chain of water molecules in a well of a few kcal/mol, using a favorable electrostatic binding pocket as an anchoring point. Moreover, the structure of the water chain matches precisely the electronic density maxima observed in x-ray diffraction experiments. This result questions the general assumption that the AmtB pore only contains ammonia (NH) molecules diffusing in a single file fashion. The probable presence of water molecules in the pore would influence the relative stability of NH and , and thus calls for a reassessment of the overall permeation mechanism in ammonium transporters. Abstract | Full Text | PDF (246 kb)

• …more

Copyright © 2007 The Biophysical Society. All rights reserved.
Biophysical Journal, Volume 92, Issue 12, L103-L105, 15 June 2007

doi:10.1529/biophysj.107.109165

Biophysical Letters

On the Equivalence Point for Ammonium (De)protonation during Its Transport through the AmtB Channel

David L. Bostick and Charles L. Brooks^,  

Department of Molecular Biology and Center for Theoretical Biological Physics, The Scripps Research Institute, La Jolla, California 92037

Address reprint requests and inquiries to Charles Brooks.

Abstract

Structural characterization of the bacterial channel, AmtB, provides a glimpse of how members of its family might control the protonated state of permeant ammonium to allow for its selective passage across the membrane. In a recent study, we employed a combination of simulation techniques that suggested ammonium is deprotonated and reprotonated near dehydrative phenylalanine landmarks (F107 and F31, respectively) during its passage from the periplasm to the cytoplasm. At these landmarks, ammonium is forced to maintain a critical number (∼3) of hydrogen bonds, suggesting that the channel controls ammonium (de)protonation by controlling its coordination/hydration. In the work presented here, a free energy-based analysis of ammonium hydration in dilute aqueous solution indicates, explicitly, that at biological pH, the transition from ammonium () to ammonia (NH₃) occurs when these species are constrained to donate three hydrogen bonds or less. This result demonstrates the viability of the proposal that AmtB indirectly controls ammonium (de)protonation by directly controlling its hydration.

AmtB exists in the membrane as a homotrimer. Each monomer of this protein forms a channel that passively transports ammonium () in the form of its “gas” ammonia (NH₃) intermediate across the membranes of bacteria; for conciseness we will henceforth refer to both and NH₃ species, together, as Am. Structural models of AmtB resulting from x-ray diffraction 1,2 have provided initial configurations for a plethora of computational 3,4,5,6,7,8,9,10,13 studies aimed at understanding this channel's mechanistic aspects and implications for homologous human counterparts.

The center of an AmtB monomer forms a narrow hydrophobic pore (lumen) connecting cytoplasmic and periplasmic vestibules, both accessible to aqueous solution. Diffraction studies revealed an binding site in the cytoplasmic vestibule (site Am1 1,2) where the cation donates hydrogen bonds to the backbone carbonyl group of A162, the side-chain hydroxyl oxygen of S219, and ∼2–3 water molecules 3,5,7. Aromatic groups (F107 and F215) form a floor for site Am1, capping the hydrophobic lumen to help prevent entrance of water from the periplasm (see Fig. 1). These aromatic groups rotate at low free energy cost to allow translocation of Am 3,5,7 under the influence of an electrochemical gradient.

Display large version of this figure

Figure 1

(A) Synopsis of binding sites and deprotonation regions. Some helices have been removed from the protein (orange) to allow visibility of the binding sites. The positions of the x-ray sites are shown as red spheres, and Am binding sites derived from simulation 3 are shown as either blue spheres (for NH₃) or blue spheres with white hydrogen atoms (for ). (De)protonation regions are marked in green (see arrows). The periplasmic and cytoplasmic (de)protonation regions coincide with the phenyl groups of F107 and F31, respectively. (B) Hydration of NH₃ (diamonds with dashed line) and (circles with solid line) as a function of the transport axis (Z). The origin coincides with the center of mass of the AmtB trimer (Z>0 is the periplasmic membrane leaflet, and Z<0 is the cytoplasmic). Dotted lines denote lipid phosphate positions, gray bars mark sites Am1 and Am5, red bars mark sites Am2-4, and green bars (see arrows) mark (de)protonation regions (Am equivalence points). (C) Proposed equivalence point for deprotonation near the periplasmic end of the lumen. Here, is stripped to three hydrogen bonds (donated to carbonyl groups of A162 and F215 and one water), and has full access to vestibular water, continuously connected to bulk solution, to allow the escape of a proton to the bulk in the form of hydronium ion.

In the presence of AmSO₄, the x-ray structure 1 displayed three luminal binding sites (Am2, Am3, and Am4—see Figure 1A), where Am interacts closely with His residues (H168 and H318). Calculations of the apparent pK_a of luminal Am 3,10 indicate that these sites may only be occupied by neutral NH₃. As such, it would appear that the disallowance of permanently charged species in the lumen is the most Am-selective feature of AmtB. An aromatic group (F31) just below site Am4 helps to prevent hydration of the lumen, and provides a low free energy barrier for NH₃ passage to the cytoplasmic vestibule (Figure 1AB). Just below the lumen, a fifth site (Am5) was revealed by a molecular dynamics (MD) study 3. At this site, calculations of the apparent pK_a3 suggest Am must exist in its protonated form, where it donates hydrogen bonds to a carboxyl oxygen of D313, the hydroxyl oxygen of S263, and surrounding water (Figure 1AB).

Combining knowledge of experimental and computational results 1,2,3,10, it appears that AmtB deprotonates between sites Am1–2, and reprotonates NH₃ between sites Am4–5 to allow Am flux toward the cytoplasm. However, it is difficult to determine, experimentally, how the channel controls these (de)protonation events. Computational studies, though they should help clarify the (de)protonation mechanism, have proposed disparate explanations 3,4,5,7. Lin et al. 5 and Nygaard et al. 7 both proposed that a highly conserved Asp residue (D160), whose mutation is known to destroy AmtB's transport capability 11, plays a key role in deprotonation. Lin et al. 5 observed that water forms a hydrogen bonded network between at Am1 and the carboxylate of D160. This led them to suggest that the charged carboxylate drives deprotonation at site Am1, and accepts a proton donated by using hydronium as an intermediate. On the other hand, Nygaard et al. 7 proposed that deprotonation occurs near site Am2, after moves from Am1 across the stacked (F107/F215) aromatic moieties. In this configuration, it was suggested that donates a proton to D160 via the backbone carbonyl group of A162 and the amide N–H of G163 using an imidic acid mechanism.

Luzhkov et al. 10 presented results that would suggest that D160 does not function as a proton acceptor. Rather, their calculations showed that the apparent pK_a of D160's carboxylate is downshifted (from its standard value of ∼3.9) by 0.3–5.1 units when site Am1 is unoccupied. When occupies Am1, the apparent pK_a of D160 shifts even further downward by 9.2 units, making its protonation effectively impossible. Our own results 3, as well as those of Luzhkov et al., showed that D160 is engaged in persistent hydrogen bonds with the protein, and that the negative charge of D160 stabilizes Am in its protonated form, shifting its apparent pK_a upward by ∼4 units. Taken together, these results indicate that the importance of D160, as evidenced by mutational studies 11, is more likely due to recruitment of from the periplasm and stabilizing its binding at site Am1 rather than accepting a proton as suggested by Lin et al. and Nygaard et al.

Recently we utilized a combination of MD simulation techniques 3, showing that the equivalence points for Am (de)protonation coincide with the periplasmic and cytoplasmic phenyl groups of F107 and F31, respectively (Figure 1A). Near these specific regions, Am was seen to be stripped to ∼3 or fewer hydrogen bonds (Figure 1B). At the periplasmic (de)protonation site (Figure 1BC), near F107, may donate two hydrogen bonds to protein and ∼1 to water. At the cytoplasmic (de)protonation site, it appears that water provides all ∼3 hydrogen bonds (Figure 1B). Given that, a), water ionizes more easily than a carbonyl group, b), the carboxylate of D160 is persistently engaged in hydrogen bonds with the protein that shift its apparent pK_a downward 10, and c), at both equivalence points, Am has full access to vestibular water, we proposed that water is the only plausible proton acceptor for After accepting this proton, it is most likely that the proton escapes to the periplasm in the form of hydronium.

Our previous study showed a clear correlation between the protonated form of Am and the number of available hydrogen bonds. However, we did not directly demonstrate that the channel need only constrain Am to ∼3 or fewer hydrogen bonds to deprotonate Such a demonstration would require a computational experiment that rules out other effects, such as the local electric field imposed by the channel at the (de)protonation regions. To address this issue in a most general way, we performed simulations of both and NH₃ in dilute aqueous solution (see Supplementary Material ). These simple simulations allowed us to isolate and directly probe the dependence of Am's apparent pK_a on coordination (hydrogen-bond) number.

In the spirit of previous work 12, we performed a free energy characterization of solute ( and NH₃) hydration preferences based upon population analysis from MD trajectories. From this analysis, we derived the apparent pK_a of Am as a function of its coordination number as follows (see Supplementary Material ):

The resulting apparent pK_a profile is shown in Fig. 2. This analysis indicates that if a local environment provides only ∼3 or fewer hydrogen bonds, Am will be favored in its deprotonated form, NH₃. The equivalence point between and NH₃, itself, appears to occur near the midpoint between three and four available acceptors. In tandem with previous results 3,10, this indicates that the functional role of D160 is to allow for AmtB's structural and electrostatic ability to recruit from the periplasm, and not to drive the deprotonation of Since Fig. 2 describes the local pK_a of Am in an isotropic medium—that of pure water—the analysis shows that loss of a proton from occurs with ∼3 hydrogen bonds regardless of any external field provided by AmtB at the calculated equivalence points shown in Fig. 1. The result we show here appears to be independent of the force field chosen to describe Am or water (Supplementary Material, Fig. S1 ), and suggests that AmtB's control over Am hydration, or equivalently, the number of hydrogen bonds, is the sole control over (de)protonation provided by the protein, as we have suggested 3.

Display large version of this figure

Figure 2

(A) Apparent pK_a of Am in a local environment that provides N_Am-H hydrogen bonds for and N_Am hydrogen bonds for NH₃ derived from population analysis of MD trajectories for Am in dilute aqueous solution. The green contour line indicates the equivalence point (pK_a=7). The diagonal dashed line indicates where N_Am-H=N_Am. (B) Apparent pK_a of Am for a local hydration environment where N_Am-H=N_Am. Note that the equivalence point occurs in an environment providing 3.6 hydrogen bonds or less (see the dotted lines).

It is also interesting to consider this result in light of a recent study suggesting a stable water chain can enter the lumen 13 from the cytoplasm. If, indeed, this occurs, might occupy the lumen more favorably than previously thought. Our previous work 3 suggested that if enters the lumen, it may be hydrated by as many as 2–3 water molecules (Figure 1B). However, the data presented here (Fig. 2) indicate Am will exist as NH₃ if only three hydrogen bond partners are provided. Thus, the (de)protonation sites we suggest (Fig. 1 and Bostick and Brooks 3) can still hold true despite a hydrated pore. Also, the preference of NH₃ for threefold (or less) coordination forces us to consider the possibility that NH₃ and H₂O might coexist in a confined luminal environment. We suggest that future computational study aimed at determining Am's protonation state in a hydrated lumen may shed light on this issue.