Magnesium ions mitigate metastable states in the regulatory landscape of mRNA elements.

Residing in the 5' untranslated region of the mRNA, the 2'-deoxyguanosine (2'-dG) riboswitch mRNA element adopts an alternative structure upon binding of the 2'-dG molecule, which terminates transcription. RNA conformations are generally strongly affected by positively charged metal ions (especially Mg2+). We have quantitatively explored the combined effect of ligand (2'-dG) and Mg2+ binding on the energy landscape of the aptamer domain of the 2'-dG riboswitch with both explicit solvent all-atom molecular dynamics simulations (99 µsec aggregate sampling for the study) and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) experiments. We show that both ligand and Mg2+ are required for the stabilization of the aptamer domain; however, the two factors act with different modalities. The addition of Mg2+ remodels the energy landscape and reduces its frustration by the formation of additional contacts. In contrast, the binding of 2'-dG eliminates the metastable states by nucleating a compact core for the aptamer domain. Mg2+ ions and ligand binding are required to stabilize the least stable helix, P1 (which needs to unfold to activate the transcription platform), and the riboswitch core formed by the backbone of the P2 and P3 helices. Mg2+ and ligand also facilitate a more compact structure in the three-way junction region.

Conformational Dynamics of Loop L3 in OmpF: Implications toward Antibiotic Translocation and Voltage Gating.

In the present work, we delineate the molecular mechanism of a bulky antibiotic permeating through a bacterial channel and uncover the role of conformational dynamics of the constriction loop in this process. Using the temperature accelerated sliced sampling approach, we shed light onto the dynamics of the L3 loop, in particular the F118 to S125 segment, at the constriction regions of the OmpF porin. We complement the findings with single channel electrophysiology experiments and applied-field simulations, and we demonstrate the role of hydrogen-bond stabilization in the conformational dynamics of the L3 loop. A molecular mechanism of permeation is put forward wherein charged antibiotics perturb the network of stabilizing hydrogen-bond interactions and induce conformational changes in the L3 segment, thereby aiding the accommodation and permeation of bulky antibiotic molecules across the constriction region. We complement the findings with single channel electrophysiology experiments and demonstrate the importance of the hydrogen-bond stabilization in the conformational dynamics of the L3 loop. The generality of the present observations and experimental results regarding the L3 dynamics enables us to identify this L3 segment as the source of gating. We propose a mechanism of OmpF gating that is in agreement with previous experimental data that showed the noninfluence of cysteine double mutants that tethered the L3 tip to the barrel wall on the OmpF gating behavior. The presence of similar loop stabilization networks in porins of other clinically relevant pathogens suggests that the conformational dynamics of the constriction loop is possibly of general importance in the context of antibiotic permeation through porins.

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.

Identification of Single Amino Acid Chiral and Positional Isomers Using an Electrostatically Asymmetric Nanopore.

Chirality is essential in nearly all biological organizations and chemical reactions but is rarely considered due to technical limitations in identifying L/D isomerization. Using OmpF, a membrane channel from Escherichia coli with an electrostatically asymmetric constriction zone, allows discriminating chiral amino acids in a single peptide. The heterogeneous distribution of charged residues in OmpF causes a strong lateral electrostatic field at the constriction. This laterally asymmetric constriction zone forces the sidechains of the peptides to specific orientations within OmpF, causing distinct ionic current fluctuations. Using statistical analysis of the respective ionic current variations allows distinguishing the presence and position of a single amino acid with different chiralities. To explore potential applications, the disease-related peptide ß-Amyloid and its d-Asp1 isoform and a mixture of the icatibant peptide drug (HOE 140) and its d-Ser7 mutant have been discriminated. Both chiral isomers were not applicable to be distinguished by mass spectroscopy approaches. These findings highlight a novel sensing mechanism for identifying single amino acids in single peptides and even for achieving single-molecule protein sequencing.

Large-Peptide Permeation Through a Membrane Channel: Understanding Protamine Translocation Through CymA from Klebsiella Oxytoca*.

Quantifying the passage of the large peptide protamine (Ptm) across CymA, a passive channel for cyclodextrin uptake, is in the focus of this study. Using a reporter-pair-based fluorescence membrane assay we detected the entry of Ptm into liposomes containing CymA. The kinetics of the Ptm entry was independent of its concentration suggesting that the permeation through CymA is the rate-limiting factor. Furthermore, we reconstituted single CymA channels into planar lipid bilayers and recorded the ion current fluctuations in the presence of Ptm. To this end, we were able to resolve the voltage-dependent entry of single Ptm peptide molecules into the channel. Extrapolation to zero voltage revealed about 1-2 events per second and long dwell times, in agreement with the liposome study. Applied-field and steered molecular dynamics simulations added an atomistic view of the permeation events. It can be concluded that a concentration gradient of 1 µm Ptm leads to a translocation rate of about one molecule per second and per channel.

Computational Modeling of Ion Transport in Bulk and through a Nanopore Using the Drude Polarizable Force Field.

In the past two decades, molecular dynamics simulations have become the method of choice for elucidating the transport mechanisms of ions through various membrane channels. Often, these simulations heavily rely on classical nonpolarizable force fields (FFs), which lack electronic polarizability in the treatment of the electrostatics. The recent advancements in the Drude polarizable FF lead to a complete set of parameters for water, ions, protein, and lipids, allowing for a more realistic modeling of membrane proteins. However, the quality of these Drude FFs remains untested for such systems. Here, we examine the quality of this FF set in two ways, i.e., (i) in simple ionic aqueous solution simulations and (ii) in more complex membrane channel simulations. First, the aqueous solutions of KCl, NaCl, MgCl2, and CaCl2 salts are simulated using the polarizable Drude and the nonpolarizable CHARMM36 FFs. The bulk conductivity has been estimated for both FF sets using applied-field simulations for several concentrations and temperatures in the case of all investigated salts and compared to experimental findings. An excellent improvement in the ability of the Drude FF to reproduce the experimental bulk conductivities for KCl, NaCl, and MgCl2 solutions can be observed but not in the case of CaCl2. Moreover, the outer membrane channel OmpC from the bacterium Escherichia coli has been employed to examine the ability of the polarizable and nonpolarizable FFs to reproduce ion transport-related quantities known from experiment. Unbiased and applied-field simulations have been performed in the presence of KCl using both FF sets. Unlike for the bulk systems of aqueous salt solutions, it has been found that the Drude FF is not accurate in modeling KCl transport properties across the OmpC porin.

Role of Electroosmosis in the Permeation of Neutral Molecules: CymA and Cyclodextrin as an Example.

To quantify the flow of small uncharged molecules into and across nanopores, one often uses ion currents. The respective ion-current fluctuations caused by the presence of the analyte make it possible to draw some conclusions about the direction and magnitude of the analyte flow. However, often this flow appears to be asymmetric with respect to the applied voltage. As a possible reason for this asymmetry, we identified the electroosmotic flow (EOF), which is the water transport associated with ions driven by the external transmembrane voltage. As an example, we quantify the contribution of the EOF through a nanopore by investigating the permeation of α-cyclodextrin through CymA, a cyclodextrin-specific channel from Klebsiella oxytoca. To understand the results from electrophysiology on a molecular level, all-atom molecular dynamics simulations are used to detail the effect of the EOF on substrate entry to and exit from a CymA channel in which the N-terminus has been deleted. The combined experimental and computational results strongly suggest that one needs to account for the significant contribution of the EOF when analyzing the penetration of cyclodextrins through the CymA pore. This example study at the same time points to the more general finding that the EOF needs to be considered in translocation studies of neutral molecules and, at least in many cases, should be able to help in discriminating between translocation and binding events.

Atomistic Simulation of Molecules Interacting with Biological Nanopores: From Current Understanding to Future Directions.

Biological nanopores have been at the focus of numerous studies due to their role in many biological processes as well as their (prospective) technological applications. Among many other topics, recent studies on nanopores have addressed two key areas: antibiotic permeation through bacterial channels and sensing of analytes. Although the two areas are quite far apart in terms of their objectives, in both cases atomistic simulations attempt to understand the solute dynamics and the solute-protein interactions within the channel lumen. While decades of studies on various channels have culminated in an improved understanding of the key molecular factors and led to practical applications in some cases, successful utilization is limited. In this Perspective we summarize recent progress in understanding key issues in molecular simulations of antibiotic translocation and in the development of nanopore sensors. Moreover, we comment on possible advancements in computational algorithms that can potentially resolve some of the issues.

Exploring the Energy Landscape of Riboswitches Using Collective Variables Based on Tertiary Contacts.

Messenger RNA regulatory elements, such as riboswitches, can display a high degree of flexibility. By characterizing their energy landscapes, and corresponding distributions of 3D configurations, structure-function relationships can be elucidated. Molecular dynamics simulation with enhanced sampling is an important strategy used to computationally access free energy landscapes characterizing the accessible 3D conformations of RNAs. While tertiary contacts are thought to play important roles in RNA dynamics, it is difficult, in explicit solvent, to sample the formation and breakage of tertiary contacts, such as helix-helix interactions, pseudoknot interactions, and junction interactions, while maintaining intact secondary structure elements. To this end, we extend previously developed collective variables and metadynamics efforts, to establish a simple metadynamics protocol, which utilizes only one collective variable, based on multiple tertiary contacts, to characterize the underlying free energy landscape of any RNA molecule. We develop a modified collective variable, the tertiary contacts distance (QTC), which can probe the formation and breakage of all or selectively chosen tertiary contacts of the RNA. The SAM-I riboswitch in the presence of three ionic and substrate conditions was investigated and validated against the structure ensemble previously generated using SAXS experiments. This efficient and easy to implement all-atom MD simulation based approach incorporating metadynamics to study RNA conformational dynamics can also be transferred to any other type of biomolecule.

Changes in Salt Concentration Modify the Translocation of Neutral Molecules through a ΔCymA Nanopore in a Non-monotonic Manner.

The voltage-dependent transport through biological and artificial nanopores is being used in many applications such as DNA or protein sequencing and sensing. The primary approach to determine the transport has been to measure the temporal ion current fluctuations caused by solutes when applying external voltages. Crossing the nanoscale confinement in the presence of an applied electric field primarily relies on two factors, i.e., the electrophoretic drag and electroosmosis. The electroosmotic flow (EOF) is a voltage-dependent ion-associated flow of solvent molecules, i.e., usually water, and depends on many factors, such as pH, temperature, pore diameter, and also the concentration of ions. The exact interplay between these factors is so far poorly understood. In this joint experimental and computational study, we have investigated the dependence of the EOF on the concentration of the buffer salt by probing the transport of α-cyclodextrin molecules through the ΔCymA channel. For five different KCl concentrations in the range between 0.125 and 2 M, we performed applied-field molecular dynamics simulations and analyzed the ionic flow and the EOF across the ΔCymA pore. To our surprise, the concentration-dependent net ionic flux changes non-monotonically and nonlinearly and the EOF is seen to follow the same pattern. On the basis of these findings, we were able to correlate the concentration-dependent EOF with experimental kinetic constants for the translocation of α-cyclodextrin through the ΔCymA nanopore. Overall, the results further improve our understanding of the EOF-mediated transport through nanopores and show that the EOF needs to seriously be taken into consideration when analyzing the permeation of (neutral) substrates through nanopores.

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.

Improved Sampling and Free Energy Estimates for Antibiotic Permeation through Bacterial Porins.

Antibiotics enter into bacterial cells via protein channels that serve as low-energy pathways through the outer membrane, which is otherwise impenetrable. Insights into the molecular mechanisms underlying the transport processes are vital for the development of effective antibacterials. A much-desired prerequisite is an accurate and reproducible determination of free energy surfaces for antibiotic translocation, enabling quantitative and meaningful comparisons of permeation mechanisms for different classes of antibiotics. Inefficient sampling along the orthogonal degrees of freedom, for example, in umbrella sampling and metadynamics approaches, is however a key limitation affecting the accuracy and the convergence of free energy estimates. To overcome this limitation, two sampling methods have been employed in the present study that, respectively, combine umbrella sampling and metadynamics-style biasing schemes with temperature acceleration for improved sampling along orthogonal degrees of freedom. As a model for the transport of bulky solutes, the ciprofloxacin-OmpF system has been selected. The well-tempered metadynamics approach with multiple walkers is compared with its "temperature-accelerated" variant in terms of improvements in sampling and convergence of free energy estimates. We find that the inclusion of collective variables governing solute degrees of freedom and solute-water interactions within the sampling scheme largely alleviates sampling issues. Concerning improved sampling and convergence of free energy estimates from independent simulations, the temperature-accelerated sliced sampling approach that combines umbrella sampling with temperature-accelerated molecular dynamics performs even better as shown for the ciprofloxacin-OmpF system.

Voltage-Dependent Transport of Neutral Solutes through Nanopores: A Molecular View.

The permeation of (neutral) molecules through nanopores in the presence of external voltages depends on several factors including pore electrostatics, electrophoretic force, and electro-osmotic drag. In earlier single-channel electrophysiology experiments, voltage-dependent asymmetric transport of neutral α-cyclodextrin (α-CD) molecules through the biological nanopore ΔCymA was observed. The voltage-dependent ion-associated flow of water, the so-called electro-osmotic flow, has been suggested to be the key factor behind the observed asymmetric behavior. The influence of pore electrostatics and electrophoretic force and their interplay with the electro-osmotic drag with varying buffers and voltages has not yet been analyzed at the molecular level. Hence, the detailed physical mechanism behind this intriguing permeation process is in part still unclear. In the present study, we have performed 36 µs all-atom free energy calculations by combining applied-field molecular dynamics simulations with metadynamics techniques. The influence of several ionic conditions as well as external voltages on the permeation of α-CD molecules across the ΔCymA pore has been investigated. To decipher the thermodynamic and kinetic details, the lowest energy paths and the permeation times for α-CD translocation have been estimated. In the presence of KCl or MgCl2 salts, the charge of the cations is found to control the direction and magnitude of the electro-osmotic flow, which in turn strongly affects α-CD permeation. Overall, the present findings significantly improve the fundamental understanding of the voltage-dependent transport of neutral solutes across nanopores.

Dynamic interaction of fluoroquinolones with magnesium ions monitored using bacterial outer membrane nanopores.

Divalent ions are known to have a severe effect on the translocation of several antibiotic molecules into (pathogenic) bacteria. In the present study we have investigated the effect of divalent ions on the permeability of norfloxacin across the major outer membrane channels from E. coli (OmpF and OmpC) and E. aerogenes (Omp35 and Omp36) at the single channel level. To understand the rate limiting steps in permeation, we reconstituted single porins into planar lipid bilayers and analyzed the ion current fluctuations caused in the presence of norfloxacin. Moreover, to obtain an atomistic view, we complemented the experiments with millisecond-long free energy calculations based on temperature-accelerated Brownian dynamics simulations to identify the most probable permeation pathways of the antibiotics through the respective pores. Both, the experimental analysis and the computational modelling, suggest that norfloxacin is able to permeate through the larger porins, i.e., OmpF, OmpC, and Omp35, whereas it only binds to the slightly narrower porin Omp36. Moreover, divalent ions can bind to negatively charged residues inside the porin, reversing the ion selectivity of the pore. In addition, the divalent ions can chelate with the fluoroquinolone molecules and alter their physicochemical properties. The results suggest that the conjugation with either pores or molecules must break when the antibiotic molecules pass the lumen of the porin, with the conjugation to the antibiotic being more stable than that to the respective pore. In general, the permeation or binding process of fluoroquinolones in porins occurs irrespective of the presence of divalent ions, but the presence of divalent ions can vary the kinetics significantly. Thus, a detailed investigation of the interplay of divalent ions with antibiotics and pores is of key importance in developing new antimicrobial drugs.

Exploration of Free Energy Surfaces Across a Membrane Channel Using Metadynamics and Umbrella Sampling.

To reach their site of action, it is essential for antibiotic molecules to cross the bacterial outer membrane. The progress of enhanced sampling techniques in molecular dynamics simulations enables us to understand these translocations at an atomic level. To this end, calculations of free energy surfaces for these permeation processes are of key importance. Herein, we investigate the translocation of a variety of anionic solutes through the outer membrane pore OprO of the Gram-negative bacterium Pseudomonas aeruginosa using the metadynamics and umbrella sampling techniques at the all-atom level. Free energy calculations have been performed employing these two distinct methods in order to illustrate the difference in computed free energies, if any. The investigated solutes range from a single atomic chloride ion over a multiatomic monophosphate ion to a more bulky fosmidomycin antibiotic. The role of complexity of the permeating solutes in estimating accurate free energy profiles is demonstrated by performing extensive convergence analysis. For simple monatomic ions, good agreement between the well-tempered metadynamics and the umbrella sampling approaches is achieved, while for the permeation of the monophosphate ion differences start to appear. In the case of larger molecules such as fosmidomycin it is a tough challenge to achieve converged free energy profiles. This issue is mainly due to neglecting orthogonal degrees of freedom during the free energy calculations. Nevertheless, the freely driven metadynamics approach leads to clearly advantageous results. Additionally, atomistic insights of the translocation mechanisms of all three solutes are discussed.

DFTB/MM Molecular Dynamics Simulations of the FMO Light-Harvesting Complex.

Because of the size of light-harvesting complexes and the involvement of electronic degrees of freedom, computationally these systems need to be treated with a combined quantum-classical description. To this end, Born-Oppenheimer molecular dynamics simulations have been employed in a quantum mechanics/molecular mechanics (QM/MM) fashion for the ground state followed by excitation energy calculations again in a QM/MM scheme for the Fenna-Matthews-Olson (FMO) complex. The self-consistent-charge density functional tight-binding (DFTB) method electrostatically coupled to a classical description of the environment was applied to perform the ground-state dynamics. Subsequently, long-range-corrected time-dependent DFTB calculations were performed to determine the excitation energy fluctuations of the individual bacteriochlorophyll a molecules. The spectral densities obtained using this approach show an excellent agreement with experimental findings. In addition, the fluctuating site energies and couplings were used to estimate the exciton transfer dynamics.

Molecular Interactions of Cephalosporins with the Deep Binding Pocket of the RND Transporter AcrB.

The drug/proton antiporter AcrB, part of the major efflux pump AcrABZ-TolC in Escherichia coli, is characterized by its impressive ability to transport chemically diverse compounds, conferring a multidrug resistance phenotype. However, the molecular features differentiating between good and poor substrates of the pump have yet to be identified. In this work, we combined molecular docking with molecular dynamics simulations to study the interactions between AcrB and two representative cephalosporins, cefepime and ceftazidime (a good and poor substrate of AcrB, respectively). Our analysis revealed different binding preferences of the two compounds toward the subsites of the large deep binding pocket of AcrB. Cefepime, although less hydrophobic than ceftazidime, showed a higher affinity than ceftazidime for the so-called hydrophobic trap, a region known for binding inhibitors and substrates. This supports the hypothesis that surface complementarity between the molecule and AcrB, more than the intrinsic hydrophobicity of the antibiotic, is a feature required for the interaction within this region. Oppositely, the preference of ceftazidime for binding outside the hydrophobic trap might not be optimal for triggering allosteric conformational changes needed to the transporter to accomplish its function. Altogether, our findings could provide valuable information for the design of new antibiotics less susceptible to the efflux mechanism.

A Multidisciplinary Approach toward Identification of Antibiotic Scaffolds for Acinetobacter baumannii.

Research efforts to discover potential new antibiotics for Gram-negative bacteria suffer from high attrition rates due to the synergistic action of efflux systems and the limited permeability of the outer membrane (OM). One strategy to overcome the OM permeability barrier is to identify small molecules that are natural substrates for abundant OM channels and use such compounds as scaffolds for the design of efficiently permeating antibacterials. Here we present a multidisciplinary approach to identify such potential small-molecule scaffolds. Focusing on the pathogenic bacterium Acinetobacter baumannii, we use OM proteomics to identify DcaP as the most abundant channel during infection in rodents. The X-ray crystal structure of DcaP reveals a trimeric, porin-like structure and suggests that dicarboxylic acids are potential transport substrates. Electrophysiological experiments and all-atom molecular dynamics simulations confirm this notion and provide atomistic information on likely permeation pathways and energy barriers for several small molecules, including a clinically relevant ß-lactamase inhibitor.

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.

Characterization of Ciprofloxacin Permeation Pathways across the Porin OmpC Using Metadynamics and a String Method.

The rapid spreading of antimicrobial resistance in Gram-negative bacteria has become a major threat for humans as well as animals. As one of the main factors involved, the permeability of the outer membrane has attracted a great deal of attention recently. However, the knowledge regarding the translocation mechanisms for most available antibiotics is so far rather limited. Here, a theoretical study concerning the diffusion route of ciprofloxacin across the outer membrane porin OmpC from E. coli is presented. To this end, we establish a protocol to characterize meaningful permeation pathways by combining metadynamics with the zero-temperature string method. It was found that the lowest-energy pathway requires a reorientation of ciprofloxacin in the extracellular side of the porin before reaching the constriction region with its carboxyl group ahead. Several affinity sites have been identified, and their metastability has been evaluated using unbiased simulations. Such a detailed understanding is potentially very helpful in guiding the development of next generation antibiotics.