Computer simulations suggest direct and stable tip to tip interaction between the outer membrane channel TolC and the isolated docking domain of the multidrug RND efflux transporter AcrB.

One way by which bacteria achieve antibiotics resistance is preventing drug access to its target molecule for example through an overproduction of multi-drug efflux pumps of the resistance nodulation division (RND) protein super family of which AcrAB-TolC in Escherichia coli is a prominent example. Although representing one of the best studied efflux systems, the question of how AcrB and TolC interact is still unclear as the available experimental data suggest that either both proteins interact in a tip to tip manner or do not interact at all but are instead connected by a hexamer of AcrA molecules. Addressing the question of TolC-AcrB interaction, we performed a series of 100 ns - 1 µs-molecular dynamics simulations of membrane-embedded TolC in presence of the isolated AcrB docking domain (AcrB(DD)). In 5/6 simulations we observe direct TolC-AcrB(DD) interaction that is only stable on the simulated time scale when both proteins engage in a tip to tip manner. At the same time we find TolC opening and closing freely on extracellular side while remaining closed at the inner periplasmic bottleneck region, suggesting that either the simulated time is too short or additional components are required to unlock TolC.

Efflux pump-mediated antibiotics resistance: insights from computational structural biology.

The continuous rise of bacterial resistance against formerly effective pharmaceuticals is a major challenge for biomedical research. Since the first computational studies published seven years ago the simulation-based investigation of antibiotics resistance mediated by multidrug efflux pumps of the resistance nodulation division (RND) protein super family has grown into a vivid field of research. Here we review the employment of molecular dynamics computer simulations to investigate RND efflux pumps focusing on our group's recent contributions to this field studying questions of energy conversion and substrate transport in the inner membrane antiporter AcrB in Escherichia coli as well as access regulation and gating mechanism in the outer membrane efflux ducts TolC and OprM in E. coli and Pseudomonas aeruginosa.

Unilateral access regulation: ground state dynamics of the Pseudomonas aeruginosa outer membrane efflux duct OprM.

Acting as an efflux duct in the MexA-MexB-OprM multidrug efflux pump, OprM plays a major role in the antibiotic resistance capability of Pseudomonas aeruginosa, trafficking substrates through the outer cell membrane. Whereas the available crystal structures showed restricted OprM access on both ends, the underlying gating mechanism is not yet fully understood. To gain insight into the functional mechanism of OprM access regulation, we conducted a series of five independent, unbiased molecular dynamics simulations, computing 200 ns dynamics samples of the wild-type protein in a phospholipid membrane/150 mM NaCl water environment. On the extracellular side, OprM opens and closes freely under the simulated conditions, suggesting the absence of a gating mechanism on this side of the isolated protein. On the periplasmic side, we observe an opening of the tip regions at Val408 and to a lesser degree Asp416 located 1.5 nm further into the channel, leading to OprM end conformations being up to 3 and 1.4 times, respectively, more open than the asymmetric crystal structure. If our simulations are correct, our findings imply that periplasmic gating involves only the Asp416 region and that in vivo additional components, absent in our simulation, might be required for periplasmic gating if the observed opening trend near Asp416 is not negligible. In addition to that ,we identified in each monomer a previously unreported sodium binding site in the channel interior coordinated by Asp171 and Asp230 whose functional role remains to be investigated.

Locked on one side only: ground state dynamics of the outer membrane efflux duct TolC.

Playing a major role in the expulsion of antibiotics and the secretion of cell toxins in conjunction with inner membrane transporters of three protein superfamilies, the outer membrane channel TolC occurs in at least two states blocking or permitting the passage of substrates. The details of the underlying gating mechanism are not fully understood. Addressing the questions of extracellular access control and periplasmic gating mechanism, we conducted a series of independent, unbiased 150-300 ns molecular dynamics simulations of wild-type TolC in a phospholipid membrane/150 mM NaCl water environment. We find that TolC opens and closes freely on the extracellular side, suggesting the absence of a gating mechanism on this side in the isolated protein. On the periplasmic side, we observe the outer periplasmic bottleneck region adopting in all simulations a conformation more open than the TolC wild-type crystal structures until in one run the successive binding of two sodium ions induces the transition to a conformation more closed than any of the available TolC X-ray structures. Concurrent with a heightened sodium residence probability near Asp374, the inner periplasmic bottleneck region at Asp374 remains closed throughout the simulations unless all NaCl is removed from the system, inducing a reopening of the outer and inner bottleneck. Our findings suggest that TolC is locked only on the periplasmic side in a sodium-dependent manner.

dxTuber: detecting protein cavities, tunnels and clefts based on protein and solvent dynamics.

Empty space in a protein structure can provide valuable insight into protein properties such as internal hydration, structure stabilization, substrate translocation, storage compartments or binding sites. This information can be visualized by means of cavity analysis. Numerous tools are available depicting cavities directly or identifying lining residues. So far, all available techniques base on a single conformation neglecting any form of protein and cavity dynamics. Here we report a novel, grid-based cavity detection method that uses protein and solvent residence probabilities derived from molecular dynamics simulations to identify (I) internal cavities, (II) tunnels or (III) clefts on the protein surface. Driven by a graphical user interface, output can be exported in PDB format where cavities are described as individually selectable groups of adjacent voxels representing regions of high solvent residence probability. Cavities can be analyzed in terms of solvent density, cavity volume and cross-sectional area along a principal axis. To assess dxTuber performance we performed test runs on a set of six example proteins representing the three main classes of protein cavities and compared our findings to results obtained with SURFNET, CAVER and PyMol.