Hostname: page-component-78c5997874-lj6df Total loading time: 0 Render date: 2024-11-13T06:59:22.648Z Has data issue: false hasContentIssue false

Titratable residues that drive RND efflux: Insights from molecular simulations

Published online by Cambridge University Press:  01 April 2024

Robert Clark
Affiliation:
Department of Biochemistry, University of Oxford, Oxford, UK
Kahlan E. Newman
Affiliation:
School of Chemistry, University of Southampton, Southampton, UK
Syma Khalid*
Affiliation:
Department of Biochemistry, University of Oxford, Oxford, UK School of Chemistry, University of Southampton, Southampton, UK
*
Corresponding author: Syma Khalid; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

The resistance–nodulation–division efflux machinery confers antimicrobial resistance to Gram-negative bacteria by actively pumping antibiotics out of the cell. The protein complex is powered by proton motive force; however, the proton transfer mechanism itself and indeed even its stoichiometry is still unclear. Here we review computational studies from the last decade that focus on elucidating the number of protons transferred per conformational cycle of the pump. Given the difficulties in studying proton movement using even state-of-the-art structural biology methods, the contributions from computational studies have been invaluable from a mechanistic perspective.

Type
Perspective
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2024. Published by Cambridge University Press

Introduction

Gram-negative bacteria, which comprise two-thirds of ESKAPE pathogens, enhance their antimicrobial resistance profiles using resistance–nodulation–division (RND) efflux machinery. These tripartite protein complexes span the cell envelope and confer multidrug resistance by expelling molecules from within the periplasm/proximal to the inner membrane into the extracellular milieu, thereby reducing intracellular accumulation of antibiotics (Figure 1a). Such machinery has thus been the focus of intense research in an attempt to better understand the process of efflux (Alav et al., Reference Alav, Kobylka, Kuth, Pos, Picard, Blair and Bavro2021). The most widely studied of these systems is AcrAB-TolC from the archetypal Gram-negative bacterium, E. coli. AcrB, the homotrimeric RND protein of this assembly, drives substrate recognition and extrusion (Figure 1b). AcrB protomers cycle the Access/Loose, Binding/Tight, and Extrusion/Open states to bind substrates and expel them into the periplasmic adaptor protein (AcrA) before diffusing out of the cell via the outer membrane factor (TolC). In this review we will refer to these RND protomer states only as Access (A), Binding (B), and Extrusion (E).

Figure 1. AcrB forms part of the tripartite multidrug efflux pump, AcrAB-TolC. (a) Surface representation of AcrAB-TolC (PDB ID: 5O66; Du et al., Reference Du, Wang, James, Voss, Klimont, Ohene-Agyei, Venter, Chiu and Luisi2014) and its situation within the cell envelope of Gram-negative bacteria. AcrB sits within the inner membrane; TolC sits within the outer membrane; and AcrA is the periplasmic adaptor protein. Substrates, such as certain classes of antibiotics, are thought to enter AcrB either in the periplasm near the inner membrane, or from the outer leaflet of the inner membrane directly. Protons enter AcrB from a different entry point, and transit through the TMD to the cytoplasm. The precise mechanism of the proton relay is unclear. TolC is coloured in purple, AcrA is coloured in orange, AcrB is coloured according to state. The Access state is coloured pink, the Binding state is coloured yellow, and the Extrusion state is coloured in blue. (b) Cartoon representation of AcrB (PDB ID: 4DX5; Eicher et al., Reference Eicher, Cha, Seeger, Brandstätter, El-Delik, Bohnert, Kern, Verrey, Grütter, Diederichs and Pos2012). AcrB has three main domains: the Docking domain which interacts with AcrA; the Porter domain contains the drug entry site; and the Transmembrane domain uses the proton motive force to induce drug extrusion. The titratable residues thought to orchestrate proton transfer are shown in stick representation and coloured according to element.

During conformational cycling, the periplasmic cleft, defined by the porter subdomains PC1 and PC2, is open in the Access and Binding states (Figures 2 and 3a), while the exit gate in the docking domain is closed. It is thought that drugs are bound in a proximal binding pocket in the Access state before moving to the distal pocket to induce transition of the protomer to the Binding state (Figures 2 and 3b). Transformation to the Extrusion state results in the closing of the periplasmic cleft and the opening of the exit gate through which the drug is extruded into the periplasmic adaptor protein. The B to E transition is putatively an energy-dependent step (Seeger et al., Reference Seeger, Diederichs, Eicher, Brandstätter, Schiefner, Verrey and Pos2008), powered by proton motive force: protons are transferred from the periplasm to the cytoplasm via the transmembrane domain (TMD). Buried within the TMD are essential/highly conserved polar/charged residues that constitute a proton relay network, facilitating proton transfer across this otherwise hydrophobic region. D407, D408, K940, R971, and T973 have been identified as essential to function (Murakami and Yamaguchi, Reference Murakami and Yamaguchi2003; Takatsuka and Nikaido, Reference Takatsuka and Nikaido2006), with nearby N941, G403, and S481 also highly conserved. The exact mechanism of proton transfer, proton stoichiometry for each full cycle, and which of the acidic residues is protonated at each stage of the cycle remains uncertain.

Figure 2. Cartoon schematic of functional rotation and substrate transport through inner membrane RND transporter proteins. (a) Side-view of the RND protein in the inner membrane. (b) Top-down (from the periplasm) view of the porter domain. For visual clarity, protonation and drug extrusion are considered for a single protomer only. Substrates enter the proximal binding pocket of the Access protomer from the periplasm/periplasmic leaflet of the inner membrane. Substrate binding induces a conformational change in the protomer to the Binding state, and the substrate moves to the distal binding pocket. Protonation then occurs in the relay in the transmembrane domain, inducing a conformational change to the Extrusion state. The periplasmic cleft closes and the exit gate opens, allowing the substrate to exit into the periplasmic adaptor protein (not shown).

Figure 3. AcrB undergoes a conformational cycle that functionally connects the movement of protons through the protein to drug extrusion. (a) Porter domains of AcrB viewed from the periplasm, showing the different conformational changes that correspond to drug access, binding and extrusion. The Access state is coloured pink, the Binding state is coloured yellow, and the Extrusion state is coloured in blue. Adapted from Matsunaga et al. (Reference Matsunaga, Yamane, Terada, Moritsugu, Fujisaki, Murakami, Ikeguchi and Kidera2018). (b) Drugs are thought to move through the proximal binding pocket (purple) in the Access state, bind to the distal pocket (dark blue) in the Binding state, and leave through the exit gate (cyan) in the Extrusion state. A drug (minocycline), is shown moving through these different regions. The drug was placed manually in the Access and Extrusion states for illustrative purposes, but the position of the drug in the binding state comes from PDB ID:4DX5 (Eicher et al., Reference Eicher, Cha, Seeger, Brandstätter, El-Delik, Bohnert, Kern, Verrey, Grütter, Diederichs and Pos2012). The definition of the different regions comes from Vargiu and Nikaido (Reference Vargiu and Nikaido2012).

Structural data and early simulation studies

Computational studies have been widely used to investigate efflux machinery (Athar et al., Reference Athar, Gervasoni, Catte, Basciu, Malloci, Ruggerone and Vargiu2023). Early simulations did not focus explicitly on the proton relay network, instead utilising postulated protonation states (Seeger et al., Reference Seeger, Diederichs, Eicher, Brandstätter, Schiefner, Verrey and Pos2008). Simulations by Fischer and Kandt (Reference Fischer and Kandt2011) demonstrated alternating access for water (important for mediating proton transfer) entering the TMD from the periplasm and cytoplasm in different protomer states. Further simulations by the same group displayed flexibility in the porter domain, with opening and closing of the periplasmic cleft (Fischer and Kandt, Reference Fischer and Kandt2013). For an extensive discussion of RND simulation studies, we direct the reader to a recent review (Athar et al., Reference Athar, Gervasoni, Catte, Basciu, Malloci, Ruggerone and Vargiu2023). Here we turn our attention only to simulation studies where the protonation states of the proton relay network residues have been a primary focus.

As there are two titratable residues in the relay network (D407 and D408), it is assumed that the number of protons transferred will be one or two (Alav et al., Reference Alav, Kobylka, Kuth, Pos, Picard, Blair and Bavro2021). Based on structural data, the current consensus is that both D407 and D408 are deprotonated in the A and B protomers: K940 sits between these two residues, within salt-bridging distance of both (Figure 4a,b). This relay network arrangement has been denoted the ‘engaged’ state. In the Extrusion protomer, K940 orients away from the aspartates and towards T978/N941 (‘disengaged’ state, Figure 4c). The protonation states of D407 and D408 in the Extrusion protomer remain unclear and cannot be elucidated explicitly via state-of-the-art experimental structure-determination techniques as hydrogens are not well-resolved. Protonation may be inferred from side-chain orientations, but this can yield ambiguous or conflicting results. In the Extrusion state the distance between the G403 backbone oxygen and D407 carboxylate (<4 Å, Figure 4b) suggests D407 is protonated to participate in hydrogen bonding. Simultaneously, the D407 carboxylate is also within salt-bridging distance (~4 Å) of the K940 amine nitrogen, implying D407 is deprotonated! Experiments such as carbodiimide labelling to test protonation states (Seeger et al., Reference Seeger, von Ballmoos, Verrey and Pos2009) can be a powerful investigative tool, but the bulky reagents are restricted in the narrow water wires of the TMD. Modelling and molecular dynamics (MD) simulations have been used to address this knowledge gap. These studies generally fall into two categories, those that support: single proton transfer; or transfer of two protons per conformational cycle (Figure 4c).

Figure 4. The titratable residues of AcrB adopt different configurations in different states. (a) Close-up view of key residues of the proton relay. Each image is rotated, so viewed from the same orientation. In the Access and Binding states the relay network is considered ‘engaged’, and in the Extrusion states the network is considered ‘disengaged’. Adapted from Matsunaga et al. (Reference Matsunaga, Yamane, Terada, Moritsugu, Fujisaki, Murakami, Ikeguchi and Kidera2018). (b) Inter-residue distances for residues in the proton relay network of Gram-negative bacterial RND protein structures deposited in the PDB (raw data and relevant references available in the Supplementary material). Distances were calculated as the shortest distance between side chain oxygen/nitrogen atoms, except for G403 where the backbone carbonyl oxygen was used. (c) Conformational transition from the Binding to the Extrusion state is proposed to proceed via one protonation event (D407/D408H) or two (D407H/D408H). Structures shown in (a) and (c) are of PDB ID: 4DX5.

Two protons

Computational support for protonation of both D407 and D408 is based largely on structures and simulations presented by Pos and colleagues (Eicher et al., Reference Eicher, Seeger, Anselmi, Zhou, Brandstätter, Verrey, Diederichs, Faraldo-Gómez and Pos2014). Eicher et al. extended the analysis of their previous crystal structure (4DX5 (Eicher et al., Reference Eicher, Cha, Seeger, Brandstätter, El-Delik, Bohnert, Kern, Verrey, Grütter, Diederichs and Pos2012)) to focus on different conformations adopted by key residues of the AcrB TMD, including D407 and D408. They compared the different structures and calculated likely protonation states of key residues using two separate approaches.

Ten independent Metropolis Monte Carlo simulations (Brooks et al., Reference Brooks, Brooks, Mackerell, Nilsson, Petrella, Roux, Won, Archontis, Bartels, Boresch, Caflisch, Caves, Cui, Dinner, Feig, Fischer, Gao, Hodoscek, Im, Kuczera, Lazaridis, Ma, Ovchinnikov, Paci, Pastor, Post, Pu, Schaefer, Tidor, Venable, Woodcock, Wu, Yang, York and Karplus2009) were conducted to evaluate protonation states of K940, R971, and all aspartate, histidine and glutamate residues. Energies were evaluated using the Poisson-Boltzmann equation. Each simulation comprised five cycles: 100,000 steps for annealing, starting at 400 K to encourage greater exploration, followed by 500,000 steps at 300 K for equilibration. Proposed changes in protonation states were accepted or rejected based on the Metropolis criterion. By averaging the protonation state over the final 500,000 steps across all 10 simulations, a global average was obtained. The periplasmic pH was assumed to be 5.5 and the cytoplasmic pH 7.5, with a dielectric constant (εp) of 4. It was noted a pH of 7.5 across the protein led to analogous results. These simulations predicted that in the E state, D408 would be protonated 100% of the time and D407 90% of the time. The study also evaluated the likelihood of protonation for the mutants K940A, R971A, D407N and D408N. When present, both D407 and D408 were predicted to be protonated in the E states. These mutants were shown previously to be functionally inactive (Guan and Nakae, Reference Guan and Nakae2001; Takatsuka and Nikaido, Reference Takatsuka and Nikaido2006): Pos and colleagues related the change in activity of these mutants to alterations in the proton relay network.

Acknowledging the dependence of this electrostatics approach on dielectric constants, free energy perturbation (FEP) was also used to re-evaluate the protonation states of key residues. All-atom simulations of the relevant AcrB protomers were used, alchemically transforming the protonation state of specific residues to calculate the change in potential energy of the system. pKas were calculated using known side chain analogues, yielding values of 8.2 and 10.8 in the E state for D407 and D408, respectively. This agreed with their protonation probability prediction that both D407 and D408 were protonated in the E state, and matched crystal structure side chain placements (Figure 4a). It should be noted there were large errors for D407, indicating that the calculations may not be fully converged (Table 1): while FEP can be accurate, it can also be limited by sampling issues which would bias results towards the simulation starting states (Coskun et al., Reference Coskun, Chen, Clark, Lu, Harder, Wang, Friesner and Miller2022).

Table 1. Varied pKa values were calculated for D407, D408 and K940 in the studies presented in the text

Methodology and choice of starting structure can lead to varied results. The highlighted residues are those of the key E state. Also included are values calculated for the AcrB structures 2DHH and 4DX5 using PROPKA v3.5 (Olsson et al., Reference Olsson2011; Søndergaard et al., Reference Søndergaard, Olsson, Rostkowski and Jensen2011).

Eicher et al. confirmed two protonation events would create an alternating access water wire, facilitating proton transfer to and from solution. Using unbiased all-atom MD, membrane-bound AcrB was simulated using the predicted protonation states. Access to the periplasmic and cytoplasmic solvent was found to alternate with conformational state (Fischer and Kandt, Reference Fischer and Kandt2011). Overall, the evidence presented in the study is consistent with a two-proton model.

The two-proton model has also been tested on a homologue of AcrB. A recent paper by Fairweather et al. investigated substrate binding and the effect of protonation on the RND protein MtrD (Neisseria gonorrhoeae, 4MT1; Bolla et al., Reference Bolla, Su, Do, Radhakrishnan, Kumar, Long, Chou, Delmar, Lei, Rajashankar, Shafer and Yu2014; Fairweather et al., Reference Fairweather, Gupta, Chitsaz, Booth, Brown and O’Mara2021). Studying the role of protonation, united atom (UA; Schmid et al., Reference Schmid, Eichenberger, Choutko, Riniker, Winger, Mark and van Gunsteren2011) conventional MD simulations were used, changing the protonation states of D405 and D406 (AcrB D407/D408 analogues) and the presence of progesterone substrate molecules. In the UA representation, five systems were simulated with varying protonation states and number of substrates. In the absence of substrates and protonation, MtrD adopted a symmetric A/A/A state – the putative resting state (Su et al., Reference Su, Li, Gu, Takatsuka, McDermott, Nikaido and Yu2006). Protonation of both aspartates in the absence of a substrate was sufficient to begin the opening of the exit gate (distance distributions with peaks at ~8.6 and ~ 12 Å, where the closed gate is ~5–6 Å) in the same monomer but caused instability at the relay site (K948-T985 distance distribution displays three peaks) and did not induce conformational cycling in neighbouring protomers.

Substrate binding accelerated transitions towards the asymmetric structure. Protonation in the extrusion protomer and binding of a single progesterone in the proximal binding pocket of an adjacent protomer stabilised an asymmetric state (reduced fluctuations in K948-T985 and exit site distances) but did not allow complete transition on the sampled timescales. The presence of progesterone in the binding pockets of both adjacent protomers allowed further asymmetry to develop: binding of substrates in two protomers, and protonation of the final protomer was required for the functional rotation of MtrD to be observed. Binding of multiple substrates has been shown to accelerate conformational cycling in other simulation studies (Wang et al., Reference Wang, Weng and Wang2015; Matsunaga et al., Reference Matsunaga, Yamane, Terada, Moritsugu, Fujisaki, Murakami, Ikeguchi and Kidera2018). Consistent with studies on AcrB, protonation of both aspartates induced rotation within the TMD, and the most extensive water pathways were observed in the binding protomer (Fischer and Kandt, Reference Fischer and Kandt2011; Eicher et al., Reference Eicher, Seeger, Anselmi, Zhou, Brandstätter, Verrey, Diederichs, Faraldo-Gómez and Pos2014). While this study yielded results suggesting protonation of both residues (and the presence of >1 substrate) can induce conformational cycling, the authors noted the proton relay network was never fully disengaged, despite protonation of D405 and D406. Two explanations were offered: the timescale of this transition may be beyond that sampled; or disengagement of the relay network may be driven by a different proton stoichiometry, that is, the two-proton model may not be appropriate for MtrD.

One proton

In 2013, Ikeguchi and colleagues used conventional all-atom MD simulations to investigate the possible protonation states of the titratable residues of AcrB (2DHH; (Murakami et al., Reference Murakami, Nakashima, Yamashita, Matsumoto and Yamaguchi2006; Yamane et al., Reference Yamane, Murakami and Ikeguchi2013). All possible combinations of protonation states for D407/D408 were tested in the E state protomer of this asymmetric structure. While each simulation was only 100 ns, they were able to observe differences in stability of the E-state depending on the protonation pattern. Their results suggested only D408H D407 stabilised the E-state with a disengaged relay network, as determined through inter-residue distances defined by K940, D407, D408, and T978. Protonation of only D407 led to immediate rearrangement of K940: the K940-T978 hydrogen bond was disrupted, and a salt bridge was formed between K940 and E408 which became unstable after 50 ns. D408H D407H was also classified as unstable as K940 became dynamic, switching interactions between D408 and T978 frequently, in addition to the salt bridge between K940 and E407 becoming unstable after 90 ns. The rearrangements and subsequent instability in both systems implies that protonation of D407 is incompatible with the E-state of this crystal structure.

They further probed these findings with two D408 D407 simulations from snapshots taken from the D407 D408H simulation. Deprotonation resulted in immediate re-engagement of the relay network, consistent with a transition to the A state. Principal component analysis of porter domain dynamics supported these results: D407D408H remained stable in the E state, whereas D407 D408 displayed a structural transition from E towards A. Together these results suggested that protonation of D408 alone stabilised the E state, and its deprotonation was sufficient to induce the transition to the A state. As stated by the authors, one weakness of their approach was the limited simulation time; the MD was unbiased, and consequently the timescale sampled was insufficient to explore full conformational transitions. We note that while all possible combinations of states were simulated, only one replicate of each system was reported.

In support of their conclusions, pKa values for the two aspartates were calculated using PROPKA (Li et al., Reference Li, Robertson and Jensen2005) and multiconformation continuum electrostatics (MCCE; Song et al., Reference Song, Mao and Gunner2009). Both methods predicted D408 to be more readily protonated than D407, and at pH 7 only D408 would be protonated. There was a disparity in the values predicted, dependent on the method and selected parameters (Table 1). The difference in values highlights the large effect dielectric constants can have on resultant pKas (Chan et al., Reference Chan, Wilbanks, Makhatadze and Wong2012).

Jewel et al. (Reference Jewel, Liu and Dutta2017) simulated the conformational changes of apo-AcrB utilising a hybrid coarse-grain regime. AcrB (asymmetric, 2DHH (Murakami et al., Reference Murakami, Nakashima, Yamashita, Matsumoto and Yamaguchi2006) was modelled in the UA PACE forcefield (Han et al., Reference Han, Wan, Jiang and Wu2010), and surrounding lipids and water molecules were represented in the MARTINI 2 forcefield (Marrink et al., Reference Marrink, Risselada, Yefimov, Tieleman and de Vries2007). Due to its coarse-grained nature, this regime can sample longer timescales than atomistic simulations, with the caveat of reduced resolution. Similar to the methodology of Yamane et al. (Reference Yamane, Murakami and Ikeguchi2013), this study explored the different combinations of permanent protonation states for D407/D408 in the extrusion protomer: both deprotonated, one protonated, or both protonated, for a total of four systems. When both residues were deprotonated in the E protomer, the relay network approached an engaged (Access) state: K940 moves away from T978 to form salt bridges with D407 and D408 within 100 ns. Extending the simulations to the microsecond timescale, this system also displayed closing of the exit channel (Y758–Q124 decreases from ~12 Å to ~5 Å) and opening of the periplasmic cleft, again consistent with an E → A transition.

Protonation of D408 was found to maintain the orientation of the residues observed in the extrusion state of the crystal structure (K940 was always closer to T978 than D408); the exit channel remained open; and the periplasmic cleft remained closed. This contrasts with systems where D407 was protonated: in both cases, the residues of the relay network were described as “rather dynamic”, and the exit channel opens to a greater extent than that observed in the crystal structure (>15 Å vs ~ 12 Å in the crystal structure). It was noted that this may be due to the coarse-graining scheme or the relatively large time step (5 fs), but may also indicate instability of the E-state with D407H. Taken together, these results are consistent with the atomistic study by Yamane et al. (Reference Yamane, Murakami and Ikeguchi2013): protonation of D407 destabilises the E state; protonation of D408 stabilises the E state; deprotonation of D408 drives the Extrusion to Access transition.

Following on from this, Jewel et al. (Reference Jewel, Dinh, Liu and Dutta2020) used the same hybrid coarse-graining regime to explore conformational changes induced by changes in protonation state and the presence of an indole substrate. Four systems were simulated, except this time protonation states were chosen to induce the B → E transition (rather than E → A). Three protonation schemes for the binding monomer were simulated with indole present in the binding site of that same monomer: D407 D408; D407 D408H; D407H D408. The fourth system had all D407/D408 residues deprotonated, and no substrate. Both systems with protonated aspartates displayed closing of the periplasmic cleft (measured by a large increase in the T676-F563 distance) and an opening of the exit channel (Q124 and Y758 distance increase to >12 Å), suggesting a B → E transition. Interestingly, the system with bound indole but no protonation of D407/D408 also displayed opening at the exit channel, though the periplasmic cleft remained open. In the absence of substrates or protonation, the final system tended towards the A/A/A state with a closed exit channel and open periplasmic cleft: the putative resting state of RND proteins (Su et al., Reference Su, Li, Gu, Takatsuka, McDermott, Nikaido and Yu2006). These results suggest that protonation of either D407 or D408 or the binding of a substrate could induce conformational cycling in the efflux pump.

Constant pH MD (CpHMD) has also been used to investigate these residues. Yue et al. (Reference Yue, Chen, Zgurskaya and Shen2017) used hybrid-solvent continuous CpHMD to investigate the protonation states and conformational dynamics of the TMD of AcrB. Only the TMD of AcrB (4DX5 (Eicher et al., Reference Eicher, Cha, Seeger, Brandstätter, El-Delik, Bohnert, Kern, Verrey, Grütter, Diederichs and Pos2012) was simulated to reduce computational cost and ensure adequate sampling. The pH replica exchange protocol (Wallace and Shen, Reference Wallace and Shen2011) was used, with 24 replicas over the pH range 1–8.5. While conventional MD assumes fixed protonation states which must be selected a priori, CpHMD allows titratable residues to be protonated/deprotonated over the course of the simulations according to the solution pH and local environment. CpHMD also enables the calculation of the pKa values of all titratable residues simultaneously, in contrast to traditional FEP methods which yield only one pKa at a time while fixing the protonation states of the other residues. In this case, CpHMD was coupled with replica exchange to achieve convergence of pKa values: by exchanging systems at varying pH values, each utilising different random walks, increased sampling of the potential energy surface of the system is achieved by accelerating the crossing of energy barriers.

Using this method, several observations were made on the pH-dependent behaviour of the titratable residues. Firstly, at pH < 4.0 the native E-state is maintained, but as pH increases a new conformational state (denoted O*/E*) is sampled: at low pH the proton relay network is disengaged (K940-T978 hydrogen bonding), but at higher pH, where both D407 and D408 are deprotonated, K940 adopts an orientation more consistent with the A-state. A large-scale motion accompanies this change: a lateral rotation within the TMD resembling the E → A transition, also observed by Eicher et al. (Reference Eicher, Seeger, Anselmi, Zhou, Brandstätter, Verrey, Diederichs, Faraldo-Gómez and Pos2014).

pKa calculations were consistent with D407/D408 in the A and B states (pKas of 2.2/3.0 in A; 2.6/2.8 in B). Due to their spatial proximity and similar pKas, protonation of these residues is strongly coupled and D407/D408 were considered as a dyad. Titration was therefore described by a stepwise model with two macroscopic pKas: the first/second protonation events in the B state had pKas of 3.4/2.1. Separation of 1.3 units suggests that binding of one proton is more likely than two in the B → E transition. In the E state, the pKas of the two protonation events increase to 5.0/3.4, with D408 having the greater pKa (D407 3.4 vs D408 5.0) and the greater increase in pKa (+0.8 for D407 vs + 2.2 for D408). While the authors note that their pKa values were systematically underestimated, the relative order of the values and the difference between them were considered robust: taken together the results suggest that only D408 accepts a proton in the B → E transition.

Free energy calculation methodologies were used by Matsunaga et al. (Reference Matsunaga, Yamane, Terada, Moritsugu, Fujisaki, Murakami, Ikeguchi and Kidera2018) to investigate the pathway of AcrB functional rotation. They hypothesised that D408 protonation/deprotonation drives the conformational cycling. Two systems were generated, both starting in the BEA state (4DX5 (Eicher et al., Reference Eicher, Cha, Seeger, Brandstätter, El-Delik, Bohnert, Kern, Verrey, Grütter, Diederichs and Pos2012)) with minocycline bound in the B protomer. These systems were driven towards the EAB state using targeted all-atom MD (Klauda et al., Reference Klauda, Venable, Freites, O’Connor, Tobias, Mondragon-Ramirez, Vorobyov, MacKerell and Pastor2010; Best et al., Reference Best, Zhu, Shim, Lopes, Mittal, Feig and Mackerell2012). In the first system the D408 of protomer II (E) was protonated, predicted to stabilise the initial BEA state. In the second system, D408 of protomer I (B) was protonated, predicted to induce functional rotation to the EAB state (Figure 3a) and hence stabilise the end state. Making use of the string method (Branduardi and Faraldo-Gómez, Reference Branduardi and Faraldo-Gómez2013) and umbrella sampling (Torrie and Valleau, Reference Torrie and Valleau1977), the minimum free energy pathway between the initial and final states was calculated. It was demonstrated that the local energy minimum for system 1 lay close to the initial BEA state while for system 2 the local minimum was closer to the EAB state, concluding that protonation of D408 of the B state could drive functional rotation to the E state.

To assess how the protonation state affected the free energy, system 1 was alchemically transformed towards system 2 by protonating D408 in protomer I and deprotonating D408 in protomer II. Due to restraints on collective variables (Cartesian coordinates of Cα atoms of the porter domain, TMHs accommodating the relay network, and a flexible loop between TMHs 5 and 6; side chain atoms of D407 and D408), free energy calculations were stated to evaluate the contribution of protonation to any free energy differences. By protonating protomer I (B), the free energy of the minimum was increased by 25.9 ± 0.5 kcal/mol. Electrostatic potential maps showed that the TMD of protomer I was destabilised upon protonation due to repulsion between D408H and proximal cationic residues (K940, R971). The free energy change on transforming the protonation state of system 2 into that of system 1 was also evaluated, but this had a reduced effect (0.4 ± 6.9 kcal/mol), suggesting protonation may no longer affect energetics at this stage in the cycle.

This study further related protonation to conformational changes in the TMD. Comparing the water distribution at the local energy minima of system 1 and system 2 showed the sheer vertical motion had created alternating access to the periplasm and cytoplasm: a possible proton release pathway consistent with previous findings (Fischer and Kandt, Reference Fischer and Kandt2011). The sheer motion and TM helix tilts were also correlated to porter domain opening/closing using mutual information analysis (McClendon et al., Reference McClendon, Friedland, Mobley, Amirkhani and Jacobson2009). This study does not discount the possibility that D407 may be protonated during conformational cycling. However, it does demonstrate that singular protonation of D408 in the B state can cause functional rotation of AcrB to a stable protonated E state, and that this functional rotation is translated to the porter domain.

Conclusions & outlook

Here we have presented a series of computational studies aiming to further elucidate details of the proton transfer mechanism in RND transporters. The literature from the last decade generally favours a one-proton model in which only D408 becomes protonated on transition from the Binding to the Extrusion state, and its deprotonation is sufficient to induce cycling back to the Access state. However, inconsistencies across these studies and conflicting experimental data remain, which cannot currently be reconciled. For example, we compared the pKas estimated by PROPKA v3.5 for 2DHH (Murakami et al., Reference Murakami, Nakashima, Yamashita, Matsumoto and Yamaguchi2006) and 4DX5 (Eicher et al., Reference Eicher, Cha, Seeger, Brandstätter, El-Delik, Bohnert, Kern, Verrey, Grütter, Diederichs and Pos2012), finding the E state D407 in particular had differing results (6.71 vs 7.86, respectively). Due to different crystallisation conditions and ligands present, it is hard to say whether the structures are within the same native ensemble or perhaps reflect different states (Best et al., Reference Best, Lindorff-Larsen, DePristo and Vendruscolo2006). In general, the pKas shown in Table 1 highlight the difficulty in determining protonation states, as each method carries with it its own biases (Coskun et al., Reference Coskun, Chen, Clark, Lu, Harder, Wang, Friesner and Miller2022; Wilson et al., Reference Wilson, Karttunen, de Groot and Gapsys2023). It was recently demonstrated that all commonly used high throughput pKa estimation approaches, including PROPKA and electrostatics-based methods, are prone to inaccuracies (Wei et al., Reference Wei, Hogues and Sulea2023), so it is not trivial to tease apart these details. Incorporating well-sampled ensemble states with accurate pKas will be essential to accurately assess protonation states. Additionally, some simulation studies have shown substrate binding to accelerate or even be required for full conformational cycling (Wang et al., Reference Wang, Weng and Wang2015; Fairweather et al., Reference Fairweather, Gupta, Chitsaz, Booth, Brown and O’Mara2021). The interdependence of protonation events and substrate recognition is yet to be fully understood. Finally, it is important to note that simulated systems will inevitably simplify true biological system complexity: studies presented here have simulated the RND transporter isolated in simple symmetric phospholipid bilayers. Of particular importance is the local environment. Not only will AcrB dynamics be affected by coupling to substrates, AcrA, and TolC in vivo, but the local proton gradient across the inner membrane will affect the likelihood of protonation for these essential residues. With the maturation of increasingly powerful experimental techniques and simulation studies able to capture greater system complexity, we anticipate further exploration of this family of proteins and their functional mechanism in the future.

Open peer review

To view the open peer review materials for this article, please visit http://doi.org/10.1017/qrd.2024.6.

Supplementary material

The supplementary material for this article can be found at http://doi.org/10.1017/qrd.2024.6.

Acknowledgements

R.C. is funded by the Medical Research Council (grant numbers MR/N013468/1 and MR/W006731/1); Magdalen College; and the Department of Biochemistry, Oxford. K.E.N. is supported by a Ph.D. Studentship from the Engineering and Physical Sciences Research Council (Project Number: 2446840) and travel funding from St Anne’s College, University of Oxford. S.K. is funded by EPSRC grant numbers EP/V030779 and EP/X035603.

Competing interest

The authors declare none.

Footnotes

R.C. and K.E.N. contributed equally to this study.

References

Alav, I, Kobylka, J, Kuth, MS, Pos, KM, Picard, M, Blair, JMA and Bavro, VN (2021) Structure, assembly, and function of tripartite efflux and type 1 secretion Systems in Gram-Negative Bacteria. Chemical Reviews 121(9), 54795596.CrossRefGoogle ScholarPubMed
Athar, M, Gervasoni, S, Catte, A, Basciu, A, Malloci, G, Ruggerone, P and Vargiu, AV (2023) Tripartite efflux pumps of the RND superfamily: What did we learn from computational studies? Microbiology 169(3), 001307.CrossRefGoogle ScholarPubMed
Best, RB, Lindorff-Larsen, K, DePristo, MA and Vendruscolo, M (2006) Relation between native ensembles and experimental structures of proteins. Proceedings of the National Academy of Sciences of the United States of America 103(29), 1090110906.CrossRefGoogle ScholarPubMed
Best, RB., Zhu, X, Shim, J, Lopes, PEM, Mittal, J, Feig, M and Mackerell, AD Jr. (2012) Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral angles. Journal of Chemical Theory and Computation 8(9), 32573273.CrossRefGoogle Scholar
Bolla, JR, Su, CC, Do, SV, Radhakrishnan, A, Kumar, N, Long, F, Chou, TH, Delmar, JA, Lei, HT, Rajashankar, KR, Shafer, WM and Yu, EW (2014) Crystal structure of the Neisseria gonorrhoeae MtrD inner membrane multidrug efflux pump. PLoS One 9(6), e97903.CrossRefGoogle ScholarPubMed
Branduardi, D and Faraldo-Gómez, JD (2013) String method for calculation of minimum free-energy paths in Cartesian space in freely tumbling systems. Journal of Chemical Theory and Computation 9(9), 41404154.CrossRefGoogle ScholarPubMed
Brooks, BR, Brooks, CL, Mackerell, AD, Nilsson, L, Petrella, RJ, Roux, B, Won, Y, Archontis, G, Bartels, C, Boresch, S, Caflisch, A, Caves, L, Cui, Q, Dinner, AR, Feig, M, Fischer, S, Gao, J, Hodoscek, M, Im, W, Kuczera, K, Lazaridis, T, Ma, J, Ovchinnikov, V, Paci, E, Pastor, RW, Post, CB, Pu, JZ, Schaefer, M, Tidor, B, Venable, RM, Woodcock, HL, Wu, X, Yang, W, York, DM and Karplus, M et al. (2009) CHARMM: The biomolecular simulation program. Journal of Computational Chemistry 30(10), 15451614.CrossRefGoogle ScholarPubMed
Chan, CH, Wilbanks, CC, Makhatadze, GI and Wong, KB (2012) Electrostatic contribution of surface charge residues to the stability of a thermophilic protein: Benchmarking experimental and predicted pKa values. PLoS One 7(1), e30296.CrossRefGoogle Scholar
Coskun, D, Chen, W, Clark, AJ, Lu, C, Harder, ED, Wang, L, Friesner, RA and Miller, EB (2022) Reliable and accurate prediction of single-residue pKa values through free energy perturbation calculations. Journal of Chemical Theory and Computation 18(12), 71937204.CrossRefGoogle ScholarPubMed
Du, D, Wang, Z, James, NR, Voss, JE, Klimont, E, Ohene-Agyei, T, Venter, H, Chiu, W and Luisi, BF (2014) Structure of the AcrAB–TolC multidrug efflux pump. Nature 509(7501), 512515.CrossRefGoogle ScholarPubMed
Eicher, T, Cha, HJ, Seeger, MA, Brandstätter, L, El-Delik, J, Bohnert, JA, Kern, WV, Verrey, F, Grütter, MG, Diederichs, K and Pos, KM (2012). Transport of drugs by the multidrug transporter AcrB involves an access and a deep binding pocket that are separated by a switch-loop. Proceedings of the National Academy of Sciences of the United States of America 109(15), 56875692.CrossRefGoogle Scholar
Eicher, T, Seeger, MA, Anselmi, C, Zhou, W, Brandstätter, L, Verrey, F, Diederichs, K, Faraldo-Gómez, JD and Pos, KM (2014) Coupling of remote alternating-access transport mechanisms for protons and substrates in the multidrug efflux pump AcrB. eLife 3, e03145.CrossRefGoogle ScholarPubMed
Fairweather, SJ, Gupta, V, Chitsaz, M, Booth, L, Brown, MH and O’Mara, ML. (2021). Coordination of substrate binding and protonation in the MtrD efflux pump controls the functionally rotating transport mechanism. ACS Infectious Diseases 7(6), 18331847.CrossRefGoogle ScholarPubMed
Fischer, N and Kandt, C. (2011). Three ways in, one way out: Water dynamics in the trans-membrane domains of the inner membrane translocase AcrB. Proteins 79(10), 28712885.CrossRefGoogle ScholarPubMed
Fischer, N and Kandt, C.. (2013). Porter domain opening and closing motions in the multi-drug efflux transporter AcrB. Biochimica et Biophysica Acta 1828(2), 632641.CrossRefGoogle ScholarPubMed
Guan, L and Nakae, T. (2001). Identification of essential charged residues in transmembrane segments of the multidrug transporter MexB of pseudomonas aeruginosa. Journal of Bacteriology 183(5), 17341739.CrossRefGoogle ScholarPubMed
Han, W, Wan, CK, Jiang, F and Wu, YD (2010) PACE force field for protein simulations. 1. Full parameterization of version 1 and verification. Journal of Chemical Theory and Computation 6(11), 33733389.CrossRefGoogle ScholarPubMed
Jewel, Y, Liu, J and Dutta, P (2017) Coarse-grained simulations of conformational changes in the multidrug efflux transporter AcrB. Molecular BioSystems 13(10), 20062014.CrossRefGoogle ScholarPubMed
Jewel, Y, Dinh, QV, Liu, J and Dutta, P. (2020). Substrate-dependent transport mechanism in AcrB of multidrug resistant bacteria. Proteins 88(7), 853864.CrossRefGoogle ScholarPubMed
Klauda, JB, Venable, RM, Freites, JA, O’Connor, JW, Tobias, DJ, Mondragon-Ramirez, C, Vorobyov, I, MacKerell, AD Jr and Pastor, RW. (2010) Update of the CHARMM all-atom additive force field for lipids: Validation on six lipid types. The Journal of Physical Chemistry. B 114(23), 78307843.CrossRefGoogle ScholarPubMed
Li, H, Robertson, AD and Jensen, JH. (2005). Very fast empirical prediction and rationalization of protein pKa values. Proteins 61(4), 704721.CrossRefGoogle ScholarPubMed
Marrink, SJ, Risselada, HJ, Yefimov, S, Tieleman, DP and de Vries, AH (2007) The MARTINI force field: Coarse grained model for biomolecular simulations. The Journal of Physical Chemistry B 11(27), 78127824.CrossRefGoogle Scholar
Matsunaga, Y, Yamane, T, Terada, T, Moritsugu, K, Fujisaki, H, Murakami, S, Ikeguchi, M and Kidera, A (2018) Energetics and conformational pathways of functional rotation in the multidrug transporter AcrB. eLife 7, e31715.CrossRefGoogle ScholarPubMed
McClendon, CL, Friedland, G, Mobley, DL, Amirkhani, H and Jacobson, MP (2009) Quantifying correlations between allosteric sites in thermodynamic ensembles. Journal of Chemical Theory and Computation 5(9), 24862502.CrossRefGoogle ScholarPubMed
Murakami, S, Nakashima, R, Yamashita, E, Matsumoto, T and Yamaguchi, A (2006) Crystal structures of a multidrug transporter reveal a functionally rotating mechanism. Nature 443(7108), 173179.CrossRefGoogle ScholarPubMed
Murakami, S and Yamaguchi, A (2003) Multidrug-exporting secondary transporters. Current Opinion in Structural Biology 13(4), 443452.CrossRefGoogle ScholarPubMed
Olsson, MHM (2011) Protein electrostatics and pKa blind predictions; contribution from empirical predictions of internal ionizable residues. Proteins: Structure, Function, and Bioinformatics 79(12), 33333345. https://doi.org/10.1002/prot.23113CrossRefGoogle ScholarPubMed
Schmid, N, Eichenberger, AP, Choutko, A, Riniker, S, Winger, M, Mark, AE and van Gunsteren, WF (2011) Definition and testing of the GROMOS force-field versions 54A7 and 54B7. European Biophysics Journal: EBJ 40(7), 843856.CrossRefGoogle ScholarPubMed
Seeger, MA., von Ballmoos, C, Verrey, F and Pos, KM (2009). Crucial role of Asp408 in the proton translocation pathway of multidrug transporter AcrB: Evidence from site-directed mutagenesis and Carbodiimide labeling. Biochemistry 48(25), 58015812.CrossRefGoogle ScholarPubMed
Seeger, MA., Diederichs, K, Eicher, T, Brandstätter, L, Schiefner, A, Verrey, F and Pos, KM (2008) The AcrB efflux pump: Conformational cycling and peristalsis lead to multidrug resistance. Current Drug Targets 9(9), 729749.CrossRefGoogle ScholarPubMed
Søndergaard, CR, Olsson, MHM, Rostkowski, M and Jensen, JH (2011) Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. Journal of Chemical Theory and Computation 7(7), 22842295.CrossRefGoogle ScholarPubMed
Song, Y, Mao, J and Gunner, MR (2009) MCCE2: Improving protein pKa calculations with extensive side chain Rotamer sampling. Journal of Computational Chemistry 30(14), 22312247.CrossRefGoogle ScholarPubMed
Su, CC, Li, M, Gu, R, Takatsuka, Y, McDermott, G, Nikaido, H and Yu, EW (2006) Conformation of the AcrB multidrug efflux pump in mutants of the putative proton relay pathway. Journal of Bacteriology 188(20), 72907296.CrossRefGoogle ScholarPubMed
Takatsuka, Y and Nikaido, H (2006) Threonine-978 in the transmembrane segment of the multidrug efflux pump AcrB of Escherichia Coli is crucial for drug transport as a probable component of the proton relay network. Journal of Bacteriology 188(20), 72847289.CrossRefGoogle ScholarPubMed
Torrie, GM and Valleau, JP (1977) Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. Journal of Computational Physics 23(2), 187199.CrossRefGoogle Scholar
Vargiu, AV and Nikaido, H (2012) Multidrug binding properties of the AcrB efflux pump characterized by molecular dynamics simulations. Proceedings of the National Academy of Sciences of the United States of America 109(50), 2063720642.CrossRefGoogle ScholarPubMed
Wallace, JA and Shen, JK (2011) Continuous constant pH molecular dynamics in explicit solvent with pH-based replica exchange. Journal of Chemical Theory and Computation 7(8), 26172629.CrossRefGoogle ScholarPubMed
Wang, B, Weng, J and Wang, W (2015) Substrate binding accelerates the conformational transitions and substrate dissociation in multidrug efflux transporter AcrB. Frontiers in Microbiology 6, 302.Google ScholarPubMed
Wei, W, Hogues, H and Sulea, T (2023) Comparative performance of high-throughput methods for protein P predictions. Journal of Chemical Information and Modeling 63(16), 51695181.CrossRefGoogle Scholar
Wilson, CJ, Karttunen, M, de Groot, BL and Gapsys, V (2023) Accurately predicting protein pKa values using nonequilibrium alchemy. Journal of Chemical Theory and Computation 19(21), 78337845.CrossRefGoogle ScholarPubMed
Yamane, T, Murakami, S and Ikeguchi, M (2013) Functional rotation induced by alternating protonation states in the multidrug transporter AcrB: All-atom molecular dynamics simulations. Biochemistry 52(43), 76487658.CrossRefGoogle ScholarPubMed
Yue, Z, Chen, W, Zgurskaya, HI and Shen, J (2017) Constant pH molecular dynamics reveals how proton release drives the conformational transition of a transmembrane efflux pump. Journal of Chemical Theory and Computation 13(12), 64056414.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. AcrB forms part of the tripartite multidrug efflux pump, AcrAB-TolC. (a) Surface representation of AcrAB-TolC (PDB ID: 5O66; Du et al.,2014) and its situation within the cell envelope of Gram-negative bacteria. AcrB sits within the inner membrane; TolC sits within the outer membrane; and AcrA is the periplasmic adaptor protein. Substrates, such as certain classes of antibiotics, are thought to enter AcrB either in the periplasm near the inner membrane, or from the outer leaflet of the inner membrane directly. Protons enter AcrB from a different entry point, and transit through the TMD to the cytoplasm. The precise mechanism of the proton relay is unclear. TolC is coloured in purple, AcrA is coloured in orange, AcrB is coloured according to state. The Access state is coloured pink, the Binding state is coloured yellow, and the Extrusion state is coloured in blue. (b) Cartoon representation of AcrB (PDB ID: 4DX5; Eicher et al., 2012). AcrB has three main domains: the Docking domain which interacts with AcrA; the Porter domain contains the drug entry site; and the Transmembrane domain uses the proton motive force to induce drug extrusion. The titratable residues thought to orchestrate proton transfer are shown in stick representation and coloured according to element.

Figure 1

Figure 2. Cartoon schematic of functional rotation and substrate transport through inner membrane RND transporter proteins. (a) Side-view of the RND protein in the inner membrane. (b) Top-down (from the periplasm) view of the porter domain. For visual clarity, protonation and drug extrusion are considered for a single protomer only. Substrates enter the proximal binding pocket of the Access protomer from the periplasm/periplasmic leaflet of the inner membrane. Substrate binding induces a conformational change in the protomer to the Binding state, and the substrate moves to the distal binding pocket. Protonation then occurs in the relay in the transmembrane domain, inducing a conformational change to the Extrusion state. The periplasmic cleft closes and the exit gate opens, allowing the substrate to exit into the periplasmic adaptor protein (not shown).

Figure 2

Figure 3. AcrB undergoes a conformational cycle that functionally connects the movement of protons through the protein to drug extrusion. (a) Porter domains of AcrB viewed from the periplasm, showing the different conformational changes that correspond to drug access, binding and extrusion. The Access state is coloured pink, the Binding state is coloured yellow, and the Extrusion state is coloured in blue. Adapted from Matsunaga et al. (2018). (b) Drugs are thought to move through the proximal binding pocket (purple) in the Access state, bind to the distal pocket (dark blue) in the Binding state, and leave through the exit gate (cyan) in the Extrusion state. A drug (minocycline), is shown moving through these different regions. The drug was placed manually in the Access and Extrusion states for illustrative purposes, but the position of the drug in the binding state comes from PDB ID:4DX5 (Eicher et al., 2012). The definition of the different regions comes from Vargiu and Nikaido (2012).

Figure 3

Figure 4. The titratable residues of AcrB adopt different configurations in different states. (a) Close-up view of key residues of the proton relay. Each image is rotated, so viewed from the same orientation. In the Access and Binding states the relay network is considered ‘engaged’, and in the Extrusion states the network is considered ‘disengaged’. Adapted from Matsunaga et al. (2018). (b) Inter-residue distances for residues in the proton relay network of Gram-negative bacterial RND protein structures deposited in the PDB (raw data and relevant references available in the Supplementary material). Distances were calculated as the shortest distance between side chain oxygen/nitrogen atoms, except for G403 where the backbone carbonyl oxygen was used. (c) Conformational transition from the Binding to the Extrusion state is proposed to proceed via one protonation event (D407/D408H) or two (D407H/D408H). Structures shown in (a) and (c) are of PDB ID: 4DX5.

Figure 4

Table 1. Varied pKa values were calculated for D407, D408 and K940 in the studies presented in the text

Supplementary material: File

Clark et al. supplementary material

Clark et al. supplementary material
Download Clark et al. supplementary material(File)
File 40 KB

Author comment: Titratable residues that drive RND efflux: Insights from molecular simulations — R0/PR1

Comments

Dear Editors,

Please find attached, a perspective/review on the insights gained from computational studies on the role of the ionisation states of titratable residues in bacterial RND efflux pumps. While the mechanism of action of these pumps has been investigated using many experimental and simulation methods over the years- we feel for this special edition of QRB Discovery, the aspect we have chosen is particularly pertinent as the calculations have provided insights that really would not have been possible with other methods.

Thank you for considering our manuscript.

Yours faithfully,

Syma Khalid

Oxford

Decision: Titratable residues that drive RND efflux: Insights from molecular simulations — R0/PR2

Comments

Minor revisions

Author comment: Titratable residues that drive RND efflux: Insights from molecular simulations — R1/PR3

Comments

We thank the reviewers for their constructive comments. We have modified the manuscript in accordance with their suggestions, and below we have responded directly to each point they raised.

Reviewer: 1

Comments to the Author

The manuscript provides an insightful review and commentary on the topic of protonation linked with drug transport in the RND family of transmembrane transporters. The focus is mostly on the AcrB homotrimer, for which structural data and simulation results are most advanced. This work presents a very nice overview of the most relevant computational research papers and development in the field. The paper is generally well written and structured and would be of interest to people in the community and the wider field.

There are only a few minor points which will hopefully be useful for the authors to consider:

Minor comments

1. Lines 103-15: “predicted that in the E state, D408 would be protonated 100% of the time and D407 90%”. Can it be briefly explained how this was devised from the simulations in the cited paper?

A brief description of the methodology used has been added (lines 100 to 106).

2. Lines 159-161: Please paraphrase instead of directly quoting (not very common for scientific papers).

We have changed the wording of this section to paraphrase the results rather than quoting the paper.

3. Overall, the content and information is nicely summarized. One minor consideration would be the writing style preference in some sections. The authors addressed frequently the findings from the papers by “e.g. line 117 their findings’’, “e.g. Line 161 the authors’’ , “the authors noted”. It may be worth rephrasing some of these in passive third voice as it will read easier.

We have changed the wording to the passive third voice in the majority of these instances.

Not a critical point, but something for the authors to reflect on is the puzzle of the apparent broad substrates recognition of the AcrB system, although it does seem to have preferences. Could this be something that simulation might explain in the future? Is there a process like recognition in a transport transition state that might be linked with proton movement?

We thank the reviewer for this thought-provoking comment. We note that there have been studies that support the idea that protonation and substrate presence/recognition are linked in some way. For example, the study of MtrD (Fairweather et al, 2021) showed that coordinated substrate binding (in two protomers) and protonation at the relay site were required for the conformation to change towards an asymmetric state. Additionally, Wang, Weng, and Wang (2015) presented simulations that suggest substrate presence accelerates the conformational cycling process. While it is not currently clear exactly how substrate recognition and protonation are linked, simulation studies have shown that there is some dependence between the two. An additional comment has been made in the conclusions to address this point (lines 342-346).

_____________________________________________________________________________________

Reviewer: 2

Comments to the Author

The mini-review by Clark, Newman, and Khalid summarizes several past studies applying molecular simulations to establish a mechanistic connection between proton transfer and conformational cycling in the AcrB efflux pump. The authors mostly focus on proton stoichiometry, critically reviewing studies that favor either a one-proton or two-proton pathway, and highlighting the former is better but not conclusively supported by the reports.

The article is well-written and reasonably comprehensive within its defined scope. My only major complaint concerns the somewhat scarce figures, leaving much of the mechanical discussion unvisualized, and perhaps hard to digest for a reader with little experience with the system in question.

Below I’m listing specific issues or questions for the Authors:

- The proximal/distal binding pockets, as well as the exit gate, should be visualized or labeled if the discussion of the mechanism is to give the reader a proper structural context.

We thank the reviewer for their feedback on the figures. In response, we have split up existing figures, generated additional panels, and added a new figure. What is now Figure 3 addresses this particular point. Panel B displays a single protomer of AcrB, with the residues lining the two pockets and the exit gate shown to add context to discussion around domain motions, binding pockets, and functional rotation.

- The notation D408H/D407H in Fig 2 should be revised so as not to be confused with aspartate-histidine mutations.

The notation in this figure (now Figure 4, panel C) has been changed: the H that indicates protonation is now superscript.

- For the comparison of pKas predicted by PROPKA v3.5 for 2DHH and 4DX5, it would be interesting to know the sensitivity of the prediction to equilibrium fluctuations, i.e. whether the two numbers reasonably correspond to samples from the same equilibrium ensemble, or indicate some meaningful structural difference between the experimentally resolved states.

In the conclusion we have made additional comments to address these points (lines 334-336, 338-341). 2DHH and 4DX5 did use different crystallisation conditions and have different ligands present, and the resolution of the two structures is different. We also note that pKa estimation approaches discussed in this review, including PROPKA and electrostatics-based methods, are prone to inaccuracies; it is challenging to know exactly from where the differences originate.

- When the phrase “functional rotation” is introduced, some background, perhaps visual, should be first given on how the conformational transitions entail a (partial) rotation of the subunits; this is to avoid readers thinking about truly rotating elements in the complex.

We have added an additional figure (now Figure 2) to summarise the functional rotation mechanism in cartoon form, with protomers labelled to make it clear that each individual protomer cycles the A, B, and E states rather than the whole complex rotating.

- When mentioning “restraints on collective variables”, given how unphysiologically high the reported free energy difference is, it would be interesting to briefly mention which CVs were constrained (rotamers? distances?).

We have added a brief summary of the collective variables selected in that paper on lines 304 to 306.

- Overall, it would be nice to include a separate figure with a summary of proposed mechanisms, perhaps accommodating all above requests for additional visualizations or labels.

As discussed, addressing the above comments, we have added an additional figure and panels to existing figures. Figure 3B displays the binding pockets and exit gates for drug extrusion. Figure 2 summarises the functional rotation mechanism in cartoon form, with proton transfer indicated for the protomer of interest. Figure 4C displays the proposed protonation states.

Decision: Titratable residues that drive RND efflux: Insights from molecular simulations — R1/PR4

Comments

No accompanying comment.