Coarse-grained simulations of conformational changes in the multidrug efflux transporter AcrB

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.

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.

Molecular Dynamics Investigation of MFS Efflux Pump MdfA Reveals an Intermediate State between Its Inward and Outward Conformations

Multidrug resistance poses a major challenge to antibiotic therapy. A principal cause of antibiotic resistance is through active export by efflux pumps embedded in the bacterial membrane. Major facilitator superfamily (MFS) efflux pumps constitute a major group of transporters, which are often related to quinolone resistance in clinical settings. Although a rocker-switch model is proposed for description of their conformational transitions, detailed changes in this process remain poorly understood. Here we used MdfA from E. coli as a representative MFS efflux pump to investigate factors that can affect its conformational transition in silico. Molecular dynamics (MD) simulations of MdfA's inward and outward conformations revealed an intermediate state between these two conformations. By comparison of the subtle differences between the intermediate state and the average state, we indicated that conformational transition from outward to inward was initiated by protonation of the periplasmic side. Subsequently, hydrophilic interaction of the periplasmic side with water was promoted and the regional structure of helix 1 was altered to favor this process. As the hydrophobic interaction between MdfA and membrane was also increased, energy was concentrated and stored for the opposite transition. In parallel, salt bridges at the cytoplasmic side were altered to lower probabilities to facilitate the entrance of substrate. In summary, we described the total and local changes during MdfA's conformational transition, providing insights for the development of potential inhibitors.

Role of TM3 in claudin-15 strand flexibility: A molecular dynamics study

Claudins are cell-cell adhesion proteins within tight junctions that connect epithelial cells together. Claudins polymerize into a network of strand-like structures within the membrane of adjoining cells and create ion channels that control paracellular permeability to water and small molecules. Tight junction morphology and barrier function is tissue specific and regulated by claudin subtypes. Here, we present a molecular dynamics study of claudin-15 strands within lipid membranes and the role of a single-point mutation (A134P) on the third transmembrane helix (TM3) of claudin-15 in determining the morphology of the strand. Our results indicate that the A134P mutation significantly affects the lateral flexibility of the strands, increasing the persistence length of claudin-15 strands by a factor of three. Analyses of claudin-claudin contact in our μsecond-long trajectories show that the mutation does not alter the intermolecular contacts (interfaces) between claudins. However, the dynamics and frequency of interfacial contacts are significantly affected. The A134P mutation introduces a kink in TM3 of claudin-15 similar to the one observed in claudin-3 crystal structure. The kink on TM3 skews the rotational flexibility of the claudins in the strands and limits their fluctuation in one direction. This asymmetric movement in the context of the double rows reduces the lateral flexibility of the strand and leads to higher persistence lengths of the mutant.

What have molecular simulations contributed to understanding of Gram-negative bacterial cell envelopes?

Bacterial cell envelopes are compositionally complex and crowded and while highly dynamic in some areas, their molecular motion is very limited, to the point of being almost static in others. Therefore, it is no real surprise that studying them at high resolution across a range of temporal and spatial scales requires a number of different techniques. Details at atomistic to molecular scales for up to tens of microseconds are now within range for molecular dynamics simulations. Here we review how such simulations have contributed to our current understanding of the cell envelopes of Gram-negative bacteria.

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.

Effect of Calcium ion on synaptotagmin-like protein during pre-fusion of vesicle for exocytosis in blood-brain barrier

Background

Calcium signaling and membrane fusion play key roles in exocytosis of drug-containing vesicles through the blood-brain barrier (BBB). Identifying the role of synaptotagmin-like protein4-a (Slp4-a) in the presence of Ca²⁺ ions, at the pre-fusion stage of a vesicle with the basolateral membrane of endothelial cell, can reveal brain drug transportation across BBB.

Methods

We utilized molecular dynamics (MD) simulations with a coarse-grained PACE force field to investigate the behaviors of Slp4-a with vesicular and endothelial membranes at the pre-fusion stage of exocytosis since all-atom MD simulation or experiments are more time-consuming and expensive to capture these behaviors.

Results

The Slp4-a pulls lipid membranes (vesicular and endothelial) into close proximity and disorganizes lipid arrangement at contact points, which are predictors for initiation of fusion. Our MD results also indicate that Slp4-a needs Ca²⁺ to bind with weakly-charged POPE lipids (phosphatidylethanolamine).

Conclusions

Slp4-a is an important trigger for membrane fusion in BBB exocytosis. It binds to lipid membranes at multiple binding sites and triggers membrane disruption for fusion in calcium-dependent case.

General significance

Understanding the prefusion process of the vesicle will help to design better drug delivery mechanisms to the brain through formidable BBB.

•

Role of Ca²⁺ on Slp4-a is studied for vesicle pre-fusion in EC to initiate exocytosis.

•

Coarse-grained MD is used to study large scale conformation change of Slp-4a.

•

Interaction between C2A domain and lipids is much stronger than that of C2B.

•

Slp4-a can bind to bilayer membrane in Ca²⁺-bound case to close membrane gap.

Studying dynamics without explicit dynamics: A structure-based study of the export mechanism by AcrB

Resistance-nodulation-cell division family proteins are transmembrane proteins identified as large spectrum drug transporters involved in multidrug resistance. A prototypical case in this superfamily, responsible for antibiotic resistance in selected gram-negative bacteria, is AcrB. AcrB forms a trimer using the proton motive force to efflux drugs, implementing a functional rotation mechanism. Unfortunately, the size of the system (1049 amino acid per monomer and membrane) has prevented a systematic dynamical exploration, so that the mild understanding of this coupled transport jeopardizes our ability to counter it. The large number of crystal structures of AcrB prompts studies to further our understanding of the mechanism. To this end, we present a novel strategy based on two key ingredients, which are to study dynamics by exploiting information embodied in the numerous crystal structures obtained to date, and to systematically consider subdomains, their dynamics, and their interactions. Along the way, we identify the subdomains responsible for dynamic events, refine the states (A, B, E) of the functional rotation mechanism, and analyze the evolution of intramonomer and intermonomer interfaces along the functional cycle. Our analysis shows the relevance of AcrB's efflux mechanism as a template within the HAE1 family but not beyond. It also paves the way to targeted simulations exploiting the most relevant degrees of freedom at certain steps, and to a targeting of specific interfaces to block the drug efflux. Our work shows that complex dynamics can be unveiled from static snapshots, a strategy that may be used on a variety of molecular machines of large size.

Chlorpromazine and Amitriptyline Are Substrates and Inhibitors of the AcrB Multidrug Efflux Pump

Efflux is an important mechanism in Gram-negative bacteria conferring multidrug resistance. Inhibition of efflux is an encouraging strategy to restore the antibacterial activity of antibiotics. Chlorpromazine and amitriptyline have been shown to behave as efflux inhibitors. However, their mode of action is poorly understood. Exposure of Salmonella enterica serovar Typhimurium and Escherichia coli to chlorpromazine selected for mutations within genes encoding RamR and MarR, regulators of the multidrug tripartite efflux pump AcrAB-TolC. Further experiments with S. Typhimurium containing AcrB D408A (a nonfunctional efflux pump) and chlorpromazine or amitriptyline resulted in the reversion of the mutant acrB allele to the wild type. Together, this suggests these drugs are AcrB efflux substrates. Subsequent docking studies with AcrB from S. Typhimurium and E. coli, followed by molecular dynamics simulations and free energy calculations showed that chlorpromazine and amitriptyline bind at the hydrophobic trap, a preferred binding site for substrates and inhibitors within the distal binding pocket of AcrB. Based on these simulations, we suggest that chlorpromazine and amitriptyline inhibit AcrB-mediated efflux by interfering with substrate binding. Our findings provide evidence that these drugs are substrates and inhibitors of AcrB, yielding molecular details of their mechanism of action and informing drug discovery of new efflux inhibitors.IMPORTANCE Efflux pumps of the resistance nodulation-cell division (RND) superfamily are major contributors to multidrug resistance for most of the Gram-negative ESKAPE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) pathogens. The development of inhibitors of these pumps would be highly desirable; however, several issues have thus far hindered all efforts at designing new efflux inhibitory compounds devoid of adverse effects. An alternative route to de novo design relies on the use of marketed drugs, for which side effects on human health have been already assessed. In this work, we provide experimental evidence that the antipsychotic drugs chlorpromazine and amitriptyline are inhibitors of the AcrB transporter, the engine of the major RND efflux pumps in Escherichia coli and Salmonella enterica serovar Typhimurium. Furthermore, in silico calculations have provided a molecular-level picture of the inhibition mechanism, allowing rationalization of experimental data and paving the way for similar studies with other classes of marketed compounds.

Substrate-dependent transport mechanism in AcrB of multidrug resistant bacteria

The multidrug resistance (MDR) system effectively expels antibiotics out of bacteria causing serious issues during bacterial infection. In addition to drug, indole, a common metabolic waste of bacteria, is expelled by MDR system of gram-negative bacteria for their survival. Experimental results suggest that AcrB, one of the key components of MDR system, undergoes large scale conformation changes during the pumping due to proton-motive process. However, due to extremely short time scale, it is difficult to observe (experimentally) those changes in the AcrB, which might facilitate the pumping process. Molecular simulations can shed light to understand the conformational changes for transport of indole in AcrB. Examination of conformational changes using all-atom simulation is, however, impractical. Here, we develop a hybrid coarse-grained force field to study the conformational changes of AcrB in presence of indole in the porter domain of monomer II. Using the coarse-grained force field, we investigated the conformational changes of AcrB for a number of model systems considering the effect of protonation in aspartic acid (Asp) residues Asp407 and Asp408 in the transmembrane domain of monomer II. Our results show that in the presence of indole, protonation of Asp408 or Asp407 residue causes conformational changes from binding state to extrusion state in monomer II, while remaining two monomers (I and III) approach access state in AcrB protein. We also observed that all three AcrB monomers prefer to go back to access state in the absence of indole. Steered molecular dynamics simulations were performed to demonstrate the feasibility of indole transport mechanism for protonated systems. Identification of indole transport pathway through AcrB can be very helpful in understanding the drug efflux mechanism used by the MDR bacteria.

Computer simulations of the activity of RND efflux pumps

The putative mechanism by which bacterial RND-type multidrug efflux pumps recognize and transport their substrates is a complex and fascinating enigma of structural biology. How a single protein can recognize a huge number of unrelated compounds and transport them through one or just a few mechanisms is an amazing feature not yet completely unveiled. The appearance of cooperativity further complicates the understanding of structure-dynamics-activity relationships in these complex machineries. Experimental techniques may have limited access to the molecular determinants and to the energetics of key processes regulating the activity of these pumps. Computer simulations are a complementary approach that can help unveil these features and inspire new experiments. Here we review recent computational studies that addressed the various molecular processes regulating the activity of RND efflux pumps.

Constant pH Molecular Dynamics Reveals How Proton Release Drives the Conformational Transition of a Transmembrane Efflux Pump

AcrB is the inner-membrane transporter of an E. coli AcrAB-TolC tripartite efflux complex, which plays a major role in the intrinsic resistance to clinically important antibiotics. AcrB pumps a wide range of toxic substrates by utilizing the proton gradient between periplasm and cytoplasm. Crystal structures of AcrB revealed three distinct conformational states of the transport cycle, substrate access, binding, and extrusion or loose (L), tight (T), and open (O) states. However, the specific residue(s) responsible for proton binding/release and the mechanism of proton-coupled conformational cycling remain controversial. Here we use the newly developed membrane hybrid-solvent continuous constant pH molecular dynamics technique to explore the protonation states and conformational dynamics of the transmembrane domain of AcrB. Simulations show that both Asp407 and Asp408 are deprotonated in the L/T states, while only Asp408 is protonated in the O state. Remarkably, release of a proton from Asp408 in the O state results in large conformational changes, such as the lateral and vertical movement of transmembrane helices as well as the salt-bridge formation between Asp408 and Lys940 and other side chain rearrangements among essential residues. Consistent with the crystallographic differences between the O and L protomers, simulations offer dynamic details of how proton release drives the O-to-L transition in AcrB and address the controversy regarding the proton/drug stoichiometry. This work offers a significant step toward characterizing the complete cycle of proton-coupled drug transport in AcrB and further validates the membrane hybrid-solvent CpHMD technique for studies of proton-coupled transmembrane proteins which are currently poorly understood.