Investigating reaction pathways in rare events simulations of antibiotics diffusion through protein channels

In Gram-negative bacteria, outer-membrane protein channels, such as OmpF of Escherichia coli, constitute the entry point of various classes of antibiotics. While antibacterial research and development is declining, bacterial resistance to antibiotics is rising and there is an emergency call for a new way to develop potent antibacterial agents and to bring them to the market faster and at reduced cost. An emerging strategy is to follow a bottom-up approach based on microscopically founded computational based screening, however such strategy needs better-tuned methods. Here we propose to use molecular dynamics (MD) simulations combined with the metadynamics algorithm, to study antibiotic translocation through OmpF at a molecular scale. This recently designed algorithm overcomes the time scale problem of classical MD by accelerating some reaction coordinates. It is expected that the initial assumption of the reaction coordinates is a key determinant for the efficiency and accuracy of the simulations. Previous studies using different computational schemes for a similar process only used one reaction coordinate, which is the directionality. Here we go further and see how it is possible to include more informative reaction coordinates, accounting explicitly for: (i) the antibiotic flexibility and (ii) interactions with the channel. As model systems, we select two compounds covering the main classes of antibiotics, ampicillin and moxifloxacine. We decipher the molecular mechanism of translocation of each antibiotic and highlight the important parameters that should be taken into account for improving further simulations. This will benefit the screening and design for antibiotics with better permeation properties.     

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.     

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.     

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.     

Brownian dynamics simulation of ion flow through porin channels

Bacterial porins, which allow the passage of solutes across the outer bacterial membrane, are structurally well characterized. They therefore lend themselves to detailed studies of the determinants of ion flow through transmembraneous channels. In a comparative study, we have performed Brownian dynamics simulations to obtain statistically significant transfer efficiencies for cations and anions through matrix porin OmpF, osmoporin OmpK36, phosphoporin PhoE and two OmpF charge mutants.The simulations show that the electrostatic potential at the highly charged channel constriction serves to enhance ion permeability of either cations or anions, dependent on the type of porin. At the same time translocation of counterions is not severely impeded. At the constriction, cations and anions follow distinct trajectories, due to the segregation of basic and acidic protein residues.Simulated ion selectivity and relative conductance agree well with experimental values, and are dependent crucially on the charge constellation at the pore constriction. The experimentally observed decrease in ion selectivity and single channel conductance with increasing ionic strength is well reproduced and can be attributed to electrostatic shielding of the pore lining.     

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.     

Insertion of short amino-functionalized single-walled carbon nanotubes into phospholipid bilayer occurs by passive diffusion

Carbon nanotubes have been proposed to be efficient nanovectors able to deliver genetic or therapeutic cargo into living cells. However, a direct evidence of the molecular mechanism of their translocation across cell membranes is still needed. Here, we report on an extensive computational study of short (5 nm length) pristine and functionalized single-walled carbon nanotubes uptake by phospholipid bilayer models using all-atom molecular dynamics simulations. Our data support the hypothesis of a direct translocation of the nanotubes through the phospholipid membrane. We find that insertion of neat nanotubes within the bilayer is a "nanoneedle" like process, which can often be divided in three consecutive steps: landing and floating, penetration of the lipid headgroup area and finally sliding into the membrane core. The presence of functional groups at moderate concentrations does not modify the overall scheme of diffusion mechanism, provided that their deprotonated state favors translocation through the lipid bilayer.     

Mechanism for translocation of fluoroquinolones across lipid membranes

Classical atom-scale molecular dynamics simulations, constrained free energy calculations, and quantum mechanical (QM) calculations are employed to study the diffusive translocation of ciprofloxacin (CPFX) across lipid membranes. CPFX is considered here as a representative of the fluoroquinolone antibiotics class. Neutral and zwitterionic CPFX coexist at physiological pH, with the latter being predominant. Simulations reveal that only the neutral form permeates the bilayer, and it does so through a novel mechanism that involves dissolution of concerted stacks of zwitterionic ciprofloxacins. Subsequent QM analysis of the observed molecular stacking shows the important role of partial charge neutralization in the stacks, highlighting how the zwitterionic form of the drug is neutralized for translocation. The findings propose a translocation mechanism in which zwitterionic CPFX molecules approach the membrane in stacks, but they diffuse through the membrane as neutral CPFX monomers due to intermolecular transfer of protons favored by partial solvation loss. The mechanism is expected to be of importance in the permeation and translocation of a variety of ampholitic drugs with stacking tendencies.     

Electroosmosis Dominates Electrophoresis of Antibiotic Transport Across the Outer Membrane Porin F

We report that the dynamics of antibiotic capture and transport across a voltage-biased OmpF nanopore is dominated by the electroosmotic flow rather than the electrophoretic force. By reconstituting an OmpF porin in an artificial lipid bilayer and applying an electric field across it, we are able to elucidate the permeation of molecules and their mechanism of transport. This field gives rise to an electrophoretic force acting directly on a charged substrate but also indirectly via coupling to all other mobile ions, causing an electroosmotic flow. The directionality and magnitude of this flow depends on the selectivity of the channel. Modifying the charge state of three different substrates (norfloxacin, ciprofloxacin, and enoxacin) by varying the pH between 6 and 9 while the charge and selectivity of OmpF is conserved allows us to work under conditions in which electroosmotic flow and electrophoretic forces add or oppose. This configuration allows us to identify and distinguish the contributions of the electroosmotic flow and the electrophoretic force on translocation. Statistical analysis of the resolvable dwell times reveals rich kinetic details regarding the direction and the stochastic movement of antibiotics inside the nanopore. We quantitatively describe the electroosmotic velocity component experienced by the substrates and their diffusion coefficients inside the porin with an estimate of the energy barrier experienced by the molecules caused by the interaction with the channel wall, which slows down the permeation by several orders of magnitude.     

Crystal structures of the OmpF porin: function in a colicin translocon

The OmpF porin in the Escherichia coli outer membrane (OM) is required for the cytotoxic action of group A colicins, which are proposed to insert their translocation and active domains through OmpF pores. A crystal structure was sought of OmpF with an inserted colicin segment. A 1.6 A OmpF structure, obtained from crystals formed in 1 M Mg2+, has one Mg2+ bound in the selectivity filter between Asp113 and Glu117 of loop 3. Co-crystallization of OmpF with the unfolded 83 residue glycine-rich N-terminal segment of colicin E3 (T83) that occludes OmpF ion channels yielded a 3.0 A structure with inserted T83, which was obtained without Mg2+ as was T83 binding to OmpF. The incremental electron density could be modelled as an extended poly-glycine peptide of at least seven residues. It overlapped the Mg2+ binding site obtained without T83, explaining the absence of peptide binding in the presence of Mg2+. Involvement of OmpF in colicin passage through the OM was further documented by immuno-extraction of an OM complex, the colicin translocon, consisting of colicin E3, BtuB and OmpF.     

Molecular dynamics simulations of DNA within a nanopore: arginine-phosphate tethering and a binding/sliding mechanism for translocation

Protein nanopores show great potential as low-cost detectors in DNA sequencing devices. To date, research has largely focused on the staphylococcal pore α-hemolysin (αHL). In the present study, we have developed simplified models of the wild-type αHL pore and various mutants in order to study the translocation dynamics of single-stranded DNA under the influence of an applied electric field. The model nanopores reflect the experimentally measured conductance values in planar lipid bilayers. We show that interactions between rings of cationic amino acids and DNA backbone phosphates result in metastable tethering of nucleic acid molecules within the pore, leading us to propose a "binding and sliding" mechanism for translocation. We also observe folding of DNA into nonlinear conformational intermediates during passage through the confined nanopore environment. Despite adopting nonlinear conformations, the DNA hexamer always exits the pore in the same orientation as it enters (3' to 5') in our simulations. The observations from our simulations help to rationalize experimentally determined trends in residual current and translocation efficiency for αHL and its mutants.     

Structure-Based Simulations of the Translocation Mechanism of the Hepatitis C Virus NS3 Helicase along Single-Stranded Nucleic Acid

The NS3 helicase of Hepatitis C virus is an ATP-fueled molecular motor that can translocate along single-stranded (ss) nucleic acid, and unwind double-stranded nucleic acids. It makes a promising antiviral target and an important prototype system for helicase research. Despite recent progress, the detailed mechanism of NS3 helicase remains unknown. In this study, we have combined coarse-grained (CG) and atomistic simulations to probe the translocation mechanism of NS3 helicase along ssDNA. At the residue level of detail, our CG simulations have captured functionally important interdomain motions of NS3 helicase and reproduced single-base translocation of NS3 helicase along ssDNA in the 3′–5′ direction, which is in good agreement with experimental data and the inchworm model. By combining the CG simulations with residue-specific perturbations to protein-DNA interactions, we have identified a number of key residues important to the translocation machinery that agree with previous structural and mutational studies. Additionally, our atomistic simulations with targeted molecular dynamics have corroborated the findings of CG simulations and further revealed key protein-DNA hydrogen bonds that break/form during the transitions. This study offers, to our knowledge, the most detailed and realistic simulations of translocation mechanism of NS3 helicase. The simulation protocol established in this study will be useful for designing inhibitors that target the translocation machinery of NS3 helicase, and for simulations of a variety of nucleic-acid-based molecular motors.     

Dynamics and Interactions of OmpF and LPS: Influence on Pore Accessibility and Ion Permeability

The asymmetric outer membrane of Gram-negative bacteria is formed of the inner leaflet with phospholipids and the outer leaflet with lipopolysaccharides (LPS). Outer membrane protein F (OmpF) is a trimeric porin responsible for the passive transport of small molecules across the outer membrane of Escherichia coli. Here, we report the impact of different levels of heterogeneity in LPS environments on the structure and dynamics of OmpF using all-atom molecular dynamics simulations. The simulations provide insight into the flexibility and dynamics of LPS components that are highly dependent on local environments, with lipid A being the most rigid and O-antigen being the most flexible. Increased flexibility of O-antigen polysaccharides is observed in heterogeneous LPS systems, where the adjacent O-antigen repeating units are weakly interacting and thus more dynamic, compared to homogeneous LPS systems in which LPS interacts strongly with each other with limited overall flexibility due to dense packing. The model systems were validated by comparing molecular-level details of interactions between OmpF surface residues and LPS core sugars with experimental data, establishing the importance of LPS core oligosaccharides in shielding OmpF surface epitopes recognized by monoclonal antibodies. There are LPS environmental influences on the movement of bulk ions (K⁺ and Cl⁻), but the ion selectivity of OmpF is mainly affected by bulk ion concentration.     

Altered antibiotic transport in OmpC mutants isolated from a series of clinical strains of multi-drug resistant E. coli

Antibiotic-resistant bacteria, particularly gram negative species, present significant health care challenges. The permeation of antibiotics through the outer membrane is largely effected by the porin superfamily, changes in which contribute to antibiotic resistance. A series of antibiotic resistant E. coli isolates were obtained from a patient during serial treatment with various antibiotics. The sequence of OmpC changed at three positions during treatment giving rise to a total of four OmpC variants (denoted OmpC20, OmpC26, OmpC28 and OmpC33, in which OmpC20 was derived from the first clinical isolate). We demonstrate that expression of the OmpC K12 porin in the clinical isolates lowers the MIC, consistent with modified porin function contributing to drug resistance. By a range of assays we have established that the three mutations that occur between OmpC20 and OmpC33 modify transport of both small molecules and antibiotics across the outer membrane. This results in the modulation of resistance to antibiotics, particularly cefotaxime. Small ion unitary conductance measurements of the isolated porins do not show significant differences between isolates. Thus, resistance does not appear to arise from major changes in pore size. Crystal structures of all four OmpC clinical mutants and molecular dynamics simulations also show that the pore size is essentially unchanged. Molecular dynamics simulations suggest that perturbation of the transverse electrostatic field at the constriction zone reduces cefotaxime passage through the pore, consistent with laboratory and clinical data. This subtle modification of the transverse electric field is a very different source of resistance than occlusion of the pore or wholesale destruction of the transverse field and points to a new mechanism by which porins may modulate antibiotic passage through the outer membrane.     

Antibiotic translocation through membrane channels: temperature-dependent ion current fluctuation for catching the fast events

Temperature-dependent facilitated permeation of antibiotics through membrane channels was investigated. Here we reconstituted single OmpF trimers from the outer membrane of Escherichia coli (E. coli) into a planar lipid bilayer. The penetration of ampicillin through OmpF causes fluctuation in the ion current, and analysis of the fluctuations at different temperatures allows us to determine the mode of permeation. The residence time of the drug inside the channel decays strongly with temperature, reaching the resolution limit of the instrument at 30°C. The number of events increases exponentially with temperature up to 30°C and then gradually decreases as temperature increases. At room temperature, we observe about 25 events per second per monomer of the trimeric channel and an extrapolation to 37°C gives roughly 50 events. The activation energy for ampicillin translocation through OmpF is estimated to be around 13 kT. Temperature-dependent study gives new insights into the faster translocation of small substrates through biological nanopores.     

Microscopic Kinetics of DNA Translocation through synthetic nanopores

We have previously demonstrated that a nanometer-diameter pore in a nanometer-thick metal-oxide-semiconductor-compatible membrane can be used as a molecular sensor for detecting DNA. The prospects for using this type of device for sequencing DNA are avidly being pursued. The key attribute of the sensor is the electric field-induced (voltage-driven) translocation of the DNA molecule in an electrolytic solution across the membrane through the nanopore. To complement ongoing experimental studies developing such pores and measuring signals in response to the presence of DNA, we conducted molecular dynamics simulations of DNA translocation through the nanopore. A typical simulated system included a patch of a silicon nitride membrane dividing water solution of potassium chloride into two compartments connected by the nanopore. External electrical fields induced capturing of the DNA molecules by the pore from the solution and subsequent translocation. Molecular dynamics simulations suggest that 20-basepair segments of double-stranded DNA can transit a nanopore of 2.2 x 2.6 nm(2) cross section in a few microseconds at typical electrical fields. Hydrophobic interactions between DNA bases and the pore surface can slow down translocation of single-stranded DNA and might favor unzipping of double-stranded DNA inside the pore. DNA occluding the pore mouth blocks the electrolytic current through the pore; these current blockades were found to have the same magnitude as the blockade observed when DNA transits the pore. The feasibility of using molecular dynamics simulations to relate the level of the blocked ionic current to the sequence of DNA was investigated.     

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.     

Metastable network model of protein transport through nuclear pores

To reconcile the observed selectivity and the high rate of translocation of cargo-importin complexes through nuclear pores, we propose that the core of the nuclear pore complex is blocked by a metastable network of phenylalanine and glycine nucleoporins. Although the network arrests the unfacilitated passage of objects larger than its mesh size, cargo-importin complexes act as catalysts that reduce the free energy barrier between the cross-linked and the dissociated states of the Nups, and open the network. Using Brownian dynamics simulations we calculate the distribution of passage times through the network for inert particles and cargo-importin complexes of different sizes and discuss the implications of our results for experiments on translocation of proteins through the nuclear pore complex.     

Free energy profiles for H+ conduction along hydrogen-bonded chains of water molecules.

The molecular mechanism for proton conduction along hydrogen-bonded chains, or "proton wires," is studied with free energy simulations. The complete transport of a charge along a proton wire requires two complementary processes: 1) translocation of an excess proton (propagation of an ionic defect), and 2) reorientation of the hydrogen-bonded chain (propagation of a bonding defect). The potential of mean force profile for these two steps is computed in model systems comprising a single-file chain of nine dissociable and polarizable water molecules represented by the PM6 model of Stillinger and co-workers. Results of molecular dynamics simulations with umbrella sampling indicate that the unprotonated chain is preferably polarized, and that the inversion of its total dipole moment involves an activation free energy of 8 kcal/mol. In contrast, the rapid translocation of an excess H+ across a chain extending between two spherical solvent droplets is an activationless process. These results suggest that the propagation of a bonding defect constitutes a limiting step for the passage of several protons along single-file chains of water molecules, whereas the ionic translocation may be fast enough to occur within the lifetime of transient hydrogen-bonded water chains in biological membranes.