Generating Concentration Gradients across Membranes for Molecular Dynamics Simulations of Periodic Systems.

The plasma membrane forms the boundary between a living entity and its environment and acts as a barrier to permeation and flow of substances. Several computational means of calculating permeability have been implemented for molecular dynamics (MD) simulations-based approaches. Except for double bilayer systems, most permeability studies have been performed under equilibrium conditions, in large part due to the challenges associated with creating concentration gradients in simulations utilizing periodic boundary conditions. To enhance the scientific understanding of permeation and complement the existing computational means of characterizing membrane permeability, we developed a non-equilibrium method that enables the generation and maintenance of steady-state gradients in MD simulations. We utilize PBCs advantageously by imposing a directional bias to the motion of permeants so that their crossing of the boundary replenishes the gradient, like a previous study on ions. Under these conditions, a net flow of permeants across membranes may be observed to determine bulk permeability by a direct application of J=PΔc. In the present study, we explore the results of its application to an exemplary O2 and POPC bilayer system, demonstrating accurate and precise permeability measurements. In addition, we illustrate the impact of permeant concentration and the choice of thermostat on the permeability. Moreover, we demonstrate that energetics of permeation can be closely examined by the dissipation of the gradient across the membrane to gain nuanced insights into the thermodynamics of permeability.

Assembly and Analysis of Cell-Scale Membrane Envelopes.

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms─a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

Carbon dioxide transport across membranes.

Carbon dioxide (CO2) movement across cellular membranes is passive and governed by Fick's law of diffusion. Until recently, we believed that gases cross biological membranes exclusively by dissolving in and then diffusing through membrane lipid. However, the observation that some membranes are CO2 impermeable led to the discovery of a gas molecule moving through a channel; namely, CO2 diffusion through aquaporin-1 (AQP1). Later work demonstrated CO2 diffusion through rhesus (Rh) proteins and NH3 diffusion through both AQPs and Rh proteins. The tetrameric AQPs exhibit differential selectivity for CO2 versus NH3 versus H2O, reflecting physico-chemical differences among the small molecules as well as among the hydrophilic monomeric pores and hydrophobic central pores of various AQPs. Preliminary work suggests that NH3 moves through the monomeric pores of AQP1, whereas CO2 moves through both monomeric and central pores. Initial work on AQP5 indicates that it is possible to create a metal-binding site on the central pore's extracellular face, thereby blocking CO2 movement. The trimeric Rh proteins have monomers with hydrophilic pores surrounding a hydrophobic central pore. Preliminary work on the bacterial Rh homologue AmtB suggests that gas can diffuse through the central pore and three sets of interfacial clefts between monomers. Finally, initial work indicates that CO2 diffuses through the electrogenic Na/HCO3 cotransporter NBCe1. At least in some cells, CO2-permeable proteins could provide important pathways for transmembrane CO2 movements. Such pathways could be amenable to cellular regulation and could become valuable drug targets.

Microscopic view of lipids and their diverse biological functions.

Biological membranes and their diverse lipid constituents play key roles in a broad spectrum of cellular and physiological processes. Characterization of membrane-associated phenomena at a microscopic level is therefore essential to our fundamental understanding of such processes. Due to the semi-fluid and dynamic nature of lipid bilayers, and their complex compositions, detailed characterization of biological membranes at an atomic scale has been refractory to experimental approaches. Computational modeling and simulation offer a highly complementary toolset with sufficient spatial and temporal resolutions to fill this gap. Here, we review recent molecular dynamics studies focusing on the diversity of lipid composition of biological membranes, or aiming at the characterization of lipid-protein interaction, with the overall goal of dissecting how lipids impact biological roles of the cellular membranes.