Effect of the F610A mutation on substrate extrusion in the AcrB transporter: explanation and rationale by molecular dynamics simulations

Statistical thermodynamics for functionally rotating mechanism of the multidrug efflux transporter AcrB

AcrB, a homotrimer, is the pivotal part of a multidrug efflux pump. A "functionally rotating" picture has been proposed for the drug transport by AcrB, but its mechanism remains unresolved. Here, we investigate the energetics of the whole functional rotation cycle using our theoretical methods. We find that the packing efficiency of AcrB is ununiform, and this ununiformity plays imperative roles primarily through the solvent-entropy effect. When a proton binds to or dissociates from a protomer, the packing properties of this protomer and its two interfaces are perturbed overall in the direction that the solvent translational entropy is lowered. The packing properties of the other two protomers are then reorganized with the recovery or maintenance of closely packed interfaces, so that the solvent-entropy loss can be compensated. The functional structural change by an isolated protomer would cause a seriously large free-energy increase. By forming a trimer, any free-energy increase caused by a protomer is always canceled out by the free-energy decrease brought by the other two protomers via the mechanism mentioned above. The functional structural rotation is thus accomplished using the free-energy decrease arising from the transfer of only a single proton per cycle. The similarities to F1-ATPase are also discussed.

Improving olefin tolerance and production in E. coli using native and evolved AcrB

Microorganisms can be engineered for the production of chemicals utilized in the polymer industry. However many such target compounds inhibit microbial growth and might correspondingly limit production levels. Here, we focus on compounds that are precursors to bioplastics, specifically styrene and representative alpha-olefins; 1-hexene, 1-octene, and 1-nonene. We evaluated the role of the Escherichia coli efflux pump, AcrAB-TolC, in enhancing tolerance towards these olefin compounds. AcrAB-TolC is involved in the tolerance towards all four compounds in E. coli. Both styrene and 1-hexene are highly toxic to E. coli. Styrene is a model plastics precursor with an established route for production in E. coli (McKenna and Nielsen, 2011). Though our data indicates that AcrAB-TolC is important for its optimal production, we observed a strong negative selection against the production of styrene in E. coli. Thus we used 1-hexene as a model compound to implement a directed evolution strategy to further improve the tolerance phenotype towards this alpha-olefin. We focused on optimization of AcrB, the inner membrane domain known to be responsible for substrate binding, and found several mutations (A279T, Q584R, F617L, L822P, F927S, and F1033Y) that resulted in improved tolerance. Several of these mutations could also be combined in a synergistic manner. Our study shows efflux pumps to be an important mechanism in host engineering for olefins, and one that can be further improved using strategies such as directed evolution, to increase tolerance and potentially production. Biotechnol. Bioeng. 2015;112: 879–888. © 2015 The Authors. Biotechnology and Bioengineering Published by John Wiley & Periodicals, Inc.

Computational study of correlated domain motions in the AcrB efflux transporter

As active part of the major efflux system in E. coli bacteria, AcrB is responsible for the uptake and pumping of toxic substrates from the periplasm toward the extracellular space. In combination with the channel protein TolC and membrane fusion protein AcrA, this efflux pump is able to help the bacterium to survive different kinds of noxious compounds. With the present study we intend to enhance the understanding of the interactions between the domains and monomers, for example, the transduction of mechanical energy from the transmembrane domain into the porter domain, correlated motions of different subdomains within monomers, and cooperative effects between monomers. To this end, targeted molecular dynamics simulations have been employed either steering the whole protein complex or specific parts thereof. By forcing only parts of the complex towards specific conformational states, the risk for transient artificial conformations during the simulations is reduced. Distinct cooperative effects between the monomers in AcrB have been observed. Possible allosteric couplings have been identified providing microscopic insights that might be exploited to design more efficient inhibitors of efflux systems.

Fisetin inhibits Listeria monocytogenes virulence by interfering with the oligomerization of listeriolysin O

Listeriolysin O (LLO), an essential virulence determinant of Listeria monocytogenes, is a pore-forming toxin whose primary function is to facilitate cytosolic bacterial replication by breaching the phagosomal membranes, which is critical for the pathogen to evade host immune recognition. The critical role of LLO in the virulence of L. monocytogenes renders it an ideal target for designing novel antivirulence therapeutics. We found that fisetin, a natural flavonoid without antimicrobial activity, is a potent antagonist of LLO-mediated hemolysis. Fisetin effectively inhibits L. monocytogenes infection in both tissue culture and animal infection models. Molecular modeling and mutational analysis revealed that fisetin directly engages loop 2 and loop 3 of LLO, leading to the blockage of cholesterol binding and the reduction of its oligomerization, thus inhibiting its hemolytic activity. Our results establish fisetin as an effective antitoxin agent for LLO, which can be further developed into novel therapeutics against infections caused by L. monocytogenes.

Random mutagenesis of the multidrug transporter AcrB from Escherichia coli for identification of putative target residues of efflux pump inhibitors

Efflux is an important mechanism of bacterial multidrug resistance (MDR), and the inhibition of MDR pumps by efflux pump inhibitors (EPIs) could be a promising strategy to overcome MDR. 1-(1-Naphthylmethyl)-piperazine (NMP) and phenylalanine-arginine-β-naphthylamide (PAβN) are model EPIs with activity in various Gram-negative bacteria expressing AcrB, the major efflux pump of Escherichia coli, or similar homologous pumps of the resistance-nodulation-cell division class. The aim of the present study was to generate E. coli AcrB mutants resistant to the inhibitory action of the two model EPIs and to identify putative EPI target residues in order to better understand mechanisms of pump inhibition. Using an in vitro random mutagenesis approach focusing on the periplasmic domain of AcrB, we identified the double mutation G141D N282Y, which substantially compromised the synergistic activity of NMP with linezolid, was associated with similar intracellular linezolid concentrations in the presence and absence of NMP, and did not impair the intrinsic MICs of various pump substrates and dye accumulation. We propose that these mutations near the outer face of the distal substrate binding pocket reduce NMP trapping. Other residues found to be relevant for efflux inhibition by NMP were G288 and A279, but mutations at these sites also changed the susceptibility to several pump substrates. Unlike with NMP, we were unable to generate AcrB periplasmic domain mutants with resistance or partial resistance to the EPI activity of PAβN, which is consistent with the modes of action of PAβN differing from those of NMP.

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.

Electrostatics and flexibility drive membrane recognition and early penetration by the antimicrobial peptide dendrimer bH1

Molecular dynamics simulations of the polycationic antimicrobial peptide dendrimer (Leu)8(DapLeu)4(DapPhe)2DapLys-NH2 binding to membranes suggest that electrostatic interactions with the polyanionic lipopolysaccharide (LPS) and conformational flexibility of the 2,3-diaminopropanoic acid (Dap) branching units drive its selective insertion into microbial membranes.

Molecular modeling reveals the novel inhibition mechanism and binding mode of three natural compounds to staphylococcal α-hemolysin

α-Hemolysin (α-HL) is a self-assembling, channel-forming toxin that is produced as a soluble monomer by Staphylococcus aureus strains. Until now, α-HL has been a significant virulence target for the treatment of S. aureus infection. In our previous report, we demonstrated that some natural compounds could bind to α-HL. Due to the binding of those compounds, the conformational transition of α-HL from the monomer to the oligomer was blocked, which resulted in inhibition of the hemolytic activity of α-HL. However, these results have not indicated how the binding of the α-HL inhibitors influence the conformational transition of the whole protein during the oligomerization process. In this study, we found that three natural compounds, Oroxylin A 7-O-glucuronide (OLG), Oroxin A (ORA), and Oroxin B (ORB), when inhibiting the hemolytic activity of α-HL, could bind to the “stem” region of α-HL. This was completed using conventional Molecular Dynamics (MD) simulations. By interacting with the novel binding sites of α-HL, the ligands could form strong interactions with both sides of the binding cavity. The results of the principal component analysis (PCA) indicated that because of the inhibitors that bind to the “stem” region of α-HL, the conformational transition of α-HL from the monomer to the oligomer was restricted. This caused the inhibition of the hemolytic activity of α-HL. This novel inhibition mechanism has been confirmed by both the steered MD simulations and the experimental data obtained from a deoxycholate-induced oligomerization assay. This study can facilitate the design of new antibacterial drugs against S. aureus.

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.

Functional rotation induced by alternating protonation states in the multidrug transporter AcrB: all-atom molecular dynamics simulations

The multidrug transporter AcrB actively exports a wide variety of noxious compounds using proton-motive force as an energy source in Gram-negative bacteria. AcrB adopts an asymmetric structure comprising three protomers with different conformations that are sequentially converted during drug export; these cyclic conformational changes during drug export are referred to as functional rotation. To investigate functional rotation driven by proton-motive force, all-atom molecular dynamics simulations were performed. Using different protonation states for the titratable residues in the middle of the transmembrane domain, our simulations revealed the correlation between the specific protonation states and the side-chain configurations. Changing the protonation state for Asp408 induced a spontaneous structural transition, which suggests that the proton translocation stoichiometry may be one proton per functional rotation cycle. Furthermore, our simulations demonstrate that alternating the protonation states in the transmembrane domain induces functional rotation in the porter domain, which is primarily responsible for drug transport.

Drug uptake pathways of multidrug transporter AcrB studied by molecular simulations and site-directed mutagenesis experiments

Multidrug resistance has been a critical issue in current chemotherapy. In Escherichia coli , a major efflux pump responsible for the multidrug resistance contains a transporter AcrB. Crystallographic studies and mutational assays of AcrB provided much of structural and overall functional insights, which led to the functionally rotating mechanism. However, the drug uptake pathways are somewhat controversial because at least two possible pathways, the vestibule and the cleft paths, were suggested. Here, combining molecular simulations and site-directed mutagenesis experiments, we addressed the uptake mechanism finding that the drug uptake pathways can be significantly different depending on the properties of drugs. First, in the computational free energy analysis of drug movements along AcrB tunnels, we found a ligand-dependent drug uptake mechanism. With the same molecular sizes, drugs that are both strongly hydrophobic and lipophilic were preferentially taken in via the vestibule path, while other drugs favored the cleft path. Second, direct simulations realized totally about 3500 events of drug uptake by AcrB for a broad range of drug property. These simulations confirmed the ligand-dependent drug uptake and further suggested that a smaller drug favors the vestibule path, while a larger one is taken in via the cleft path. Moreover, the direct simulations identified an alternative uptake path which is not visible in the crystal structure. Third, site-directed mutagenesis of AcrB in E. coli verified that mutations of residues located along the newly identified path significantly reduced the efflux efficiency, supporting its relevance in in vivo function.

Molecular Dynamics Computer Simulations of Multidrug RND Efflux Pumps

Over-expression of multidrug efflux pumps of the Resistance Nodulation Division (RND) protein super family counts among the main causes for microbial resistance against pharmaceuticals. Understanding the molecular basis of this process is one of the major challenges of modern biomedical research, involving a broad range of experimental and computational techniques. Here we review the current state of RND transporter investigation employing molecular dynamics simulations providing conformational samples of transporter components to obtain insights into the functional mechanism underlying efflux pump-mediated antibiotics resistance in Escherichia coli and Pseudomonas aeruginosa.

Self-association and complexation of the anti-cancer drug doxorubicin with PEGylated hyperbranched polyesters in an aqueous environment

Fully atomistic molecular dynamics simulations were employed in order to examine in detail the self-assembly characteristics and the complexation behavior of the anticancer drug doxorubicin with PEGylated hyperbranched polyesters in an aqueous environment. We have examined two variants of the polymeric compound by altering the length of the hydrophilic poly(ethylene glycol) arms attached to the hydrophobic hyperbranched core. By comparing the clustering properties of the drug molecules in a polymer-free system to those in the polymer-containing models, we were able to assess the effects related to the presence and to the structural features of the polymer moiety. In addition, we have distinguished the effects associated with the neutral and protonated drug molecules separately. It was found that, in the presence of the polymeric material, the drug molecules formed clusters preferentially close to the polymer's periphery, the characteristics of which depended on the structural details of the polymeric host and on the charge of the drug molecules. Hydrogen bonding was found to contribute to the polymer/drug complexation, with the nature of the prevailing donor/acceptor pairs depending on the charge of the drug. Dynamic analysis of the drugs' motion revealed that in the polymer-containing systems the drug molecules experienced a larger degree of confinement within the formed clusters compared to that describing their polymer-free analogues, while the structural coherence of the clusters was found to be more persistent in the system with the larger poly(ethylene glycol) arms. The results described in this work, through the monitoring of both static and dynamic aspects of the self-association and the complexation behavior of the neutral and charged molecules of doxorubicin with the polymeric host, may help toward the elucidation of the key parameters that are involved in the formation of effective polymer-based carriers for drug molecules of the anthracycline family used in cancer chemotherapy.

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.

Recognition of imipenem and meropenem by the RND-transporter MexB studied by computer simulations

Basic understanding of the means by which multidrug efflux systems can efficiently recognize and transport drugs constitutes a fundamental step toward development of compounds able to tackle the continuous outbreak of new bacterial strains resistant to traditional antibiotics. We applied a series of computational techniques, from molecular docking to molecular dynamics simulations and free energy estimate methods, to determine the differences in the binding properties of imipenem and meropenem, two potent antibiotics of the carbapenem family, to MexB, the RND transporter of the major efflux system of Pseudomonas aeruginosa. We identified and characterized two affinity sites in the periplasmic domain of the transporter, sharing strong similarities with the distal and proximal binding pockets identified in AcrB, the homologue of MexB in Escherichia coli. According to our results, meropenem has a higher affinity to the distal binding pocket than imipenem while both compounds are weakly bound to the proximal pocket. This different behavior is mainly due to the hydration properties of the nonpharmacophore part of the two compounds, being that of imipenem less bulky and hydrophobic. Our data provide for the first time a rationale at molecular level for the experimental evidence indicating meropenem as a compound strongly affected by MexB contrary to imipenem, which is apparently poorly transported by the same pump.

Porter domain opening and closing motions in the multi-drug efflux transporter AcrB

Acriflavine resistance protein B acts as the active transporter in the multi-drug efflux pump Acriflavine resistance proteins A / B - Tolerance to colicins protein in Escherichia coli. Within the same reaction cycle intermediate all Acriflavine resistance protein B X-ray structures display highly similar conformations of the substrate-recruiting and transporting porter domain. To assess if this structural homogeneity is an intrinsic feature of Acriflavine resistance protein B or stems from other causes we performed a series of six independent, unbiased 100 ns molecular dynamics simulations of membrane-embedded, asymmetric, substrate-free wild type Acriflavine resistance protein B in a 150 mM NaCl solution. We find the porter domain more flexible than previously assumed displaying clear opening and closing motions of the proximal binding pocket (L and T-state) and the exit of the drug transport channels (O-intermediate). Concurrently the hydrophobic binding pocket favors a closed conformation in all three protomers. Our findings suggest that the conformational homogeneity seen in the crystal structures is likely an effect of bound but structurally unresolved substrate. Our simulations further imply that each of the known three reaction cycle intermediates occurs in at least two variants, the Thr676 loop independently regulates porter domain access and likely plays a key role in substrate transport. On a 100 ns time scale we find no evidence supporting the proposed LLL resting state in the absence of substrate. If the proximal binding pocket dynamics have an inhibiting effect on Acriflavine resistance protein B pump activity lowering the life time of substrate-accessible conformations, the observed dynamics could provide a structural explanation for the Acriflavine resistance protein B activity-enhancing effect of the adaptor protein Acriflavine resistance protein A stabilizing PC1 and PC2 subdomain orientations.

Mutations in MexB that affect the efflux of antibiotics with cytoplasmic targets

Drug efflux pumps such as MexAB-OprM from Pseudomonas aeruginosa confer resistance to a wide range of chemically different compounds. Within the tripartite assembly, the inner membrane protein MexB is mainly responsible for substrate recognition. Recently, considerable advances have been made in elucidating the drug efflux pathway through the large periplasmic domains of resistance-nodulation-division (RND) transporters. However, little is known about the role of amino acids in other parts of the protein. We have investigated the role of two conserved phenylalanine residues that are aligned around the cytoplasmic side of the central cavity of MexB. The two conserved phenylalanine residues have been mutated to alanine residues (FAFA MexB). The interaction of the wild-type and mutant proteins with a variety of drugs from different classes was investigated by assays of cytotoxicity and drug transport. The FAFA mutation affected the efflux of compounds that have targets inside the cell, but antibiotics that act on cell wall synthesis and membrane probes were unaffected. Combined, our results indicate the presence of a hitherto unidentified cytoplasmic-binding site in RND drug transporters and enhance our understanding of the molecular mechanisms that govern drug resistance in Gram-negative pathogens.

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.