Tuning the affinity of anion binding sites in porin channels with negatively charged residues: molecular details for OprP

Importance of the lysine cluster in the translocation of anions through the pyrophosphate specific channel OprO

Pseudomonas aeruginosa is a Gram-negative bacterium with an intrinsic resistance towards antibiotics due to the lack of a large diffusion pores. Exchange of substances with the environment is done mainly through a set of narrow and substrate-specific porins in its outer membrane that filter molecules according to their size and chemical composition. Among these proteins are OprP and OprO involved in the selective uptake of mono- and pyrophosphates, respectively. Both proteins are homotrimers and each monomer features an hourglass-shaped channel structure including a periplasmic cavity with a lysine cluster. In this study, we focus on the characterization of this lysine cluster in OprO. The importance of these lysine residues was shown with alanine substitutions in single channel conductance experiments, by titration of mono- and pyrophosphate in multi-channel analysis and by molecular dynamics simulations. All obtained data demonstrated that the closer the mutated lysine residues are to arginine 133, the lower gets the single channel conductance. It was found that the ion flow through each monomer can follow two different lysine paths indicating that phosphate ions have a larger freedom on the periplasmic side of the constriction region. Our results emphasize the important role of the lysine residue 121 in the binding site together with arginine 133 and aspartic acid 94. An improved understanding of the ion mobility across these channels can potentially lead to an optimized permeation of (phosphonic acid containing) antibiotics through the outer membrane of P. aeruginosa and the development of new drug molecules.

Permeation of Fosfomycin through the Phosphate-Specific Channels OprP and OprO of Pseudomonas aeruginosa

Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen responsible for many nosocomial infections. It is quite resistant to various antibiotics, caused by the absence of general diffusion pores in the outer membrane. Instead, it contains many substrate-specific channels. Among them are the two phosphate- and pyrophosphate-specific porins OprP and OprO. Phosphonic acid antibiotics such as fosfomycin and fosmidomycin seem to be good candidates for using these channels to enter P. aeruginosa bacteria. Here, we investigated the permeation of fosfomycin through OprP and OprO using electrophysiology and molecular dynamics (MD) simulations. The results were compared to those of the fosmidomycin translocation, for which additional MD simulations were performed. In the electrophysiological approach, we noticed a higher binding affinity of fosfomycin than of fosmidomycin to OprP and OprO. In MD simulations, the ladder of arginine residues and the cluster of lysine residues play an important role in the permeation of fosfomycin through the OprP and OprO channels. Molecular details on the permeation of fosfomycin through OprP and OprO channels were derived from MD simulations and compared to those of fosmidomycin translocation. In summary, this study demonstrates that the selectivity of membrane channels can be employed to improve the permeation of antibiotics into Gram-negative bacteria and especially into resistant P. aeruginosa strains.

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.

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.

Multisystem combined uranium resistance mechanisms and bioremediation potential of Stenotrophomonas bentonitica BII-R7: Transcriptomics and microscopic study

The potential use of microorganisms in the bioremediation of U pollution has been extensively described. However, a lack of knowledge on molecular resistance mechanisms has become a challenge for the use of these technologies. We reported on the transcriptomic and microscopic response of Stenotrophomonas bentonitica BII-R7 exposed to 100 and 250 μM of U. Results showed that exposure to 100 μM displayed up-regulation of 185 and 148 genes during the lag and exponential phases, respectively, whereas 143 and 194 were down-regulated, out of 3786 genes (>1.5-fold change). Exposure to 250 μM of U showed up-regulation of 68 genes and down-regulation of 290 during the lag phase. Genes involved in cell wall and membrane protein synthesis, efflux systems and phosphatases were up-regulated under all conditions tested. Microscopic observations evidenced the formation of U-phosphate minerals at membrane and extracellular levels. Thus, a biphasic process is likely to occur: the increased cell wall would promote the biosorption of U to the cell surface and its precipitation as U-phosphate minerals enhanced by phosphatases. Transport systems would prevent U accumulation in the cytoplasm. These findings contribute to an understanding of how microbes cope with U toxicity, thus allowing for the development of efficient bioremediation strategies.

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 MgCl₂ 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.

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.

Computational and Experimental Studies of Inhibitor Design for Aldolase A

Glycolytic enzyme fructose-bisphosphate aldolase A is an emerging therapeutic target in cancer. Recently, we have solved the crystal structure of murine aldolase in complex with naphthalene-2,6-diyl bisphosphate (ND1) that served as a template of the design of bisphosphate-based inhibitors. In this work, a series of ND1 analogues containing difluoromethylene (-CF₂), methylene (-CH₂), or aldehyde substitutions were designed. All designed compounds were studied using molecular dynamics (MD) simulations with the AMOEBA force field. Both energetics and structural analyses have been done to understand the calculated binding free energies. The average distances between ligand and protein atoms for ND1 were very similar to those for the ND1 crystal structure, which indicates that our MD simulation is sampling the correct conformation well. CF₂ insertion lowers the binding free energy by 10-15 kcal/mol, while CF₂ substitution slightly increases the binding free energy, which matches the experimental measurement. In addition, we found that NDB with two CF₂ insertions, the strongest binder, is entropically driven, while others including NDA with one CF₂ insertion are all enthalpically driven. This work provides insights into the mechanisms underlying protein-phosphate binding and enhances the capability of applying computational and theoretical frameworks to model, predict, and design diagnostic strategies targeting cancer.

Elucidating the Phosphate Binding Mode of Phosphate-Binding Protein: The Critical Effect of Buffer Solution

Phosphate is an essential component of cell functions, and the specific transport of phosphorus into a cell is mediated by phosphate-binding protein (PBP). The mechanism of PBP-phosphate recognition remains controversial: on the basis of similar binding affinities at acidic and basic pHs, it is believed that the hydrogen network in the binding site is flexible to adapt to different protonation states of phosphates. However, only hydrogen (1H) phosphate was observed in the sub-angstrom X-ray structures. To address this inconsistency, we performed molecular dynamics simulations using the AMOEBA polarizable force field. Structural and free energy data from simulations suggested that 1H phosphate was the preferred bound form at both pHs. The binding of dihydrogen (2H) phosphate disrupted the hydrogen-bond network in the PBP pocket, and the computed affinity was much weaker than that of 1H phosphate. Furthermore, we showed that the discrepancy in the studies described above is resolved if the interaction between phosphate and the buffer agent is taken into account. The calculated apparent binding affinities are in excellent agreement with experimental measurements. Our results suggest the high specificity of PBP for 1H phosphate and highlight the importance of the buffer solution for the binding of highly charged ligands.

Understanding antibiotic resistance via outer membrane permeability

Collective antibiotic drug resistance is a global threat, especially with respect to Gram-negative bacteria. The low permeability of the bacterial outer cell wall has been identified as a challenging barrier that prevents a sufficient antibiotic effect to be attained at low doses of the antibiotic. The Gram-negative bacterial cell envelope comprises an outer membrane that delimits the periplasm from the exterior milieu. The crucial mechanisms of antibiotic entry via outer membrane includes general diffusion porins (Omps) responsible for hydrophilic antibiotics and lipid-mediated pathway for hydrophobic antibiotics. The protein and lipid arrangements of the outer membrane have had a strong impact on the understanding of bacteria and their resistance to many types of antibiotics. Thus, one of the current challenges is effective interpretation at the molecular basis of the outer membrane permeability. This review attempts to develop a state of knowledge pertinent to Omps and their effective role in solute influx. Moreover, it aims toward further understanding and exploration of prospects to improve our knowledge of physicochemical limitations that direct the translocation of antibiotics via bacterial outer membrane.

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.

Structure, function and regulation of Pseudomonas aeruginosa porins

Pseudomonas aeruginosa is a Gram-negative bacterium belonging to the γ-proteobacteria. Like other members of the Pseudomonas genus, it is known for its metabolic versatility and its ability to colonize a wide range of ecological niches, such as rhizosphere, water environments and animal hosts, including humans where it can cause severe infections. Another particularity of P. aeruginosa is its high intrinsic resistance to antiseptics and antibiotics, which is partly due to its low outer membrane permeability. In contrast to Enterobacteria, pseudomonads do not possess general diffusion porins in their outer membrane, but rather express specific channel proteins for the uptake of different nutrients. The major outer membrane 'porin', OprF, has been extensively investigated, and displays structural, adhesion and signaling functions while its role in the diffusion of nutrients is still under discussion. Other porins include OprB and OprB2 for the diffusion of glucose, the two small outer membrane proteins OprG and OprH, and the two porins involved in phosphate/pyrophosphate uptake, OprP and OprO. The remaining nineteen porins belong to the so-called OprD (Occ) family, which is further split into two subfamilies termed OccD (8 members) and OccK (11 members). In the past years, a large amount of information concerning the structure, function and regulation of these porins has been published, justifying why an updated review is timely.

Conversion of OprO into an OprP-like Channel by Exchanging Key Residues in the Channel Constriction

Under phosphate-limiting conditions, the channels OprP and OprO are induced and expressed in the outer membrane of Pseudomonas aeruginosa. Despite their large homology, the phosphate-specific OprP and the diphosphate-specific OprO pores show structural differences in their binding sites situated in the constriction region. Previously, it was shown that the mutation of amino acids in OprP (Y62F and Y114D) led to an exchange in substrate specificity similar to OprO. To support the role of these key amino acids in the substrate sorting of these specific channels, the reverse mutants for OprO (F62Y, D114Y, and F62Y/D114Y) were created in this study. The phosphate and diphosphate binding of the generated channels was studied in planar lipid bilayers. Our results show that mutations of key residues indeed reverse the substrate specificity of OprO to OprP and support the view that just a few strategically positioned amino acids are mainly responsible for its substrate specificity.

Exploring bacterial outer membrane barrier to combat bad bugs

One of the main fundamental mechanisms of antibiotic resistance in Gram-negative bacteria comprises an effective change in the membrane permeability to antibiotics. The Gram-negative bacterial complex cell envelope comprises an outer membrane that delimits the periplasm from the exterior environment. The outer membrane contains numerous protein channels, termed as porins or nanopores, which are mainly involved in the influx of hydrophilic compounds, including antibiotics. Bacterial adaptation to reduce influx through these outer membrane proteins (Omps) is one of the crucial mechanisms behind antibiotic resistance. Thus to interpret the molecular basis of the outer membrane permeability is the current challenge. This review attempts to develop a state of knowledge pertinent to Omps and their effective role in antibiotic influx. Further, it aims to study the bacterial response to antibiotic membrane permeability and hopefully provoke a discussion toward understanding and further exploration of prospects to improve our knowledge on physicochemical parameters that direct the translocation of antibiotics through the bacterial membrane protein channels.

Structure, Dynamics, and Substrate Specificity of the OprO Porin from Pseudomonas aeruginosa

The outer membrane (OM) of Gram-negative bacteria functions as a selective permeability barrier between cell and environment. For nutrient acquisition, the OM contains a number of channels that mediate uptake of small molecules by diffusion. Many of these channels are specific, i.e., they prefer certain substrates over others. In electrophysiological experiments, the OM channels OprP and OprO from Pseudomonas aeruginosa show a specificity for phosphate and diphosphate, respectively. In this study we use x-ray crystallography, free-energy molecular dynamics (MD) simulations, and electrophysiology to uncover the atomic basis for the different substrate specificity of these highly similar channels. A structural analysis of OprP and OprO revealed two crucial differences in the central constriction region. In OprP there are two tyrosine residues, Y62 and Y114, whereas the corresponding residues in OprO are phenylalanine F62 and aspartate D114. To probe the importance of these two residues in generating the different substrate specificities, the double mutants were generated in silico and in vitro. Applied-field MD simulations and electrophysiological experiments demonstrated that the double mutations interchange the phosphate and diphosphate specificities of OprP and OprO. Our findings outline a possible strategy to rationally design channel specificity by modification of a small number of residues that may be applicable to other pores as well.

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.

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.

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.