Chapter 12 Gas Conduction of Lipid Bilayers and Membrane Channels

Elucidating the complex membrane binding of a protein with multiple anchoring domains using extHMMM

Membrane binding is a crucial mechanism for many proteins, but understanding the specific interactions between proteins and membranes remains a challenging endeavor. Coagulation factor Va (FVa) is a large protein whose membrane interactions are complicated due to the presence of multiple anchoring domains that individually can bind to lipid membranes. Using molecular dynamics simulations, we investigate the membrane binding of FVa and identify the key mechanisms that govern its interaction with membranes. Our results reveal that FVa can either adopt an upright or a tilted molecular orientation upon membrane binding. We further find that the domain organization of FVa deviates (sometimes significantly) from its crystallographic reference structure, and that the molecular orientation of the protein matches with domain reorganization to align the C2 domain toward its favored membrane-normal orientation. We identify specific amino acid residues that exhibit contact preference with phosphatidylserine lipids over phosphatidylcholine lipids, and we observe that mostly electrostatic effects contribute to this preference. The observed lipid-binding process and characteristics, specific to FVa or common among other membrane proteins, in concert with domain reorganization and molecular tilt, elucidate the complex membrane binding dynamics of FVa and provide important insights into the molecular mechanisms of protein-membrane interactions. An updated version of the HMMM model, termed extHMMM, is successfully employed for efficiently observing membrane bindings of systems containing the whole FVa molecule.

Membrane Permeability: Characteristic Times and Lengths for Oxygen and a Simulation-Based Test of the Inhomogeneous Solubility-Diffusion Model

The balance of normal and radial (lateral) diffusion of oxygen in phospholipid membranes is critical for biological function. Based on the Smoluchowski equation for the inhomogeneous solubility-diffusion model, Bayesian analysis (BA) can be applied to molecular dynamics trajectories of oxygen to extract the free energy and the normal and radial diffusion profiles. This paper derives a theoretical formalism to convert these profiles into characteristic times and lengths associated with entering, escaping, or completely crossing the membrane. The formalism computes mean first passage times and holds for any process described by rate equations between discrete states. BA of simulations of eight model membranes with varying lipid composition and temperature indicate that oxygen travels 3 to 5 times further in the radial than in the normal direction when crossing the membrane in a time of 15 to 32 ns, thereby confirming the anisotropy of passive oxygen transport in membranes. Moreover, the preceding times and distances estimated from the BA are compared to the aggregate of 280 membrane exits explicitly observed in the trajectories. BA predictions for the distances of oxygen radial diffusion within the membrane are statistically indistinguishable from the corresponding simulation values, yet BA oxygen exit times from the membrane interior are approximately 20% shorter than the simulation values, averaged over seven systems. The comparison supports the BA approach and, therefore, the applicability of the Smoluchowski equation to membrane diffusion. Given the shorter trajectories required for the BA, these results validate the BA as a computationally attractive alternative to direct observation of exits when estimating characteristic times and radial distances. The effect of collective membrane undulations on the BA is also discussed.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

A liposome-encapsulated spin trap for the detection of nitric oxide

Electron paramagnetic resonance (EPR) is one of the few methods that allows for the unambiguous detection of nitric oxide (NO). However, the dithiocarbamate-iron spin traps employed with this method inhibit the activity of nitric oxide synthase and catalyze NO production from nitrite. These disadvantages limit EPR's application to biological NO detection. We present a liposome-encapsulated spin-trap (LEST) method for the capture and in situ detection of NO by EPR. The method shows a linear response for [NO]≥4µM and can detect [NO]≥40nM in a 500µL sample (≥20 pmol). The kinetics of NO production can be followed in real time over minutes to hours. LEST does not inhibit the activity of inducible nitric oxide synthase or nitrate reductase and shows minimal abiotic NO production in the presence of nitrite and NADH. Nitrate reductase-like activity is detected in cell lysates of the coccolithophore Emiliania huxleyi and is elevated in virus-infected culture. This method shows particular promise for NO detection in cell lysates and crude preparations of NO-producing tissues.

Movement of NH₃ through the human urea transporter B: a new gas channel

Aquaporins and Rh proteins can function as gas (CO₂ and NH₃) channels. The present study explores the urea, H₂O, CO₂, and NH₃ permeability of the human urea transporter B (UT-B) (SLC14A1), expressed in Xenopus oocytes. We monitored urea uptake using [¹⁴C]urea and measured osmotic water permeability (Pf) using video microscopy. To obtain a semiquantitative measure of gas permeability, we used microelectrodes to record the maximum transient change in surface pH (ΔpHS) caused by exposing oocytes to 5% CO₂/33 mM HCO₃⁻ (pHS increase) or 0.5 mM NH₃/NH₄⁺ (pHS decrease). UT-B expression increased oocyte permeability to urea by >20-fold, and Pf by 8-fold vs. H₂O-injected control oocytes. UT-B expression had no effect on the CO₂-induced ΔpHS but doubled the NH₃-induced ΔpHS. Phloretin reduced UT-B-dependent urea uptake (Jurea*) by 45%, Pf* by 50%, and (- ΔpHS*)NH₃ by 70%. p-Chloromercuribenzene sulfonate reduced Jurea* by 25%, Pf* by 30%, and (ΔpHS*)NH₃ by 100%. Molecular dynamics (MD) simulations of membrane-embedded models of UT-B identified the monomeric UT-B pores as the main conduction pathway for both H₂O and NH₃ and characterized the energetics associated with permeation of these species through the channel. Mutating each of two conserved threonines lining the monomeric urea pores reduced H₂O and NH₃ permeability. Our data confirm that UT-B has significant H₂O permeability and for the first time demonstrate significant NH₃ permeability. Thus the UTs become the third family of gas channels. Inhibitor and mutagenesis studies and results of MD simulations suggest that NH₃ and H₂O pass through the three monomeric urea channels in UT-B.

Visualizing functional motions of membrane transporters with molecular dynamics simulations

Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins.

Accelerating membrane insertion of peripheral proteins with a novel membrane mimetic model

Characterizing atomic details of membrane binding of peripheral membrane proteins by molecular dynamics (MD) has been significantly hindered by the slow dynamics of membrane reorganization associated with the phenomena. To expedite lateral diffusion of lipid molecules without sacrificing the atomic details of such interactions, we have developed a novel membrane representation, to our knowledge, termed the highly mobile membrane-mimetic (HMMM) model to study binding and insertion of various molecular species into the membrane. The HMMM model takes advantage of an organic solvent layer to represent the hydrophobic core of the membrane and short-tailed phospholipids for the headgroup region. We demonstrate that using these components, bilayer structures are formed spontaneously and rapidly, regardless of the initial position and orientation of the lipids. In the HMMM membrane, lipid molecules exhibit one to two orders of magnitude enhancement in lateral diffusion. At the same time, the membrane atomic density profile of the headgroup region produced by the HMMM model is essentially identical to those obtained for full-membrane models, indicating the faithful representation of the membrane surface by the model. We demonstrate the efficiency of the model in capturing spontaneous binding and insertion of peripheral proteins by using the membrane anchor (γ-carboxyglutamic-acid-rich domain; GLA domain) of human coagulation factor VII as a test model. Achieving full insertion of the GLA domain consistently in 10 independent unbiased simulations within short simulation times clearly indicates the robustness of the HMMM model in capturing membrane association of peripheral proteins very efficiently and reproducibly. The HMMM model will provide significant improvements to the current all-atom models by accelerating lipid dynamics to examine protein-membrane interactions more efficiently.

Insight into the Properties of Cardiolipin Containing Bilayers from Molecular Dynamics Simulations, Using a Hybrid All-Atom/United-Atom Force Field

Simulation of three models of cardiolipin (CL) containing membranes using a new set of parameters for tetramyristoyl and tetraoleoyl CLs has been developed in the framework of the united-atom CHARMM27-UA and the all-atom CHARMM36 force fields with the aim of performing molecular dynamics (MD) simulations of cardiolipin-containing mixed-lipid membranes. The new parameters use a hybrid representation of all-atom head groups in conjunction with implicit-hydrogen united-atom (UA) to describe the oleoyl and myristoyl chains of the CLs, in lieu of the fully atomistic description, thereby allowing longer simulations to be undertaken. The physicochemical properties of the bilayers were determined and compared with previously reported data. Furthermore, using tetramyristoyl CL mixed with POPG and POPE lipids, a mitochondrial membrane was simulated. The results presented here show the different behavior of the bilayers as a result of the lipid composition, where the length of the acyl chain and the conformation of the headgroup can be associated with the mitochondrial membrane properties. The new hybrid CL parameters prove to be well suited for the simulation of the molecular structure of CL-containing bilayers and can be extended to other lipid bilayers composed of CLs with different acyl chains or alternate head groups.

Atomic view of calcium-induced clustering of phosphatidylserine in mixed lipid bilayers

Membranes play key regulatory roles in biological processes, with bilayer composition exerting marked effects on binding affinities and catalytic activities of a number of membrane-associated proteins. In particular, proteins involved in diverse processes such as vesicle fusion, intracellular signaling cascades, and blood coagulation interact specifically with anionic lipids such as phosphatidylserine (PS) in the presence of Ca(2+) ions. While Ca(2+) is suspected to induce PS clustering in mixed phospholipid bilayers, the detailed structural effects of this ion on anionic lipids are not established. In this study, combining magic angle spinning (MAS) solid-state NMR (SSNMR) measurements of isotopically labeled serine headgroups in mixed lipid bilayers with molecular dynamics (MD) simulations of PS lipid bilayers in the presence of different counterions, we provide site-resolved insights into the effects of Ca(2+) on the structure and dynamics of lipid bilayers. Ca(2+)-induced conformational changes of PS in mixed bilayers are observed in both liposomes and Nanodiscs, a nanoscale membrane mimetic of bilayer patches. Site-resolved multidimensional correlation SSNMR spectra of bilayers containing (13)C,(15)N-labeled PS demonstrate that Ca(2+) ions promote two major PS headgroup conformations, which are well resolved in two-dimensional (13)C-(13)C, (15)N-(13)C, and (31)P-(13)C spectra. The results of MD simulations performed on PS lipid bilayers in the presence or absence of Ca(2+) provide an atomic view of the conformational effects underlying the observed spectra.

Exploring transmembrane diffusion pathways with molecular dynamics

Transmembrane exchange of materials is a fundamental process in biology. Molecular dynamics provides a powerful method to investigate in great detail various aspects of the phenomenon, particularly the permeation of small uncharged molecules, which continues to pose a challenge to experimental studies. We will discuss some of the recent simulation studies investigating the role of lipid-mediated and protein-mediated mechanisms in permeation of water and gas molecules across the membrane.

Nitric oxide conduction by the brain aquaporin AQP4

Involvement of aquaporins in gas conduction across the membrane and the physiological significance of this process have attracted marked attention from both experimental and theoretical studies. Previous work demonstrated that AQP1 is permeable to both CO(2) and O(2). Here we employ various simulation techniques to examine the permeability of the brain aquaporin AQP4 to NO and O(2) and to describe energetics and pathways associated with these phenomena. The energy barrier to NO and O(2) permeation through AQP4 central pore is found to be only approximately 3 kcal mol(-1). The results suggest that the central pore of AQP4, similar to that of AQP1, can indeed conduct gas molecules. Interestingly, despite a longer and narrower central pore, AQP4 appears to provide an energetically more favorable permeation pathway for gas molecules than AQP1, mainly due to the different orientation of its charged residues near the pore entrance. Although the low barrier against gas permeation through AQP4 indicates that it can participate in gas conduction across the cellular membrane, physiological relevance of the phenomenon remains to be established experimentally, particularly since pure lipid bilayers appear to present a more favorable pathway for gas conduction across the membrane. With an energy well of -1.8 kcal mol(-1), the central pore of AQP4 may also act as a reservoir for NO molecules to accumulate in the membrane.