Ions and counterions in a biological channel: a molecular dynamics simulation of OmpF porin from Escherichia coli in an explicit membrane with 1 M KCl aqueous salt solution

Orphan receptor GPR88 as a potential therapeutic target for CNS disorders - an in silico approach

The G-protein-coupled receptors are a part of the largest and most physiologically relevant family of membrane proteins. One-third of the medications, now on the market, target the GPCR receptor family, which is one of the most important therapeutic targets for many disorders. In the reported work, we have focussed on orphan GPR88 receptor which is a part of the GPCR protein family and a potential target for central nervous system disorders. GPR88 is known to show the highest expression in the striatum, which is a key region in motor control and cognitive functions. Recent studies have reported that GPR88 is activated by two agonists, 2-PCCA and RTI-13951-33. In this study, we have predicted the three-dimensional protein structure for the orphan GPR88 using the homology modeling approach. We then used shape-based screening techniques based on known agonists and structure-based virtual screening methods employing docking to uncover novel GPR88 ligands. The screened GPR88-ligand complexes were further subjected to molecular dynamics simulation studies. The selected ligands could fasten the development of novel treatments for the vast list of movement and central nervous system disorders.Communicated by Ramaswamy H. Sarma.

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.

Structural basis of pH-dependent activation in a CLC transporter

CLCs are dimeric chloride channels and anion/proton exchangers that regulate processes such as muscle contraction and endo-lysosome acidification. Common gating controls their activity; its closure simultaneously silences both protomers, and its opening allows them to independently transport ions. Mutations affecting common gating in human CLCs cause dominant genetic disorders. The structural rearrangements underlying common gating are unknown. Here, using single-particle cryo-electron microscopy, we show that the prototypical Escherichia coli CLC-ec1 undergoes large-scale rearrangements in activating conditions. The slow, pH-dependent remodeling of the dimer interface leads to the concerted opening of the intracellular H+ pathways and is required for transport. The more frequent formation of short water wires in the open H+ pathway enables Cl- pore openings. Mutations at disease-causing sites favor CLC-ec1 activation and accelerate common gate opening in the human CLC-7 exchanger. We suggest that the pH activation mechanism of CLC-ec1 is related to the common gating of CLC-7.

How is Membrane Permeation of Small Ionizable Molecules Affected by Protonation Kinetics?

According to the pH-partition hypothesis, the aqueous solution adjacent to a membrane is a mixture of the ionization states of the permeating molecule at fixed Henderson-Hasselbalch concentrations, such that each state passes through the membrane in parallel with its own specific permeability. An alternative view, based on the assumption that the rate of switching ionization states is instantaneous, represents the permeation of ionizable molecules via an effective Boltzmann-weighted average potential (BWAP). Such an assumption is used in constant-pH molecular dynamics simulations. The inhomogeneous solubility-diffusion framework can be used to compute the pH-dependent membrane permeability for each of these two limiting treatments. With biased WTM-eABF molecular dynamics simulations, we computed the potential of mean force and diffusivity of each ionization state of two weakly basic small molecules: nicotine, an addictive drug, and varenicline, a therapeutic for treating nicotine addiction. At pH = 7, the BWAP effective permeability is greater than that determined by pH-partitioning by a factor of 2.5 for nicotine and 5 for varenicline. To assess the importance of ionization kinetics, we present a Smoluchowski master equation that includes explicitly the protonation and deprotonation processes coupled with the diffusive motion across the membrane. At pH = 7, the increase in permeability due to the explicit ionization kinetics is negligible for both nicotine and varenicline. This finding is reaffirmed by combined Brownian dynamics and Markov state model simulations for estimating the permeability of nicotine while allowing changes in its ionization state. We conclude that for these molecules the pH-partition hypothesis correctly captures the physics of the permeation process. The small free energy barriers for the permeation of nicotine and varenicline in their deprotonated neutral forms play a crucial role in establishing the validity of the pH-partitioning mechanism. Essentially, BWAP fails because ionization kinetics are too slow on the time scale of membrane crossing to affect the permeation of small ionizable molecules such as nicotine and varenicline. For the singly protonated state of nicotine, the computational results agree well with experimental measurements (P1 = 1.29 × 10-7 cm/s), but the agreement for neutral (P0 = 6.12 cm/s) and doubly protonated nicotine (P2 = 3.70 × 10-13 cm/s) is slightly worse, likely due to factors associated with the aqueous boundary layer (neutral form) or leaks through paracellular pathways (doubly protonated form).

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

Influence of Membrane Asymmetry on OmpF Insertion, Orientation and Function

The effect of asymmetric membranes containing lipopolysaccharides (LPS) on the outer membrane protein F (OmpF) reconstitution, channel orientation, and antibiotic permeation across the outer membrane was investigated. After forming an asymmetric planar lipid bilayer composed of LPS on one and phospholipids on the other side, the membrane channel OmpF was added. The ion current recordings demonstrate that LPS has a strong influence on the OmpF membrane insertion, orientation, and gating. Enrofloxacin was used as an example of an antibiotic interacting with the asymmetric membrane and with OmpF. The enrofloxacin caused the blockage of the ion current through the OmpF, depending on the side of addition, the transmembrane voltage applied, and the composition of the buffer. Furthermore, the enrofloxacin changed the phase behavior of the LPS-containing membranes, demonstrating that its membrane activity influences the function of OmpF and potentially the membrane permeability.

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.

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.

Charged Nanochannels in Covalent Organic Framework Membranes Enabling Efficient Ion Exclusion

Controllable ion transport through nanochannels is crucial for biological and artificial membrane systems. Covalent organic frameworks (COFs) with regular and tunable nanochannels are emerging as an ideal material platform to develop synthetic membranes for ion transport. However, ion exclusion by COF membranes remains challenging because most COF materials have large-sized nanochannels leading to nonselective transport of small ions. Here we develop ionic COF membranes (iCOFMs) to control ion transport through charged framework nanochannels, the interior surfaces of which are covered with arrayed sulfonate groups to render superior charge density. The overlap of an electrical double layer in charged nanochannels blocks the entry of co-ions, narrows their passageways, and concomitantly restrains the permeation of counterions via the charge balance. These highly charged large-sized nanochannels within the iCOFM enable ion exclusion while maintaining intrinsically high water permeability. Our results reveal possibilities for controllable ion transport based on COF membranes for water purification, ionic separation, sensing, and energy conversion.

Role of internal loop dynamics in antibiotic permeability of outer membrane porins

Antibiotic resistance in Gram-negative pathogens has been identified as an urgent threat to human health by the World Health Organization. The major challenge with treating infections by these pathogens is developing antibiotics that can traverse the dense bacterial outer membrane (OM) formed by a mesh of lipopolysaccharides. Effective antibiotics permeate through OM porins, which have evolved for nutrient diffusion; however, the conformational states of these porins regulating permeation are still unclear. Here, we used molecular dynamics simulations, free energy calculations, Markov-state modeling, and whole-cell accumulation assays to provide mechanistic insight on how a porin shifts between open and closed states. We provide a mechanism of how Gram-negative bacteria confer resistance to antibiotics.

Gram-negative bacteria pose a serious public health concern due to resistance to many antibiotics, caused by the low permeability of their outer membrane (OM). Effective antibiotics use porins in the OM to reach the interior of the cell; thus, understanding permeation properties of OM porins is instrumental to rationally develop broad-spectrum antibiotics. A functionally important feature of OM porins is undergoing open–closed transitions that modulate their transport properties. To characterize the molecular basis of these transitions, we performed an extensive set of molecular dynamics (MD) simulations of Escherichia coli OM porin OmpF. Markov-state analysis revealed that large-scale motion of an internal loop, L3, underlies the transition between energetically stable open and closed states. The conformation of L3 is controlled by H bonds between highly conserved acidic residues on the loop and basic residues on the OmpF β-barrel. Mutation of key residues important for the loop’s conformation shifts the equilibrium between open and closed states and regulates translocation of permeants (ions and antibiotics), as observed in the simulations and validated by our whole-cell accumulation assay. Notably, one mutant system G119D, which we find to favor the closed state, has been reported in clinically resistant bacterial strains. Overall, our accumulated ∼200 µs of simulation data (the wild type and mutants) along with experimental assays suggest the involvement of internal loop dynamics in permeability of OM porins and antibiotic resistance in Gram-negative bacteria.

Rationalizing the generation of broad spectrum antibiotics with the addition of a positive charge

Antibiotic resistance of Gram-negative bacteria is largely attributed to the low permeability of their outer membrane (OM). Recently, we disclosed the eNTRy rules, a key lesson of which is that the introduction of a primary amine enhances OM permeation in certain contexts. To understand the molecular basis for this finding, we perform an extensive set of molecular dynamics (MD) simulations and free energy calculations comparing the permeation of aminated and amine-free antibiotic derivatives through the most abundant OM porin of E. coli, OmpF. To improve sampling of conformationally flexible drugs in MD simulations, we developed a novel, Monte Carlo and graph theory based algorithm to probe more efficiently the rotational and translational degrees of freedom visited during the permeation of the antibiotic molecule through OmpF. The resulting pathways were then used for free-energy calculations, revealing a lower barrier against the permeation of the aminated compound, substantiating its greater OM permeability. Further analysis revealed that the amine facilitates permeation by enabling the antibiotic to align its dipole to the luminal electric field of the porin and form favorable electrostatic interactions with specific, highly-conserved charged residues. The importance of these interactions in permeation was further validated with experimental mutagenesis and whole cell accumulation assays. Overall, this study provides insights on the importance of the primary amine for antibiotic permeation into Gram-negative pathogens that could help the design of future antibiotics. We also offer a new computational approach for calculating free-energy of processes where relevant molecular conformations cannot be efficiently captured.

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.

Specific adsorption of trivalent cations in biological nanopores determines conductance dynamics and reverses ionic selectivity

Adsorption processes are central to ionic transport in industrial and biological membrane systems. Multivalent cations modulate the conductive properties of nanofluidic devices through interactions with charged surfaces that depend principally on the ion charge number. Considering that ion channels are specialized valves that demand a sharp specificity in ion discrimination, we investigate the adsorption dynamics of trace amounts of different salts of trivalent cations in biological nanopores. We consider here OmpF from Escherichia coli, an archetypical protein nanopore, to probe the specificity of biological nanopores to multivalent cations. We systematically compare the effect of three trivalent electrolytes on OmpF current-voltage relationships and characterize the degree of rectification induced by each ion. We also analyze the open channel current noise to determine the existence of equilibrium/non-equilibrium mechanisms of ion adsorption and evaluate the extent of charge inversion through selectivity measurements. We show that the interaction of trivalent electrolytes with biological nanopores occurs via ion-specific adsorption yielding differential modulation of ion conduction and selectivity inversion. We also demonstrate the existence of non-equilibrium fluctuations likely related to ion-dependent trapping-detrapping processes. Our study provides fundamental information relevant to different biological and electrochemical systems where transport phenomena involve ion adsorption in charged surfaces under nanoscale confinement.

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.

Antibiotic uptake through porins located in the outer membrane of Gram-negative bacteria

Introduction: Making selective inhibitors of novel Gram-negative targets is not a substantial challenge - getting them into Gram-negative bacteria to reach their lethal target is the bottleneck. Poor permeability of the antibiotic requires high concentration causing off target activity. The lack of simple experimental techniques to measure antibiotic uptake as well as the local concentration at the target site creates a particular bottleneck in understanding and in improving the antibiotic activity.Areas covered: Here we recall current approaches to quantify the uptake. For a few antibiotics with known evidence for channel-limited permeation, the flux across a single OmpF or OmpC channel has been measured. For a typical concentration gradient of 1 µM of antibiotics the uptake varies between one up to few hundred molecules per second and per channel.Expert opinion: The current research effort is on quantifying the flux for a larger list of compounds on a cellular (mass spectra, fluorescence) or at single channel level (electrophysiology). A larger dataset of single channel permeabilities under various condition will be a powerful tool for understanding and improving the activity of antibiotics.

A constricted opening in Kir channels does not impede potassium conduction

The canonical mechanistic model explaining potassium channel gating is of a conformational change that alternately dilates and constricts a collar-like intracellular entrance to the pore. It is based on the premise that K⁺ ions maintain a complete hydration shell while passing between the transmembrane cavity and cytosol, which must be accommodated. To put the canonical model to the test, we locked the conformation of a Kir K⁺ channel to prevent widening of the narrow collar. Unexpectedly, conduction was unimpaired in the locked channels. In parallel, we employed all-atom molecular dynamics to simulate K⁺ ions moving along the conduction pathway between the lower cavity and cytosol. During simulations, the constriction did not significantly widen. Instead, transient loss of some water molecules facilitated K⁺ permeation through the collar. The low free energy barrier to partial dehydration in the absence of conformational change indicates Kir channels are not gated by the canonical mechanism.

Manipulation of charge distribution in the arginine and glutamate clusters of the OmpG pore alters sugar specificity and ion selectivity

OmpG is a general diffusion pore in the E. coli outer membrane with a molecular architecture comprising a 14-stranded β-barrel scaffold and unique structural features. In contrast to other non-specific porins, OmpG lacks a central constriction zone and has an exceptionally wide pore diameter of about 13 Å. The equatorial plane of OmpG harbors an annulus of four alternating basic and acidic patches whose function is only poorly characterized. We have investigated the role of charge distribution for ion selectivity and sugar transport with the help of OmpG variants mutated in the annulus. Substituting the glutamate residues of the annulus for histidines or alanines led to a strong reduction in cation selectivity. Replacement of the glutamates in the annulus by histidine residues also disfavored the passage of pentoses and hexoses relative to disaccharides. Our results demonstrate that despite the wide pore diameter, an annulus only consisting of two opposing basic patches confers reduced cation and monosaccharide transport compared to OmpG wild type. Furthermore, randomization of charged residues in the annulus had the potential to abolish pH-dependency of sugar transport. Our results indicate that E¹⁵, E³¹, R⁹², R¹¹¹ and R²¹¹ in the annulus form electrostatic interactions with R²²⁸, E²²⁹ and D²³² in loop L6 that influence pH-dependency of sugar transport.

Essential and Detrimental - an Update on Intracellular Iron Trafficking and Homeostasis

Chloroplasts, mitochondria and vacuoles represent characteristic organelles of the plant cell, with a predominant function in cellular metabolism. Chloroplasts are the site of photosynthesis and therefore basic and essential for photoautotrophic growth of plants. Mitochondria produce energy during respiration and vacuoles act as internal waste and storage compartments. Moreover, chloroplasts and mitochondria are sites for the biosynthesis of various compounds of primary and secondary metabolism. For photosynthesis and energy generation, the internal membranes of chloroplasts and mitochondria are equipped with electron transport chains. To perform proper electron transfer and several biosynthetic functions, both organelles contain transition metals and here iron is by far the most abundant. Although iron is thus essential for plant growth and development, it becomes toxic when present in excess and/or in its free, ionic form. The harmful effect of the latter is caused by the generation of oxidative stress. As a consequence, iron transport and homeostasis have to be tightly controlled during plant growth and development. In addition to the corresponding transport and homeostasis proteins, the vacuole plays an important role as an intracellular iron storage and release compartment at certain developmental stages. In this review, we will summarize current knowledge on iron transport and homeostasis in chloroplasts, mitochondria and vacuoles. In addition, we aim to integrate the physiological impact of intracellular iron homeostasis on cellular and developmental processes.

Molecular Dynamics Simulations of Membrane Permeability

This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.

Rapid and Accurate Determination of Nanopore Ionic Current Using a Steric Exclusion Model

Nanopore sensing has emerged as a versatile approach to detection and identification of biomolecules. Presently, researchers rely on experience and intuition for choosing or modifying the nanopores to detect a target analyte. The field would greatly benefit from a computational method that could relate the atomic-scale geometry of the nanopores and analytes to the blockade nanopore currents they produce. Existing computational methods are either computationally too expensive to be used routinely in experimental laboratories or not sensitive enough to account for the atomic structure of the pore and the analytes. Here, we demonstrate a robust and inexpensive computational approach-the steric exclusion model (SEM) of nanopore conductance-that is orders of magnitude more efficient than all-atom MD and yet is sensitive enough to account for the atomic structure of the nanopore and the analyte. The method combines the computational efficiency of a finite element solver with the atomic precision of a nanopore conductance map to yield unprecedented speed and accuracy of ionic current prediction. We validate our SEM approach through comparison with the current blockades computed using the all-atom molecular dynamics method for a range of proteins confined to a solid-state nanopore, biological channels embedded in a lipid bilayer membranes, and blockade currents produced by DNA homopolymers in MspA. We illustrate potential applications of SEM by computing blockade currents produced by nucleosome proteins in a solid-state nanopore, individual amino acids in MspA, and by testing the effect of point mutations on amino acid distinguishability. We expect our SEM approach to become an integral part of future development of the nanopore sensing field.

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.

Outer Membrane Porins

The transport of small molecules across membranes is essential for the import of nutrients and other energy sources into the cell and, for the export of waste and other potentially harmful byproducts out of the cell. While hydrophobic molecules are permeable to membranes, ions and other small polar molecules require transport via specialized membrane transport proteins . The two major classes of membrane transport proteins are transporters and channels. With our focus here on porins-major class of non-specific diffusion channel proteins , we will highlight some recent structural biology reports and functional assays that have substantially contributed to our understanding of the mechanism that mediates uptake of small molecules, including antibiotics, across the outer membrane of Enterobacteriaceae . We will also review advances in the regulation of porin expression and porin biogenesis and discuss these pathways as new therapeutic targets.

Molecular Dynamics Simulation for All

The impact of molecular dynamics (MD) simulations in molecular biology and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomolecules in full atomic detail and at very fine temporal resolution. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of experimental structural data, have increased the appeal of biomolecular simulation to experimentalists-a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomolecules, in uncovering the structural basis for disease, and in the design and optimization of small molecules, peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further experimental work.

Nutrient transport suggests an evolutionary basis for charged archaeal surface layer proteins

Surface layers (S-layers) are two-dimensional, proteinaceous, porous lattices that form the outermost cell envelope component of virtually all archaea and many bacteria. Despite exceptional sequence diversity, S-layer proteins (SLPs) share important characteristics such as their ability to form crystalline sheets punctuated with nano-scale pores, and their propensity for charged amino acids, leading to acidic or basic isoelectric points. However, the precise function of S-layers, or the role of charged SLPs and how they relate to cellular metabolism is unknown. Nano-scale lattices affect the diffusion behavior of low-concentration solutes, even if they are significantly smaller than the pore size. Here, we offer a rationale for charged S-layer proteins in the context of the structural evolution of S-layers. Using the ammonia-oxidizing archaea (AOA) as a model for S-layer geometry, and a 2D electrodiffusion reaction computational framework to simulate diffusion and consumption of the charged solute ammonium (NH₄⁺), we find that the characteristic length scales of nanoporous S-layers elevate the concentration of NH₄⁺ in the pseudo-periplasmic space. Our simulations suggest an evolutionary, mechanistic basis for S-layer charge and shed light on the unique ability of some AOA to oxidize ammonia in environments with nanomolar NH₄⁺ availability, with broad implications for comparisons of ecologically distinct populations.

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.

Ampicillin permeation across OmpF, the major outer-membrane channel in Escherichia coli

The outer cell wall of the Gram-negative bacteria is a crucial barrier for antibiotics to reach their target. Here, we show that the chemical stability of the widely used antibiotic ampicillin is a major factor in the permeation across OmpF to reach the target in the periplasm. Using planar lipid bilayers we investigated the interactions and permeation of OmpF with ampicillin, its basic pH-induced primary degradation product (penicilloic acid), and the chemically more stable benzylpenicillin. We found that the solute-induced ion current fluctuation is 10 times higher with penicilloic acid than with ampicillin. Furthermore, we also found that ampicillin can easily permeate through OmpF, at an ampicillin gradient of 10 μm and a conductance of G_amp ≅ 3.8 fS, with a flux rate of roughly 237 molecules/s of ampicillin at V_m = 10 mV. The structurally related benzylpenicillin yields a lower conductance of G_amp ≅ 2 fS, corresponding to a flux rate of ≈120 molecules/s. In contrast, the similar sized penicilloic acid was nearly unable to permeate through OmpF. MD calculations show that, besides their charge difference, the main differences between ampicillin and penicilloic acid are the shape of the molecules, and the strength and direction of the dipole vector. Our results show that OmpF can impose selective permeation on similar sized molecules based on their structure and their dipolar properties.

Communication between the leaflets of asymmetric membranes revealed from coarse-grain molecular dynamics simulations

We use coarse-grain molecular simulations to investigate the structural and dynamics differences between an asymmetric and a symmetrical membrane, both containing beta barrel transmembrane proteins. We find in where the dynamics of the two leaflets differ greatly, the slowest leaflet dominates the structural effects and importance of protein-lipid interactions.

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.

Multiscale modeling of a rectifying bipolar nanopore: explicit-water versus implicit-water simulations

In a multiscale modeling approach, we present computer simulation results for a rectifying bipolar nanopore at two modeling levels. In an all-atom model, we use explicit water to simulate ion transport directly with the molecular dynamics technique. In a reduced model, we use implicit water and apply the Local Equilibrium Monte Carlo method together with the Nernst-Planck transport equation. This hybrid method makes the fast calculation of ion transport possible at the price of lost details. We show that the implicit-water model is an appropriate representation of the explicit-water model when we look at the system at the device (i.e., input vs. output) level. The two models produce qualitatively similar behavior of the electrical current for different voltages and model parameters. Looking at the details of concentration and potential profiles, we find profound differences between the two models. These differences, however, do not influence the basic behavior of the model as a device because they do not influence the z-dependence of the concentration profiles which are the main determinants of current. These results then address an old paradox: how do reduced models, whose assumptions should break down in a nanoscale device, predict experimental data? Our simulations show that reduced models can still capture the overall device physics correctly, even though they get some important aspects of the molecular-scale physics quite wrong; reduced models work because they include the physics that is necessary from the point of view of device function. Therefore, reduced models can suffice for general device understanding and device design, but more detailed models might be needed for molecular level understanding.

Refinement of OprH-LPS Interactions by Molecular Simulations

The outer membrane (OM) of Gram-negative bacteria is composed of lipopolysaccharide (LPS) in the outer leaflet and phospholipids in the inner leaflet. The outer membrane protein H (OprH) of Pseudomonas aeruginosa provides an increased stability to the OMs by directly interacting with LPS. Here we report the influence of various P. aeruginosa and, for comparison, Escherichia coli LPS environments on the physical properties of the OMs and OprH using all-atom molecular dynamics simulations. The simulations reveal that although the P. aeruginosa OMs are thinner hydrophobic bilayers than the E. coli OMs, which is expected from the difference in the acyl chain length of their lipid A, this effect is almost imperceptible around OprH due to a dynamically adjusted hydrophobic match between OprH and the OM. The structure and dynamics of the extracellular loops of OprH show distinct behaviors in different LPS environments. Including the O-antigen greatly reduces the flexibility of the OprH loops and increases the interactions between these loops and LPS. Furthermore, our study shows that the interactions between OprH and LPS mainly depend on the secondary structure of OprH and the chemical structure of LPS, resulting in distinctive patterns in different LPS environments.

CHARMM-GUI 10 years for biomolecular modeling and simulation

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

Predictive compound accumulation rules yield a broad-spectrum antibiotic

Most small molecules are unable to rapidly traverse the outer membrane of Gram-negative bacteria and accumulate inside these cells, making the discovery of much-needed drugs against these pathogens challenging. Current understanding of the physicochemical properties that dictate small-molecule accumulation in Gram-negative bacteria is largely based on retrospective analyses of antibacterial agents, which suggest that polarity and molecular weight are key factors. Here we assess the ability of over 180 diverse compounds to accumulate in Escherichia coli. Computational analysis of the results reveals major differences from the retrospective studies, namely that the small molecules that are most likely to accumulate contain an amine, are amphiphilic and rigid, and have low globularity. These guidelines were then applied to convert deoxynybomycin, a natural product that is active only against Gram-positive organisms, into an antibiotic with activity against a diverse panel of multi-drug-resistant Gram-negative pathogens. We anticipate that these findings will aid in the discovery and development of antibiotics against Gram-negative bacteria.

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.

Charged Substrate and Product Together Contribute Like a Nonreactive Species to the Overall Electrostatic Steering in Diffusion-Reaction Processes

The Debye-Hückel limiting law is used to study the binding kinetics of substrate-enzyme system as well as to estimate the reaction rate of a electrostatically steered diffusion-controlled reaction process. It is based on a linearized Poisson-Boltzmann model and known for its accurate predictions in dilute solutions. However, the substrate and product particles are in nonequilibrium states and are possibly charged, and their contributions to the total electrostatic field cannot be explicitly studied in the Poisson-Boltzmann model. Hence the influences of substrate and product on reaction rate coefficient were not known. In this work, we consider all the charged species, including the charged substrate, product, and mobile salt ions in a Poisson-Nernst-Planck model, and then compare the results with previous work. The results indicate that both the charged substrate and product can significantly influence the reaction rate coefficient with different behaviors under different setups of computational conditions. It is interesting to find that when substrate and product are both considered, under an overall neutral boundary condition for all the bulk charged species, the computed reaction rate kinetics recovers a similar Debye-Hückel limiting law again. This phenomenon implies that the charged product counteracts the influence of charged substrate on reaction rate coefficient. Our analysis discloses the fact that the total charge concentration of substrate and product, though in a nonequilibrium state individually, obeys an equilibrium Boltzmann distribution, and therefore contributes as a normal charged ion species to ionic strength. This explains why the Debye-Hückel limiting law still works in a considerable range of conditions even though the effects of charged substrate and product particles are not specifically and explicitly considered in the theory.

Challenges in structural approaches to cell modeling

Computational modeling is essential for structural characterization of biomolecular mechanisms across the broad spectrum of scales. Adequate understanding of biomolecular mechanisms inherently involves our ability to model them. Structural modeling of individual biomolecules and their interactions has been rapidly progressing. However, in terms of the broader picture, the focus is shifting toward larger systems, up to the level of a cell. Such modeling involves a more dynamic and realistic representation of the interactomes in vivo, in a crowded cellular environment, as well as membranes and membrane proteins, and other cellular components. Structural modeling of a cell complements computational approaches to cellular mechanisms based on differential equations, graph models, and other techniques to model biological networks, imaging data, etc. Structural modeling along with other computational and experimental approaches will provide a fundamental understanding of life at the molecular level and lead to important applications to biology and medicine. A cross section of diverse approaches presented in this review illustrates the developing shift from the structural modeling of individual molecules to that of cell biology. Studies in several related areas are covered: biological networks; automated construction of three-dimensional cell models using experimental data; modeling of protein complexes; prediction of non-specific and transient protein interactions; thermodynamic and kinetic effects of crowding; cellular membrane modeling; and modeling of chromosomes. The review presents an expert opinion on the current state-of-the-art in these various aspects of structural modeling in cellular biology, and the prospects of future developments in this emerging field.

How Tolerant are Membrane Simulations with Mismatch in Area per Lipid between Leaflets?

Difficulties in estimating the correct number of lipids in each leaflet of complex bilayer membrane simulation systems make it inevitable to introduce a mismatch in lipid packing (i.e., area per lipid) and thus alter the lateral pressure of each leaflet. To investigate potential impacts of such mismatch on simulation results, we performed molecular dynamics simulations of saturated and monounsaturated lipid bilayers with and without gramicidin A or WALP23 at various mismatches by adjusting the number of lipids in the lower leaflet from no mismatch to a 25% reduction compared to that in the upper leaflet. All simulations were stable under the constant pressure barostat, but the mismatch induces asymmetric lipid packing between the leaflets, so that the upper leaflet becomes more ordered, and the lower leaflet becomes less ordered. The mismatch impacts on various bilayer properties are mild up to 5–10% mismatch, and bilayers with fully saturated chains appear to be more prone to these impacts than those with unsaturated tails. The nonvanishing leaflet surface tensions and the free energy derivatives with respect to the bilayer curvature indicate that the bilayer would be energetically unstable in the presence of mismatch. We propose a quantitative criterion for allowable mismatch based on the energetics derived from a continuum elastic model, which grows as a square root of the number of the lipids in the system. On the basis of this criterion, we infer that the area per lipid mismatch up to 5% would be tolerable in various membrane simulations of reasonable all-atom system sizes (40–160 lipids per leaflet).

Living on the edge: Simulations of bacterial outer-membrane proteins

Gram-negative bacteria are distinguished in part by a second, outer membrane surrounding them. This membrane is distinct from others, possessing an outer leaflet composed not of typical phospholipids but rather large, highly charged molecules known as lipopolysaccharides. Therefore, modeling the structure and dynamics of proteins embedded in the outer membrane requires careful consideration of their native environment. In this review, we examine how simulations of such outer-membrane proteins have evolved over the last two decades, culminating most recently in detailed, highly accurate atomistic models of the outer membrane. We also draw attention to how the simulations have coupled with experiments to produce novel insights unattainable through a single approach. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.