Folding and self-assembly of the TatA translocation pore based on a charge zipper mechanism.

We propose a concept for the folding and self-assembly of the pore-forming TatA complex from the Twin-arginine translocase and of other membrane proteins based on electrostatic "charge zippers." Each subunit of TatA consists of a transmembrane segment, an amphiphilic helix (APH), and a C-terminal densely charged region (DCR). The sequence of charges in the DCR is complementary to the charge pattern on the APH, suggesting that the protein can be "zipped up" by a ladder of seven salt bridges. The length of the resulting hairpin matches the lipid bilayer thickness, hence a transmembrane pore could self-assemble via intra- and intermolecular salt bridges. The steric feasibility was rationalized by molecular dynamics simulations, and experimental evidence was obtained by monitoring the monomer-oligomer equilibrium of specific charge mutants. Similar "charge zippers" are proposed for other membrane-associated proteins, e.g., the biofilm-inducing peptide TisB, the human antimicrobial peptide dermcidin, and the pestiviral E(RNS) protein.

Supramolecular Hydrogels and Water Channels of Differing Diameters from Dipeptide Isomers.

Dipeptides stereoisomers and regioisomers composed of norleucine (Nle) and phenylalanine (Phe) self-assemble into hydrogels under physiological conditions that are suitable for cell culture. The supramolecular behavior, however, differs as the packing modes comprise amphipathic layers or water channels, whose diameter is defined by either four or six dipeptide molecules. A variety of spectroscopy, microscopy, and synchrotron-radiation-based techniques unveil fine details of intermolecular interactions that pinpoint the relationship between the chemical structure and ability to form supramolecular architectures that define soft biomaterials.

Single-atom substitution enables supramolecular diversity from dipeptide building blocks.

Dipeptides are popular building blocks for supramolecular gels that do not persist in the environment and may find various applications. In this work, we show that a simple substitution on the aromatic side-chain of phenylalanine with either fluorine or iodine enables supramolecular diversity upon self-assembly at neutral pH, leading to hydrogels or crystals. Each building block is characterized by 1H- and 13C-NMR spectroscopy, LC-MS, circular dichroism, and molecular models. The supramolecular behaviour is monitored with a variety of techniques, including circular dichroism, oscillatory rheology, transmission electron microscopy, attenuated total reflectance Fourier-transformed infrared spectroscopy, visible Raman spectroscopy, synchrotron-radiation single-crystal X-ray diffraction and UV Resonance Raman spectroscopy, allowing key differences to be pinpointed amongst the halogenated analogues.

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.

Coupling enhanced sampling of the apo-receptor with template-based ligand conformers selection: performance in pose prediction in the D3R Grand Challenge 4.

We report the performance of our newly introduced Ensemble Docking with Enhanced sampling of pocket Shape (EDES) protocol coupled to a template-based algorithm to generate near-native ligand conformations in the 2019 iteration of the Grand Challenge (GC4) organized by the D3R consortium. Using either AutoDock4.2 or HADDOCK2.2 docking programs (each software in two variants of the protocol) our method generated native-like poses among the top 5 submitted for evaluation for most of the 20 targets with similar performances. The protein selected for GC4 was the human beta-site amyloid precursor protein cleaving enzyme 1 (BACE-1), a transmembrane aspartic-acid protease. We identified at least one pose whose heavy-atoms RMSD was less than 2.5 Å from the native conformation for 16 (80%) and 17 (85%) of the 20 targets using AutoDock and HADDOCK, respectively. Dissecting the possible sources of errors revealed that: (i) our EDES protocol (with minor modifications) was able to sample sub-ångstrom conformations for all 20 protein targets, reproducing the correct conformation of the binding site within ~ 1 Å RMSD; (ii) as already shown by some of us in GC3, even in the presence of near-native protein structures, a proper selection of ligand conformers is crucial for the success of ensemble-docking calculations. Importantly, our approach performed best among the protocols exploiting only structural information of the apo protein to generate conformations of the receptor for ensemble-docking calculations.

Holo-like and Druggable Protein Conformations from Enhanced Sampling of Binding Pocket Volume and Shape.

Understanding molecular recognition of small molecules by proteins in atomistic detail is key for drug design. Molecular docking is a widely used computational method to mimic ligand-protein association in silico. However, predicting conformational changes occurring in proteins upon ligand binding is still a major challenge. Ensemble docking approaches address this issue by considering a set of different conformations of the protein obtained either experimentally or from computer simulations, e.g., molecular dynamics. However, holo structures prone to host (the correct) ligands are generally poorly sampled by standard molecular dynamics simulations of the apo protein. In order to address this limitation, we introduce a computational approach based on metadynamics simulations called ensemble docking with enhanced sampling of pocket shape (EDES) that allows holo-like conformations of proteins to be generated by exploiting only their apo structures. This is achieved by defining a set of collective variables that effectively sample different shapes of the binding site, ultimately mimicking the steric effect due to the ligand. We assessed the method on three challenging proteins undergoing different extents of conformational changes upon ligand binding. In all cases our protocol generates a significant fraction of structures featuring a low RMSD from the experimental holo geometry. Moreover, ensemble docking calculations using those conformations yielded in all cases native-like poses among the top-ranked ones.

Aminoacyl ß-naphthylamides as substrates and modulators of AcrB multidrug efflux pump.

Efflux pumps of the resistance-nodulation division superfamily, such as AcrB, make a major contribution to multidrug resistance in Gram-negative bacteria. Inhibitors of such pumps would improve the efficacy of antibiotics, and ameliorate the crisis in health care caused by the prevalence of multidrug resistant Gram-negative pathogens. Phenylalanyl-arginine ß-naphthylamide (PAßN), is a well-known inhibitor of AcrB and its homologs. However, its mechanism of inhibition is not clear. Because the hydrolysis of PAßN in Escherichia coli was nearly entirely dependent on an aminopeptidase, PepN, expression of PepN in periplasm allowed us to carry out a quantitative determination of PAßN efflux kinetics through the determination of its periplasmic concentrations by quantitation of the first hydrolysis product, phenylalanine, after a short period of treatment. We found that PAßN is efficiently pumped out by AcrB, with a sigmoidal kinetics. We also examined the behavior of PAßN homologs, Ala ß-naphthylamide, Arg ß-naphthylamide, and Phe ß-naphthylamide, as substrates of AcrB and as modulators of nitrocefin efflux through AcrB. Furthermore, molecular dynamics simulations indicated that the mode of binding of these compounds to AcrB affects the modulatory activity on the efflux of other substrates. These results, and the finding that PAßN changes the nitrocefin kinetics into a sigmoidal one, suggested that PAßN inhibited the efflux of other drugs by binding to the bottom of the distal binding pocket, the so-called hydrophobic trap, and also by interfering with the binding of other drug substrates to the upper part of the binding pocket.

Molecular basis for inhibition of AcrB multidrug efflux pump by novel and powerful pyranopyridine derivatives.

The Escherichia coli AcrAB-TolC efflux pump is the archetype of the resistance nodulation cell division (RND) exporters from Gram-negative bacteria. Overexpression of RND-type efflux pumps is a major factor in multidrug resistance (MDR), which makes these pumps important antibacterial drug discovery targets. We have recently developed novel pyranopyridine-based inhibitors of AcrB, which are orders of magnitude more powerful than the previously known inhibitors. However, further development of such inhibitors has been hindered by the lack of structural information for rational drug design. Although only the soluble, periplasmic part of AcrB binds and exports the ligands, the presence of the membrane-embedded domain in AcrB and its polyspecific binding behavior have made cocrystallization with drugs challenging. To overcome this obstacle, we have engineered and produced a soluble version of AcrB [AcrB periplasmic domain (AcrBper)], which is highly congruent in structure with the periplasmic part of the full-length protein, and is capable of binding substrates and potent inhibitors. Here, we describe the molecular basis for pyranopyridine-based inhibition of AcrB using a combination of cellular, X-ray crystallographic, and molecular dynamics (MD) simulations studies. The pyranopyridines bind within a phenylalanine-rich cage that branches from the deep binding pocket of AcrB, where they form extensive hydrophobic interactions. Moreover, the increasing potency of improved inhibitors correlates with the formation of a delicate protein- and water-mediated hydrogen bond network. These detailed insights provide a molecular platform for the development of novel combinational therapies using efflux pump inhibitors for combating multidrug resistant Gram-negative pathogens.

AcrB drug-binding pocket substitution confers clinically relevant resistance and altered substrate specificity.

The incidence of multidrug-resistant bacterial infections is increasing globally and the need to understand the underlying mechanisms is paramount to discover new therapeutics. The efflux pumps of Gram-negative bacteria have a broad substrate range and transport antibiotics out of the bacterium, conferring intrinsic multidrug resistance (MDR). The genomes of pre- and posttherapy MDR clinical isolates of Salmonella Typhimurium from a patient that failed antibacterial therapy and died were sequenced. In the posttherapy isolate we identified a novel G288D substitution in AcrB, the resistance-nodulation division transporter in the AcrAB-TolC tripartite MDR efflux pump system. Computational structural analysis suggested that G288D in AcrB heavily affects the structure, dynamics, and hydration properties of the distal binding pocket altering specificity for antibacterial drugs. Consistent with this hypothesis, recreation of the mutation in standard Escherichia coli and Salmonella strains showed that G288D AcrB altered substrate specificity, conferring decreased susceptibility to the fluoroquinolone antibiotic ciprofloxacin by increased efflux. At the same time, the substitution increased susceptibility to other drugs by decreased efflux. Information about drug transport is vital for the discovery of new antibacterials; the finding that one amino acid change can cause resistance to some drugs, while conferring increased susceptibility to others, could provide a basis for new drug development and treatment strategies.

Effect of site-directed mutations in multidrug efflux pump AcrB examined by quantitative efflux assays.

BACKGROUND: The Resistance-Nodulation-Division (RND) family transporter AcrB plays a major role in the intrinsic and increased resistance of Escherichia coli to a large number of antibiotics. The distal binding pocket within this multidrug efflux transporter is very large, but the effort to define the roles of various residues facing this pocket through site-directed mutagenesis so far involved only the determination of minimal inhibitory concentrations of drugs in mutants. METHODS: We measured in intact E. coli cells the kinetics of efflux of two substrates, nitrocefin (a cephalosporin) that is predicted mainly to bind to the upper, "groove" domain of the pocket, and L-alanyl-ß-naphthylamide (Ala-Naph) that is likely to bind to the lower, "cave" domain, in a number of site-directed mutants of AcrB, where a hydrophobic or aromatic residue was changed into alanine. RESULTS: The efflux of nitrocefin became attenuated by some mutations in the groove domain, such as I278A and F178A, but in some experiments a mutation in the cave domain, F628A produced a similar result. In some cases an increased value of KM was detected. The efflux of Ala-Naph was increased by mutations in the cave domain, such as F136A and I626A, but also by those in the groove domain (I277A, I278A, F178A). In most cases the increased Vmax values appeared to be responsible. F610A mutation had a profound effect on the efflux of both substrates, as reported earlier. CONCLUSIONS: Our data show for the first time effects of various substrate-binding pocket mutations on the kinetics of efflux of two substrates by the AcrB pump. They also confirm interactions between substrates and drugs predicted by MD simulation studies, and also reveal areas that need future research.

Multidrug binding properties of the AcrB efflux pump characterized by molecular dynamics simulations.

Multidrug resistance in Gram-negative bacteria, to which multidrug efflux pumps such as the AcrB transporter makes a major contribution, is becoming a major public health problem. Unfortunately only a few compounds have been cocrystallized with AcrB, and thus computational approaches are essential in elucidating the interaction between diverse ligands and the pump protein. We used molecular dynamics simulation to examine the binding of nine substrates, two inhibitors, and two nonsubstrates to the distal binding pocket of AcrB, identified earlier by X-ray crystallography. This approach gave us more realistic views of the binding than the previously used docking approach, as the explicit water molecules contributed to the process and the flexible binding site was often seen to undergo large structural changes. We analyzed the interaction in detail in terms of the binding energy, hydrophobic surface-matching, and the residues involved in the process. We found that all substrates tested bound to the pocket, whereas the binding to this site was not preferred for the nonsubstrates. Interestingly, both inhibitors [Phe-Arg-ß-naphthylamide and 1-(1-naphtylmethyl)-piperazine] tended to move out of the pocket at least partially, getting into contact with a glycine-rich loop that separates the distal pocket from the more proximal region of the protein and is thought to control the access of substrates to the distal pocket.

The Phe-Phe Motif for Peptide Self-Assembly in Nanomedicine.

Since its discovery, the Phe-Phe motif has gained in popularity as a minimalist building block to drive the self-assembly of short peptides and their analogues into nanostructures and hydrogels. Molecules based on the Phe-Phe motif have found a range of applications in nanomedicine, from drug delivery and biomaterials to new therapeutic paradigms. Here we discuss the various production methods for this class of compounds, and the characterization, nanomorphologies, and application of their self-assembled nanostructures. We include the most recent findings on their remarkable properties, which hold substantial promise for the creation of the next generation nanomedicines.

Molecular mechanism of viral resistance to a potent non-nucleoside inhibitor unveiled by molecular simulations.

Recently, we reported on a potent benzimidazole derivative (227G) that inhibits the growth of the bovine viral diarrhea virus (BVDV) in cell-based and enzyme assays at nanomolar concentrations. The target of 227G is the viral RNA-dependent RNA polymerase (RdRp), and the I261M mutation located in motif I of the RdRp finger domain was found to induce drug resistance. Here we propose a molecular mechanism for the retained functionality of the enzyme in the presence of the inhibitor, on the basis of a thorough computational study of the apo and holo forms of the BVDV RdRp either in the wild type (wt) or in the form carrying the I261M mutation. Our study shows that although the mutation affects to some extent the structure of the apoenzyme, the functional dynamics of the protein appear to be largely maintained, which is consistent with the retained functionality of this natural mutant. Despite the binding site of 227G not collapsing or undergoing drastic structural changes upon introduction of the I261M substitution, these alterations reflect crucially on the binding mode of 227G, which is significantly different from that found in wt RdRp. In particular, while in the wt system the four loops lining the template entrance site embrace 227G and close the template passageway, in the I261M variant the template entrance is only marginally occluded, allowing in principle the translocation of the template to the interior of the enzyme. In addition, the mutated enzyme in the presence of 227G retains several characteristics of the wt apoprotein. Our work provides an original molecular picture of a resistance mechanism that is consistent with published experimental data.

Molecular mechanism of MBX2319 inhibition of Escherichia coli AcrB multidrug efflux pump and comparison with other inhibitors.

Efflux pumps of the resistance nodulation division (RND) superfamily, such as AcrB, make a major contribution to multidrug resistance in Gram-negative bacteria. The development of inhibitors of the RND pumps would improve the efficacy of current and next-generation antibiotics. To date, however, only one inhibitor has been cocrystallized with AcrB. Thus, in silico structure-based analysis is essential for elucidating the interaction between other inhibitors and the efflux pumps. In this work, we used computer docking and molecular dynamics simulations to study the interaction between AcrB and the compound MBX2319, a novel pyranopyridine efflux pump inhibitor with potent activity against RND efflux pumps of Enterobacteriaceae species, as well as other known inhibitors (D13-9001, 1-[1-naphthylmethyl]-piperazine, and phenylalanylarginine-ß-naphthylamide) and the binding of doxorubicin to the efflux-defective F610A variant of AcrB. We also analyzed the binding of a substrate, minocycline, for comparison. Our results show that MBX2319 binds very tightly to the lower part of the distal pocket in the B protomer of AcrB, strongly interacting with the phenylalanines lining the hydrophobic trap, where the hydrophobic portion of D13-9001 was found to bind by X-ray crystallography. Additionally, MBX2319 binds to AcrB in a manner that is similar to the way in which doxorubicin binds to the F610A variant of AcrB. In contrast, 1-(1-naphthylmethyl)-piperazine and phenylalanylarginine-ß-naphthylamide appear to bind to somewhat different areas of the distal pocket in the B protomer of AcrB than does MBX2319. However, all inhibitors (except D13-9001) appear to distort the structure of the distal pocket, impairing the proper binding of substrates.

Ru[(bpy)2(dppz)]²âº and Rh[(bpy)2(chrysi)]³âº targeting double strand DNA: the shape of the intercalating ligand tunes the free energy landscape of deintercalation.

Octahedral metal complexes can bind to double strand (ds) DNA either by intercalation or by insertion, this latter mechanism being observed in the case of mismatched base pairs (bps). In this work we modeled the process of deintercalation from the major groove for Δ-Ru[(bpy)2(dppz)](2+) (1) and Δ-Rh[(bpy)2(chrysi)](3+) (2), prototypical examples of metallo-intercalators and metallo-insertors, respectively. By using advanced sampling techniques, we show that the two complexes have comparable deintercalation barriers and that in both systems the main cost of deintercalation is due to disruption of π-π stacking interactions between the intercalating moiety and the bps flanking the binding site. A striking difference between dppz and chrysi is found in their intercalation modes, being their longest axes, respectively, perpendicular and parallel to the P-P direction between opposite DNA strands. This leads the two ligands to deintercalate from the DNA through different mechanisms. Compound 1 goes through the formation of a metastable short-lived intermediate, with an overall free energy barrier of ~14.5 kcal/mol, in line with experimental findings. Due to the length of the dppz intercalating moiety, an extended plateau appears in the free energy landscape at ~3 kcal/mol above the most stable minimum. Compound 2 must cross a similar barrier (~15.5 kcal/mol), but does not form intermediates along the deintercalation path, and the deintercalation profile is steeper than that found for 1. Thus, the shape of the intercalating moiety affects the deintercalation mechanism of these inorganic molecules. This work is a first step to rationalize from a computational perspective the factors tuning the preferential binding mode of inorganic molecules (such as diagnostic probes, therapeutic agents, or regulators of DNA expression) to ds DNA.

Predicting binding events in very flexible, allosteric, multi-domain proteins.

Knowledge of the structures formed by proteins and small ligands is of fundamental importance for understanding molecular principles of chemotherapy and for designing new and more effective drugs. Due to the still high costs and to the several limitations of experimental techniques, it is most often desirable to predict these ligand-protein complexes in silico, particularly when screening for new putative drugs from databases of millions of compounds. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology aimed at generating bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites. Validation was performed on the enzyme adenylate kinase (ADK), a paradigmatic example of proteins that undergo very large conformational changes upon ligand binding. By only exploiting the unbound structure and the putative binding sites of the protein, we generated a significant fraction of bound-like structures, which employed in ensemble-docking calculations allowed to find native-like poses of substrates, inhibitors, and catalytically incompetent binders. Our protocol provides a general framework for the generation of bound-like conformations of flexible proteins that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening for difficult targets, prediction of the impact of amino acid mutations on structure and dynamics, and protein engineering.

Insulin amyloid fibril formation reduction by tripeptide stereoisomers.

Insulin fibrillation is a problem for diabetic patients that can occur during storage and transport, as well as at the subcutaneous injection site, with loss of bioactivity, inflammation, and various adverse effects. Tripeptides are ideal additives to stabilise insulin formulations, thanks to their low cost of production and inherent cytocompatibility. In this work, we analysed the ability of eight tripeptide stereoisomers to inhibit the fibrillation of human insulin in vitro. The sequences contain proline as ß-breaker and Phe-Phe as binding motif for the amyloid-prone aromatic triplet found in insulin. Experimental data based on spectroscopy, fluorescence, microscopy, and calorimetric techniques reveal that one stereoisomer is a more effective inhibitor than the others, and cell live/dead assays confirmed its high cytocompatibility. Importantly, in silico data revealed the key regions of insulin engaged in the interaction with this tripeptide, rationalising the molecular mechanism behind insulin fibril formation reduction.

Peptide Stereochemistry Effects from pKa-Shift to Gold Nanoparticle Templating in a Supramolecular Hydrogel.

The divergent supramolecular behavior of a series of tripeptide stereoisomers was elucidated through spectroscopic, microscopic, crystallographic, and computational techniques. Only two epimers were able to effectively self-organize into amphipathic structures, leading to supramolecular hydrogels or crystals, respectively. Despite the similarity between the two peptides' turn conformations, stereoconfiguration led to different abilities to engage in intramolecular hydrogen bonding. Self-assembly further shifted the pKa value of the C-terminal side chain. As a result, across the pH range 4-6, only one epimer predominated sufficiently as a zwitterion to reach the critical molar fraction, allowing gelation. By contrast, the differing pKa values and higher dipole moment of the other epimer favored crystallization. The four stereoisomers were further tested for gold nanoparticle (AuNP) formation, with the supramolecular hydrogel being the key to control and stabilize AuNPs, yielding a nanocomposite that catalyzed the photodegradation of a dye. Importantly, the AuNP formation occurred without the use of reductants other than the peptide, and the redox chemistry was investigated by LC-MS, NMR, and infrared scattering-type near field optical microscopy (IR s-SNOM). This study provides important insights for the rational design of simple peptides as minimalistic and green building blocks for functional nanocomposites.

Predicting permeation of compounds across the outer membrane of P. aeruginosa using molecular descriptors.

The ability Gram-negative pathogens have at adapting and protecting themselves against antibiotics has increasingly become a public health threat. Data-driven models identifying molecular properties that correlate with outer membrane (OM) permeation and growth inhibition while avoiding efflux could guide the discovery of novel classes of antibiotics. Here we evaluate 174 molecular descriptors in 1260 antimicrobial compounds and study their correlations with antibacterial activity in Gram-negative Pseudomonas aeruginosa. The descriptors are derived from traditional approaches quantifying the compounds' intrinsic physicochemical properties, together with, bacterium-specific from ensemble docking of compounds targeting specific MexB binding pockets, and all-atom molecular dynamics simulations in different subregions of the OM model. Using these descriptors and the measured inhibitory concentrations, we design a statistical protocol to identify predictors of OM permeation/inhibition. We find consistent rules across most of our data highlighting the role of the interaction between the compounds and the OM. An implementation of the rules uncovered in our study is shown, and it demonstrates the accuracy of our approach in a set of previously unseen compounds. Our analysis sheds new light on the key properties drug candidates need to effectively permeate/inhibit P. aeruginosa, and opens the gate to similar data-driven studies in other Gram-negative pathogens.

Some ligands enhance the efflux of other ligands by the Escherichia coli multidrug pump AcrB.

By measuring quantitatively the active efflux of cephalosporins by the RND (resistance-nodulation-division) family efflux pump AcrB in intact cells of Escherichia coli, we found that the simultaneous presence of another substrate, such as chloramphenicol, benzene, cyclohexane, or Arg ß-naphthilamide, significantly enhanced the extrusion of cephalosporins. The stimulation occurred also in a strain expressing the covalently linked trimer of AcrB, and thus cannot be ascribed to the enhanced assembly of the trimer from AcrB monomers. When Val139 of AcrB was changed into Phe, the stimulation by benzene was found to occur at much lower concentration of the solvent. A plausible explanation of these observations is that the AcrB pump is constructed to pump out very rapidly the solvent or chloramphenicol molecules, and thus the efflux of cephalosporins, which presumably bind to a different subsite within the large binding pocket of AcrB, can become facilitated. Computer simulations of ligand binding to AcrB, both by docking and by molecular dynamics simulations, produced results supporting and extending this hypothesis. Benzene and the cephalosporin nitrocefin can bind simultaneously to the distal binding pocket of AcrB, both in the wild type and in the V139F variant. Interestingly, while the binding position and strength of benzene are almost unaffected by the presence of nitrocefin, this latter substrate is significantly displaced toward the exit gate in both wild type and mutant transporter in the presence of benzene. Additionally, the cephalosporin efflux may be enhanced by the binding of solvents (sometimes to the cephalosporin-free protomer), which could accelerate AcrB conformational changes necessary for substrate extrusion.