Energetics and conformational pathways of functional rotation in the multidrug transporter AcrB

Impact of exporter proteins and their engineering on the productivity of Corynebacterium

Enhancing product efflux is crucial in improving fermentative bioproduction. Despite advancements in metabolic engineering guided by the design-build-test-learn cycle, membrane transport engineering of product efflux remains underdeveloped, limiting the efficient production of target chemicals. This review explores the historical findings on product efflux, regardless of passive or active transport, in fermentation engineering, focusing on Corynebacterium species, and highlights the potential of multidrug transporters as valuable screening sources for efflux improvement. Furthermore, the review emphasizes the importance of understanding the machinery of efflux transporters to optimize their functionality. Molecular dynamics simulations are a promising tool for exploring novel strategies to advance fermentation-related processes. These insights provide a framework for overcoming current challenges in membrane transport engineering of product efflux and improving industrial-scale bioproduction. KEY POINTS: • Review of strategies to enhance product efflux in Corynebacterium species. • Multidrug transporters are key tools for optimizing metabolite efflux. • Efflux transporter mechanisms analyzed to improve microbial productivity. • Molecular dynamics simulations employed for understanding transporter mechanisms.

Dissecting Large-Scale Structural Transitions in Membrane Transporters Using Advanced Simulation Technologies

Membrane transporters are integral membrane proteins that act as gatekeepers of the cell, controlling fundamental processes such as recruitment of nutrients and expulsion of waste material. At a basic level, transporters operate using the "alternating access model," in which transported substances are accessible from only one side of the membrane at a time. This model usually involves large-scale structural changes in the transporter, which often cannot be captured using unbiased, conventional molecular simulation techniques. In this article, we provide an overview of some of the major simulation techniques that have been applied to characterize the structural dynamics and energetics involved in the transition of membrane transporters between their functional states. After briefly introducing each technique, we discuss some of their advantages and limitations and provide some recent examples of their application to membrane transporters.

Specificity and tunability of efflux pumps: A new role for the proton gradient?

Efflux pumps that transport antibacterial drugs out of bacterial cells have broad specificity, commonly leading to broad spectrum resistance and limiting treatment strategies for infections. It remains unclear how efflux pumps can maintain this broad spectrum specificity to diverse drug molecules while limiting the efflux of other cytoplasmic content. We have investigated the origins of this broad specificity using theoretical models informed by the experimentally determined structural and kinetic properties of efflux pumps. We developed a set of mathematical models describing operation of efflux pumps as a discrete cyclic stochastic process across a network of states characterizing pump conformations and the presence/absence of bound ligands and protons. These include a minimal three-state model that lends itself to clear analytic calculations as well as a five-state model that relaxes some of the simpler model's most strict assumptions. We found that the pump specificity is determined not solely by the drug affinity to the pump-as is commonly assumed-but it is also directly affected by the periplasmic pH and the transmembrane potential. Therefore, changes to the proton concentration gradient and voltage drop across the membrane can influence how effective the pump is at extruding a particular drug molecule. Furthermore, we found that while both the proton concentration gradient across the membrane and the transmembrane potential contribute to the thermodynamic force driving the pump, their effects on the efflux enter not strictly in a combined proton motive force. Rather, they have two distinguishable effects on the overall throughput. These results highlight the unexpected effects of thermodynamic driving forces out of equilibrium and illustrate how efflux pump structure and function are conducive to the emergence of multidrug resistance.

Structural and functional diversity of Resistance-Nodulation-Division (RND) efflux pump transporters with implications for antimicrobial resistance

SUMMARYThe discovery of bacterial efflux pumps significantly advanced our understanding of how bacteria can resist cytotoxic compounds that they encounter. Within the structurally and functionally distinct families of efflux pumps, those of the Resistance-Nodulation-Division (RND) superfamily are noteworthy for their ability to reduce the intracellular concentration of structurally diverse antimicrobials. RND systems are possessed by many Gram-negative bacteria, including those causing serious human disease, and frequently contribute to resistance to multiple antibiotics. Herein, we review the current literature on the structure-function relationships of representative transporter proteins of tripartite RND efflux pumps of clinically important pathogens. We emphasize their contribution to bacterial resistance to clinically used antibiotics, host defense antimicrobials and other biocides, as well as highlighting structural similarities and differences among efflux transporters that help bacteria survive in the face of antimicrobials. Furthermore, we discuss technical advances that have facilitated and advanced efflux pump research and suggest future areas of investigation that will advance antimicrobial development efforts.

Insights into Allosteric Inhibition of the AcrB Efflux Pump: Role of Distinct Binding Pockets, Protomer Preferences, and Crosstalk Disruption

AcrB, a key component in bacterial efflux processes, exhibits distinct binding pockets that influence inhibitor interactions. In addition to the well-known distal binding pocket within the periplasmic domain, a noteworthy pocket amidst the transmembrane (TM) helices serves as an alternate binding site for inhibitors. The bacterial efflux mechanism involves a pivotal functional rotation of the TM protein, inducing conformational changes in each protomer and propelling drugs toward the outer membrane domain. Surprisingly, inhibitors binding to the TM domain display a preference for L protomers over T protomers. Metadynamics simulations elucidate that Lys940 in the TM domain of AcrB can adopt two conformations in L protomers, whereas the energy barrier for such transitions is higher in T protomers. This phenomenon results in stable inhibitor binding in l protomers. Upon a detailed analysis of unbinding pathways using random accelerated molecular dynamics and umbrella sampling, we have identified three distinct routes for ligand exit from the allosteric site, specifically involving regions within the TM domains─TM4, TM5, and TM10. To explore allosteric crosstalk, we focused on the following key residues: Val452 from the TM domain and Ala831 from the porter domain. Surprisingly, our findings reveal that inhibitor binding disrupts this communication. The shortest path connecting Val452 and Ala831 increases upon inhibitor binding, suggesting sabotage of the natural interdomain communication dynamics. This result highlights the intricate interplay between inhibitor binding and allosteric signaling within our studied system.

Titratable residues that drive RND efflux: Insights from molecular simulations

The resistance-nodulation-division efflux machinery confers antimicrobial resistance to Gram-negative bacteria by actively pumping antibiotics out of the cell. The protein complex is powered by proton motive force; however, the proton transfer mechanism itself and indeed even its stoichiometry is still unclear. Here we review computational studies from the last decade that focus on elucidating the number of protons transferred per conformational cycle of the pump. Given the difficulties in studying proton movement using even state-of-the-art structural biology methods, the contributions from computational studies have been invaluable from a mechanistic perspective.

Crystal structures of multidrug efflux transporters from Burkholderia pseudomallei suggest details of transport mechanism

BpeB and BpeF are multidrug efflux transporters from Burkholderia pseudomallei that enable multidrug resistance. Here, we report the crystal structures of BpeB and BpeF at 2.94 Å and 3.0 Å resolution, respectively. BpeB was found as an asymmetric trimer, consistent with the widely-accepted functional rotation mechanism for this type of transporter. One of the monomers has a distinct structure that we interpret as an intermediate along this functional cycle. Additionally, a detergent molecule bound in a previously undescribed binding site provides insights into substrate translocation through the pathway. BpeF shares structural similarities with the crystal structure of OqxB from Klebsiella pneumoniae, where both are symmetric trimers composed of three "binding"-state monomers. The structures of BpeB and BpeF further our understanding of the functional mechanisms of transporters belonging to the HAE1-RND superfamily.

Tripartite efflux pumps of the RND superfamily: what did we learn from computational studies?

Bacterial resistance to antibiotics has been long recognized as a priority to address for human health. Among all micro-organisms, the so-called multi-drug resistant (MDR) bacteria, which are resistant to most, if not all drugs in our current arsenal, are particularly worrisome. The World Health Organization has prioritized the ESKAPE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa and Enterobacter species) pathogens, which include four Gram-negative bacterial species. In these bacteria, active extrusion of antimicrobial compounds out of the cell by means of 'molecular guns' known as efflux pumps is a main determinant of MDR phenotypes. The resistance-nodulation-cell division (RND) superfamily of efflux pumps connecting the inner and outer membrane in Gram-negative bacteria is crucial to the onset of MDR and virulence, as well as biofilm formation. Thus, understanding the molecular basis of the interaction of antibiotics and inhibitors with these pumps is key to the design of more effective therapeutics. With the aim to contribute to this challenge, and complement and inspire experimental research, in silico studies on RND efflux pumps have flourished in recent decades. Here, we review a selection of such investigations addressing the main determinants behind the polyspecificity of these pumps, the mechanisms of substrate recognition, transport and inhibition, as well as the relevance of their assembly for proper functioning, and the role of protein-lipid interactions. The journey will end with a perspective on the role of computer simulations in addressing the challenges posed by these beautifully complex machineries and in supporting the fight against the spread of MDR bacteria.

Analysis of AcrB in Klebsiella pneumoniae reveals natural variants promoting enhanced multidrug resistance

Infections caused by Klebsiella pneumoniae are often difficult to manage due to the high frequency of multidrug resistance, often conferred by efflux pumps. In this study, we analyzed sequence variations of the major RND family multidrug efflux pump AcrB from 387 assembled K. pneumoniae genomes. We confirm that AcrB is a highly-conserved efflux pump in K. pneumoniae, and identified several variants that were prevalent in clinical isolates. Molecular dynamics analyses on two of these variants (L118M and S966A) suggested conformational changes that may correlate with increased drug efflux capabilities. The L118M change resulted in enhanced protein rigidity while the flexibility of drug binding pockets was stable or increased, and the interactions between the proximal pockets and water molecules were stronger. For S966A, the significantly enlarged proximal pocket suggested higher drug accommodation ability. Consistent with these predictions, the L118M and S966A variants conferred a slightly increased ability to grow in the presence of tetracycline and to survive cefoxitin exposure when overexpressed. In summary, our results suggest that the emergence of enhanced-function AcrB variants may be a potential risk for increased antibiotic resistance in clinical K. pneumoniae isolates.

Transport Dynamics of MtrD: An RND Multidrug Efflux Pump from Neisseria gonorrhoeae

The MtrCDE system confers multidrug resistance to Neisseria gonorrhoeae, the causative agent of gonorrhea. Using free and directed molecular dynamics (MD) simulations, we analyzed the interactions between MtrD and azithromycin, a transport substrate of MtrD, and a last-resort clinical treatment for multidrug-resistant gonorrhea. We then simulated the interactions between MtrD and streptomycin, an apparent nonsubstrate of MtrD. Using known conformations of MtrD homologues, we simulated a potential dynamic transport cycle of MtrD using targeted MD techniques (TMD), and we noted that forces were not applied to ligands of interest. In these TMD simulations, we observed the transport of azithromycin and the rejection of streptomycin. In an unbiased, long-time scale simulation of AZY-bound MtrD, we observed the spontaneous diffusion of azithromycin through the periplasmic cleft. Our simulations show how the peristaltic motions of the periplasmic cleft facilitate the transport of substrates by MtrD. Our data also suggest that multiple transport pathways for macrolides may exist within the periplasmic cleft of MtrD.

Structural and energetic analysis of metastable intermediate states in the E1P-E2P transition of Ca2+-ATPase

Ion pumps (or P-type ATPases) are membrane proteins, which transport ions through biological membranes against a concentration gradient, a function essential for many biological processes, such as muscle contraction, neurotransmission, and metabolism. Molecular mechanisms underlying active ion transport by ion pumps have been investigated by biochemical experiments and high-resolution structure analyses. Here, the transition of sarcoplasmic reticulum Ca²⁺-ATPase upon dissociation of Ca²⁺ is investigated using atomistic molecular dynamics simulations. We find intermediate structures along the pathway are stabilized by transient interactions between A- and P-domains as well as lipid molecules in the transmembrane helices.

Sarcoplasmic reticulum (SR) Ca²⁺-ATPase transports two Ca²⁺ ions from the cytoplasm to the SR lumen against a large concentration gradient. X-ray crystallography has revealed the atomic structures of the protein before and after the dissociation of Ca²⁺, while biochemical studies have suggested the existence of intermediate states in the transition between E1P⋅ADP⋅2Ca²⁺ and E2P. Here, we explore the pathway and free energy profile of the transition using atomistic molecular dynamics simulations with the mean-force string method and umbrella sampling. The simulations suggest that a series of structural changes accompany the ordered dissociation of ADP, the A-domain rotation, and the rearrangement of the transmembrane (TM) helices. The luminal gate then opens to release Ca²⁺ ions toward the SR lumen. Intermediate structures on the pathway are stabilized by transient sidechain interactions between the A- and P-domains. Lipid molecules between TM helices play a key role in the stabilization. Free energy profiles of the transition assuming different protonation states suggest rapid exchanges between Ca²⁺ ions and protons when the Ca²⁺ ions are released toward the SR lumen.

Multiscale aspects of molecular motions: From molecular vibrations, conformational changes of biomolecules to cellular dynamics

Molecular aspects of living systems are important because it is the most basic aspects of life as exemplified in biochemistry and structural biology. Since molecules move due to interactive forces between atoms, physics plays an important role to understand the dynamic phenomena of living systems. Here we review our multiscale approaches for computationally treating different levels of molecular motions: vibrational dynamics of molecules, conformational change of biomolecules, and cellular dynamics using statistical-mechanics-based models.

Structure and function relationship of OqxB efflux pump from Klebsiella pneumoniae

OqxB is an RND (Resistance-Nodulation-Division) efflux pump that has emerged as a factor contributing to the antibiotic resistance in Klebsiella pneumoniae. OqxB underwent horizontal gene transfer and is now seen in other Gram-negative bacterial pathogens including Escherichia coli, Enterobacter cloacae and Salmonella spp., further disseminating multi-drug resistance. In this study, we describe crystal structure of OqxB with n-dodecyl-β-D-maltoside (DDM) molecules bound in its substrate-binding pocket, at 1.85 Å resolution. We utilize this structure in computational studies to predict the key amino acids contributing to the efflux of fluoroquinolones by OqxB, distinct from analogous residues in related transporters AcrB and MexB. Finally, our complementation assays with mutated OqxB and minimum inhibitory concentration (MIC) experiments with clinical isolates of E. coli provide further evidence that the predicted structural features are indeed involved in ciprofloxacin efflux.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

Insight into the AcrAB-TolC Complex Assembly Process Learned from Competition Studies

The RND family efflux pump AcrAB-TolC in E. coli and its homologs in other Gram-negative bacteria are major players in conferring multidrug resistance to the cells. While the structure of the pump complex has been elucidated with ever-increasing resolution through crystallography and Cryo-EM efforts, the dynamic assembly process remains poorly understood. Here, we tested the effect of overexpressing functionally defective pump components in wild type E. coli cells to probe the pump assembly process. Incorporation of a defective component is expected to reduce the efflux efficiency of the complex, leading to the so called "dominant negative" effect. Being one of the most intensively studied bacterial multidrug efflux pumps, many AcrA and AcrB mutations have been reported that disrupt efflux through different mechanisms. We examined five groups of AcrB and AcrA mutants, defective in different aspects of assembly and substrate efflux. We found that none of them demonstrated the expected dominant negative effect, even when expressed at concentrations many folds higher than their genomic counterpart. The assembly of the AcrAB-TolC complex appears to have a proof-read mechanism that effectively eliminated the formation of futile pump complex.

Coordination of Substrate Binding and Protonation in the N. gonorrhoeae MtrD Efflux Pump Controls the Functionally Rotating Transport Mechanism

Multidrug resistance is a serious problem that threatens the effective treatment of the widespread sexually transmitted disease gonorrhea, caused by the Gram-negative bacterium Neisseria gonorrhoeae. The drug efflux pump primarily implicated in N. gonorrhoeae antimicrobial resistance is the inner membrane transporter MtrD, which forms part of the tripartite multiple transferable resistance (Mtr) CDE efflux system. A structure of MtrD was first solved in 2014 as a symmetrical homotrimer, and then, recently, as an asymmetrical homotrimer. Through a series of molecular dynamics simulations and mutagenesis experiments, we identify the combination of substrate binding and protonation states of the proton relay network that drives the transition from the symmetric to the asymmetric conformation of MtrD. We characterize the allosteric coupling between the functionally important local regions that control conformational changes between the access, binding, and extrusion states and allow for transition to the asymmetric MtrD conformation. We also highlight a significant rotation of the transmembrane helices caused by protonation of the proton relay network, which widens the intermonomeric gap that is a hallmark of the rotational transporter mechanism. This is the first analysis and description of the transport mechanism for the N. gonorrhoeae MtrD transporter and provides evidence that antimicrobial efflux in MtrD follows the functionally rotating transport mechanism seen in protein homologues from the same transport protein superfamily.

Structure, Assembly, and Function of Tripartite Efflux and Type 1 Secretion Systems in Gram-Negative Bacteria

Tripartite efflux pumps and the related type 1 secretion systems (T1SSs) in Gram-negative organisms are diverse in function, energization, and structural organization. They form continuous conduits spanning both the inner and the outer membrane and are composed of three principal components-the energized inner membrane transporters (belonging to ABC, RND, and MFS families), the outer membrane factor channel-like proteins, and linking the two, the periplasmic adaptor proteins (PAPs), also known as the membrane fusion proteins (MFPs). In this review we summarize the recent advances in understanding of structural biology, function, and regulation of these systems, highlighting the previously undescribed role of PAPs in providing a common architectural scaffold across diverse families of transporters. Despite being built from a limited number of basic structural domains, these complexes present a staggering variety of architectures. While key insights have been derived from the RND transporter systems, a closer inspection of the operation and structural organization of different tripartite systems reveals unexpected analogies between them, including those formed around MFS- and ATP-driven transporters, suggesting that they operate around basic common principles. Based on that we are proposing a new integrated model of PAP-mediated communication within the conformational cycling of tripartite systems, which could be expanded to other types of assemblies.

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

Perturbed structural dynamics underlie inhibition and altered efflux of the multidrug resistance pump AcrB

Resistance–nodulation–division efflux pumps play a key role in inherent and evolved multidrug resistance in bacteria. AcrB, a prototypical member of this protein family, extrudes a wide range of antimicrobial agents out of bacteria. Although high-resolution structures exist for AcrB, its conformational fluctuations and their putative role in function are largely unknown. Here, we determine these structural dynamics in the presence of substrates using hydrogen/deuterium exchange mass spectrometry, complemented by molecular dynamics simulations, and bacterial susceptibility studies. We show that an efflux pump inhibitor potentiates antibiotic activity by restraining drug-binding pocket dynamics, rather than preventing antibiotic binding. We also reveal that a drug-binding pocket substitution discovered within a multidrug resistant clinical isolate modifies the plasticity of the transport pathway, which could explain its altered substrate efflux. Our results provide insight into the molecular mechanism of drug export and inhibition of a major multidrug efflux pump and the directive role of its dynamics.

AcrB is a prototypical resistance–nodulation–division (RND) bacterial transporter, conferring resistance to a variety of antibiotics. HDX-MS and other, complementary approaches offer insight into AcrB structural dynamics and suggest the molecular mechanisms underlying drug export and inhibition of this multidrug-resistance conferring pump.

Antibiotic export by MexB multidrug efflux transporter is allosterically controlled by a MexA-OprM chaperone-like complex

The tripartite multidrug efflux system MexAB-OprM is a major actor in Pseudomonas aeruginosa antibiotic resistance by exporting a large variety of antimicrobial compounds. Crystal structures of MexB and of its Escherichia coli homolog AcrB had revealed asymmetric trimers depicting a directional drug pathway by a conformational interconversion (from Loose and Tight binding pockets to Open gate (LTO) for drug exit). It remains unclear how MexB acquires its LTO form. Here by performing functional and cryo-EM structural investigations of MexB at various stages of the assembly process, we unveil that MexB inserted in lipid membrane is not set for active transport because it displays an inactive LTC form with a Closed exit gate. In the tripartite complex, OprM and MexA form a corset-like platform that converts MexB into the active form. Our findings shed new light on the resistance nodulation cell division (RND) cognate partners which act as allosteric factors eliciting the functional drug extrusion.

Cryo-EM analysis of a membrane protein embedded in the liposome

Membrane proteins (MPs) used to be the most difficult targets for structural biology when X-ray crystallography was the mainstream approach. With the resolution revolution of single-particle electron cryo-microscopy (cryo-EM), rapid progress has been made for structural elucidation of isolated MPs. The next challenge is to preserve the electrochemical gradients and membrane curvature for a comprehensive structural elucidation of MPs that rely on these chemical and physical properties for their biological functions. Toward this goal, here we present a convenient workflow for cryo-EM structural analysis of MPs embedded in liposomes, using the well-characterized AcrB as a prototype. Combining optimized proteoliposome isolation, cryo-sample preparation on graphene grids, and an efficient particle selection strategy, the three-dimensional (3D) reconstruction of AcrB embedded in liposomes was obtained at 3.9 Å resolution. The conformation of the homotrimeric AcrB remains the same when the surrounding membranes display different curvatures. Our approach, which can be widely applied to cryo-EM analysis of MPs with distinctive soluble domains, lays out the foundation for cryo-EM analysis of integral or peripheral MPs whose functions are affected by transmembrane electrochemical gradients or/and membrane curvatures.

Cryo-EM Structure and Molecular Dynamics Analysis of the Fluoroquinolone Resistant Mutant of the AcrB Transporter from Salmonella

Salmonella is an important genus of Gram-negative pathogens, treatment of which has become problematic due to increases in antimicrobial resistance. This is partly attributable to the overexpression of tripartite efflux pumps, particularly the constitutively expressed AcrAB-TolC. Despite its clinical importance, the structure of the Salmonella AcrB transporter remained unknown to-date, with much of our structural understanding coming from the Escherichia coli orthologue. Here, by taking advantage of the styrene maleic acid (SMA) technology to isolate membrane proteins with closely associated lipids, we report the very first experimental structure of Salmonella AcrB transporter. Furthermore, this novel structure provides additional insight into mechanisms of drug efflux as it bears the mutation (G288D), originating from a clinical isolate of Salmonella Typhimurium presenting an increased resistance to fluoroquinolones. Experimental data are complemented by state-of-the-art molecular dynamics (MD) simulations on both the wild type and G288D variant of Salmonella AcrB. Together, these reveal several important differences with respect to the E. coli protein, providing insights into the role of the G288D mutation in increasing drug efflux and extending our understanding of the mechanisms underlying antibiotic resistance.

Free Energy Analysis of a Conformational Change of Heme ABC Transporter BhuUV-T

The heme ATP-binding cassette (ABC) transporter BhuUV-T of bacterial pathogen Burkholderia cenocepacia is required to transport heme across the inner cell membrane. The current hypothesis is that the binding of two ATPs to the nucleotide-binding domains of the transporter drives the initial steps of the transport cycle in which the empty transport sites are reoriented from the cytosol to the periplasm. Molecular details are missing because the structure of a key occluded intermediate remains hypothetical. Here we perform molecular simulations to analyze the free energy surface (FES) of the first step of the reorientation, namely the transition from an open inward-facing (IF) transport site to an occluded (Occ) conformation. We have modeled the latter structure in silico in a previous study. A simple annealing procedure removes residual bias originating from non-equilibrium targeted molecular dynamics. The calculated FES reveals the role of the ATPs in inducing the IF → Occ conformational change and validates the modeled Occ conformation.

Computational Dissection of Membrane Transport at a Microscopic Level

Membrane transporters are key gatekeeper proteins at cellular membranes that closely control the traffic of materials. Their function relies on structural rearrangements of varying degrees that facilitate substrate translocation across the membrane. Characterizing these functionally important molecular events at a microscopic level is key to our understanding of membrane transport, yet challenging to achieve experimentally. Recent advances in simulation technology and computing power have rendered molecular dynamics (MD) simulation a powerful biophysical tool to investigate a wide range of dynamical events spanning multiple spatial and temporal scales. Here, we review recent studies of diverse membrane transporters using computational methods, with an emphasis on highlighting the technical challenges, key lessons learned, and new opportunities to illuminate transporter structure and function.

Molecular Interactions of Carbapenem Antibiotics with the Multidrug Efflux Transporter AcrB of Escherichia coli

The drug/proton antiporter AcrB, engine of the major efflux pump AcrAB(Z)-TolC of Escherichia coli and other bacteria, is characterized by its impressive ability to transport chemically diverse compounds, conferring a multi-drug resistance (MDR) phenotype. Although hundreds of small molecules are known to be AcrB substrates, only a few co-crystal structures are available to date. Computational methods have been therefore intensively employed to provide structural and dynamical fingerprints related to transport and inhibition of AcrB. In this work, we performed a systematic computational investigation to study the interaction between representative carbapenem antibiotics and AcrB. We focused on the interaction of carbapenems with the so-called distal pocket, a region known for its importance in binding inhibitors and substrates of AcrB. Our findings reveal how the different physico-chemical nature of these antibiotics is reflected on their binding preference for AcrB. The molecular-level information provided here could help design new antibiotics less susceptible to the efflux mechanism.

Proton Control of Transitions in an Amino Acid Transporter

Amino acid transport into the cell is often coupled to the proton electrochemical gradient, as found in the solute carrier 36 family of proton-coupled amino acid transporters. Although no structure of a human proton-coupled amino acid transporter exists, the crystal structure of a related homolog from bacteria, GkApcT, has recently been solved in an inward-occluded state and allows an opportunity to examine how protons are coupled to amino acid transport. Our working hypothesis is that release of the amino acid substrate is facilitated by the deprotonation of a key glutamate residue (E115) located at the bottom of the binding pocket, which forms part of the intracellular gate, allowing the protein to transition from an inward-occluded to an inward-open conformation. During unbiased molecular dynamics simulations, we observed a transition from the inward-occluded state captured in the crystal structure to a much more open state, which we consider likely to be representative of the inward-open state associated with substrate release. To explore this and the role of protons in these transitions, we have used umbrella sampling to demonstrate that the transition from inward occluded to inward open is more energetically favorable when E115 is deprotonated. That E115 is likely to be protonated in the inward-occluded state and deprotonated in the inward-open state is further confirmed via the use of absolute binding free energies. Finally, we also show, via the use of absolute binding free energy calculations, that the affinity of the protein for alanine is very similar regardless of either the conformational state or the protonation of E115, presumably reflecting the fact that all the key interactions are deep within the binding cavity. Together, our results give a detailed picture of the role of protons in driving one of the major transitions in this transporter.

Molecular Interactions of Cephalosporins with the Deep Binding Pocket of the RND Transporter AcrB

The drug/proton antiporter AcrB, part of the major efflux pump AcrABZ-TolC in Escherichia coli, is characterized by its impressive ability to transport chemically diverse compounds, conferring a multidrug resistance phenotype. However, the molecular features differentiating between good and poor substrates of the pump have yet to be identified. In this work, we combined molecular docking with molecular dynamics simulations to study the interactions between AcrB and two representative cephalosporins, cefepime and ceftazidime (a good and poor substrate of AcrB, respectively). Our analysis revealed different binding preferences of the two compounds toward the subsites of the large deep binding pocket of AcrB. Cefepime, although less hydrophobic than ceftazidime, showed a higher affinity than ceftazidime for the so-called hydrophobic trap, a region known for binding inhibitors and substrates. This supports the hypothesis that surface complementarity between the molecule and AcrB, more than the intrinsic hydrophobicity of the antibiotic, is a feature required for the interaction within this region. Oppositely, the preference of ceftazidime for binding outside the hydrophobic trap might not be optimal for triggering allosteric conformational changes needed to the transporter to accomplish its function. Altogether, our findings could provide valuable information for the design of new antibiotics less susceptible to the efflux mechanism.

Identification and characterization of carbapenem binding sites within the RND-transporter AcrB

Understanding the molecular determinants for recognition, binding and transport of antibiotics by multidrug efflux systems is important for basic research and useful for the design of more effective antimicrobial compounds. Imipenem and meropenem are two carbapenems whose antibacterial activity is known to be poorly and strongly affected by MexAB-OprM, the major efflux pump transporter in Pseudomonas aeruginosa. However, not much is known regarding recognition and transport of these compounds by AcrAB-TolC, which is the MexAB-OprM homologue in Escherichia coli and by definition the paradigm model for structural studies on efflux pumps. Prompted by this motivation, we unveiled the molecular details of the interaction of imipenem and meropenem with the transporter AcrB by combining computer simulations with biophysical experiments. Regarding the interaction with the two main substrate binding regions of AcrB, the so-called access and deep binding pockets, molecular dynamics simulations revealed imipenem to be more mobile than meropenem in the former, while comparable mobilities were observed in the latter. This result is in line with isothermal titration calorimetry, differential scanning experiments, and binding free energy calculations, indicating a higher affinity for meropenem than imipenem at the deep binding pocket, while both sharing similar affinities at the access pocket. Our findings rationalize how different physico-chemical properties of compounds reflect on their interactions with AcrB. As such, they constitute precious information to be exploited for the rational design of antibiotics able to evade efflux pumps.

Exploring Configuration Space and Path Space of Biomolecules Using Enhanced Sampling Techniques-Searching for Mechanism and Kinetics of Biomolecular Functions

To understand functions of biomolecules such as proteins, not only structures but their conformational change and kinetics need to be characterized, but its atomistic details are hard to obtain both experimentally and computationally. Here, we review our recent computational studies using novel enhanced sampling techniques for conformational sampling of biomolecules and calculations of their kinetics. For efficiently characterizing the free energy landscape of a biomolecule, we introduce the multiscale enhanced sampling method, which uses a combined system of atomistic and coarse-grained models. Based on the idea of Hamiltonian replica exchange, we can recover the statistical properties of the atomistic model without any biases. We next introduce the string method as a path search method to calculate the minimum free energy pathways along a multidimensional curve in high dimensional space. Finally we introduce novel methods to calculate kinetics of biomolecules based on the ideas of path sampling: one is the Onsager⁻Machlup action method, and the other is the weighted ensemble method. Some applications of the above methods to biomolecular systems are also discussed and illustrated.

Molecular Determinants of the Promiscuity of MexB and MexY Multidrug Transporters of Pseudomonas aeruginosa

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.

Molecular mechanisms of AcrB-mediated multidrug export

The over-expression of multidrug efflux pumps belonging to the Resistance-Nodulation-Division (RND) superfamily is one of the main causes of multidrug-resistance (MDR) in Gram-negative pathogenic bacteria. AcrB is the most thoroughly studied RND transporter and has functioned as a model for our understanding of efflux-mediated MDR. This multidrug-exporter can recognize and transport a wide range of structurally unrelated compounds (including antibiotics, dyes, bile salts and detergents), while it shows a strict inhibitor specificity. Here we discuss our current knowledge of AcrB, and include recent advances, regarding its structure, mechanism of drug transport, substrate recognition, different intramolecular entry pathways and the drug export driven by remote conformational coupling.