Towards an understanding of complex biological membranes from atomistic molecular dynamics simulations

Computer simulation has emerged as a powerful tool for studying the structural and functional properties of complex biological membranes. In the last few years, the use of recently developed simulation methodologies and current generation force fields has permitted novel applications of molecular dynamics simulations, which have enhanced our understanding of the different physical processes governing biomembrane structure and dynamics. This review focuses on frontier areas of research with important biomedical applications. We have paid special attention to polyunsaturated lipids, membrane proteins and ion channels, surfactant additives in membranes, and lipid–DNA gene transfer complexes.     

K and Na Conduction in Selective and Nonselective Ion Channels Via Molecular Dynamics Simulations

Generations of scientists have been captivated by ion channels and how they control the workings of the cell by admitting ions from one side of the cell membrane to the other. Elucidating the molecular determinants of ion conduction and selectivity are two of the most fundamental issues in the field of biophysics. Combined with ongoing progress in structural studies, modeling and simulation have been an integral part of the development of the field. As of this writing, the relentless growth in computational power, the development of new algorithms to tackle the so-called rare events, improved force-field parameters, and the concomitant increasing availability of membrane protein structures, allow simulations to contribute even further, providing more-complete models of ion conduction and selectivity in ion channels. In this report, we give an overview of the recent progress made by simulation studies on the understanding of ion permeation in selective and nonselective ion channels.     

K+ and Na+ Conduction in Selective and Nonselective Ion Channels Via Molecular Dynamics Simulations

Generations of scientists have been captivated by ion channels and how they control the workings of the cell by admitting ions from one side of the cell membrane to the other. Elucidating the molecular determinants of ion conduction and selectivity are two of the most fundamental issues in the field of biophysics. Combined with ongoing progress in structural studies, modeling and simulation have been an integral part of the development of the field. As of this writing, the relentless growth in computational power, the development of new algorithms to tackle the so-called rare events, improved force-field parameters, and the concomitant increasing availability of membrane protein structures, allow simulations to contribute even further, providing more-complete models of ion conduction and selectivity in ion channels. In this report, we give an overview of the recent progress made by simulation studies on the understanding of ion permeation in selective and nonselective ion channels.     

Interactions of alpha-helices with lipid bilayers: a review of simulation studies

Membrane proteins, of which the majority seem to contain one or more alpha-helix, constitute approx. 30% of most genomes. A complete understanding of the nature of helix/bilayer interactions is necessary for an understanding of the structural principles underlying membrane proteins. This review describes computer simulation studies of helix/bilayer interactions. Key experimental studies of the interactions of alpha-helices and lipid bilayers are briefly reviewed. Surface associated helices are found in some membrane-bound enzymes (e.g. prostaglandin synthase), and as stages in the mechanisms of antimicrobial peptides and of pore-forming bacterial toxins. Transmembrane alpha-helices are found in most integral membrane proteins, and also in channels formed by amphipathic peptides or by bacterial toxins. Mean field simulations, in which the lipid bilayer is approximated as a hydrophobic continuum, have been used in studies of membrane-active peptides (e.g. alamethicin, melittin, magainin and dermaseptin) and of simple membrane proteins (e.g. phage Pf1 coat protein). All atom molecular dynamics simulations of fully solvated bilayers with transmembrane helices have been applied to: the constituent helices of bacteriorhodopsin; peptide-16 (a simple model TM helix); and a number of pore-lining helices from ion channels. Surface associated helices (e.g. melittin and dermaseptin) have been simulated, as have alpha-helical bundles such as bacteriorhodopsin and alamethicin. From comparison of the results from the two classes of simulation, it emerges that a major theoretical challenge is to exploit the results of all atom simulations in order to improve the mean field approach.     

Molecular dynamics simulations of membrane proteins

Membrane proteins control the traffic across cell membranes and thereby play an essential role in cell function from transport of various solutes to immune response via molecular recognition. Because it is very difficult to determine the structures of membrane proteins experimentally, computational methods have been increasingly used to study their structure and function. Here we focus on two classes of membrane proteins—ion channels and transporters—which are responsible for the generation of action potentials in nerves, muscles, and other excitable cells. We describe how computational methods have been used to construct models for these proteins and to study the transport mechanism. The main computational tool is the molecular dynamics (MD) simulation, which can be used for everything from refinement of protein structures to free energy calculations of transport processes. We illustrate with specific examples from gramicidin and potassium channels and aspartate transporters how the function of these membrane proteins can be investigated using MD simulations.     

Molecular dynamics simulations and membrane protein structure quality

Despite a growing repertoire of membrane protein structures (currently ∼120 unique structures), considerations of low resolution and crystallization in the absence of a lipid bilayer require the development of techniques to assess the global quality of membrane protein folds. This is also the case for assessment of, e.g. homology models of human membrane proteins based on structures of (distant) bacterial homologues. Molecular dynamics (MD) simulations may be used to help evaluate the quality of a membrane protein structure or model. We have used a structure of the bacterial ABC transporter MsbA which has the correct transmembrane helices but an incorrect handedness and topology of their packing to test simulation methods of quality assessment. An MD simulation of the MsbA model in a lipid bilayer is compared to a simulation of another bacterial ABC transporter, BtuCD. The latter structure has demonstrated good conformational stability in the same bilayer environment and over the same timescale (20 ns) as for the MsbA model simulation. A number of comparative analyses of the two simulations were performed to assess changes in the structural integrity of each protein. The results show a significant difference between the two simulations, chiefly due to the dramatic structural deformations of MsbA. We therefore propose that MD could become a useful quality control tool for membrane protein structural biology. In particular, it provides a way in which to explore the global conformational stability of a model membrane protein fold.     

Interaction of local anesthetics with the K+ channel pore domain

Local anesthetics and related drugs block ionic currents of Na⁺, K⁺ and Ca²⁺ conducted across the cell membrane by voltage-dependent ion channels. Many of these drugs bind in the permeation pathway, occlude the pore and stop ion movement. However channel-blocking drugs have also been associated with decreased membrane stability of certain tetrameric K⁺ channels, similar to the destabilization of channel function observed at low extracellular K⁺ concentration. Such drug-dependent stability may result from electrostatic repulsion of K⁺ from the selectivity filter by a cationic drug molecule bound in the central cavity of the channel. In this study we used the pore domain of the KcsA K⁺ channel protein to test this hypothesis experimentally with a biochemical assay of tetramer stability and theoretically by computational simulation of local anesthetic docking to the central cavity. We find that two common local anesthetics, lidocaine and tetracaine, promote thermal dissociation of the KcsA tetramer in a K⁺-dependent fashion. Docking simulations of these drugs with open, open-inactivated and closed crystal structures of KcsA yield many energetically favorable drug-channel complexes characterized by nonbonded attraction to pore-lining residues and electrostatic repulsion of K⁺. The results suggest that binding of cationic drugs to the inner cavity can reduce tetramer stability of K⁺ channels.     

Molecular dynamics simulations of proteins in lipid bilayers

With recent advances in X-ray crystallography of membrane proteins promising many new high-resolution structures, molecular dynamics simulations will become increasingly valuable for understanding membrane protein function, as they can reveal the dynamic behavior concealed in the static structures. Dramatic increases in computational power, in synergy with more efficient computational methodologies, now allow us to carry out molecular dynamics simulations of any structurally known membrane protein in its native environment, covering timescales of up to 0.1 micros. At the frontiers of membrane protein simulations are ion channels, aquaporins, passive and active transporters, and bioenergetic proteins.     

Structural determinants of lateral gate opening in the protein translocon

The heterotrimeric SecY/Sec61 complex is a protein-conducting channel that provides a passage for proteins across the membrane as well as a means to integrate nascent proteins into the membrane. While the first function is common among membrane protein channels and transporters, the latter is unique. Insertion of nascent membrane proteins, one transmembrane segment at a time, by SecY likely occurs through a lateral gate in the channel. Molecular dynamics simulations have been used to investigate the mechanism of gate opening. Opening and closing the gate under different conditions allowed us to identify structural elements that resist opening as well as those that aid closure. SecE, considered to act as a clamp keeping the lateral gate closed, was found to play no such role. Loosening of the plug by lateral gate opening, a potential step in channel gating, was also observed. The simulations revealed that lipids on time scales of up to 1 micros do not flood channels with an open lateral gate.     

Simulation of MscL gating in a bilayer under stress

The initial stages of the gating of the mechanoselective channel of large conductance from Mycobacterium tuberculosis have been studied in atomic detail using molecular dynamics simulation techniques. A truncated form of the protein complex embedded in a palmitoyloleoylphosphatidylcholine lipid bilayer and surrounded by explicit water was simulated under different pressure conditions to mimic the effects of tension and compression within the membrane on the protein. As a direct result of lateral tension being applied to the membrane, an increase in the tilt of a subset of the transmembrane helices was observed. This in turn led to the enlargement of the pore and the disruption of the hydrophobic gate consisting of residues Ile-14 and Val-21. The simulations suggest that opening occurs in a sequential staged process. Such a mechanism could explain the partial opening or staged conductance observed in patch-clamp experiments using related large conductance mechanosensitive channel complexes.     

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.     

Dynamics and function in a bacterial ABC transporter: simulation studies of the BtuCDF system and its components

ABC transporters are integral membrane proteins which couple the energy of ATP hydrolysis to the translocation of solutes across cell membranes. BtuCD is a approximately 1100-residue protein found in the inner membrane of Gram-negative bacteria which transports vitamin B12. Vitamin B12 is bound in the periplasm by BtuF, which delivers the solute to the periplasmic entrance of the transporter protein complex BtuCD. Molecular dynamics simulations of the BtuCD and BtuCDF complexes (in a lipid bilayer) and of the isolated BtuD and BtuF proteins (in water) have been used to explore the conformational dynamics of this complex transport system. Overall, seven simulations have been performed, with and without bound ATP, corresponding to a total simulation time of 0.1 micros. Binding of ATP drives closure of the nucleotide-binding domains (NBDs) in BtuD in a symmetrical fashion, but not in BtuCD. It seems that ATP constrains the flexibility of the NBDs in BtuCD such that their closure may only occur upon binding of BtuF to the complex. Upon introduction of BtuF, and concomitant with NBD association, one ATP-binding site displays a closure, while the opposite site remains relatively unchanged. This asymmetry may reflect an initial step in the "alternating hydrolysis" mechanism and is consistent with measurements of nucleotide-binding stoichiometries. Principal components analysis of the simulation of BtuCD reveals motions that are comparable to those suggested in current transport models.     

Filter flexibility in a mammalian K channel: models and simulations of Kir6.2 mutants

The single-channel conductance varies significantly between different members of the inward rectifier (Kir) family of potassium channels. Mutations at three sites in Kir6.2 have been shown to produce channels with reduced single-channel conductance, the largest reduction (to 40% of wild-type) being for V127T. We have used homology modeling (based on a KcsA template) combined with molecular dynamics simulations in a phosphatidycholine bilayer to explore whether changes in structural dynamics of the filter were induced by three such mutations: V127T, M137C, and G135F. Overall, 12 simulations of Kir6.2 models, corresponding to a total simulation time of 27 ns, have been performed. In these simulations we focused on distortions of the selectivity filter, and on the presence/absence of water molecules lying behind the filter, which form interactions with the filter and the remainder of the protein. Relative to the wild-type simulation, the V127T mutant showed significant distortion of the filter such that approximately 50% of the simulation time was spent in a closed conformation. While in this conformation, translocation of K(+) ions between sites S1 and S2 was blocked. The distorted filter conformation resembles that of the bacterial channel KcsA when crystallized in the presence of a low [K(+)]. This suggests filter distortion may be a possible general model for determining the conductance of K channels.     

Coarse-grained simulations of protein-protein association: an energy landscape perspective

Understanding protein-protein association is crucial in revealing the molecular basis of many biological processes. Here, we describe a theoretical simulation pipeline to study protein-protein association from an energy landscape perspective. First, a coarse-grained model is implemented and its applications are demonstrated via molecular dynamics simulations for several protein complexes. Second, an enhanced search method is used to efficiently sample a broad range of protein conformations. Third, multiple conformations are identified and clustered from simulation data and further projected on a three-dimensional globe specifying protein orientations and interacting energies. Results from several complexes indicate that the crystal-like conformation is favorable on the energy landscape even if the landscape is relatively rugged with metastable conformations. A closer examination on molecular forces shows that the formation of associated protein complexes can be primarily electrostatics-driven, hydrophobics-driven, or a combination of both in stabilizing specific binding interfaces. Taken together, these results suggest that the coarse-grained simulations and analyses provide an alternative toolset to study protein-protein association occurring in functional biomolecular complexes.     

Membrane-delimited cell signaling complexes: Direct ion channel regulation by G Proteins

Summary Ion channels are signaling molecules and by them-selves perform no work. In this regard they are un like the usual membrane enzyme effectors for G proteins. The pathways of G protein receptor, G protein and ion channels are, therefore, purely infor mational in function. Because a single G protein may have several ion channels as effectors, the effects should be coordinated and this seems to be the case. Inhibition of Ca²⁺ current and stimulation of K⁺ currents would have a greater impact than either alone. Additional flexibility is provided by spontane ous noise in the complexes of G protein receptor, G protein, and ion channel. By having a non-zero setpoint, the range of control is extended and the responses become bi-directional.     

Structure and dynamics of one-dimensional ionic solutions in biological transmembrane channels.

The structure and dynamics of solvated alkali metal cations in transmembrane channels are treated using the molecular dynamics simulation technique. The simulations are based on a modified Fischer-Brickmann model (Fischer, W., and J. Brickmann, 1983, Biophys. Chem., 18:323-337) for gramicidin A-type channels. The trajectories of all particles in the channel as well as two-dimensional pair correlation functions are analyzed. It is found from the analysis of the stationary simulation state that one-dimensional solvation complexes are formed and that the number of water molecules in the channel varies for different alkali metal cations.     

Membrane proteins: molecular dynamics simulations

Molecular dynamics simulations of membrane proteins are making rapid progress, because of new high-resolution structures, advances in computer hardware and atomistic simulation algorithms, and the recent introduction of coarse-grained models for membranes and proteins. In addition to several large ion channel simulations, recent studies have explored how individual amino acids interact with the bilayer or snorkel/anchor to the headgroup region, and it has been possible to calculate water/membrane partition free energies. This has resulted in a view of bilayers as being adaptive rather than purely hydrophobic solvents, with important implications, for example, for interaction between lipids and arginines in the charged S4 helix of voltage-gated ion channels. However, several studies indicate that the typical current simulations fall short of exhaustive sampling, and that even simple protein-membrane interactions require at least ca. 1mus to fully sample their dynamics. One new way this is being addressed is coarse-grained models that enable mesoscopic simulations on multi-mus scale. These have been used to model interactions, self-assembly and membrane perturbations induced by proteins. While they cannot replace all-atom simulations, they are a potentially useful technique for initial insertion, placement, and low-resolution refinement.     

Regulation of channel function due to physical energetic coupling with a lipid bilayer

Regulation of membrane protein functions due to hydrophobic coupling with a lipid bilayer has been investigated. An energy formula describing interactions between lipid bilayer and integral ion channels with different structures, which is based on the screened Coulomb interaction approximation, has been developed. Here the interaction energy is represented as being due to charge-based interactions between channel and lipid bilayer. The hydrophobic bilayer thickness channel length mismatch is found to induce channel destabilization exponentially while negative lipid curvature linearly. Experimental parameters related to channel dynamics are consistent with theoretical predictions. To measure comparable energy parameters directly in the system and to elucidate the mechanism at an atomistic level we performed molecular dynamics (MD) simulations of the ion channel forming peptide–lipid complexes. MD simulations indicate that peptides and lipids experience electrostatic and van der Waals interactions for short period of time when found within each other’s proximity. The energies from these two interactions are found to be similar to the energies derived theoretically using the screened Coulomb and the van der Waals interactions between peptides (in ion channel) and lipids (in lipid bilayer) due to mainly their charge properties. The results of in silico MD studies taken together with experimental observable parameters and theoretical energetic predictions suggest that the peptides induce ion channels inside lipid membranes due to peptide–lipid physical interactions. This study provides a new insight helping better understand of the underlying mechanisms of membrane protein functions in cell membrane leading to important biological implications.     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.