Curvature sensing lipid dynamics in a mitochondrial inner membrane model

Membrane curvature is essential for many cellular structures and processes, and factors such as leaflet asymmetry, lipid composition, and proteins all play important roles. Cardiolipin is the signature lipid of mitochondrial membranes and is essential for maintaining the highly curved shapes of the inner mitochondrial membrane (IMM) and the spatial arrangement of membrane proteins. In this study, we investigate the partitioning behavior of various lipids present in the IMM using coarse-grained molecular dynamics simulations. This study explores curved bilayer systems containing phosphatidylcholine (PC), phosphatidylethanolamine (PE), and cardiolipin (CDL) in binary and ternary component mixtures. Curvature properties such as mean and Gaussian curvatures, as well as the distribution of lipids into the various curved regions of the cristae models, are quantified. Overall, this work represents an advance beyond previous studies on lipid curvature sensing by simulating these systems in a geometry that has the morphological features and scales of curvature consistent with regions of the IMM. We find that CDL has a stronger preference for accumulating in regions of negative curvature than PE lipids, in agreement with previous results. Furthermore, we find lipid partitioning propensity is dominated by sensitivity to mean curvature, while there is a weaker correlation with Gaussian curvature.

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.

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.

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.

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.