Resistance-Nodulation-cell Division (RND) transporters AcrB and AcrD of Escherichia coli expel a wide range of substrates out of the cell in conjunction with AcrA and TolC, contributing to the onset of bacterial multidrug resistance. Despite sharing an overall sequence identity of ~66% (similarity ~80%), these RND transporters feature distinct substrate specificity patterns whose underlying basis remains elusive. We performed exhaustive comparative analyses of the putative substrate binding pockets considering crystal structures, homology models and conformations extracted from multi-copy μs-long molecular dynamics simulations of both AcrB and AcrD. The impact of physicochemical and topographical properties (volume, shape, lipophilicity, electrostatic potential, hydration and distribution of multi-functional sites) within the pockets on their substrate specificities was quantitatively assessed. Differences in the lipophilic and electrostatic potentials among the pockets were identified. In particular, the deep pocket of AcrB showed the largest lipophilicity convincingly pointing out its possible role as a lipophilicity-based selectivity filter. Furthermore, we identified dynamic features (not inferable from sequence analysis or static structures) such as different flexibilities of specific protein loops that could potentially influence the substrate recognition and transport profile. Our findings can be valuable for drawing structure (dynamics)-activity relationship to be employed in drug design.

Secondary multidrug transporters of the resistance-nodulation-cell division (RND) superfamily contribute crucially to antibiotic resistance in Gram-negative bacteria. Compared to the most studied transporter AcrB of Escherichia coli, little is known about the molecular determinants of distinct polyspecificities of the most important RND transporters MexB and MexY of Pseudomonas aeruginosa. In an effort to add knowledge on this topic, we performed an exhaustive atomic-level comparison of the main putative recognition sites (access and deep binding pockets) in these two Mex transporters. We identified an underlying link between some structural, chemical and dynamical features of the binding pockets and the physicochemical nature of the corresponding substrates recognized by either one or both pumps. In particular, mosaic-like lipophilic and electrostatic surfaces of the binding pockets provide for both proteins several multifunctional sites for diffuse binding of diverse substrates. Specific lipophilicity signatures of the weakly conserved deep pocket suggest a key role of this site as a selectivity filter as in Acr transporters. Finally, the different dynamics of the bottom-loop in MexB and MexY support its possible role in binding of large substrates. Our work represents the first comparative study of the major RND transporters in P. aeruginosa and also the first structure-based study of MexY, for which no experimental structure is available yet.

We provided a mechanistic understanding of a fundamental process related to multi-drug transport. Our results can help rationalizing the behavior of other polyspecific transporters and designing compounds avoiding extrusion or inhibitors of efflux pumps.

Secondary multidrug (Mdr) transporters utilize ion concentration gradients to actively remove antibiotics and other toxic compounds from cells. The model Mdr transporter MdfA from Escherichia coli exchanges dissimilar drugs for protons. The transporter should open at the cytoplasmic side to enable access of drugs into the Mdr recognition pocket. Here we show that the cytoplasmic rim around the Mdr recognition pocket represents a previously overlooked important regulatory determinant in MdfA. We demonstrate that increasing the positive charge of the electrically asymmetric rim dramatically inhibits MdfA activity and sometimes even leads to influx of planar, positively charged compounds, resulting in drug sensitivity. Our results suggest that unlike the mutants with the electrically modified rim, the membrane-embedded wild-type MdfA exhibits a significant probability of an inward-closed conformation, which is further increased by drug binding. Since MdfA binds drugs from its inward-facing environment, these results are intriguing and raise the possibility that the transporter has a sensitive, drug-induced conformational switch, which favors an inward-closed state.

The putative mechanism by which bacterial RND-type multidrug efflux pumps recognize and transport their substrates is a complex and fascinating enigma of structural biology. How a single protein can recognize a huge number of unrelated compounds and transport them through one or just a few mechanisms is an amazing feature not yet completely unveiled. The appearance of cooperativity further complicates the understanding of structure-dynamics-activity relationships in these complex machineries. Experimental techniques may have limited access to the molecular determinants and to the energetics of key processes regulating the activity of these pumps. Computer simulations are a complementary approach that can help unveil these features and inspire new experiments. Here we review recent computational studies that addressed the various molecular processes regulating the activity of RND efflux pumps.