Molecular Mechanics Study of Flow and Surface Influence in Ligand-Protein Association

Recognition of quinolone antibiotics by the multidrug efflux transporter MexB of Pseudomonas aeruginosa

The drug/proton antiporter MexB is the engine of the major efflux pump MexAB-OprM in Pseudomonas aeruginosa. This protein is known to transport a large variety of compounds, including antibiotics, thus conferring a multi-drug resistance phenotype. Due to the difficulty of producing co-crystals, only two X-ray structures of MexB in a complex with ligands are available to date, and mechanistic aspects are largely hypothesized based on the body of data collected for the homologous protein AcrB of Escherichia coli. In particular, a recent study (Ornik-Cha, Wilhelm, Kobylka et al., Nat. Commun., 2021, 12, 6919) reported a co-crystal structure of AcrB in a complex with levofloxacin, an antibiotic belonging to the important class of (fluoro)-quinolones. In this work, we performed a systematic ensemble docking campaign coupled to the cluster analysis and molecular-mechanics optimization of docking poses to study the interaction between 36 quinolone antibiotics and MexB. We additionally investigated surface complementarity between each molecule and the transporter and thoroughly assessed the computational protocol adopted against the known experimental data. Our study reveals different binding preferences of the investigated compounds towards the sub-sites of the large deep binding pocket of MexB, supporting the hypothesis that MexB substrates oscillate between different binding modes with similar affinity. Interestingly, small changes in the molecular structure translate into significant differences in MexB–quinolone interactions. All the predicted binding modes are available for download and visualization at the following link: https://www.dsf.unica.it/dock/mexb/quinolones.

Putative binding modes (BMs) of quinolones to the bacterial efflux transporter MexB were identified. Multiple interaction patterns are possible, supporting the hypothesis that substrates oscillate between different BMs with similar affinity.

GeomBD3: Brownian Dynamics Simulation Software for Biological and Engineered Systems

GeomBD3 is a robust Brownian dynamics simulation package designed to easily handle natural or engineered systems in diverse environments and arrangements. The software package described herein allows users to design, execute, and analyze BD simulations. The simulations use all-atom, rigid molecular models that diffuse according to overdamped Langevin dynamics and interact through electrostatic, Lennard-Jones, and ligand desolvation potentials. The program automatically calculates molecular association rates, surface residence times, and association statistics for any number of user-defined association criteria. Users can also extract molecular association pathways, diffusion coefficients, intermolecular interaction energies, intermolecular contact probability maps, and more using the provided supplementary analysis scripts. We detail the use of the package from start to finish and apply it to a protein-ligand system and a large nucleic acid biosensor. GeomBD3 provides a versatile tool for researchers from various disciplines that can aid in rational design of engineered systems or play an explanatory role as a complement to experiments. GeomBD version 3 is available on our website at http://chemcha-gpu0.ucr.edu/geombd3/ and KBbox at https://kbbox.h-its.org/toolbox/methods/molecular-simulation/geombd/.