Drug export and allosteric coupling in a multidrug transporter revealed by molecular simulations

Multiple-Basin Go̅-Martini for Investigating Conformational Transitions and Environmental Interactions of Proteins

Proteins are inherently dynamic molecules, and their conformational transitions among various states are essential for numerous biological processes, which are often modulated by their interactions with surrounding environments. Although molecular dynamics (MD) simulations are widely used to investigate these transitions, all-atom (AA) methods are often limited by short time scales and high computational costs, and coarse-grained (CG) implicit-solvent Go̅-like models are usually incapable of studying the interactions between proteins and their environments. Here, we present an approach called Multiple-basin Go̅-Martini, which combines the recent Go̅-Martini model with an exponential mixing scheme to facilitate the simulation of spontaneous protein conformational transitions in explicit environments. We demonstrate the versatility of our method through five diverse case studies: GlnBP, Arc, Hinge, SemiSWEET, and TRAAK, representing ligand-binding proteins, fold-switching proteins, de novo designed proteins, transporters, and mechanosensitive ion channels, respectively. Multiple-basin Go̅-Martini offers a new computational tool for investigating protein conformational transitions, identifying key intermediate states, and elucidating essential interactions between proteins and their environments, particularly protein-membrane interactions. In addition, this approach can efficiently generate thermodynamically meaningful data sets of protein conformational space, which may enhance deep learning-based models for predicting protein conformation distributions.

Allosteric Regulation of Enzymatic Catalysis through Mechanical Force

Mechanical force has been increasingly recognized to play crucial roles in regulating various cellular processes, which has inspired wide interest in elucidating the biophysical mechanism underlying these mechanobiological processes. In this work, we investigate the mechanical regulation of enzyme catalysis by developing a residue-resolved computational model capable of describing the full catalytic cycle of enzymes under mechanical force. Intriguingly, for a model enzyme, adenylate kinase, we showed that applying tensile forces with biologically relevant strength can increase the enzymatic activity. Further analysis showed that mechanical tensile force allosterically modifies the global free energy landscape and conformational dynamics of the protein, which then promotes the rate-limiting product release step of the enzymatic cycle. The effect of mechanical allostery on enzyme catalysis depends on the intrinsic conformational propensity of the enzymes. The crucial role of mechanical allostery in enzymatic catalysis elucidated in this work sheds important insights into the biophysical principle of enzymatic regulation and suggests a possible strategy for fine-tuning the functioning dynamics of biological enzymes.

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.

Exploring imidazo[4,5-g]quinoline-4,9-dione derivatives as Acinetobacter baumannii efflux pump inhibitor: an in silico approach

Antimicrobial resistance (AMR) is fast becoming a medical crisis affecting the entire global population. World Health Organization (WHO) statistics show that globally 0.7 million people are dying yearly due to the emergence of AMR. By 2050, the expected number of lives lost will be 10 million per year. Acinetobacter baumannii is a dreadful nosocomial pathogen that has developed multidrug resistance (MDR) to several currently prescribed antibiotics worldwide. Overexpression of drug efflux transporters (DETs) is one of the mechanisms of multidrug resistance (MDR) in Acinetobacter baumannii. Therefore, blocking the DET can raise the efficacy of the existing antibiotics by increasing their residence time inside the bacteria. In silico screening of five synthetic compounds against three drug efflux pump from A. baumannii has identified KSA5, a novel imidazo[4,5-g]quinoline-4,9-dione derivative, to block the efflux of antibiotics. Molecular docking and simulation results showed that KSA5 could bind to adeB, adeG, and adeJ by consistently interacting with ligand-binding site residues. KSA5 has a higher binding free energy and a lower HOMO-LUMO energy gap than PAβN, suggesting a better ability to interact and inhibit DETs. Further analysis showed that KSA5 is a drug-like molecule with optimal physicochemical and ADME properties. Hence, KSA5 could be combined with antibiotics to overcome antimicrobial resistance.Communicated by Ramaswamy H. Sarma.

Switching Go̅-Martini for Investigating Protein Conformational Transitions and Associated Protein-Lipid Interactions

Proteins are dynamic biomolecules that can transform between different conformational states when exerting physiological functions, which is difficult to simulate using all-atom methods. Coarse-grained (CG) Go̅-like models are widely used to investigate large-scale conformational transitions, which usually adopt implicit solvent models and therefore cannot explicitly capture the interaction between proteins and surrounding molecules, such as water and lipid molecules. Here, we present a new method, named Switching Go̅-Martini, to simulate large-scale protein conformational transitions between different states, based on the switching Go̅ method and the CG Martini 3 force field. The method is straightforward and efficient, as demonstrated by the benchmarking applications for multiple protein systems, including glutamine binding protein (GlnBP), adenylate kinase (AdK), and β2-adrenergic receptor (β2AR). Moreover, by employing the Switching Go̅-Martini method, we can not only unveil the conformational transition from the E2Pi-PL state to E1 state of the type 4 P-type ATPase (P4-ATPase) flippase ATP8A1-CDC50 but also provide insights into the intricate details of lipid transport.

Molecular insights into the determinants of substrate specificity and efflux inhibition of the RND efflux pumps AcrB and AdeB

Gram-negative bacterial members of the Resistance Nodulation and cell Division (RND) superfamily form tripartite efflux pump systems that span the cell envelope. One of the intriguing features of the multiple drug efflux members of this superfamily is their ability to recognize different classes of antibiotics, dyes, solvents, bile salts, and detergents. This review provides an overview of the molecular mechanisms of multiple drug efflux catalysed by the tripartite RND efflux system AcrAB-TolC from Eschericha coli. The determinants for sequential or simultaneous multiple substrate binding and efflux pump inhibitor binding are discussed. A comparison is made with the determinants for substrate binding of AdeB from Acinetobacter baumannii, which acts within the AdeABC multidrug efflux system. There is an apparent general similarity between the structures of AcrB and AdeB and their substrate specificity. However, the presence of distinct conformational states and different drug efflux capacities as revealed by single-particle cryo-EM and mutational analysis suggest that the drug binding and transport features exhibited by AcrB may not be directly extrapolated to the homolog AdeB efflux pump.

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.

Role of Repeated Conformational Transitions in Substrate Binding of Adenylate Kinase

The catalytic cycle of the enzyme adenylate kinase involves large conformational motions between open and closed states. A previous single-molecule experiment showed that substrate binding tends to accelerate both the opening and the closing rates and that a single turnover event often involves multiple rounds of conformational switching. In this work, we showed that the repeated conformational transitions of adenylate kinase are essential for the relaxation of incorrectly bound substrates into the catalytically competent conformation by combining all-atom and coarse-grained molecular simulations. In addition, free energy calculations based on all-atom and coarse-grained models demonstrated that the enzyme with incorrectly bound substrates has much a lower free energy barrier for domain opening compared to that with the correct substrate conformation, which may explain the the acceleration of the domain opening rate by substrate binding. The results of this work provide mechanistic understanding to previous experimental observations and shed light onto the interplay between conformational dynamics and enzyme catalysis.

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.

A second shell residue modulates a conserved ATP-binding site with radically different affinities for ATP

BACKGROUND

Prediction of ligand binding and design of new function in enzymes is a time-consuming and expensive process. Crystallography gives the impression that proteins adopt a fixed shape, yet enzymes are functionally dynamic. Molecular dynamics offers the possibility of probing protein movement while predicting ligand binding. Accordingly, we choose the bacterial F₁F_o ATP synthase ε subunit to unravel why ATP affinity by ε subunits from Bacillus subtilis and Bacillus PS3 differs ~500-fold, despite sharing identical sequences at the ATP-binding site.

METHODS

We first used the Bacillus PS3 ε subunit structure to model the B. subtilis ε subunit structure and used this to explore the utility of molecular dynamics (MD) simulations to predict the influence of residues outside the ATP binding site. To verify the MD predictions, point mutants were made and ATP binding studies were employed.

RESULTS

MD simulations predicted that E102 in the B. subtilis ε subunit, outside of the ATP binding site, influences ATP binding affinity. Engineering E102 to alanine or arginine revealed a ~10 or ~54 fold increase in ATP binding, respectively, confirming the MD prediction that E102 drastically influences ATP binding affinity.

CONCLUSIONS

These findings reveal how MD can predict how changes in the "second shell" residues around substrate binding sites influence affinity in simple protein structures. Our results reveal why seemingly identical ε subunits in different ATP synthases have radically different ATP binding affinities.

GENERAL SIGNIFICANCE

This study may lead to greater utility of molecular dynamics as a tool for protein design and exploration of protein design and function.

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.

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.

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.

Exploring Successful Parameter Region for Coarse-Grained Simulation of Biomolecules by Bayesian Optimization and Active Learning

Accompanied with an increase of revealed biomolecular structures owing to advancements in structural biology, the molecular dynamics (MD) approach, especially coarse-grained (CG) MD suitable for macromolecules, is becoming increasingly important for elucidating their dynamics and behavior. In fact, CG-MD simulation has succeeded in qualitatively reproducing numerous biological processes for various biomolecules such as conformational changes and protein folding with reasonable calculation costs. However, CG-MD simulations strongly depend on various parameters, and selecting an appropriate parameter set is necessary to reproduce a particular biological process. Because exhaustive examination of all candidate parameters is inefficient, it is important to identify successful parameters. Furthermore, the successful region, in which the desired process is reproducible, is essential for describing the detailed mechanics of functional processes and environmental sensitivity and robustness. We propose an efficient search method for identifying the successful region by using two machine learning techniques, Bayesian optimization and active learning. We evaluated its performance using F1-ATPase, a biological rotary motor, with CG-MD simulations. We successfully identified the successful region with lower computational costs (12.3% in the best case) without sacrificing accuracy compared to exhaustive search. This method can accelerate not only parameter search but also biological discussion of the detailed mechanics of functional processes and environmental sensitivity based on MD simulation studies.

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.

Building a macro-mixing dual-basin Gō model using the Multistate Bennett Acceptance Ratio

The dual-basin Gō-model is a structural-based coarsegrained model for simulating a conformational transition between two known structures of a protein. Two parameters are required to produce a dual-basin potential mixed using two single-basin potentials, although the determination of mixing parameters is usually not straightforward. Here, we have developed an efficient scheme to determine the mixing parameters using the Multistate Bennett Acceptance Ratio (MBAR) method after short simulations with a set of parameters. In the scheme, MBAR allows us to predict observables at various unsimulated conditions, which are useful to improve the mixing parameters in the next round of iterative simulations. The number of iterations that are necessary for obtaining the converged mixing parameters are significantly reduced in the scheme. We applied the scheme to two proteins, the glutamine binding protein and the ribose binding protein, for showing the effectiveness in the parameter determination. After obtaining the converged parameters, both proteins show frequent conformational transitions between open and closed states, providing the theoretical basis to investigate structure-dynamics-function relationships of the proteins.

Rotational Mechanism Model of the Bacterial V1 Motor Based on Structural and Computational Analyses

V₁-ATPase exemplifies the ubiquitous rotary motor, in which a central shaft DF complex rotates inside a hexagonally arranged catalytic A₃B₃ complex, powered by the energy from ATP hydrolysis. We have recently reported a number of crystal structures of the Enterococcus hirae A₃B₃DF (V₁) complex corresponding to its nucleotide-bound intermediate states, namely the forms waiting for ATP hydrolysis (denoted as catalytic dwell), ATP binding (ATP-binding dwell), and ADP release (ADP-release dwell) along the rotatory catalytic cycle of ATPase. Furthermore, we have performed microsecond-scale molecular dynamics simulations and free-energy calculations to investigate the conformational transitions between these intermediate states and to probe the long-time dynamics of the molecular motor. In this article, the molecular structure and dynamics of the V₁-ATPase are reviewed to bring forth a unified model of the motor’s remarkable rotational mechanism.

Chromatin remodelers couple inchworm motion with twist-defect formation to slide nucleosomal DNA

ATP-dependent chromatin remodelers are molecular machines that control genome organization by repositioning, ejecting, or editing nucleosomes, activities that confer them essential regulatory roles on gene expression and DNA replication. Here, we investigate the molecular mechanism of active nucleosome sliding by means of molecular dynamics simulations of the Snf2 remodeler translocase in complex with a nucleosome. During its inchworm motion driven by ATP consumption, the translocase overwrites the original nucleosome energy landscape via steric and electrostatic interactions to induce sliding of nucleosomal DNA unidirectionally. The sliding is initiated at the remodeler binding location via the generation of a pair of twist defects, which then spontaneously propagate to complete sliding throughout the entire nucleosome. We also reveal how remodeler mutations and DNA sequence control active nucleosome repositioning, explaining several past experimental observations. These results offer a detailed mechanistic picture of remodeling important for the complete understanding of these key biological processes.

Nucleosomes are the protein-DNA complexes underlying Eukaryotic genome organization, and serve as regulators of gene expression by occluding DNA to other proteins. This regulation requires the precise positioning of nucleosomes along DNA. Chromatin remodelers are the molecular machines that consume ATP to slide nucleosome at their correct locations, but the mechanisms of remodeling are still unclear. Based on the static structural information of a remodeler bound on nucleosome, we performed molecular dynamics computer simulations revealing the details of how remodelers slide nucleosomal DNA: the inchworm-like motion of remodelers create small DNA deformations called twist defects, which then spontaneously propagate throughout the nucleosome to induce sliding. These simulations explain several past experimental findings and are important for our understanding of genome organization.

Mapping the Dynamic Functions and Structural Features of AcrB Efflux Pump Transporter Using Accelerated Molecular Dynamics Simulations

Multidrug efflux pumps confer resistance to their bacterial hosts by pumping out a diverse range of compounds, including most antibiotics. Being more familiar with the details of functional dynamics and conformations of these types of pumps could help in discovering approaches to stop them functioning properly. Computational approaches, particularly conventional molecular dynamics simulations followed by diverse post simulation analysis, are powerful methods that help researchers by opening a new window to study phenomena that are not detectable in as much detail in vitro or in vivo as they are in silico. In this study, accelerated molecular dynamics simulations were applied to study the dynamics of AcrB efflux pump transporters in interaction with PAβN and tetracycline as an inhibitor and a substrate, respectively, to compare the differences in the dynamics and consequently the mechanism of action of the pump. The different dynamics for PAβN -bound form of AcrB compared to the TET-bound form is likely to affect the rotating mechanism typically observed for AcrB transporter. This shows the dynamics of the active AcrB transporter is different in a substrate-bound state compared to an inhibitor-bound state. This advances our knowledge and helps to unravel the mechanism of tripartite efflux pumps.

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

The multidrug transporter AcrB transports a broad range of drugs out of the cell by means of the proton-motive force. The asymmetric crystal structure of trimeric AcrB suggests a functionally rotating mechanism for drug transport. Despite various supportive forms of evidence from biochemical and simulation studies for this mechanism, the link between the functional rotation and proton translocation across the membrane remains elusive. Here, calculating the minimum free energy pathway of the functional rotation for the complete AcrB trimer, we describe the structural and energetic basis behind the coupling between the functional rotation and the proton translocation at atomic resolution. Free energy calculations show that protonation of Asp408 in the transmembrane portion of the drug-bound protomer drives the functional rotation. The conformational pathway identifies vertical shear motions among several transmembrane helices, which regulate alternate access of water in the transmembrane as well as peristaltic motions that pump drugs in the periplasm.

Exploration of Multistate Conformational Dynamics upon Ligand Binding of a Monomeric Enzyme Involved in Pyrophosphoryl Transfer

HPPK (6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase) is a monomeric protein with 158 residues, which undergoes large-scale conformational changes between apo, open, and holo states responding to ligand binding for its function. It has been explored widely as an excellent target for potential antibacterial drug development. However, little is known about how conformational dynamics between the native states influences the substrate recognition and the functionality of enzymatic catalysis. Here, we report a coarse-grained triple-basin structure-based model upon ligand binding to describe such multiple-state system by the molecular dynamics simulation. With our model, we have made theoretical predictions that are in good agreement with the experimental measurements. Our results revealed the intrinsic conformational fluctuations between apo and open states without ligand binding. We found that HPPK can switch to the activated holo state upon the ordered binding of the two ligands (ATP and HP). We uncovered the underlying mechanism by which major induced fit and minor population shift pathways coexist upon ligand binding by quantitative flux analysis. Additionally, we pointed out the structural origin for the conformational changes and identified the key residues as well as contact interactions. We further explored the temperature effect on the conformational distributions and pathway weights. It gave strong support that higher temperatures promote population shift, while the induced fit pathway is always the predominant activation route of the HPPK system. These findings will provide significant insights of the mechanisms of the multistate conformational dynamics of HPPK upon ligand binding.

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.

Water-mediated interactions enable smooth substrate transport in a bacterial efflux pump

BACKGROUND

Efflux pumps of the Resistance-Nodulation-cell Division superfamily confer multi-drug resistance to Gram-negative bacteria. The most-studied polyspecific transporter belonging to this class is the inner-membrane trimeric antiporter AcrB of Escherichia coli. In previous studies, a functional rotation mechanism was proposed for its functioning, according to which the three monomers undergo concerted conformational changes facilitating the extrusion of substrates. However, the molecular determinants and the energetics of this mechanism still remain unknown, so its feasibility must be proven mechanistically.

METHODS

A computational protocol able to mimic the functional rotation mechanism in AcrB was developed. By using multi-bias molecular dynamics simulations we characterized the translocation of the substrate doxorubicin driven by conformational changes of the protein. In addition, we estimated for the first time the free energy profile associated to this process.

RESULTS

We provided a molecular view of the process in agreement with experimental data. Moreover, we showed that the conformational changes occurring in AcrB enable the formation of a layer of structured waters on the internal surface of the transport channel. This water layer, in turn, allows for a fairly constant hydration of the substrate, facilitating its diffusion over a smooth free energy profile.

CONCLUSIONS

Our findings reveal a new molecular mechanism of polyspecific transport whereby water contributes by screening potentially strong substrate-protein interactions.

GENERAL SIGNIFICANCE

We provided a mechanistic understanding of a fundamental process related to multi-drug transport. Our results can help rationalizing the behavior of other polyspecific transporters and designing compounds avoiding extrusion or inhibitors of efflux pumps.

Multiscale molecular dynamics simulations of rotary motor proteins

Protein functions require specific structures frequently coupled with conformational changes. The scale of the structural dynamics of proteins spans from the atomic to the molecular level. Theoretically, all-atom molecular dynamics (MD) simulation is a powerful tool to investigate protein dynamics because the MD simulation is capable of capturing conformational changes obeying the intrinsically structural features. However, to study long-timescale dynamics, efficient sampling techniques and coarse-grained (CG) approaches coupled with all-atom MD simulations, termed multiscale MD simulations, are required to overcome the timescale limitation in all-atom MD simulations. Here, we review two examples of rotary motor proteins examined using free energy landscape (FEL) analysis and CG-MD simulations. In the FEL analysis, FEL is calculated as a function of reaction coordinates, and the long-timescale dynamics corresponding to conformational changes is described as transitions on the FEL surface. Another approach is the utilization of the CG model, in which the CG parameters are tuned using the fluctuation matching methodology with all-atom MD simulations. The long-timespan dynamics is then elucidated straightforwardly by using CG-MD simulations.

Dynamics of Intact MexAB-OprM Efflux Pump: Focusing on the MexA-OprM Interface

Antibiotic efflux is one of the most critical mechanisms leading to bacterial multidrug resistance. Antibiotics are effluxed out of the bacterial cell by a tripartite efflux pump, a complex machinery comprised of outer membrane, periplasmic adaptor, and inner membrane protein components. Understanding the mechanism of efflux pump assembly and its dynamics could facilitate discovery of novel approaches to counteract antibiotic resistance in bacteria. We built here an intact atomistic model of the Pseudomonas aeruginosa MexAB-OprM pump in a Gram-negative membrane model that contained both inner and outer membranes separated by a periplasmic space. All-atom molecular dynamics (MD) simulations confirm that the fully assembled pump is stable in the microsecond timescale. Using a combination of all-atom and coarse-grained MD simulations and sequence covariation analysis, we characterized the interface between MexA and OprM in the context of the entire efflux pump. These analyses suggest a plausible mechanism by which OprM is activated via opening of its periplasmic aperture through a concerted interaction with MexA.

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

The multidrug resistance (MDR) system actively pumps antibiotics out of cells causing serious health problems. During the pumping, AcrB (one of the key components of MDR) undergoes a series of large-scale and proton-motive conformational changes. Capturing the conformational changes through all-atom simulations is challenging. Here, we implement a hybrid coarse-grained force field to investigate the conformational changes of AcrB in the porter domain under different protonation states of Asp407/Asp408 in the trans-membrane domain. Our results show that protonation of Asp408 in monomer III (extrusion) stabilizes the asymmetric structure of AcrB; deprotonation of Asp408 induces clear opening of the entrance and closing of the exit leading to the transition from extrusion to access state. The structural changes in the porter domain of AcrB are strongly coupled with the proton translocation stoichiometry in the trans-membrane domain. Moreover, our simulations support the postulation that AcrB should adopt the symmetric resting state in a substrate-free situation.

Rotation Mechanism of Molecular Motor V1-ATPase Studied by Multiscale Molecular Dynamics Simulation

Enterococcus hirae V₁-ATPase is a molecular motor composed of the A₃B₃ hexamer ring and the central stalk. In association with ATP hydrolysis, three catalytic AB pairs in the A₃B₃ ring undergo conformational changes, which lead to a 120° rotation of the central stalk. To understand how the conformational changes of three catalytic pairs induce the 120° rotation of the central stalk, we performed multiscale molecular dynamics (MD) simulations in which coarse-grained and all-atom MD simulations were combined using a fluctuation matching methodology. During the rotation, a catalytic AB pair spontaneously adopted an intermediate conformation, which was not included in the initial inputs of the simulations and was essentially close to the “bindable-like” structure observed in a recently solved crystal structure. Furthermore, the creation of a space between the bindable-like and tight pairs was required for the central stalk to rotate without steric hindrance. These cooperative rearrangements of the three catalytic pairs are crucial for the rotation of the central stalk.

Efflux Pump Blockers in Gram-Negative Bacteria: The New Generation of Hydantoin Based-Modulators to Improve Antibiotic Activity

Multidrug resistant (MDR) bacteria are an increasing health problem with the shortage of new active antibiotic agents. Among effective mechanisms that contribute to the spread of MDR Gram-negative bacteria are drug efflux pumps that expel clinically important antibiotic classes out of the cell. Drug pumps are attractive targets to restore the susceptibility toward the expelled antibiotics by impairing their efflux activity. Arylhydantoin derivatives were investigated for their potentiation of activities of selected antibiotics described as efflux substrates in Enterobacter aerogenes expressing or not AcrAB pump. Several compounds increased the bacterial susceptibility toward nalidixic acid, chloramphenicol and sparfloxacin and were further pharmacomodulated to obtain a better activity against the AcrAB producing bacteria.

Insights into the Inhibitory Mechanism of D13-9001 to the Multidrug Transporter AcrB through Molecular Dynamics Simulations

The resistance-nodulation-cell division transporter AcrB is responsible for energy transduction and substrate recognition in the tripartite AcrAB-TolC efflux system in Escherichia coli. Despite a broad substrate specificity, only a few compounds have been cocrystallized with AcrB inside the distal binding pocket (DBP), including doxorubicin (DOX) and D13-9001. D13-9001 is a promising efflux pump inhibitor that potentiates the efficacy of a wide variety of antibiotics. To understand its inhibition effect under the framework of functional rotating mechanism, we performed targeted and steered molecular dynamics simulations to compare the binding and extrusion processes of this inhibitor and the substrate DOX in AcrB. The results demonstrate that, with respect to DOX, the interaction of D13-9001 with the hydrophobic trap results in delayed disassociation from the DBP. Notably, the detachment of D13-9001 is tightly correlated with the side-chain reorientation of Phe628 and large-scale displacement of Tyr327. Furthermore, the inhibitor induces much more significant conformational changes at the exit gate than DOX does, thereby causing higher energy cost for extrusion and contributing to the inhibitory effect in addition to the tight binding at DBP.

Substrate-dependent dynamics of the multidrug efflux transporter AcrB of Escherichia coli

The resistance-nodulation-cell division (RND)-type xenobiotic efflux system plays a major role in the multidrug resistance of gram-negative bacteria. The only constitutively expressed RND system of Escherichia coli consists of the inner membrane transporter AcrB, the membrane fusion protein AcrA, and the outer membrane channel TolC. The latter two components are shared with another RND-type transporter AcrD, whose expression is induced by environmental stimuli. Here, we demonstrate how RND-type ternary complexes, which span two membranes and the cell wall, form in vivo. Total internal reflection fluorescence (TIRF) microscopy revealed that most fluorescent foci formed by AcrB fused to green fluorescent protein (GFP) were stationary in the presence of TolC but showed lateral displacements when tolC was deleted. The fraction of stationary AcrB-GFP foci decreased with increasing levels of AcrD. We propose that the AcrB-containing complex becomes unstable upon the induction of AcrD, which presumably replaces AcrB, a process we call “transporter exchange.” This instability is suppressed by AcrB-specific substrates, suggesting that the ternary complex is stabilised when it is in action. These results suggest that the assembly of the RND-type efflux system is dynamically regulated in response to external stimuli, shedding new light on the adaptive antibiotic resistance of bacteria.

An overview of bacterial efflux pumps and computational approaches to study efflux pump inhibitors

Micro-organisms express a wide range of transmembrane pumps known as multidrug efflux pumps that improve the micro-organism's ability to survive in severe environments and contribute to resistance against antibiotic and antimicrobial agents. There is significant interest in developing efflux inhibitors as an adjunct to treatment with current and next generation of antibiotics. A greater understanding of drug recognition and transport by multidrug efflux pumps is needed to develop clinically useful inhibitors, given the breadth of molecules that can be effluxed by these systems. We summarize some structural and functional data that could provide insights into the inhibition of transport mechanisms of these intricate molecular nanomachines with a focus on the advances in computational approaches.

Domain Motion Enhanced (DoME) Model for Efficient Conformational Sampling of Multidomain Proteins

Large conformational changes of multidomain proteins are difficult to simulate using all-atom molecular dynamics (MD) due to the slow time scale. We show that a simple modification of the structure-based coarse-grained (CG) model enables a stable and efficient MD simulation of those proteins. "Motion Tree", a tree diagram that describes conformational changes between two structures in a protein, provides information on rigid structural units (domains) and the magnitudes of domain motions. In our new CG model, which we call the DoME (domain motion enhanced) model, interdomain interactions are defined as being inversely proportional to the magnitude of the domain motions in the diagram, whereas intradomain interactions are kept constant. We applied the DoME model in combination with the Go model to simulations of adenylate kinase (AdK). The results of the DoME-Go simulation are consistent with an all-atom MD simulation for 10 μs as well as known experimental data. Unlike the conventional Go model, the DoME-Go model yields stable simulation trajectories against temperature changes and conformational transitions are easily sampled despite domain rigidity. Evidently, identification of domains and their interfaces is useful approach for CG modeling of multidomain proteins.

Stepwise substrate translocation mechanism revealed by free energy calculations of doxorubicin in the multidrug transporter AcrB

AcrB is the inner membrane transporter of the tripartite multidrug efflux pump AcrAB-TolC in E. coli, which poses a major obstacle to the treatment of bacterial infections. X-ray structures have identified two types of substrate-binding pockets in the porter domains of AcrB trimer: the proximal binding pocket (PBP) and the distal binding pocket (DBP), and suggest a functional rotating mechanism in which each protomer cycles consecutively through three distinct conformational states (access, binding and extrusion). However, the details of substrate binding and translocation between the binding pockets remain elusive. In this work, we performed atomic simulations to obtain the free energy profile of the translocation of an antibiotic drug doxorubicin (DOX) inside AcrB. Our simulation indicates that DOX binds at the PBP and DBP with comparable affinities in the binding state protomer, and overcomes a 3 kcal/mol energy barrier to transit between them. Obvious conformational changes including closing of the PC1/PC2 cleft and shrinking of the DBP were observed upon DOX binding in the PBP, resulting in an intermediate state between the access and binding states. Taken together, the simulation results reveal a detailed stepwise substrate binding and translocation process in the framework of functional rotating mechanism.

Structure, mechanism and cooperation of bacterial multidrug transporters

Cells from all domains of life encode energy-dependent trans-membrane transporters that can expel harmful substances including clinically applied therapeutic agents. As a collective body, these transporters perform as a super-system that confers tolerance to an enormous range of harmful compounds and consequently aid survival in hazardous environments. In the Gram-negative bacteria, some of these transporters serve as energy-transducing components of tripartite assemblies that actively efflux drugs and other harmful compounds, as well as deliver virulence agents across the entire cell envelope. We draw together recent structural and functional data to present the current models for the transport mechanisms for the main classes of multi-drug transporters and their higher-order assemblies.

Recent advances toward a molecular mechanism of efflux pump inhibition

Multidrug resistance (MDR) in Gram-negative pathogens, such as the Enterobacteriaceae and Pseudomonas aeruginosa, poses a significant threat to our ability to effectively treat infections caused by these organisms. A major component in the development of the MDR phenotype in Gram-negative bacteria is overexpression of Resistance-Nodulation-Division (RND)-type efflux pumps, which actively pump antibacterial agents and biocides from the periplasm to the outside of the cell. Consequently, bacterial efflux pumps are an important target for developing novel antibacterial treatments. Potent efflux pump inhibitors (EPIs) could be used as adjunctive therapies that would increase the potency of existing antibiotics and decrease the emergence of MDR bacteria. Several potent inhibitors of RND-type efflux pump have been reported in the literature, and at least three of these EPI series were optimized in a pre-clinical development program. However, none of these compounds have been tested in the clinic. One of the major hurdles to the development of EPIs has been the lack of biochemical, computational, and structural methods that could be used to guide rational drug design. Here, we review recent reports that have advanced our understanding of the mechanism of action of several potent EPIs against RND-type pumps.

The challenge of efflux-mediated antibiotic resistance in Gram-negative bacteria

The global emergence of multidrug-resistant Gram-negative bacteria is a growing threat to antibiotic therapy. The chromosomally encoded drug efflux mechanisms that are ubiquitous in these bacteria greatly contribute to antibiotic resistance and present a major challenge for antibiotic development. Multidrug pumps, particularly those represented by the clinically relevant AcrAB-TolC and Mex pumps of the resistance-nodulation-division (RND) superfamily, not only mediate intrinsic and acquired multidrug resistance (MDR) but also are involved in other functions, including the bacterial stress response and pathogenicity. Additionally, efflux pumps interact synergistically with other resistance mechanisms (e.g., with the outer membrane permeability barrier) to increase resistance levels. Since the discovery of RND pumps in the early 1990s, remarkable scientific and technological advances have allowed for an in-depth understanding of the structural and biochemical basis, substrate profiles, molecular regulation, and inhibition of MDR pumps. However, the development of clinically useful efflux pump inhibitors and/or new antibiotics that can bypass pump effects continues to be a challenge. Plasmid-borne efflux pump genes (including those for RND pumps) have increasingly been identified. This article highlights the recent progress obtained for organisms of clinical significance, together with methodological considerations for the characterization of MDR pumps.

Assembly and operation of bacterial tripartite multidrug efflux pumps

Microorganisms encode several classes of transmembrane pumps that can expel an enormous range of toxic substances, thereby improving their fitness in harsh environments and contributing to resistance against antimicrobial agents. In Gram-negative bacteria these pumps can take the form of tripartite assemblies that actively efflux drugs and other harmful compounds across the cell envelope. We describe recent structural and functional data that have provided insights into the transport mechanisms of these intricate molecular machines.

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.

Assembly and stability of Salmonella enterica ser. Typhi TolC protein in POPE and DMPE

In this work we assessed the suitability of two different lipid membranes for the simulation of a TolC protein from Salmonella enterica serovar Typhi. The TolC protein family is found in many pathogenic Gram-negative bacteria including Vibrio cholera and Pseudomonas aeruginosa and acts as an outer membrane channel for expulsion of drug and toxin from the cell. In S. typhi, the causative agent for typhoid fever, the TolC outer membrane protein is an antigen for the pathogen. The lipid environment is an important modulator of membrane protein structure and function. We evaluated the conformation of the TolC protein in the presence of DMPE and POPE bilayers using molecular dynamics simulation. The S. typhi TolC protein exhibited similar conformational dynamics to TolC and its homologues. Conformational flexibility of the protein is seen in the C-terminal, extracellular loops, and α-helical region. Despite differences in the two lipids, significant similarities in the motion of the protein in POPE and DMPE were observed, including the rotational motion of the C-terminal residues and the partially open extracellular loops. However, analysis of the trajectories demonstrated effects of hydrophobic matching of the TolC protein in the membrane, particularly in the lengthening of the lipids and subtle movements of the protein's β-barrel towards the lower leaflet in DMPE. The study exhibited the use of molecular dynamics simulation in revealing the differential effect of membrane proteins and lipids on each other. In this study, POPE is potentially a more suitable model for future simulation of the S. typhi TolC protein.

Towards understanding promiscuity in multidrug efflux pumps

Drug export from cells is a major factor in the acquisition of cellular resistance to antimicrobial and cancer chemotherapy, and poses a significant threat to future clinical management of disease. Many of the proteins that catalyse drug efflux do so with remarkably low substrate specificity, a phenomenon known as multidrug transport. For these reasons we need a greater understanding of drug recognition and transport in multidrug pumps to inform research that attempts to circumvent their action. Structural and computational studies have been heralded as being great strides towards a full elucidation of multidrug recognition and transport. In this review we summarise these advances and ask how close we are to a molecular understanding of this remarkable phenomenon.

Modulation of Bacterial Multidrug Resistance Efflux Pumps of the Major Facilitator Superfamily

Bacterial infections pose a serious public health concern, especially when an infectious disease has a multidrug resistant causative agent. Such multidrug resistant bacteria can compromise the clinical utility of major chemotherapeutic antimicrobial agents. Drug and multidrug resistant bacteria harbor several distinct molecular mechanisms for resistance. Bacterial antimicrobial agent efflux pumps represent a major mechanism of clinical resistance. The major facilitator superfamily (MFS) is one of the largest groups of solute transporters to date and includes a significant number of bacterial drug and multidrug efflux pumps. We review recent work on the modulation of multidrug efflux pumps, paying special attention to those transporters belonging primarily to the MFS.

Surface-Tension Replica-Exchange Molecular Dynamics Method for Enhanced Sampling of Biological Membrane Systems

Conformational sampling is fundamentally important for simulating complex biomolecular systems. The generalized-ensemble algorithm, especially the temperature replica-exchange molecular dynamics method (T-REMD), is one of the most powerful methods to explore structures of biomolecules such as proteins, nucleic acids, carbohydrates, and also of lipid membranes. T-REMD simulations have focused on soluble proteins rather than membrane proteins or lipid bilayers, because explicit membranes do not keep their structural integrity at high temperature. Here, we propose a new generalized-ensemble algorithm for membrane systems, which we call the surface-tension REMD method. Each replica is simulated in the NPγT ensemble, and surface tensions in a pair of replicas are exchanged at certain intervals to enhance conformational sampling of the target membrane system. We test the method on two biological membrane systems: a fully hydrated DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine) lipid bilayer and a WALP23-POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane system. During these simulations, a random walk in surface tension space is realized. Large-scale lateral deformation (shrinking and stretching) of the membranes takes place in all of the replicas without collapse of the lipid bilayer structure. There is accelerated lateral diffusion of DPPC lipid molecules compared with conventional MD simulation, and a much wider range of tilt angle of the WALP23 peptide is sampled due to large deformation of the POPC lipid bilayer and through peptide-lipid interactions. Our method could be applicable to a wide variety of biological membrane systems.

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.