A structural study of ion permeation in OmpF porin from anomalous X-ray diffraction and molecular dynamics simulations

Impact of bacterial outer membrane and general porins on cyanide diffusion and biodegradation kinetics

The present study focuses on the analysis of the diffusion process of various cyanide compounds through general porins and outer membranes of gram-negative bacteria. We demonstrate the impact of the compound-to-porin radius ratio, the charge of cyanide ion, the Donnan potential, the intrinsic porin potential, the number and length of general porins, the fraction of open channels, and the size of bacteria on the effective diffusion coefficients and permeability coefficients of cyanide compounds. Moreover, we report, for the first time, the procedure for comparison of the rate of cyanide diffusion across the outer membrane with the rate of cyanide biodegradation that allows establishing the conditions for which the biodegradation is a diffusion-limited process or the diffusion is a significantly faster process than biodegradation. We apply this procedure to several experimental studies and predict the range of extracellular cyanide concentrations for which diffusion is a significantly faster process than biodegradation. We also demonstrate how these results affect the theoretical view of the cyanide biodegradation kinetics.

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.

The Effect of Lipopolysaccharides on the Electrostatic Properties of Gram-Negative General Porins from Enterobacteriaceae

We investigated, by using all-atom molecular dynamics simulations, the effect of the outer membrane of Gram-negative bacteria, composed in the outer leaflet by polar/charged lipopolysaccharides (LPS), on the electrostatic properties of general porins from the Enterobacteriaceae family. General porins constitute the main path for the facilitated diffusion of polar antibiotics through the outer membrane. As model system we selected OmpK36 from Klebsiella pneumoniae, the ortholog of OmpC from Escherichia coli. This species presents high variability of amino acid composition of porins, with the effect to increase its resistance to the penetration of antibiotics. The various properties we analyzed seem to indicate that LPS acts as an independent layer without affecting the internal electrostatic properties of OmpK36. The only apparent effect on the microsecond time scale we sampled is the appearance of calcium ions, when present at moderate concentration in solution, inside the pore. However, we noticed increased fluctuations of the polarization density and only minor changes on its average value.

Rapid antibiotic screening based on E. coli apoptosis using a potentiometric sensor array

Phenotypic antimicrobial susceptibility testing enables reliable antibiotic screening but requires multiple strategies to identify each phenotypic change induced by different bactericidal mechanisms. Bacteria apoptosis with typical phenotypic features has never been explored for antibiotic screening. Herein, we developed an antibiotic screening method based on the measurement of antibiotic-induced phosphatidylserine (PS) exposure of apoptotic bacteria. Phosphatidylserine externalization of E. coli that can be widely used as an apoptosis marker for antibiotics with different antibacterial mechanisms was explored. A positively charged PS-binding peptide was immobilized on magnetic beads (MBs) to recognize and capture apoptotic E. coli with PS externalization. Apoptotic E. coli binding led to the charge or charge density change of MBs-peptide, resulting in a potential change on a magneto-controlled polymeric membrane potentiometric sensor. Based on the detection of apoptotic E. coli killed by antibiotics, antibiotic screening for different classes of antibiotics and silver nanoparticles was achieved within 1.5 h using a potentiometric sensor array. This approach enables sensitive, general, and time-saving antibiotic screening, and may open up a new path for antibiotic susceptibility testing.

Charging Metal-Organic Framework Membranes by Incorporating Crown Ethers to Capture Cations for Ion Sieving

Protein channels on the biofilm conditionally manipulate ion transport via regulating the distribution of charge residues, making analogous processes on artificial membranes a hot spot and challenge. Here, we employ metal-organic frameworks (MOFs) membrane with charge-adjustable subnano-channel to selectively govern ion transport. Various valent ions are binded with crown ethers embedded in the MOF cavity, which act as charged guest to regulate the channels' charge state from the negativity to positivity. Compared with the negatively charged channel, the positive counterpart obviously enhances Li+ /Mg2+ selectivity, which benefit from the reinforcement of the electrostatic repulsion between ions and the channel. Meanwhile, theoretical calculations reveal that Mg2+ transport through the more positively charged channel needed to overcome higher entrance energy barrier than that of Li+ . This work provides a subtle strategy for ion-selective transport upon regulating the charge state of insulating membrane, which paves the way for the application like seawater desalination and lithium extraction from salt lakes.

Interaction of Tryptophan- and Arginine-Rich Antimicrobial Peptide with E. coli Outer Membrane-A Molecular Simulation Approach

A short antimicrobial peptide (AMP), rich in tryptophan and arginine (P6-HRWWRWWRR-NH2), was used in molecular dynamics (MD) simulations to investigate the interaction between AMPs and lipopolysaccharides (LPS) from two E. coli outer membrane (OM) membrane models. The OM of Gram-negative bacteria is an asymmetric bilayer, with the outer layer consisting exclusively of lipopolysaccharide molecules and the lower leaflet made up of phospholipids. The mechanisms by which short AMPs permeate the OM of Gram-negative bacteria are not well understood at the moment. For this study, two types of E. coli OM membrane models were built with (i) smooth LPS composed of lipid A, K12 core and O21 O-antigen, and (ii) rough type LPS composed of lipid A and R1 core. An OmpF monomer from E. coli was embedded in both membrane models. MD trajectories revealed that AMP insertion in the LPS layer was facilitated by the OmpF-created gap and allowed AMPs to form hydrogen bonds with the phosphate groups of inner core oligosaccharides. OM proteins such as OmpF may be essential for the permeation of short AMPs such as P6 by exposing the LPS binding site or even by direct translocation of AMPs across the OM.

Design and directed evolution of noncanonical β-stereoselective metalloglycosidases

Metallohydrolases are ubiquitous in nearly all subclasses of hydrolases, utilizing metal elements to activate a water molecule and facilitate its subsequent dissociation of diverse chemical bonds. However, such a catalytic role of metal ions is rarely found with glycosidases that hydrolyze the glycosidic bonds in sugars. Herein, we design metalloglycosidases by constructing a hydrolytically active Zn-binding site within a barrel-shaped outer membrane protein OmpF. Structure- and mechanism-based redesign and directed evolution have led to the emergence of Zn-dependent glycosidases with catalytic proficiency of 2.8 × 10⁹ and high β-stereoselectivity. Biochemical characterizations suggest that the Zn-binding site constitutes a key catalytic motif along with at least one adjacent acidic residue. This work demonstrates that unprecedented metalloenzymes can be tailor-made, expanding the scope of inorganic reactivities in proteinaceous environments, resetting the structural and functional diversity of metalloenzymes, and providing the potential molecular basis of unidentified metallohydrolases and novel whole-cell biocatalysts.

What have molecular simulations contributed to understanding of Gram-negative bacterial cell envelopes?

Bacterial cell envelopes are compositionally complex and crowded and while highly dynamic in some areas, their molecular motion is very limited, to the point of being almost static in others. Therefore, it is no real surprise that studying them at high resolution across a range of temporal and spatial scales requires a number of different techniques. Details at atomistic to molecular scales for up to tens of microseconds are now within range for molecular dynamics simulations. Here we review how such simulations have contributed to our current understanding of the cell envelopes of Gram-negative bacteria.

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.

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.

Dynamics of OmpF Trimer Formation in the Bacterial Outer Membrane of Escherichia coli

The self-assembly of outer membrane protein F (OmpF) in the outer membrane of Escherichia coli Gram-negative bacteria was studied using multiscale molecular dynamics simulations. To accommodate the long time scale required for protein assembly, coarse-grained parametrization of E. coli outer membrane lipids was first developed. The OmpF monomers formed stable dimers at specific protein-protein interactions sites irrespective of the lipid membrane environment. The dimer intermediate was asymmetric but provided a template to form a symmetric trimer. Superposition analysis of the self-assembled trimer with the X-ray crystal structure of the trimer available in the protein data bank showed excellent agreement with global root-mean-square deviation of less than 2.2 Å. The free energy change associated with dimer formation was -26 ± 1 kcal mol^-1, and for a dimer to bind to a monomer and to form a trimer yielded -56 ± 4 kcal mol^-1. Based on thermodynamic data, an alternate path to trimer formation via interaction of two dimers is also presented.

Native State Organization of Outer Membrane Porins Unraveled by HDx-MS

Hydrogen-deuterium exchange (HDx) associated with mass spectrometry (MS) is emerging as a powerful tool to provide conformational information about membrane proteins. Unfortunately, as for X-ray diffraction and NMR, HDx performed on reconstituted in vitro systems might not always reflect the in vivo environment. Outer-membrane vesicles naturally released by Escherichia coli were used to carry out analysis of native OmpF through HDx-MS. A new protocol compatible with HDx analysis that avoids hindrance from the lipid contents was setup. The extent of deuterium incorporation was in good agreement with the X-ray diffraction data of OmpF as the buried β-barrels incorporated a low amount of deuterium, whereas the internal loop L3 and the external loops incorporated a higher amount of deuterium. Moreover, the kinetics of incorporation clearly highlights that peptides segregate well in two distinct groups based exclusively on a trimeric organization of OmpF in the membrane: peptides presenting fast kinetics of labeling are facing the complex surrounding environment, whereas those presenting slow kinetics are located in the buried core of the trimer. The data show that HDx-MS applied to a complex biological system is able to reveal solvent accessibility and spatial arrangement of an integral outer-membrane protein complex.

The impact of length variations in the L2 loop on the structure and thermal stability of non-specific porins: The case of OmpCs from the Yersinia pseudotuberculosis complex

Porins are integral proteins of the outer membranes of gram-negative bacteria. In membranes, they exist as homotrimers and the L2 loops contribute to their stability. Comparison of OmpC porins of the Yersinia pseudotuberculosis complex with other enterobacterial porins demonstrated L2 loop length diversity, which is caused by varying numbers of dipeptide/tripeptide repeats. The OmpC porins are highly homologous to each other, and they can be subdivided into five isoforms based on their L2 loop structure. Optical spectroscopy and SDS-PAGE experiments revealed that particularities of the L2 loops affected the structure and thermal stability of the porins. Thermal denaturation studies showed that porins with shorter loops, compared to porins with longer loops, had more stable tertiary and less stable secondary and quaternary structures. According to our comparative modeling results, the L2 loops differ in their structure by adopting different spatial positions and forming different polar bonds with a neighbor monomer. The replacement of asparagine with arginine at the C-terminus of the L2 loop shifts the loop upwards and causes the loss of contacts with the arginine clusters within the pores. The increase in the length of these loops ensures that they shift down toward the pore and restore their contacts with arginines on the channel wall, as is the case in classical nonspecific porins. Despite the fact that the surface charge density varies considerably among the OmpC porins, the L2 loops form a typical negatively charged region in the center of the trimer.

Bacterial Outer Membrane Porins as Electrostatic Nanosieves: Exploring Transport Rules of Small Polar Molecules

Transport of molecules through biological membranes is a fundamental process in biology, facilitated by selective channels and general pores. The architecture of some outer membrane pores in Gram-negative bacteria, common to other eukaryotic pores, suggests them as prototypes of electrostatically regulated nanosieve devices. In this study, we sensed the internal electrostatics of the two most abundant outer membrane channels of Escherichia coli, using norfloxacin as a dipolar probe in single molecule electrophysiology. The voltage dependence of the association rate constant of norfloxacin interacting with these nanochannels follows an exponential trend, unexpected for neutral molecules. We combined electrophysiology, channel mutagenesis, and enhanced sampling molecular dynamics simulations to explain this molecular mechanism. Voltage and temperature dependent ion current measurements allowed us to quantify the transversal electric field inside the channel as well as the distance where the applied potential drops. Finally, we proposed a general model for transport of polar molecules through these electrostatic nanosieves. Our model helps to further understand the basis for permeability in Gram-negative pathogens, contributing to fill in the innovation gap that has limited the discovery of effective antibiotics in the last 20 years.

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.

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.

Journey of poly-nucleotides through OmpF porin

OmpF is an abundant porin in many bacteria which attracts attention as a promising biological nanopore for DNA sequencing. We study the interactions of OmpF with pentameric poly-nucleotides (poly-Ns) in silico. The poly-N molecule is forced to translocate through the lumen of OmpF. Subsequently, the structural and dynamical effects of translocation steps on protein and poly-N molecules are explored in detail. The external loops of OmpF are introduced as the main region for discrimination of poly-Ns based on their organic bases. Structural network analyses of OmpF in the presence or absence of poly-Ns characterize special residues in the structural network of porin. These residues pave the way for engineering OmpF protein. The poly-N-specific pattern of OmpF's local conductance is detected in the current study. Computing the potential of mean force for translocation steps, we define the energetic barrier ahead of poly-N to move through OmpF's lumen. We suggest that fast translocation of the examined poly-N molecules through OmpF seems unattainable by small external driving forces. Our computational results suggest some abilities for OmpF porin like OmpF's potential for being used in poly-N sequencing.

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1-2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug-specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans-membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α-helical MPs as well as β-barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge-based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade.