Transitions between closed and open conformations of TolC: the effects of ions in simulations

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.

Conformational dynamics and putative substrate extrusion pathways of the N-glycosylated outer membrane factor CmeC from Campylobacter jejuni

The outer membrane factor CmeC of the efflux machinery CmeABC plays an important role in conferring antibiotic and bile resistance to Campylobacter jejuni. Curiously, the protein is N-glycosylated, with the glycans playing a key role in the effective function of this system. In this work we have employed atomistic equilibrium molecular dynamics simulations of CmeC in a representative model of the C. jejuni outer membrane to characterise the dynamics of the protein and its associated glycans. We show that the glycans are more conformationally labile than had previously been thought. The extracellular loops of CmeC visit the open and closed states freely suggesting the absence of a gating mechanism on this side, while the narrow periplasmic entrance remains tightly closed, regulated via coordination to solvated cations. We identify several cation binding sites on the interior surface of the protein. Additionally, we used steered molecular dynamics simulations to elucidate translocation pathways for a bile acid and a macrolide antibiotic. These, and additional equilibrium simulations suggest that the anionic bile acid utilises multivalent cations to climb the ladder of acidic residues that line the interior surface of the protein.

Redefining the bacterial Type I protein secretion system

Type I secretion systems (T1SS) are versatile molecular machines for protein transport across the Gram-negative cell envelope. The archetypal Type I system mediates secretion of the Escherichia coli hemolysin, HlyA. This system has remained the pre-eminent model of T1SS research since its discovery. The classic description of a T1SS is composed of three proteins: an inner membrane ABC transporter, a periplasmic adaptor protein and an outer membrane factor. According to this model, these components assemble to form a continuous channel across the cell envelope, an unfolded substrate molecule is then transported in a one-step mechanism, directly from the cytosol to the extracellular milieu. However, this model does not encapsulate the diversity of T1SS that have been characterized to date. In this review, we provide an updated definition of a T1SS, and propose the subdivision of this system into five subgroups. These subgroups are categorized as T1SSa for RTX proteins, T1SSb for non-RTX Ca²⁺-binding proteins, T1SSc for non-RTX proteins, T1SSd for class II microcins, and T1SSe for lipoprotein secretion. Although often overlooked in the literature, these alternative mechanisms of Type I protein secretion offer many avenues for biotechnological discovery and application.

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.

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.

In silico characterisation of olive phenolic compounds as potential cyclooxygenase modulators. Part 2

Non-steroidal anti-inflammatory drugs (NSAIDs) are commonly used to reduce pain, and function by targeting cyclooxygenase (COX) enzymes to inhibit the production of prostaglandins that facilitate inflammation. Since oleocanthal derived from Olea europaea is known to inhibit COX, we sought to characterise novel olive compounds with COX inhibitory activity using in silico techniques. Following on from part 1 of this study which identified 1-oleyltyrosol (1OL) and ligstroside derivative 2 (LG2) with COX inhibitory potential, the mechanisms of COX interactions by these selected compounds were further examined using molecular dynamics (MD) simulations. Classical MD simulations were carried out on COX-1 and COX-2 complexed with 1OL and LG2 to determine the stability and protein backbone fluctuation. Protein dynamics were examined using essential dynamics methods and network analysis, which identified that the N-terminal epidermal growth factor-like domain and membrane bound domains of COX-1 and -2 exhibited altered motions when ligands were bound. Distinct dynamical modules were identified, and that COX-2 inter-residue communications were more sensitive to ligand binding compared to COX-1. The use of various network metrics presents a novel approach in the characterisation of network behaviour of different ligands. It is proposed that inter-residue network metrics provide additional measures of the potential bioactivity of ligands, which may form a useful adjunct to conventional direct predictions of binding affinity, in determining the efficacy of potential small-molecule inhibitors. Overall, this two-part study characterises anti-inflammatory effects of low dosage dietary COX inhibitors, and provides a possible avenue for the development of therapeutics in inflammatory diseases.

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.

Modulation of the Neisseria gonorrhoeae drug efflux conduit MtrE

Widespread antibiotic resistance, especially of Gram-negative bacteria, has become a severe concern for human health. Tripartite efflux pumps are one of the major contributors to resistance in Gram-negative pathogens, by efficiently expelling a broad spectrum of antibiotics from the organism. In Neisseria gonorrhoeae, one of the first bacteria for which pan-resistance has been reported, the most expressed efflux complex is MtrCDE. Here we present the electrophysiological characterisation of the outer membrane component MtrE and the membrane fusion protein MtrC, obtained by a combination of planar lipid bilayer recordings and in silico techniques. Our in vitro results show that MtrE can be regulated by periplasmic binding events and that the interaction between MtrE and MtrC is sufficient to stabilize this complex in an open state. In contrast to other efflux conduits, the open complex only displays a slight preference for cations. The maximum conductance we obtain in the in vitro recordings is comparable to that seen in our computational electrophysiology simulations conducted on the MtrE crystal structure, indicating that this state may reflect a physiologically relevant open conformation of MtrE. Our results suggest that the MtrC/E binding interface is an important modulator of MtrE function, which could potentially be targeted by new efflux inhibitors.

Structural and dynamical correlations in PfHGXPRT oligomers: A molecular dynamics simulation study

PfHGXPRT is a key enzyme involved in purine nucleotide salvage pathway of the malarial parasite, Plasmodium falciparum. Atomistic molecular dynamics simulations have been performed on two types of PfHGXPRT dimers (D1 and D3) and its tetramer in their apo and ligand-bound states. A significant event in the catalytic cycle is the dynamics of a gate that provides access for the ligand molecules to the reaction center. The gate is formed by loops II and IV, the former being the most flexible. Large amplitude conformational changes have been observed in active site loop II. Upon complete occupancy of the active site, loop II gets stabilized due to specific interactions between its residues and the ligand molecules. Remote loop, X, is seen to be less fluxional in the D3 dimer than in D1 which is rationalized as due to the greater number of inter-subunit contacts in the former. The presence of ligand molecules in subunits of the tetramer further reduces the flexibility of loop X epitomizing a communication between this region and the active sites in the tetramer. These observations are in accordance with the outcomes of several experimental investigations. Participation of loop X in the oligomerization process has also been discerned. Between the two types of dimers in solution, D1 tetramerizes readily and thus would not be present as free dimers. We conjecture an equilibrium to exist between D3 and the tetramer in solution; upon binding of the ligand molecules to the D3 dimer, this equilibrium shifts toward the tetramer.

Computer simulations suggest direct and stable tip to tip interaction between the outer membrane channel TolC and the isolated docking domain of the multidrug RND efflux transporter AcrB

One way by which bacteria achieve antibiotics resistance is preventing drug access to its target molecule for example through an overproduction of multi-drug efflux pumps of the resistance nodulation division (RND) protein super family of which AcrAB-TolC in Escherichia coli is a prominent example. Although representing one of the best studied efflux systems, the question of how AcrB and TolC interact is still unclear as the available experimental data suggest that either both proteins interact in a tip to tip manner or do not interact at all but are instead connected by a hexamer of AcrA molecules. Addressing the question of TolC-AcrB interaction, we performed a series of 100 ns - 1 µs-molecular dynamics simulations of membrane-embedded TolC in presence of the isolated AcrB docking domain (AcrB(DD)). In 5/6 simulations we observe direct TolC-AcrB(DD) interaction that is only stable on the simulated time scale when both proteins engage in a tip to tip manner. At the same time we find TolC opening and closing freely on extracellular side while remaining closed at the inner periplasmic bottleneck region, suggesting that either the simulated time is too short or additional components are required to unlock TolC.

Living on the edge: Simulations of bacterial outer-membrane proteins

Gram-negative bacteria are distinguished in part by a second, outer membrane surrounding them. This membrane is distinct from others, possessing an outer leaflet composed not of typical phospholipids but rather large, highly charged molecules known as lipopolysaccharides. Therefore, modeling the structure and dynamics of proteins embedded in the outer membrane requires careful consideration of their native environment. In this review, we examine how simulations of such outer-membrane proteins have evolved over the last two decades, culminating most recently in detailed, highly accurate atomistic models of the outer membrane. We also draw attention to how the simulations have coupled with experiments to produce novel insights unattainable through a single approach. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Simulations of outer membrane channels and their permeability

Channels in the outer membrane of Gram-negative bacteria provide essential pathways for the controlled and unidirectional transport of ions, nutrients and metabolites into the cell. At the same time the outer membrane serves as a physical barrier for the penetration of noxious substances such as antibiotics into the bacteria. Most antibiotics have to pass through these membrane channels to either reach cytoplasmic bound targets or to further cross the hydrophobic inner membrane. Considering the pharmaceutical significance of antibiotics, understanding the functional role and mechanism of these channels is of fundamental importance in developing strategies to design new drugs with enhanced permeation abilities. Due to the biological complexity of membrane channels and experimental limitations, computer simulations have proven to be a powerful tool to investigate the structure, dynamics and interactions of membrane channels. Considerable progress has been made in computer simulations of membrane channels during the last decade. The goal of this review is to provide an overview of the computational techniques and their roles in modeling the transport across outer membrane channels. A special emphasis is put on all-atom molecular dynamics simulations employed to better understand the transport of molecules. Moreover, recent molecular simulations of ion, substrate and antibiotics translocation through membrane pores are briefly summarized. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Crystal structure of an antigenic outer-membrane protein from Salmonella Typhi suggests a potential antigenic loop and an efflux mechanism

ST50, an outer-membrane component of the multi-drug efflux system from Salmonella enterica serovar Typhi, is an obligatory diagnostic antigen for typhoid fever. ST50 is an excellent and unique diagnostic antigen with 95% specificity and 90% sensitivity and is used in the commercial diagnosis test kit (TYPHIDOT(TM)). The crystal structure of ST50 at a resolution of 2.98 Å reveals a trimer that forms an α-helical tunnel and a β-barrel transmembrane channel traversing the periplasmic space and outer membrane. Structural investigations suggest significant conformational variations in the extracellular loop regions, especially extracellular loop 2. This is the location of the most plausible antibody-binding domain that could be used to target the design of new antigenic epitopes for the development of better diagnostics or drugs for the treatment of typhoid fever. A molecule of the detergent n-octyl-β-D-glucoside is observed in the D-cage, which comprises three sets of Asp361 and Asp371 residues at the periplasmic entrance. These structural insights suggest a possible substrate transport mechanism in which the substrate first binds at the periplasmic entrance of ST50 and subsequently, via iris-like structural movements to open the periplasmic end, penetrates the periplasmic domain for efflux pumping of molecules, including poisonous metabolites or xenobiotics, for excretion outside the pathogen.

Focus on the Outer Membrane Factor OprM, the Forgotten Player from Efflux Pumps Assemblies

Antibiotics have been used extensively during several decades and we are now facing the emergence of multidrug resistant strains. It has become a major public concern, urging the need to discover new strategies to combat them. Among the different ways used by bacteria to resist antibiotics, the active efflux is one of the main mechanisms. In Gram-negative bacteria the efflux pumps are comprised of three components forming a long edifice crossing the complete cell wall from the inside to the outside of the cell. Blocking these pumps would permit the restoration of the effectiveness of the current antibiotherapy which is why it is important to increase our knowledge on the different proteins involved in these complexes. A tremendous number of experiments have been performed on the inner membrane protein AcrB from Escherichia coli and, to a lesser extent, the protein partners forming the AcrAB-TolC pump, but less information is available concerning the efflux pumps from other virulent Gram-negative bacteria. The present review will focus on the OprM outer membrane protein from the MexAB-OprM pump of Pseudomonas aeruginosa, highlighting similarities and differences compare to the archetypal AcrAB-TolC in terms of structure, function, and assembly properties.

Interaction of antibacterial compounds with RND efflux pumps in Pseudomonas aeruginosa

Pseudomonas aeruginosa infections are becoming increasingly difficult to treat due to intrinsic antibiotic resistance and the propensity of this pathogen to accumulate diverse resistance mechanisms. Hyperexpression of efflux pumps of the Resistance-Nodulation-Cell Division (RND)-type multidrug efflux pumps (e.g., MexAB-OprM), chromosomally encoded by mexAB-oprM, mexCD-oprJ, mexEF-oprN, and mexXY (-oprA) is often detected in clinical isolates and contributes to worrying multi-drug resistance phenotypes. Not all antibiotics are affected to the same extent by the aforementioned RND efflux pumps. The impact of efflux on antibiotic activity varies not only between different classes of antibiotics but also between members of the same family of antibiotics. Subtle differences in physicochemical features of compound-pump and compound-solvent interactions largely determine how compounds are affected by efflux activity. The combination of different high-resolution techniques helps to gain insight into the functioning of these molecular machineries. This review discusses substrate recognition patterns based on experimental evidence and computer simulations with a focus on MexB, the pump subunit of the main RND transporter in P. aeruginosa.

Architecture and roles of periplasmic adaptor proteins in tripartite efflux assemblies

Recent years have seen major advances in the structural understanding of the different components of tripartite efflux assemblies, which encompass the multidrug efflux (MDR) pumps and type I secretion systems. The majority of these investigations have focused on the role played by the inner membrane transporters and the outer membrane factor (OMF), leaving the third component of the system – the Periplasmic Adaptor Proteins (PAPs) – relatively understudied. Here we review the current state of knowledge of these versatile proteins which, far from being passive linkers between the OMF and the transporter, emerge as active architects of tripartite assemblies, and play diverse roles in the transport process. Recognition between the PAPs and OMFs is essential for pump assembly and function, and targeting this interaction may provide a novel avenue for combating multidrug resistance. With the recent advances elucidating the drug efflux and energetics of the tripartite assemblies, the understanding of the interaction between the OMFs and PAPs is the last piece remaining in the complete structure of the tripartite pump assembly puzzle.

Computational study of correlated domain motions in the AcrB efflux transporter

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

Insight derived from molecular dynamics simulations into molecular motions, thermodynamics and kinetics of HIV-1 gp120

Although the crystal structures of the HIV-1 gp120 core bound and pre-bound by CD4 are known, the details of dynamics involved in conformational equilibrium and transition in relation to gp120 function have remained elusive. The homology models of gp120 comprising the N- and C-termini and loops V3 and V4 in the CD4-bound and CD4-unbound states were built and subjected to molecular dynamics (MD) simulations to investigate the differences in dynamic properties and molecular motions between them. The results indicate that the CD4-bound gp120 adopted a more compact and stable conformation than the unbound form during simulations. For both the unbound and bound gp120, the large concerted motions derived from essential dynamics (ED) analyses can influence the size/shape of the ligand-binding channel/cavity of gp120 and, therefore, were related to its functional properties. The differences in motion direction between certain structural components of these two forms of gp120 were related to the conformational interconversion between them. The free energy calculations based on the metadynamics simulations reveal a more rugged and complex free energy landscape (FEL) for the unbound than for the bound gp120, implying that gp120 has a richer conformational diversity in the unbound form. The estimated free energy difference of ∼−6.0 kJ/mol between the global minimum free energy states of the unbound and bound gp120 indicates that gp120 can transform spontaneously from the unbound to bound states, revealing that the bound state represents a high-probability “ground state” for gp120 and explaining why the unbound state resists crystallization. Our results provide insight into the dynamics-and-function relationship of gp120, and facilitate understandings of the thermodynamics, kinetics and conformational control mechanism of HIV-1 gp120.

Free energy profiles of ion permeation and doxorubicin translocation in TolC

The outer membrane protein TolC of Escherichia coli forms a channel-tunnel pore spanning the periplasmic space and outer membrane, serving as the main exit duct for bacteria multidrug resistance and protein export. Many aspects of the transport mechanism of TolC are still unclear. Here, we have investigated the substrate permeability and gating mechanism of TolC by calculating the potential of mean forces (PMFs) for transporting sodium ion and doxorubicin through TolC using the adaptive biasing force (ABF) method. The transport mechanism is turned out to be substrate dependent. It is found that the periplasmic gate is required to open for the passage of both Na + and doxorubicin, but the conformational gating does not lead to permeation barrier for Na + at this region. The extracellular loops and K283 residues cause permeation barriers for Na + at the extracellular entrance, but not for doxorubicin due to the extensive interactions between the drug molecule and the protein. TolC exhibits high conformational flexibility during the transport of Na +, while doxorubicin seems to be able to stabilize TolC in the resting state with the periplasmic gate closed. The association of the TolC docking domain of AcrB does not lower the permeation barrier for doxorubicin at the periplasmic gate, while the gate opening induces the dissociation of the TolC–AcrB complex.

Molecular Dynamics Computer Simulations of Multidrug RND Efflux Pumps

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

Unilateral access regulation: ground state dynamics of the Pseudomonas aeruginosa outer membrane efflux duct OprM

Acting as an efflux duct in the MexA-MexB-OprM multidrug efflux pump, OprM plays a major role in the antibiotic resistance capability of Pseudomonas aeruginosa, trafficking substrates through the outer cell membrane. Whereas the available crystal structures showed restricted OprM access on both ends, the underlying gating mechanism is not yet fully understood. To gain insight into the functional mechanism of OprM access regulation, we conducted a series of five independent, unbiased molecular dynamics simulations, computing 200 ns dynamics samples of the wild-type protein in a phospholipid membrane/150 mM NaCl water environment. On the extracellular side, OprM opens and closes freely under the simulated conditions, suggesting the absence of a gating mechanism on this side of the isolated protein. On the periplasmic side, we observe an opening of the tip regions at Val408 and to a lesser degree Asp416 located 1.5 nm further into the channel, leading to OprM end conformations being up to 3 and 1.4 times, respectively, more open than the asymmetric crystal structure. If our simulations are correct, our findings imply that periplasmic gating involves only the Asp416 region and that in vivo additional components, absent in our simulation, might be required for periplasmic gating if the observed opening trend near Asp416 is not negligible. In addition to that ,we identified in each monomer a previously unreported sodium binding site in the channel interior coordinated by Asp171 and Asp230 whose functional role remains to be investigated.

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

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

Locked on one side only: ground state dynamics of the outer membrane efflux duct TolC

Playing a major role in the expulsion of antibiotics and the secretion of cell toxins in conjunction with inner membrane transporters of three protein superfamilies, the outer membrane channel TolC occurs in at least two states blocking or permitting the passage of substrates. The details of the underlying gating mechanism are not fully understood. Addressing the questions of extracellular access control and periplasmic gating mechanism, we conducted a series of independent, unbiased 150-300 ns molecular dynamics simulations of wild-type TolC in a phospholipid membrane/150 mM NaCl water environment. We find that TolC opens and closes freely on the extracellular side, suggesting the absence of a gating mechanism on this side in the isolated protein. On the periplasmic side, we observe the outer periplasmic bottleneck region adopting in all simulations a conformation more open than the TolC wild-type crystal structures until in one run the successive binding of two sodium ions induces the transition to a conformation more closed than any of the available TolC X-ray structures. Concurrent with a heightened sodium residence probability near Asp374, the inner periplasmic bottleneck region at Asp374 remains closed throughout the simulations unless all NaCl is removed from the system, inducing a reopening of the outer and inner bottleneck. Our findings suggest that TolC is locked only on the periplasmic side in a sodium-dependent manner.

Elastic network model-based normal mode analysis reveals the conformational couplings in the tripartite AcrAB-TolC multidrug efflux complex

The AcrAB-TolC drug efflux system, energized by proton movement down the transmembrane electrochemical gradient, is responsible for the resistance of the organism to a wide range of drugs. Experimental data suggest functional roles of each part of the assembly, but the detailed working mechanism of this machinery remains elusive. We used elastic network-based normal mode analysis (NMA) to explore the conformational dynamics of the AcrAB-TolC complex. The intrinsic flexibilities of the pore domain in AcrB monomer conform to the previously proposed three-step functionally rotating mechanism for asymmetric AcrB trimer. Conformational couplings across monomers in the AcrB trimer were observed, and the coupling between the transmembrane domain and the other parts of AcrB are strengthened through trimeric assembly. In the tripartite AcrAB-TolC assembly obtained through molecular docking, concerted motions were observed not only at the direct contact interfaces between various components but also between distant parts of the whole complex. The presence of AcrA was shown to significantly strengthen the motional couplings between AcrB and TolC. Overall, NMA revealed an allosteric network in the AcAB-TolC efflux system, which provides hints to our understanding of its detailed working mechanism.

Structural and functional studies of a 50 kDa antigenic protein from Salmonella enterica serovar Typhi

The high typhoid incidence rate in developing and under-developed countries emphasizes the need for a rapid, affordable and accessible diagnostic test for effective therapy and disease management. TYPHIDOT®, a rapid dot enzyme immunoassay test for typhoid, was developed from the discovery of a ∼50 kDa protein specific for Salmonella enterica serovar Typhi. However, the structure of this antigen remains unknown till today. Studies on the structure of this antigen are important to elucidate its function, which will in turn increase the efficiency of the development and improvement of the typhoid detection test. This paper described the predictive structure and function of the antigenically specific protein. The homology modeling approach was employed to construct the three-dimensional structure of the antigen. The built structure possesses the features of TolC-like outer membrane protein. Molecular docking simulation was also performed to further probe the functionality of the antigen. Docking results showed that hexamminecobalt, Co(NH(3))(6)(3+), as an inhibitor of TolC protein, formed favorable hydrogen bonds with D368 and D371 of the antigen. The single point (D368A, D371A) and double point (D368A and D371A) mutations of the antigen showed a decrease (single point mutation) and loss (double point mutations) of binding affinity towards hexamminecobalt. The architecture features of the built model and the docking simulation reinforced and supported that this antigen is indeed the variant of outer membrane protein, TolC. As channel proteins are important for the virulence and survival of bacteria, therefore this ∼50 kDa channel protein is a good specific target for typhoid detection test.

Permeation through nanochannels: revealing fast kinetics

The permeation of water soluble molecules across cell membranes is controlled by channel-forming proteins and, in particular, the channel surface determines the selectivity. An adequate method to study the properties of these channels is electrophysiology and, in particular, analyzing the ion current fluctuation in the presence of permeating solutes. Ion current fluctuation analysis provides information on possible interactions of solutes with the channel surface. Due to the limited time resolution, fast permeation events are not visible using standard techniques. Here, we demonstrate that miniaturization of the lipid bilayer; varying the temperature or changing the solvent may enhance the resolution. Although electrophysiology is considered as a single molecule technique, it does not provide atomic resolution. Molecular details of solute permeation can be revealed by combining electrophysiology and all-atom computer modeling; these methods include ion conductance, selectivity, ion pair formation, and rate limiting interactions of the solute with the channel walls during permeation.

Structural and dynamical insights into the opening mechanism of P. aeruginosa OprM channel

Originally described in bacteria, drug transporters are now recognized as major determinants in antibiotics resistance. For Gram-negative bacteria, the reversible assembly consisting of an inner membrane protein responsible for the active transport, a periplasmic protein, and an exit outer membrane channel achieves transport. The opening of the outer membrane protein OprM from Pseudomonas aeruginosa was modeled through normal mode analysis starting from a new X-ray structure solved at 2.4 A resolution in P2(1)2(1)2(1) space group. The three monomers are not linked by internal crystallographic symmetries highlighting the possible functional differences. This structure is closed at both ends, but modeling allowed for an opening that is not reduced to the classically proposed "iris-like mechanism."

Functional rotation of the transporter AcrB: insights into drug extrusion from simulations

The tripartite complex AcrAB-TolC is the major efflux system in Escherichia coli. It extrudes a wide spectrum of noxious compounds out of the bacterium, including many antibiotics. Its active part, the homotrimeric transporter AcrB, is responsible for the selective binding of substrates and energy transduction. Based on available crystal structures and biochemical data, the transport of substrates by AcrB has been proposed to take place via a functional rotation, in which each monomer assumes a particular conformation. However, there is no molecular-level description of the conformational changes associated with the rotation and their connection to drug extrusion. To obtain insights thereon, we have performed extensive targeted molecular dynamics simulations mimicking the functional rotation of AcrB containing doxorubicin, one of the two substrates that were co-crystallized so far. The simulations, including almost half a million atoms, have been used to test several hypotheses concerning the structure-dynamics-function relationship of this transporter. Our results indicate that, upon induction of conformational changes, the substrate detaches from the binding pocket and approaches the gate to the central funnel. Furthermore, we provide strong evidence for the proposed peristaltic transport involving a zipper-like closure of the binding pocket, responsible for the displacement of the drug. A concerted opening of the channel between the binding pocket and the gate further favors the displacement of the drug. This microscopically well-funded information allows one to identify the role of specific amino acids during the transitions and to shed light on the functioning of AcrB.

In nature, bacteria have to resist several toxic threats to be able to survive, from bile acids in intestines up to antibiotics. The Escherichia coli bacterium, which usually is a commensal inhabitant of human intestines, can also acquire pathogenic properties which would harm the human body. To dispose of toxic compounds, E. coli has developed a protein machinery which is called “efflux pump”. Here, we studied the dynamics of the transporter protein AcrB, a component of the E. coli major efflux system, in complex with an antibiotic (doxorubicin). We used computer simulations to complement the existing experimental data. Our purpose was to gain more detailed insights into the pumping mechanism at the molecular level. In our simulations the drug leaves the binding pocket upon induction of functional rotation in the protein, although a complete extrusion was never observed. A peristaltic motion, which starts with a zipper-like closure of the interior of the protein, is an important step for the extrusion of the drug. Interestingly, such a peristaltic mechanism of pumping has been suggested before on the basis of structural data. The molecular details obtained in this study shall deepen the understanding of the functioning of the efflux pump.

Understanding ion conductance on a molecular level: an all-atom modeling of the bacterial porin OmpF

All-atom molecular dynamics simulations of the ion current through OmpF, the major porin in the outer membrane of Escherichia coli, were performed. Starting from the crystal structure, the all-atom modeling allows us to calculate a parameter-free ion conductance in semiquantitative agreement with experiment. Discrepancies between modeling and experiment occur, e.g., at salt concentrations above 1 M KCl or at high temperatures. At lower salt concentrations, the ions have separate pathways along the channel surface. The constriction zone in the channel contains, on one side, a series of positively charges (R42, R82, R132), and on the opposite side, two negatively charged residues (D113, E117). Mutations generated in the constriction zone by removing cationic residues enhance the otherwise small cation selectivity, whereas removing the anionic residues reverses the selectivity. Reduction of the negatively charged residues decreases the conductance by half, whereas cationic residues enhance the conductance. Experiments on mutants confirm the results of the molecular-level simulations.