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

Outer membrane permeability of Pseudomonas aeruginosa through β-lactams: new evidence on the role of OprD and OpdP porins in antibiotic resistance

Gram-negative bacteria are a major concern for public health, particularly due to the continuous rise of antibiotic resistance. A major factor that helps the development of resistance is the outer membrane that is essential since it acts as a strong permeability barrier to many antibiotics that are effective against other bacteria. In this study, we determine the specific permeability coefficients for various antibiotics in Pseudomonas aeruginosa strains, which differ from each other for their porin expressions. We showed that OprD and OpdP porins contribute both to internalize meropenem and biapenem. Using qRT-PCR, we demonstrated that their expression is dependent of the various phases of cellular growth. We were able to show how the OpdP porin is less expressed in exponential growth phases, while it tends to be produced when the bacterial culture enters into the latent phase, in an inversely proportional way compared to the OprD porin. The deletion of the OpdP porin, in the presence of meropenem at concentrations equivalent to the MIC values, contributes to the selection of carbapenem-resistant strains. Therefore, the presence of mutations/deletions of the OpdP porin should receive greater consideration from a clinical point of view as the use of meropenem at nonoptimal concentrations could lead to the appearance of resistance phenotypes.IMPORTANCECarbapenem-resistant strains of Pseudomonas aeruginosa are among the major threats to public health. The permeability of the outer membrane for the β-lactam antibiotics is one of the major factors that reduce the activity of the antibiotics. In this study, we measure the low permeability coefficient of the P. aeruginosa outer membrane to β-lactams. The methodology we develop to determine the permeability can be applied to other antibiotic families and/or pathogens.

Improved Free-Energy Estimates for the Permeation of Bulky Antibiotic Molecules through Porin Channels Using Temperature-Accelerated Sliced Sampling

The estimation of accurate free energies for antibiotic permeation via the bacterial outer-membrane porins has proven to be challenging. Atomistic simulations of the process suffer from sampling issues that are typical of systems with complex and slow dynamics, even with the application of advanced sampling methods. Ultimately, the objective is to obtain accurate potential of mean force (PMF) for a large set of antibiotics and to predict permeation rates. Therefore, the computational expense becomes an important criterion as well. Simulation studies on the permeation process and similar complex processes have shown that both the sampling scheme employed and the procedure used for the generation of the initial states can critically affect the quality of the estimates obtained and the respective computational overhead. The temperature-accelerated sliced sampling method (TASS) has been shown to partly address the issues with efficient sampling of the important and slow degrees of freedom by enabling simultaneous biasing of a large number of collective variables. In this work, we investigate the effect of the procedure used for the generation of input conformations on the convergence of free-energy estimates obtained from TASS simulations. In particular, we compare the steered molecular dynamics (MD)-based procedure that has been used in previous TASS studies with the Monte Carlo pathway search method, which is used to obtain approximate permeation trajectories with minimum perturbation of the protein channel. We tested different input setups for enrofloxacin permeation through the porins OmpK35 and OmpE35. The best setup shows an improved agreement between independent PMFs in both cases at a much lower computational cost.

Specific Substrate Activity of Lotus Root Polyphenol Oxidase: Insights from Gaussian-Accelerated Molecular Dynamics and Markov State Models

Polyphenol oxidase (PPO) plays a key role in the enzymatic browning process, and this study employed Gaussian-accelerated molecular dynamics (GaMD) simulations to investigate the catalytic efficiency mechanisms of lotus root PPO with different substrates, including catechin, epicatechin, and chlorogenic acid, as well as the inhibitor oxalic acid. Key findings reveal significant conformational changes in PPO that correlate with its enzymatic activity. Upon substrate binding, the alpha-helix in the Q53-D63 region near the copper ion extends, likely stabilizing the active site and enhancing catalysis. In contrast, this helix is disrupted in the presence of the inhibitor, resulting in a decrease in enzymatic efficiency. Additionally, the F350-V378 region, which covers the substrate-binding site, forms an alpha-helix upon substrate binding, further stabilizing the substrate and promoting catalytic function. However, this alpha-helix does not form when the inhibitor is bound, destabilizing the binding site and contributing to inhibition. These findings offer new insights into the substrate-specific and inhibitor-induced structural dynamics of lotus root PPO, providing valuable information for enhancing food processing and preservation techniques.

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.

Interactions between Inhibitors and 5-Lipoxygenase: Insights from Gaussian Accelerated Molecular Dynamics and Markov State Models

Inflammation is a protective stress response triggered by external stimuli, with 5-lipoxygenase (5LOX) playing a pivotal role as a potent mediator of the leukotriene (Lts) inflammatory pathway. Nordihydroguaiaretic acid (NDGA) functions as a natural orthosteric inhibitor of 5LOX, while 3-acetyl-11-keto-β-boswellic acid (AKBA) acts as a natural allosteric inhibitor targeting 5LOX. However, the precise mechanisms of inhibition have remained unclear. In this study, Gaussian accelerated molecular dynamics (GaMD) simulation was employed to elucidate the inhibitory mechanisms of NDGA and AKBA on 5LOX. It was found that the orthosteric inhibitor NDGA was tightly bound in the protein's active pocket, occupying the active site and inhibiting the catalytic activity of the 5LOX enzyme through competitive inhibition. The binding of the allosteric inhibitor AKBA induced significant changes at the distal active site, leading to a conformational shift of residues 168-173 from a loop to an α-helix and significant negative correlated motions between residues 285-290 and 375-400, reducing the distance between these segments. In the simulation, the volume of the active cavity in the stable conformation of the protein was reduced, hindering the substrate's entry into the active cavity and, thereby, inhibiting protein activity through allosteric effects. Ultimately, Markov state models (MSM) were used to identify and classify the metastable states of proteins, revealing the transition times between different conformational states. In summary, this study provides theoretical insights into the inhibition mechanisms of 5LOX by AKBA and NDGA, offering new perspectives for the development of novel inhibitors specifically targeting 5LOX, with potential implications for anti-inflammatory drug development.

Do All Paths Lead to Rome? How Reliable is Umbrella Sampling Along a Single Path?

Molecular dynamics (MD) simulations are widely applied to estimate absolute binding free energies of protein-ligand and protein-protein complexes. A routinely used method for binding free energy calculations with MD is umbrella sampling (US), which calculates the potential of mean force (PMF) along a single reaction coordinate. Surprisingly, in spite of its widespread use, few validation studies have focused on the convergence of the free energy computed along a single path for specific cases, not addressing the reproducibility of such calculations in general. In this work, we therefore investigate the reproducibility and convergence of US along a standard distance-based reaction coordinate for various protein-protein and protein-ligand complexes, following commonly used guidelines for the setup. We show that repeating the complete US workflow can lead to differences of 2-20 kcal/mol in computed binding free energies. We attribute those discrepancies to small differences in the binding pathways. While these differences are unavoidable in the established US protocol, the popularity of the latter could hint at a lack of awareness of such reproducibility problems. To test if the convergence of PMF profiles can be improved if multiple pathways are sampled simultaneously, we performed additional simulations with an adaptive-biasing method, here the accelerated weight histogram (AWH) approach. Indeed, the PMFs obtained from AHW simulations are consistent and reproducible for the systems tested. To the best of our knowledge, our work is the first to attempt a systematic assessment of the pitfalls in one the most widely used protocols for computing binding affinities. We anticipate therefore that our results will provide an incentive for a critical reassessment of the validity of PMFs computed with US, and make a strong case to further benchmark the performance of adaptive-biasing methods for computing binding affinities.

Structure of the Bacterial Cell Envelope and Interactions with Antimicrobials: Insights from Molecular Dynamics Simulations

Bacteria have evolved over 3 billion years, shaping our intrinsic and symbiotic coexistence with these single-celled organisms. With rising populations of drug-resistant strains, the search for novel antimicrobials is an ongoing area of research. Advances in high-performance computing platforms have led to a variety of molecular dynamics simulation strategies to study the interactions of antimicrobial molecules with different compartments of the bacterial cell envelope of both Gram-positive and Gram-negative species. In this review, we begin with a detailed description of the structural aspects of the bacterial cell envelope. Simulations concerned with the transport and associated free energy of small molecules and ions through the outer membrane, peptidoglycan, inner membrane and outer membrane porins are discussed. Since surfactants are widely used as antimicrobials, a section is devoted to the interactions of surfactants with the cell wall and inner membranes. The review ends with a discussion on antimicrobial peptides and the insights gained from the molecular simulations on the free energy of translocation. Challenges involved in developing accurate molecular models and coarse-grained strategies that provide a trade-off between atomic details with a gain in sampling time are highlighted. The need for efficient sampling strategies to obtain accurate free energies of translocation is also discussed. Molecular dynamics simulations have evolved as a powerful tool that can potentially be used to design and develop novel antimicrobials and strategies to effectively treat bacterial infections.

The mechanism of an electrostatic nanofilter: overcoming entropy with electrostatics

General porins are nature's sieving machinery in the outer membrane of Gram-negative bacteria. Their unique hourglass-shaped architecture is highly conserved among different bacterial membrane proteins and other biological channels. These biological nanopores have been designed to protect the interior of the bacterial cell from leakage of toxic compounds while selectively allowing the entry of the molecules needed for cell growth and function. The mechanism of transport through porins is of utmost and direct interest for drug discovery, extending toward nanotechnology applications for blue energy, separations, and sequencing. Here we present a theoretical framework for analysing the filter of general porins in relation to translocating molecules with the aid of enhanced molecular simulations quantitatively. Using different electrostatic probes in the form of a series of related molecules, we describe the nature of this filter and how to finely tune permeability by exploiting electrostatic interactions between the pore and the translocating molecule. Eventually, we show how enhanced simulations constitute today a valid tool for characterising the mechanism and quantifying energetically the transport of molecules through nanopores.

Engineering Biological Nanopore Approaches toward Protein Sequencing

Biotechnological innovations have vastly improved the capacity to perform large-scale protein studies, while the methods we have for identifying and quantifying individual proteins are still inadequate to perform protein sequencing at the single-molecule level. Nanopore-inspired systems devoted to understanding how single molecules behave have been extensively developed for applications in genome sequencing. These nanopore systems are emerging as prominent tools for protein identification, detection, and analysis, suggesting realistic prospects for novel protein sequencing. This review summarizes recent advances in biological nanopore sensors toward protein sequencing, from the identification of individual amino acids to the controlled translocation of peptides and proteins, with attention focused on device and algorithm development and the delineation of molecular mechanisms with the aid of simulations. Specifically, the review aims to offer recommendations for the advancement of nanopore-based protein sequencing from an engineering perspective, highlighting the need for collaborative efforts across multiple disciplines. These efforts should include chemical conjugation, protein engineering, molecular simulation, machine-learning-assisted identification, and electronic device fabrication to enable practical implementation in real-world scenarios.

Converging PMF Calculations of Antibiotic Permeation across an Outer Membrane Porin with Subkilocalorie per Mole Accuracy

The emergence of multidrug-resistant pathogens led to a critical need for new antibiotics. A key property of effective antibiotics against Gram-negative bacteria is their ability to permeate through the bacterial outer membrane via transmembrane porin proteins. Molecular dynamics (MD) simulations are, in principle, capable of modeling antibiotic permeation across outer membrane porins (OMPs). However, owing to sampling problems, it has remained challenging to obtain converged potentials of mean force (PMFs) for antibiotic permeation across OMPs. Here, we investigated the convergence of PMFs along a single collective variable aimed at quantifying the permeation of the antibiotic fosmidomycin across the OprO porin. We compared standard umbrella sampling (US) with three advanced flavors of the US technique: (i) Hamiltonian replica exchange with solute tempering in combination with US, (ii) simulated tempering-enhanced US, and (iii) replica-exchange US. To quantify the PMF convergence and to reveal hysteresis problems, we computed several independent sets of US simulations starting from pulling simulations in the outward and inward permeation directions. We find that replica-exchange US in combination with well-chosen restraints is highly successful for obtaining converged PMFs of fosmidomycin permeation through OprO, reaching PMFs converged to subkilocalorie per mole accuracy.

Ionizable lipids penetrate phospholipid bilayers with high phase transition temperatures: perspectives from free energy calculations

The efficacies of modern gene-therapies strongly depend on their contents. At the same time the most potent formulations might not contain the best compounds. In this work we investigated the effect of phospholipids and their saturation on the binding ability of (6Z,9Z,28Z,31Z)-heptatriacont-6,9,28,31-tetraene-19-yl 4-(dimethylamino) butanoate (DLin-MC3-DMA) to model membranes at the neutral pH. We discovered that DLin-MC3-DMA has affinity to the most saturated monocomponent lipid bilayer 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and an aversion to the unsaturated one 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The preference to a certain membrane was also well-correlated to the phase transition temperatures of phospholipid bilayers, and to their structural and dynamical properties. Additionally, in the case of the presence of DLin-MC3-DMA in the membrane with DOPC the ionizable lipid penetrated it, which indicates possible synergistic effects. Comparisons with other ionizable lipids were performed using a model lipid bilayer of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC). Particularly, the lipids heptadecan-9-yl 8-[2-hydroxyethyl-(6-oxo-6-undecoxyhexyl)amino]octanoate (SM-102) and [(4-Hydroxybutyl) azanediyl] di(hexane-6,1-diyl) bis(2-hexyldecanoate) (ALC-0315) from modern mRNA-vaccines against COVID-19 were investigated and force fields parameters were derived for those new lipids. It was discovered that ALC-0315 binds strongest to the membrane, while DLin-MC3-DMA is not able to reside in the bilayer center. The ability to penetrate the membrane POPC by SM-102 and ALC-0315 can be related to their saturation, comparing to DLin-MC3-DMA.

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.

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.

Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

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 (Q_TC), 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.

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.

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.

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.

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.

Natural Polyphenols Inhibit the Dimerization of the SARS-CoV-2 Main Protease: The Case of Fortunellin and Its Structural Analogs

3CL-Pro is the SARS-CoV-2 main protease (MPro). It acts as a homodimer to cleave the large polyprotein 1ab transcript into proteins that are necessary for viral growth and replication. 3CL-Pro has been one of the most studied SARS-CoV-2 proteins and a main target of therapeutics. A number of drug candidates have been reported, including natural products. Here, we employ elaborate computational methods to explore the dimerization of the 3CL-Pro protein, and we formulate a computational context to identify potential inhibitors of this process. We report that fortunellin (acacetin 7-O-neohesperidoside), a natural flavonoid O-glycoside, and its structural analogs are potent inhibitors of 3CL-Pro dimerization, inhibiting viral plaque formation in vitro. We thus propose a novel basis for the search of pharmaceuticals as well as dietary supplements in the fight against SARS-CoV-2 and COVID-19.

Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators

Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.

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.

Fosmidomycin transport through the phosphate-specific porins OprO and OprP of Pseudomonas aeruginosa

The Gram-negative bacterium Pseudomonas aeruginosa is an opportunistic pathogen, responsible for many hospital-acquired infections. The bacterium is quite resistant toward many antibiotics, in particular because of the fine-tuned permeability of its outer membrane (OM). General diffusion outer membrane pores are quite rare in this organism. Instead, its OM contains many substrate-specific porins. Their expression is varying according to growth conditions and virulence. Phosphate limitations, as well as pathogenicity factors, result in the induction of the two mono- and polyphosphate-specific porins, OprP and OprO, respectively, together with an inner membrane uptake mechanism and a periplasmic binding protein. These outer membrane channels could serve as outer membrane pathways for the uptake of phosphonates. Among them are not only herbicides, but also potent antibiotics, such as fosfomycin and fosmidomycin. In this study, we investigated the interaction between OprP and OprO and fosmidomycin in detail. We could demonstrate that fosmidomycin is able to bind to the phosphate-specific binding site inside the two porins. The inhibition of chloride conductance of OprP and OprO by fosmidomycin is considerably less than that of phosphate or diphosphate, but it can be measured in titration experiments of chloride conductance and also in single-channel experiments. The results suggest that fosmidomycin transport across the OM of P. aeruginosa occurs through OprP and OprO. Our data with the ones already known in the literature show that phosphonic acid-containing antibiotics are in general good candidates to treat the infections of P. aeruginosa at the very beginning through a favorable OM transport system.

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.

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.

The conformational phase space of the photoprotective switch in the major light harvesting complex II

Extensive conformational sampling of the major light harvesting complex II defines the collective variables of the photoprotective switch.