Microscopic Mechanism of Antibiotics Translocation through a Porin

OmpF from the outer membrane of Escherichia coli is a general porin considered to be the main pathway for beta-lactam antibiotics. The availability of a high-resolution crystal structure of OmpF and new experimental techniques at the single-molecule level have opened the way to the investigation of the microscopic mechanisms that allow the passage of antibiotics through bacterial pores. We applied molecular dynamics simulations to investigate the translocation process of ampicillin (Amp) through OmpF. Using a recent algorithm capable of accelerating molecular dynamics simulations we have been able to obtain a reaction path for the translocation of Amp through OmpF. The mechanism of passage depends both on the internal degrees of freedom of Amp and on interactions of Amp with OmpF. Understanding this mechanism would help us design more efficient antibiotics and shed light on nature's way of devising channels able to enhance the transport of molecules through membranes.     

Interaction of zwitterionic penicillins with the OmpF channel facilitates their translocation

To study translocation of beta-lactam antibiotics of different size and charge across the outer bacterial membrane, we combine an analysis of ion currents through single trimeric outer membrane protein F (OmpF) porins in planar lipid bilayers with molecular dynamics simulations. Because the size of penicillin molecules is close to the size of the narrowest part of the OmpF pore, penicillins occlude the pore during their translocation. Favorably interacting penicillins cause time-resolvable transient blockages of the small-ion current through the channel and thereby provide information about their dynamics within the pore. Analyzing these random fluctuations, we find that ampicillin and amoxicillin have a relatively high affinity for OmpF. In contrast, no or only a weak interaction is detected for carbenicillin, azlocillin, and piperacillin. Molecular dynamics simulations suggest a possible pathway of these drugs through the OmpF channel and rationalize our experimental findings. For zwitterionic ampicillin and amoxicillin, we identify a region of binding sites near the narrowest part of the channel pore. Interactions with these sites partially compensate for the entropic cost of drug confinement by the channel. Whereas azlocillin and piperacillin are clearly too big to pass through the channel constriction, dianionic carbenicillin does not find an efficient binding region in the constriction zone. Carbenicillin's favorable interactions are limited to the extracellular vestibule. These observations confirm our earlier suggestion that a set of high-affinity sites at the narrowest part of the OmpF channel improves a drug's ability to cross the membrane via the pore.     

In-silico interaction studies suggest RND efflux pump mediates polymyxin resistance in Acinetobacter baumannii

Bacterial efflux pumps have emerged as antibiotic resistance determinants and confers multi-drug resistance to a broad range of antimicrobials as well as non-antibiotic substances. A study about translocation of antibiotic molecules through the efflux transporter, will contribute in determining substrate specificity. In the present study, we have explored RND family efflux pump extensively found in Acinetobacter baumannii i.e. AdeABC. Besides, another well studied RND efflux pump, AcrAB-TolC together with a non-RND efflux pump, NorM was investigated for comparative analysis. We employed a series of computational techniques ranging from molecular docking to binding free energy estimation and molecular dynamics simulations to determine the binding affinity for different classes of drugs, namely aminoglycosides, polymyxins, β-lactams, tetracyclines, glycylcyclines, quinolones and metronidazole with AdeB, AcrB, and NorM efflux proteins. Our results revealed that class polymyxins has the highest binding affinity with the RND efflux pumps i.e. AcrAB-TolC and AdeABC as well as non-RND efflux pump, NorM. The experimental validation study demonstrated bigger zone of inhibition in presence of efflux pump inhibitor than polymyxin alone thus unveiling its specificity toward efflux pump. The reported experimental data comprising of minimum inhibitory concentration of antibiotics toward these efflux pumps also support our finding based on in silico approach. To recapitulate the outcome, polymyxins shows maximum specificity toward RND as well as non-RND efflux pump and may unlatch the way to rationally develop new potential antibacterial agents as well as efflux pump inhibitors in order to combat resistance.     

Fosfomycin Permeation through the Outer Membrane Porin OmpF

Fosfomycin is a frequently prescribed drug in the treatment of acute urinary tract infections. It enters the bacterial cytoplasm and inhibits the biosynthesis of peptidoglycans by targeting the MurA enzyme. Despite extensive pharmacological studies and clinical use, the permeability of fosfomycin across the bacterial outer membrane is largely unexplored. Here, we investigate the fosfomycin permeability across the outer membrane of Gram-negative bacteria by electrophysiology experiments as well as by all-atom molecular dynamics simulations including free-energy and applied-field techniques. Notably, in an electrophysiological zero-current assay as well as in the molecular simulations, we found that fosfomycin can rapidly permeate the abundant Escherichia coli porin OmpF. Furthermore, two triple mutants in the constriction region of the porin have been investigated. The permeation rates through these mutants are slightly lower than that of the wild type but fosfomycin can still permeate. Altogether, this work unravels molecular details of fosfomycin permeation through the outer membrane porin OmpF of E. coli and moreover provides hints for understanding the translocation of phosphonic acid antibiotics through other outer membrane pores.     

Kibdelomycin A,a congener of kibdelomycin,derivatives and their antibacterial activities

Emergence of bacterial resistance has eroded the effectiveness of many life saving antibiotics leading to an urgent need for new chemical classes of antibacterial agents. We have applied a Staphylococcus aureus fitness test strategy to natural products screening to meet this challenge. In this paper we report the discovery of kibdelomycin A, a demethylated congener of kibdelomycin, the representative of a novel class of antibiotics produced by a new strain of Kibdelosporangium. Kibdelomycin A is a potent inhibitor of DNA gyrase and topoisomerase IV, inhibits DNA synthesis and shows whole cell antibiotic activity, albeit, less potently than kibdelomycin. Kibdelomycin C-33 acetate and tetrahydro-bisdechloro derivatives of kibdelomycin were prepared which helped define a basic SAR of the family.     

Experimental and simulation studies reveal mechanism of action of human defensin derivatives

Antimicrobial peptides (AMPs) are naturally occurring promising candidates which can be used as antibiotics against a wide variety of bacteria. The key component for using them as a potent antibiotic is that their mechanism of action is less prone to bacterial resistance. However, the molecular details of their mechanism of action is not yet fully understood. In this study, we try to shed light on the mode of action of AMPs, possible reason behind it, and their interaction with lipid bilayers through experimental as well as molecular dynamics (MD) simulation studies. The focal of our study was Human beta defensin 3 (hBD-3) which is a naturally occurring AMP. We chose three derivatives of hBD-3, namely CHRG01, KSR, and KLR for the detailed analysis presented in this study. These three peptides are evaluated for their antibacterial potency, secondary structure analysis and mechanism of action. The experimental results reveal that these peptides are active against gram positive as well as gram negative bacteria and kill bacteria by forming membrane pores. The MD simulation results correlate well with the antibacterial activity and shed light into the early membrane insertion dynamics. Moreover, the specific amino acids responsible for membrane disruptions are also identified from the MD simulations. Understanding the molecular level interaction of individual amino acids with the lipid bilayer will greatly help in the design of more efficient antimicrobial peptides.     

Event detection and sub‐state discovery from biomolecular simulations using higher‐order statistics: Application to enzyme adenylate kinase

Biomolecular simulations at millisecond and longer time‐scales can provide vital insights into functional mechanisms. Because post‐simulation analyses of such large trajectory datasets can be a limiting factor in obtaining biological insights, there is an emerging need to identify key dynamical events and relating these events to the biological function online, that is, as simulations are progressing. Recently, we have introduced a novel computational technique, quasi‐anharmonic analysis (QAA) (Ramanathan et al., PLoS One 2011;6:e15827), for partitioning the conformational landscape into a hierarchy of functionally relevant sub‐states. The unique capabilities of QAA are enabled by exploiting anharmonicity in the form of fourth‐order statistics for characterizing atomic fluctuations. In this article, we extend QAA for analyzing long time‐scale simulations online. In particular, we present HOST4MD—a higher‐order statistical toolbox for molecular dynamics simulations, which (1) identifies key dynamical events as simulations are in progress, (2) explores potential sub‐states, and (3) identifies conformational transitions that enable the protein to access those sub‐states. We demonstrate HOST4MD on microsecond timescale simulations of the enzyme adenylate kinase in its apo state. HOST4MD identifies several conformational events in these simulations, revealing how the intrinsic coupling between the three subdomains (LID, CORE, and NMP) changes during the simulations. Further, it also identifies an inherent asymmetry in the opening/closing of the two binding sites. We anticipate that HOST4MD will provide a powerful and extensible framework for detecting biophysically relevant conformational coordinates from long time‐scale simulations. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

An inner membrane protein N-terminal signal sequence is able to promote efficient localisation of an outer membrane protein in Escherichia coli

To test the importance of N-terminal pre-sequences in translocation of different classes of membrane proteins, we exchanged the normal signal sequence of an Escherichia coli outer membrane protein, OmpF, for the pre-sequence of the inner membrane protein, DacA. The DacA-OmpF hybrid was efficiently assembled into the outer membrane in a functionally active form. Thus the pre-sequence of DacA, despite its relatively low hydrophobicity compared with that of OmpF, contains all the essential information necessary to initiate the translocation of OmpF to the outer membrane. Since processing of DacA was also shown to be dependent upon SecA we conclude that the initiation of translocation of this inner membrane polypeptide across the envelope occurs by the same mechanism as outer membrane and periplasmic proteins. The N-terminal 11 amino acids of mature OmpF, which in the hybrid are replaced by the N-terminal nine amino acids of DacA, carry no essential assembly signals since the hybrid protein is apparently assembled with equal efficiency to OmpF.     

The Binding of Antibiotics in OmpF Porin

1. Download : Download high-res image (445KB)
2. Download : Download full-size image

Highlights? OmpF porin binds antibiotics in the extracellular and periplasmic pore vestibules ? MD simulations show that ionic current through OmpF is blocked by bound ampicillin ? Carbenicillin binding alters ion selectivity of OmpF but not total ionic current ? Disruption of binding increases the susceptibility of E. coli for antibiotics     

Membrane active antimicrobial activity and molecular dynamics study of a novel cationic antimicrobial peptide polybia-MPI,from the venom of Polybia paulista

As the frequent emergence of the resistant bacteria, the development of new agents with a new action mode attracts a great deal of interest. It is now widely accepted that antimicrobial peptides (AMPs) are promising alternatives to conventional antibiotics. In this study, antimicrobial peptide polybia-MPI and its analogs were synthesized and their antibacterial activity was studied. Our results revealed that polybia-MPI has potent antibacterial activity against both Gram-positive and Gram-negative bacteria. Its ability to make PI permeate into bacteria and lead to the leakage of calcein from model membrane LUVs, suggests a killing mechanism involving membrane perturbation. SEM and TEM microscopy experiments verified that the morphology of bacteria was changed greatly under the treatment of polybia-MPI. Compared with the conventional chemotherapy, polybia-MPI targets the cell membrane rather than entering into the cell to exert its antibacterial activity. Furthermore, molecular dynamics (MD) simulations were employed to investigate the mechanism of membrane perturbation. The results indicated that the α-helical conformation in the membrane is required for the exhibition of antibacterial activity and the membrane disturbance by polybia-MPI is a cooperative process. In conclusion, with the increasing resistance to conventional antibiotics, there is no doubt that polybia-MPI could offer a new strategy to defend the resistant bacteria.     

An assessment of optimal time scale of conformational resampling for parallel cascade selection molecular dynamics

Parallel cascade selection molecular dynamics (PaCS-MD) has been proposed as a conformational sampling method for enhancing structural transitions from a given reactant to a product by repeating cycles of short-time MD simulations. In the present paper, we assessed how the time scale of a short-time MD simulation affected the computational efficiency by changing the simulation length. In conclusion, ps-order (t_ps) PaCS-MD simulations showed a higher computational efficiency as a total simulation time over the cycles than ns-order (t_ns) PaCS-MD simulations, indicating that t_ps might be suitable for generating structural transitions efficiently.     

Molecular dynamics simulations of membrane proteins

Membrane proteins control the traffic across cell membranes and thereby play an essential role in cell function from transport of various solutes to immune response via molecular recognition. Because it is very difficult to determine the structures of membrane proteins experimentally, computational methods have been increasingly used to study their structure and function. Here we focus on two classes of membrane proteins—ion channels and transporters—which are responsible for the generation of action potentials in nerves, muscles, and other excitable cells. We describe how computational methods have been used to construct models for these proteins and to study the transport mechanism. The main computational tool is the molecular dynamics (MD) simulation, which can be used for everything from refinement of protein structures to free energy calculations of transport processes. We illustrate with specific examples from gramicidin and potassium channels and aspartate transporters how the function of these membrane proteins can be investigated using MD simulations.     

Molecular dynamic simulations of the metallo-beta-lactamase from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor

The beta-lactam-based antibiotics are among the most prescribed and effective antibacterial agents. Widespread use of these antibiotics, however, has created tremendous pressure for the emergence of resistance mechanisms in bacteria. The most common cause of antibiotic resistance is bacterial production of actamases that efficiently degrade antibiotics. The metallo-beta-lactamases are of particular clinical concern due to their transference between bacterial strains. We used molecular dynamics (MD) simulations to further study the conformational changes that occur due to binding of an inhibitor to the dicanzinc metallo-beta-lactamase from Bacteroides fragilis. Our studies confirm previous findings that the major flap is a major source of plasticity within the active site, therefore its dynamic response should be considered in drug development. However, our results also suggest the need for care in using MD simulations in evaluating loop mobility, both due to relaxation times and to the need to accurately model the zinc active site. Our study also reveals two new robust responses to ligand binding. First, there are specific localized changes in the zinc active site—a local loop flip—due to ligand intercalation that may be critical to the function of this enzyme. Second, inhibitor binding perturbs the dynamics throughout the protein, without otherwise perturbing the enzyme structure. These dynamic perturbations radiate outward from the active site and their existence suggests that long-range communication and dynamics may be important in the activity of this enzyme. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Translocation and interactions of L-arabinose in OmpF porin: A molecular dynamics study

The passage of a natural substrate, L-arabinose (L-ARA) through Escherichia coli porin embedded in an artificial bilayer, is studied by equilibrium molecular dynamics simulations. We investigate the early stage of translocation process of L-ARA from intra-cellular to extra-cellular side (Int-to-Ext) across the bilayer. The average trajectory path over all L-ARA molecules along with quantum-mechanical configuration-optimizations at PM3 level predict the existence of at least three trapping zones. The common feature within all these zones is that L-ARA remains perpendicular to the channel axis. It is remarkable how the orientation and translational-rotational motion of L-ARA molecule play a role in its transport through OmpF channel. These simulations are important for better understanding of permeation process in OmpF channel. They also provide an insight into the chiral recognition of translocation process in protein nanochannels from substrate and protein prospects and help interpret experiments on permeation process of small dipolar molecules across biological membranes.     

Asymmetric pore occupancy in crystal structure of OmpF porin from Salmonella typhi

OmpF is a major general diffusion porin of Salmonella typhi, a Gram-negative bacterium, which is an obligatory human pathogen causing typhoid. The structure of S. typhi Ty21a OmpF (PDB Id: 3NSG) determined at 2.8 ? resolution by X-ray crystallography shows a 16-stranded β-barrel with three β-barrel monomers associated to form a trimer. The packing observed in S. typhi Ty21a rfOmpF crystals has not been observed earlier in other porin structures. The variations seen in the loop regions provide a starting point for using the S. typhi OmpF for structure-based multi-valent vaccine design. Along one side of the S. typhi Ty21a OmpF pore there exists a staircase arrangement of basic residues (20R, 60R, 62K, 65R, 77R, 130R and 16K), which also contribute, to the electrostatic potential in the pore. This structure suggests the presence of asymmetric electrostatics in the porin oligomer. Moreover, antibiotic translocation, permeability and reduced uptake in the case of mutants can be understood based on the structure paving the way for designing new antibiotics.     

Non-Traditional Antibacterial Screening Approaches for the Identification of Novel Inhibitors of the Glyoxylate Shunt in Gram-Negative Pathogens

Antibacterial compounds that affect bacterial viability have traditionally been identified, confirmed, and characterized in standard laboratory media. The historical success of identifying new antibiotics via this route has justifiably established a traditional means of screening for new antimicrobials. The emergence of multi-drug-resistant (MDR) bacterial pathogens has expedited the need for new antibiotics, though many in the industry have questioned the source(s) of these new compounds. As many pharmaceutical companies'' chemical libraries have been exhaustively screened via the traditional route, we have concluded that all compounds with any antibacterial potential have been identified. While new compound libraries and platforms are being pursued, it also seems prudent to screen the libraries we currently have in hand using alternative screening approaches. One strategy involves screening under conditions that better reflect the environment pathogens experience during an infection, and identifying in vivo essential targets and pathways that are dispensable for growth in standard laboratory media in vitro. Here we describe a novel screening strategy for identifying compounds that inhibit the glyoxylate shunt in Pseudomonas aeruginosa, a pathway that is required for bacterial survival in the pulmonary environment. We demonstrate that these compounds, which were not previously identified using traditional screening approaches, have broad-spectrum antibacterial activity when they are tested under in vivo-relevant conditions. We also show that these compounds have potent activity on both enzymes that comprise the glyoxylate shunt, a feature that was supported by computational homology modeling. By dual-targeting both enzymes in this pathway, we would expect to see a reduced propensity for resistance development to these compounds. Taken together, these data suggest that understanding the in vivo environment that bacterial pathogens must tolerate, and adjusting the antibacterial screening paradigm to reflect those conditions, could identify novel antibiotics for the treatment of serious MDR pathogens.     

Exploring the permeation of fluoroquinolone metalloantibiotics across outer membrane porins by combining molecular dynamics simulations and a porin-mimetic in vitro model

The misuse and overuse of fluoroquinolones in recent years have triggered alarming levels of resistance to these antibiotics. Porin channels are crucial for the permeation of fluoroquinolones across the outer membrane of Gram-negative bacteria and modifications in porin expression are an important mechanism of bacterial resistance. One possible strategy to overcome this problem is the development of ternary copper complexes with fluoroquinolones. Compared to fluoroquinolones, these metalloantibiotics present a larger partition to the lipid bilayer and a more favorable permeation, by passive diffusion, across bacteriomimetic phospholipid-based model membranes. To rule out the porin-dependent pathway for the metalloantibiotics, we explored the permeation through OmpF (one of the most abundant porins present in the outer membrane of Gram-negative bacteria) using a multi-component approach. X-ray studies of OmpF porin crystals soaked with a ciprofloxacin ternary copper complex did not show a well-defined binding site for the compound. Molecular dynamics simulations showed that the translocation of the metalloantibiotic through this porin is less favorable than that of free fluoroquinolone, as it presented a much larger free energy barrier to cross the narrow constriction region of the pore. Lastly, permeability studies of different fluoroquinolones and their respective copper complexes using a porin-mimetic in vitro model corroborated the lower rate of permeation for the metalloantibiotics relative to the free antibiotics. Our results support a porin-independent mechanism for the influx of the metalloantibiotics into the bacterial cell. This finding brings additional support to the potential application of these metalloantibiotics in the fight against resistant infections and as an alternative to fluoroquinolones.