Brownian Dynamics Approach Including Explicit Atoms for Studying Ion Permeation and Substrate Translocation across Nanopores.

A Brownian dynamics (BD) approach including explicit atoms called BRODEA is presented to model ion permeation and molecule translocation across a nanopore confinement. This approach generalizes our previous hybrid molecular dynamics-Brownian dynamics framework ( J. Chem. Theory Comput. 2016, 12, 2401) by incorporating a widespread and enhanced set of simulation schemes based on several boundary conditions and electrostatic models, as well as a temperature accelerated method for sampling free energy surfaces and determining substrate translocation pathways. As a test case, BRODEA was applied to study the ion diffusion as well as to ciprofloxacin and enrofloxacin transport through the outer membrane porin OmpC from E. coli. The equivalence between the different simulation schemes was demonstrated and their computational efficiency evaluated. The BRODEA results are able to reproduce the main features of the ion currents and free energy surfaces determined by all-atom molecular dynamics simulations and validated by experiments. Furthermore, the BRODEA results are able to determine the ciprofloxacin and enrofloxacin permeation pathways showing a remarkable agreement with the results obtained from a computational protocol that combines metadynamics and a zero-temperature string method ( J. Chem. Theory Comput. 2017, 13, 4553; J. Phys. Chem. B 2018, 122, 1417). To our knowledge, this is the first time such antibiotic permeation pathways have been characterized by a technique based on Brownian dynamics.

BROMOCEA Code: An Improved Grand Canonical Monte Carlo/Brownian Dynamics Algorithm Including Explicit Atoms.

All-atom molecular dynamics simulations have a long history of applications studying ion and substrate permeation across biological and artificial pores. While offering unprecedented insights into the underpinning transport processes, MD simulations are limited in time-scales and ability to simulate physiological membrane potentials or asymmetric salt solutions and require substantial computational power. While several approaches to circumvent all of these limitations were developed, Brownian dynamics simulations remain an attractive option to the field. The main limitation, however, is an apparent lack of protein flexibility important for the accurate description of permeation events. In the present contribution, we report an extension of the Brownian dynamics scheme which includes conformational dynamics. To achieve this goal, the dynamics of amino-acid residues was incorporated into the many-body potential of mean force and into the Langevin equations of motion. The developed software solution, called BROMOCEA, was applied to ion transport through OmpC as a test case. Compared to fully atomistic simulations, the results show a clear improvement in the ratio of permeating anions and cations. The present tests strongly indicate that pore flexibility can enhance permeation properties which will become even more important in future applications to substrate translocation.

Simulations of outer membrane channels and their permeability.

Channels in the outer membrane of Gram-negative bacteria provide essential pathways for the controlled and unidirectional transport of ions, nutrients and metabolites into the cell. At the same time the outer membrane serves as a physical barrier for the penetration of noxious substances such as antibiotics into the bacteria. Most antibiotics have to pass through these membrane channels to either reach cytoplasmic bound targets or to further cross the hydrophobic inner membrane. Considering the pharmaceutical significance of antibiotics, understanding the functional role and mechanism of these channels is of fundamental importance in developing strategies to design new drugs with enhanced permeation abilities. Due to the biological complexity of membrane channels and experimental limitations, computer simulations have proven to be a powerful tool to investigate the structure, dynamics and interactions of membrane channels. Considerable progress has been made in computer simulations of membrane channels during the last decade. The goal of this review is to provide an overview of the computational techniques and their roles in modeling the transport across outer membrane channels. A special emphasis is put on all-atom molecular dynamics simulations employed to better understand the transport of molecules. Moreover, recent molecular simulations of ion, substrate and antibiotics translocation through membrane pores are briefly summarized. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

BROMOC-D: Brownian Dynamics/Monte-Carlo Program Suite to Study Ion and DNA Permeation in Nanopores.

A theoretical framework is presented to model ion and DNA translocation across a nanopore confinement under an applied electric field. A combined Grand Canonical Monte Carlo Brownian Dynamics (GCMC/BD) algorithm offers a general approach to study ion permeation through wide molecular pores with a direct account of ion-ion and ion-DNA correlations. This work extends previously developed theory by incorporating the recently developed coarse-grain polymer model of DNA by de Pablo and colleagues [Knotts, T. A.; Rathore, N.; Schwartz, D. C.; de Pablo, J. J. J. Chem. Phys. 2007, 126] with explicit ions for simulations of polymer dynamics. Atomistic MD simulations were used to guide model developments. The power of the developed scheme is illustrated with studies of single-stranded DNA (ss-DNA) oligomer translocation in two model cases: a cylindrical pore with a varying radius and a well-studied experimental system, the staphylococcal α-hemolysin channel. The developed model shows good agreement with experimental data for model studies of two homopolymers: ss-poly(dA)(n) and ss-poly(dC)(n). The developed protocol allows for direct evaluation of different factors (charge distribution and pore shape and size) controlling DNA translocation in a variety of nanopores.