Atomistic molecular-dynamics simulations enable prediction of the arginine permeation pathway through OccD1/OprD from Pseudomonas aeruginosa

Molecular Mechanism of Ciprofloxacin Translocation Through the Major Diffusion Channels of the ESKAPE Pathogens Klebsiella pneumoniae and Enterobacter cloacae

Experimental studies on the translocation and accumulation of antibiotics in Gram-negative bacteria have revealed details of the properties that allow efficient permeation through bacterial outer membrane porins. Among the major outer membrane diffusion channels, OmpF has been extensively studied to understand the antibiotic translocation process. In a few cases, this knowledge has also helped to improve the efficacy of existing antibacterial molecules. However, the extension of these strategies to enhance the efficacy of other existing and novel drugs require comprehensive molecular insight into the permeation process and an understanding of how antibiotic and channel properties influence the effective permeation rates. Previous studies have investigated how differences in antibiotic charge distribution can influence the observed permeation pathways through the OmpF channel, and have shown that the dynamics of the L3 loop can play a dominant role in the permeation process. Here, we perform all-atom simulations of the OmpF orthologs, OmpE35 from Enterobacter cloacae and OmpK35 from Klebsiella pneumoniae. Unbiased simulations of the porins and biased simulations of the ciprofloxacin permeation processes through these channels provide insight into the differences in the permeation pathway and energetics. In addition, we show that similar to the OmpF channel, antibiotic-induced dynamics of the L3 loop are also operative in the orthologs. However, the sequence and structural differences, influence the extent of the L3 loop fluctuations with OmpK35 showing greater stability in unbiased runs and subdued fluctuations in simulations with ciprofloxacin.

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.

Free energy and kinetics of cAMP permeation through connexin26 via applied voltage and milestoning

The connexin family is a diverse group of highly regulated wide-pore channels permeable to biological signaling molecules. Despite the critical roles of connexins in mediating selective molecular signaling in health and disease, the basis of molecular permeation through these pores remains unclear. Here, we report the thermodynamics and kinetics of binding and transport of a second messenger, adenosine-3',5'-cyclophosphate (cAMP), through a connexin26 hemichannel (Cx26). First, inward and outward fluxes of cAMP molecules solvated in KCl solution were obtained from 4 μs of ± 200 mV simulations. These fluxes data yielded a single-channel permeability of cAMP and cAMP/K⁺ permeability ratio consistent with experimentally measured values. The results from voltage simulations were then compared with the potential of mean force (PMF) and the mean first passage times (MFPTs) of a single cAMP without voltage, obtained from a total of 16.5 μs of Voronoi-tessellated Markovian milestoning simulations. Both the voltage simulations and the milestoning simulations revealed two cAMP-binding sites, for which the binding constants K_D and dissociation rates k_off were computed from PMF and MFPTs. The protein dipole inside the pore produces an asymmetric PMF, reflected in unequal cAMP MFPTs in each direction once within the pore. The free energy profiles under opposite voltages were derived from the milestoning PMF and revealed the interplay between voltage and channel polarity on the total free energy. In addition, we show how these factors influence the cAMP dipole vector during permeation, and how cAMP affects the local and nonlocal pore diameter in a position-dependent manner.

Development of the computational antibiotic screening platform (CLASP) to aid in the discovery of new antibiotics

Bacterial colonization of biotic and abiotic surfaces and antibiotic resistance are grand challenges with paramount societal impacts. However, in the face of increasing bacterial resistance to all known antibiotics, efforts to discover new classes of antibiotics have languished, creating an urgent need to accelerate the antibiotic discovery pipeline. A major deterrent in the discovering of new antibiotics is the limited permeability of molecules across the bacterial envelope. Notably, the Gram-negative bacteria have nutrient specific protein channels (or porins) that restrict the permeability of non-essential molecules, including antibiotics. Here, we have developed the Computational Antibiotic Screening Platform (CLASP) for screening of potential drug molecules through the porins. The CLASP takes advantage of coarse grain (CG) resolution, advanced sampling techniques, and a parallel computing environment to maximize its performance. The CLASP yields comprehensive thermodynamic and kinetic output data of a potential drug molecule within a few hours of wall-clock time. Its output includes the potential of mean force profile, energy barrier, the rate constant, and contact analysis of the molecule with the pore-lining residues, and the orientational analysis of the molecule in the porin channel. In our first CLASP application, we report the transport properties of six carbapenem antibiotics-biapenem, doripenem, ertapenem, imipenem, meropenem, and panipenem-through OccD3, a major channel for carbapenem uptake in Pseudomonas aeruginosa. The CLASP is designed to screen small molecule libraries with a fast turnaround time to yield structure-property relationships to discover antibiotics with high permeability. The CLASP will be freely distributed to enable accelerated antibiotic drug discovery.

How to Enter a Bacterium: Bacterial Porins and the Permeation of Antibiotics

Despite tremendous successes in the field of antibiotic discovery seen in the previous century, infectious diseases have remained a leading cause of death. More specifically, pathogenic Gram-negative bacteria have become a global threat due to their extraordinary ability to acquire resistance against any clinically available antibiotic, thus urging for the discovery of novel antibacterial agents. One major challenge is to design new antibiotics molecules able to rapidly penetrate Gram-negative bacteria in order to achieve a lethal intracellular drug accumulation. Protein channels in the outer membrane are known to form an entry route for many antibiotics into bacterial cells. Up until today, there has been a lack of simple experimental techniques to measure the antibiotic uptake and the local concentration in subcellular compartments. Hence, rules for translocation directly into the various Gram-negative bacteria via the outer membrane or via channels have remained elusive, hindering the design of new or the improvement of existing antibiotics. In this review, we will discuss the recent progress, both experimentally as well as computationally, in understanding the structure-function relationship of outer-membrane channels of Gram-negative pathogens, mainly focusing on the transport of antibiotics.

Millisecond-Long Simulations of Antibiotics Transport through Outer Membrane Channels

To reach their target site inside Gram-negative bacteria, almost all antibiotics need to cross the outer membrane. Computational modeling of such processes can be numerically demanding due to the size of the systems and especially due to the timescales involved. Recently, a hybrid Brownian and molecular dynamics approach, i.e., Brownian dynamics including explicit atoms (BRODEA), has been developed and evaluated for studying the transport of monoatomic ions through membrane channels. Later on, this numerically efficient scheme has been applied to determine the free energy surfaces of the ciprofloxacin and enrofloxacin translocation through the porin OmpC using temperature-accelerated simulations. To improve the usability and accuracy of the approach, schemes to approximate the position-dependent diffusion constant of the molecule while traversing the pore had to be established. To this end, we have studied the translocation of the charged phosphonic acid antibiotic fosfomycin through the porin OmpF from Escherichia coli devising and benchmarking several diffusion models. To test the efficiency and sensitivity of these models, the effect of OmpF mutations on the permeation of fosfomycin was analyzed. Permeation events have been recorded over millisecond-long biased and unbiased simulations, from which thermodynamics and kinetics quantities of the translocation processes were determined. As a result, the use of the BRODEA approach, together with the appropriate diffusion model, was seen to accurately reproduce the findings observed in electrophysiology experiments and all-atom molecular dynamics simulations. These results suggest that the BRODEA approach can become a valuable tool for screening numerous compounds to evaluate their outer membrane permeability, a property important in the development of new antibiotics.

Molecular Simulations of Gram-Negative Bacterial Membranes Come of Age

Gram-negative bacteria are protected by a multicompartmental molecular architecture known as the cell envelope that contains two membranes and a thin cell wall. As the cell envelope controls influx and efflux of molecular species, in recent years both experimental and computational studies of such architectures have seen a resurgence due to the implications for antibiotic development. In this article we review recent progress in molecular simulations of bacterial membranes. We show that enormous progress has been made in terms of the lipidic and protein compositions of bacterial systems. The simulations have moved away from the traditional setup of one protein surrounded by a large patch of the same lipid type toward a more bio-logically representative viewpoint. Simulations with multiple cell envelope components are also emerging. We review some of the key method developments that have facilitated recent progress, discuss some current limitations, and offer a perspective on future directions.

An Overview of Molecular Modeling for Drug Discovery with Specific Illustrative Examples of Applications

In this paper we review the current status of high-performance computing applications in the general area of drug discovery. We provide an introduction to the methodologies applied at atomic and molecular scales, followed by three specific examples of implementation of these tools. The first example describes in silico modeling of the adsorption of small molecules to organic and inorganic surfaces, which may be applied to drug delivery issues. The second example involves DNA translocation through nanopores with major significance to DNA sequencing efforts. The final example offers an overview of computer-aided drug design, with some illustrative examples of its usefulness.

Enrofloxacin Permeation Pathways across the Porin OmpC

In Gram-negative bacteria, the lack or quenching of antibiotic translocation across the outer membrane is one of the main factors for acquiring antibiotic resistance. An atomic-level comprehension of the key features governing the transport of drugs by outer-membrane protein channels would be very helpful in developing the next generation of antibiotics. In a previous study [ J. D. Prajapati et al. J. Chem. Theory Comput. 2017 , 13 , 4553 ], we characterized the diffusion pathway of a ciprofloxacin molecule through the outer membrane porin OmpC of Escherichia coli by combining metadynamics and a zero-temperature string method. Here, we evaluate the diffusion route through the OmpC porin for a similar fluoroquinolone, that is, the enrofloxacin molecule, using the previously developed protocol. As a result, it was found that the lowest-energy pathway was similar to that for ciprofloxacin; namely, a reorientation was required on the extracellular side with the carboxyl group ahead before enrofloxacin reached the constriction region. In turn, the free-energy basins for both antibiotics are located at similar positions in the space defined by selected reaction coordinates, and their affinity sites share a wide number of porin residues. However, there are some important deviations due to the chemical differences of these two drugs. On the one hand, a slower diffusion process is expected for enrofloxacin, as the permeation pathway exhibits higher overall energy barriers, mainly in the constriction region. On the other hand, enrofloxacin needs to replace some polar interactions in its affinity sites with nonpolar ones. This study demonstrates how minor chemical modifications can qualitatively affect the translocation mechanism of an antibiotic molecule.

Rationalizing the permeation of polar antibiotics into Gram-negative bacteria

The increasing level of antibiotic resistance in Gram-negative bacteria, together with the lack of new potential drug scaffolds in the pipeline, make the problem of infectious diseases a global challenge for modern medicine. The main reason that Gram-negative bacteria are particularly challenging is the presence of an outer cell-protecting membrane, which is not present in Gram-positive species. Such an asymmetric bilayer is a highly effective barrier for polar molecules. Several protein systems are expressed in the outer membrane to control the internal concentration of both nutrients and noxious species, in particular: (i) water-filled channels that modulate the permeation of polar molecules and ions according to concentration gradients, and (ii) efflux pumps to actively expel toxic compounds. Thus, besides expressing specific enzymes for drugs degradation, Gram-negative bacteria can also resist by modulating the influx and efflux of antibiotics, keeping the internal concentration low. However, there are no direct and robust experimental methods capable of measuring the permeability of small molecules, thus severely limiting our knowledge of the molecular mechanisms that ultimately control the permeation of antibiotics through the outer membrane. This is the innovation gap to be filled for Gram-negative bacteria. This review is focused on the permeation of small molecules through porins, considered the main path for the entry of polar antibiotics into Gram-negative bacteria. A fundamental understanding of how these proteins are able to filter small molecules is a prerequisite to design/optimize antibacterials with improved permeation. The level of sophistication of modern molecular modeling algorithms and the advances in new computer hardware has made the simulation of such complex processes possible at the molecular level. In this work we aim to share our experience and perspectives in the context of a multidisciplinary extended collaboration within the IMI-Translocation consortium. The synergistic combination of structural data, in vitro assays and computer simulations has proven to give new insights towards the identification and description of physico-chemical properties modulating permeation. Once similar general rules are identified, we believe that the use of virtual screening techniques will be very helpful in searching for new molecular scaffolds with enhanced permeation, and that molecular modeling will be of fundamental assistance to the optimization stage.

Molecular Epidemiology of Mutations in Antimicrobial Resistance Loci of Pseudomonas aeruginosa Isolates from Airways of Cystic Fibrosis Patients

The chronic airway infections with Pseudomonas aeruginosa in people with cystic fibrosis (CF) are treated with aerosolized antibiotics, oral fluoroquinolones, and/or intravenous combination therapy with aminoglycosides and β-lactam antibiotics. An international strain collection of 361 P. aeruginosa isolates from 258 CF patients seen at 30 CF clinics was examined for mutations in 17 antimicrobial susceptibility and resistance loci that had been identified as hot spots of mutation by genome sequencing of serial isolates from a single CF clinic. Combinatorial amplicon sequencing of pooled PCR products identified 1,112 sequence variants that were not present in the genomes of representative strains of the 20 most common clones of the global P. aeruginosa population. A high frequency of singular coding variants was seen in spuE, mexA, gyrA, rpoB, fusA1, mexZ, mexY, oprD, ampD, parR, parS, and envZ (amgS), reflecting the pressure upon P. aeruginosa in lungs of CF patients to generate novel protein variants. The proportion of nonneutral amino acid exchanges was high. Of the 17 loci, mexA, mexZ, and pagL were most frequently affected by independent stop mutations. Private and de novo mutations seem to play a pivotal role in the response of P. aeruginosa populations to the antimicrobial load and the individual CF host.

Molecular Simulations of Gram-Negative Bacterial Membranes: A Vignette of Some Recent Successes

In the following review we use recent examples from the literature to discuss progress in the area of atomistic and coarse-grained molecular dynamics simulations of selected bacterial membranes and proteins, with a particular focus on Gram-negative bacteria. As structural biology continues to provide increasingly high-resolution data on the proteins that reside within these membranes, simulations have an important role to play in linking these data with the dynamical behavior and function of these proteins. In particular, in the last few years there has been significant progress in addressing the issue of biochemical complexity of bacterial membranes such that the heterogeneity of the lipid and protein components of these membranes are now being incorporated into molecular-level models. Thus, in future we can look forward to complementary data from structural biology and molecular simulations combining to provide key details of structure-dynamics-function relationships in bacterial membranes.

Identification of unique in-frame deletions in OprD among clinical isolates of Pseudomonas aeruginosa

A large percentage of Pseudomonas aeruginosa clinical isolates have been noted to be resistant to carbapenems due to loss of function of the OprD porin, the primary mechanism of entry for carbapenems. Such modifications also substantially abolish the organism's ability to transport arginine. Here we report the identification of an in-frame deletion in oprD which confers carbapenem resistance but is expressed and retains the ability to transport arginine.

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.

Selectivity of pyoverdine recognition by the FpvA receptor of Pseudomonas aeruginosa from molecular dynamics simulations

The Gram-negative bacterium Pseudomonas aeruginosa, a ubiquitous human opportunistic pathogen, has developed resistances to multiple antibiotics. It uses its primary native siderophore, pyoverdine, to scavenge the iron essential to its growth in the outside medium and transport it back into its cytoplasm. The FpvA receptor on the bacterial outer membrane recognizes and internalizes pyoverdine bearing its iron payload, but can also bind pyoverdines from other Pseudomonads or synthetic analogues. Pyoverdine derivatives could therefore be used as vectors to deliver antibiotics into the bacterium. In this study, we use molecular dynamics and free energy calculations to characterize the mechanisms and thermodynamics of the recognition of the native pyoverdines of P. aeruginosa and P. fluorescens by FpvA. Based on these results, we delineate the features that pyoverdines with high affinity for FpvA should possess. In particular, we show that (i) the dynamics and interaction of the unbound pyoverdines with water should be optimized with equal care as the interface contacts in the complex with FpvA; (ii) the C-terminal extremity of the pyoverdine chain, which appears to play no role in the bound complex, is involved in the intermediate stages of recognition; and (iii) the length and cyclicity of the pyoverdine chain can be used to fine-tune the kinetics of the recognition mechanism.

The Nucleotide Capture Region of Alpha Hemolysin: Insights into Nanopore Design for DNA Sequencing from Molecular Dynamics Simulations

Nanopore technology for DNA sequencing is constantly being refined and improved. In strand sequencing a single strand of DNA is fed through a nanopore and subsequent fluctuations in the current are measured. A major hurdle is that the DNA is translocated through the pore at a rate that is too fast for the current measurement systems. An alternative approach is “exonuclease sequencing”, in which an exonuclease is attached to the nanopore that is able to process the strand, cleaving off one base at a time. The bases then flow through the nanopore and the current is measured. This method has the advantage of potentially solving the translocation rate problem, as the speed is controlled by the exonuclease. Here we consider the practical details of exonuclease attachment to the protein alpha hemolysin. We employ molecular dynamics simulations to determine the ideal (a) distance from alpha-hemolysin, and (b) the orientation of the monophosphate nucleotides upon release from the exonuclease such that they will enter the protein. Our results indicate an almost linear decrease in the probability of entry into the protein with increasing distance of nucleotide release. The nucleotide orientation is less significant for entry into the protein.

Theoretical analysis of ion conductance and gating transitions in the OpdK (OccK1) channel

Molecular simulations have been performed on the pore OpdK elucidating molecular details of ion conductance and a possible gating mechanism.

Molecular basis of substrate translocation through the outer membrane channel OprD of Pseudomonas aeruginosa

The dynamics and interplay of internal and external loops create two alternative paths for the permeation of substrates through the specific outer membrane channel OprD.