Explicit Lipid Research Articles

Modeling and Molecular Dynamics Simulation of the Human Gonadotropin-Releasing Hormone Receptor in a Lipid Bilayer

In the present study, a model for the human gonadotropin-releasing hormone receptor embedded in an explicit lipid bilayer was developed. The final conformation was obtained by extensive molecular dynamics simulations of a homology model based on the bovine rhodopsin crystal structure. The analysis of the receptor structure allowed us to detect a number of specific contacts between different amino acid residues, as well as water- and lipid-mediated interactions. These interactions were stable in six additional independent 35 ns long simulations at 310 and 323 K, which used the refined model as the starting structure. All loops, particularly the extracellular loop 2 and the intracellular loop 3, exhibited high fluctuations, whereas the transmembrane helices were more static. Although other models of this receptor have been previously developed, none of them have been subjected to extensive molecular dynamics simulations, and no other three-dimensional structure is publicly available. Our results suggest that the presence of ions as well as explicit solvent and lipid molecules are critical for the structure of membrane protein models, and that molecular dynamics simulations are certainly useful for their refinement.

Spatial Structure of the Dimeric Transmembrane Domain of the Growth Factor Receptor ErbB2 Presumably Corresponding to the Receptor Active State

Proper lateral dimerization of the transmembrane domains of receptor tyrosine kinases is required for biochemical signal transduction across the plasma membrane. The spatial structure of the dimeric transmembrane domain of the growth factor receptor ErbB2 embedded into lipid bicelles was obtained by solution NMR, followed by molecular dynamics relaxation in an explicit lipid bilayer. ErbB2 transmembrane segments associate in a right-handed alpha-helical bundle through the N-terminal tandem GG4-like motif Thr652-X3-Ser656-X3-Gly660, providing an explanation for the pathogenic power of some oncogenic mutations.

How to lose a kink and gain a helix: pH independent conformational changes of the fusion domains from influenza hemagglutinin in heterogeneous lipid bilayers

We have simulated two conformations of the fusion domain of influenza hemagglutinin (HA) within explicit water, salt, and heterogeneous lipid bilayers composed of POPC:POPG (4:1). Each conformation has seven different starting points in which the initial peptide structure is the same for each conformation, but the location across the membrane normal and lipid arrangement around the peptide are varied, giving a combined total simulation time of 140 ns. For the HA5 conformation (primary structure from recent NMR spectroscopy at pH = 5), the peptide exhibits a stable and less kinked structure in the lipid bilayer compared to that from the NMR studies. The relative fusogenic behavior of the different conformations has been investigated by calculation of the relative free energy of insertion into the hydrophobic region of lipid bilayer as a function of the depth of immersion. For the HA7 conformations (primary structure from recent NMR spectroscopy at pH = 7.4), while the N-terminal helix preserves its initial structure, the flexible C-terminal chain produces a transient helical motif inside the lipid bilayer. This conformational change is pH-independent, and is closely related to the peptide insertion into the lipid bilayer.

Phospholemman Transmembrane Structure Reveals Potential Interactions with Na+/K+-ATPase

Phospholemman (PLM) is a 72-residue bitopic cardiac transmembrane protein, which acts as a modulator of the Na(+)/K(+)-ATPase and the Na(+)/Ca(2+) exchanger and possibly forms taurine channels in nonheart tissue. This work presents a high resolution structural model obtained from a combination of site-specific infrared spectroscopy and experimentally constrained high throughput molecular dynamics (MD) simulations. Altogether, 37 experimental constraints, including nine long range orientational constraints, have been used during MD simulations in an explicit lipid bilayer/water system. The resulting tetrameric alpha-helical bundle has an average helix tilt of 7.3 degrees and a crossing angle close to 0 degrees . It does not reveal a hydrophilic pore, but instead strong interactions between various residues occlude any pore. The helix-helix packing is unusual, with Gly(19) and Gly(20) pointing to the outside of the helical bundle, facilitating potential interaction with other transmembrane proteins, thus providing a structural basis for the modulatory effect of PLM on the Na(+)/K(+)-ATPase. A two-stage model of interaction between PLM and the Na(+)/K(+)-ATPase is discussed involving PLM-ATPase interaction and subsequent formation of an unstable PLM trimer, which readily interacts with surrounding ATPase molecules. Further unconstrained MD simulations identified other packing models of PLM, one of which could potentially undergo a conformational transition to an open pore.

Folding Simulations of the Transmembrane Helix of Virus Protein U in an Implicit Membrane Model.

Vpu is an 81-amino-acid auxiliary membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). One of its roles is to amplify viral release by self-assembling in homo-oligomers to form functional water-filled pores enabling the flux of ions across the membrane. Various NMR and CD studies have shown that the transmembrane domain of Vpu has a helical conformation. With a recently developed implicit membrane model and an efficient Monte Carlo (MC) algorithm using concerted backbone rotations, we simulate the folding of the transmembrane domain of Vpu at atomic resolution. The implicit membrane environment is based on the generalized Born theory and enables very long time scale events, such as folding to be observed using detailed all-atom representation of the protein. Such studies are currently computationally unfeasible with fully explicit lipid bilayer molecular dynamics simulations. The correct helical transmembrane structure of Vpu is predicted from extended conformations and remains stably inserted. Tilt and kink angles agree well with experimental estimates from NMR measurements. The experimentally observed change in tilt angle in membranes of varying hydrophobic width is accurately reproduced. The extensive simulation of a pentamer of the Vpu transmembrane domain in the implicit membrane gives results similar to the ones reported previously for fully explicit bilayer simulations.

Mechanism of Activation of a G Protein-coupled Receptor, the Human Cholecystokinin-2 Receptor

G protein-coupled receptors (GPCRs) represent a major focus in functional genomics programs and drug development research, but their important potential as drug targets contrasts with the still limited data available concerning their activation mechanism. Here, we investigated the activation mechanism of the cholecystokinin-2 receptor (CCK2R). The three-dimensional structure of inactive CCK2R was homology-modeled on the basis of crystal coordinates of inactive rhodopsin. Starting from the inactive CCK2R modeled structure, active CCK2R (namely cholecystokinin-occupied CCK2R) was modeled by means of steered molecular dynamics in a lipid bilayer and by using available data from other GPCRs, including rhodopsin. By comparing the modeled structures of the inactive and active CCK2R, we identified changes in the relative position of helices and networks of interacting residues, which were expected to stabilize either the active or inactive states of CCK2R. Using targeted molecular dynamics simulations capable of converting CCK2R from the inactive to the active state, we delineated structural changes at the atomic level. The activation mechanism involved significant movements of helices VI and V, a slight movement of helices IV and VII, and changes in the position of critical residues within or near the binding site. The mutation of key amino acids yielded inactive or constitutively active CCK2R mutants, supporting this proposed mechanism. Such progress in the refinement of the CCK2R binding site structure and in knowledge of CCK2R activation mechanisms will enable target-based optimization of nonpeptide ligands.

Prediction of the 3D structure and dynamics of human DP G-protein coupled receptor bound to an agonist and an antagonist.

Prostanoids play important physiological roles in the cardiovascular and immune systems and in pain sensation in peripheral systems through their interactions with eight G-protein coupled receptors. These receptors are important drug targets, but development of subtype specific agonists and antagonists has been hampered by the lack of 3D structures for these receptors. We report here the 3D structure for the human DP G-protein coupled receptor (GPCR) predicted by the MembStruk computational method. To validate this structure, we use the HierDock computational method to predict the binding mode for the endogenous agonist (PGD2) to DP. Based on our structure, we predicted the binding of different antagonists and optimized them. We find that PGD2 binds vertically to DP in the TM1237 region with the alpha chain toward the extracellular (EC) region and the omega chain toward the middle of the membrane. This structure explains the selectivity of the DP receptor and the residues involved in the predicted binding site correlate very well with available mutation experiments on DP, IP, TP, FP, and EP subtypes. We report molecular dynamics of DP in explicit lipid and water and find that the binding of the PGD2 agonist leads to correlated rotations of helices of TM3 and TM7, whereas binding of antagonist leads to no such rotations. Thus, these motions may be related to the mechanism of activation.

Unique Dimeric Structure of BNip3 Transmembrane Domain Suggests Membrane Permeabilization as a Cell Death Trigger

BNip3 is a prominent representative of apoptotic Bcl-2 proteins with rather unique properties initiating an atypical programmed cell death pathway resembling both necrosis and apoptosis. Many Bcl-2 family proteins modulate the permeability state of the outer mitochondrial membrane by forming homo- and hetero-oligomers. The structure and dynamics of the homodimeric transmembrane domain of BNip3 were investigated with the aid of solution NMR in lipid bicelles and molecular dynamics energy relaxation in an explicit lipid bilayer. The right-handed parallel helix-helix structure of the domain with a hydrogen bond-rich His-Ser node in the middle of the membrane, accessibility of the node for water, and continuous hydrophilic track across the membrane suggest that the domain can provide an ion-conducting pathway through the membrane. Incorporation of the BNip3 transmembrane domain into an artificial lipid bilayer resulted in pH-dependent conductivity increase. A possible biological implication of the findings in relation to triggering necrosis-like cell death by BNip3 is discussed.

Conformational study of the protegrin-1 (PG-1) dimer interaction with lipid bilayers and its effect

BackgroundProtegrin-1 (PG-1) is known as a potent antibiotic peptide; it prevents infection via an attack on the membrane surface of invading microorganisms. In the membrane, the peptide forms a pore/channel through oligomerization of multiple subunits. Recent experimental and computational studies have increasingly unraveled the molecular-level mechanisms underlying the interactions of the PG-1 β-sheet motifs with the membrane. The PG-1 dimer is important for the formation of oligomers, ordered aggregates, and for membrane damaging effects. Yet, experimentally, different dimeric behavior has been observed depending on the environment: antiparallel in the micelle environment, and parallel in the POPC bilayer. The experimental structure of the PG-1 dimer is currently unavailable.ResultsAlthough the β-sheet structures of the PG-1 dimer are less stable in the bulk water environment, the dimer interface is retained by two intermolecular hydrogen bonds. The formation of the dimer in the water environment implies that the pathway of the dimer invasion into the membrane can originate from the bulk region. In the initial contact with the membrane, both the antiparallel and parallel β-sheet conformations of the PG-1 dimer are well preserved at the amphipathic interface of the lipid bilayer. These β-sheet structures illustrate the conformations of PG-1 dimer in the early stage of the membrane attack. Here we observed that the activity of PG-1 β-sheets on the bilayer surface is strongly correlated with the dimer conformation. Our long-term goal is to provide a detailed mechanism of the membrane-disrupting effects by PG-1 β-sheets which are able to attack the membrane and eventually assemble into the ordered aggregates.ConclusionIn order to understand the dimeric effects leading to membrane damage, extensive molecular dynamics (MD) simulations were performed for the β-sheets of the PG-1 dimer in explicit water, salt, and lipid bilayers composed of POPC lipids. Here, we studied PG-1 dimers when organized into a β-sheet motif with antiparallel and parallel β-sheet arrangements in an NCCN packing mode. We focus on the conformations of PG-1 dimers in the lipid bilayer, and on the correlation between the conformations and the membrane disruption effects by PG-1 dimers. We investigate equilibrium structures of the PG-1 dimers in different environments in the early stage of the dimer invasion. The dimer interface of the antiparallel β-sheets is more stable than the parallel β-sheets, similar to the experimental observation in micelle environments. However, we only observe membrane disruption effects by the parallel β-sheets of the PG-1 dimer. This indicates that the parallel β-sheets interact with the lipids with the β-sheet plane lying obliquely to the bilayer surface, increasing the surface pressure in the initial insertion into the lipid bilayer. Recent experimental observation verified that parallel PG-1 dimer is biologically more active to insert into the POPC lipid bilayer.

Effect of different treatments of long‐range interactions and sampling conditions in molecular dynamic simulations of rhodopsin embedded in a dipalmitoyl phosphatidylcholine bilayer

The present study analyzes the effect of the simulation conditions on the results of molecular dynamics simulations of G-protein coupled receptors (GPCRs) performed with an explicit lipid bilayer. Accordingly, the present work reports the analysis of different simulations of bovine rhodopsin embedded in a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer using two different sampling conditions and two different approaches for the treatment of long-range electrostatic interactions. Specifically, sampling was carried out either by using the statistical ensembles NVT or NPT (constant number of atoms, a pressure of 1 atm in all directions and fixed temperature), and the electrostatic interactions were treated either by using a twin-cutoff, or the particle mesh Ewald summation method (PME). The results of the present study suggest that the use of the NPT ensemble in combination with the PME method provide more realistic simulations. The use of NPT during the equilibration avoids the need of an a priori estimation of the box dimensions, giving the correct area per lipid. However, once the system is equilibrated, the simulations are irrespective of the sampling conditions used. The use of an electrostatic cutoff induces artifacts on both lipid thickness and the ion distribution, but has no direct effect on the protein and water molecules.

Interaction of Protegrin-1 with Lipid Bilayers: Membrane Thinning Effect

Protegrins (PG) are important in defending host tissues, preventing infection via an attack on the membrane surface of invading microorganisms. Protegrins have powerful antibiotic abilities, but the molecular-level mechanisms underlying the interactions of their β-sheet motifs with the membrane are not known. Protegrin-1 (PG-1) is composed of 18 amino acids with a high content of basic residues and two disulfide bonds. Here we focused on the stability of PG-1 at the amphipathic interface in lipid bilayers and on the details of the peptide-membrane interactions. We simulated all-atom models of the PG-1 monomer with explicit water and lipid bilayers composed of both homogeneous POPC (palmitoyl-oleyl-phosphatidylcholine) lipids and a mixture of POPC/POPG (palmitoyl-oleyl-phosphatidylglycerol) (4:1) lipids. We observed that local thinning of the lipid bilayers mediated by the peptide is enhanced in the lipid bilayer containing POPG, consistent with experimental results of selective membrane targeting. The β-hairpin motif of PG-1 is conserved in both lipid settings, whereas it is highly bent in aqueous solution. The conformational dynamics of PG-1, especially the highly charged β-hairpin turn region, are found to be mostly responsible for disturbing the membrane. Even though the eventual membrane disruption requires PG-1 oligomers, our simulations clearly show the first step of the monomeric effects. The thinning effects in the bilayer should relate to pore/channel formation in the lipid bilayer and thus be responsible for further defects in the membrane caused by oligomer.

Dynamic behavior of fully solvated β2-adrenergic receptor, embedded in the membrane with bound agonist or antagonist

Recently we predicted the 3D structure of the human beta2-adrenergic receptor (beta2AR) and of the binding site of several agonists and antagonists to beta2AR. These predictions (MembStruk and HierDock) included no explicit water and only a few lipid molecules. Here we include explicit H(2)O and an infinite lipid bilayer membrane in molecular dynamics (MD) simulations of three systems: apo-beta2AR, epinephrine-bound beta2AR, and butoxamine-bound beta2AR (epinephrine is an endogenous agonist, and butoxamine is a beta2AR selective antagonist). The predicted structures for apo-beta2AR and butoxamine-beta2AR are stable in MD, but in epinephrine-beta2AR, extracellular water trickles into the binding pocket to mediate hydrogen bonding between the catechol of epinephrine and Ser-204 on helix 5. The epinephrine-beta2AR structure shows dynamic flexibility with small, piston-like movements of helices 3 and 6 and transient interhelical hydrogen bonding between Ser-165 on transmembrane 4 and Ser-207 on transmembrane 5. These couplings and motions may play a role in protein activation. The apo-beta2AR shows less dynamic flexibility, whereas the antagonist-beta2AR structure is quite rigid. This MD validation of the structure predictions for G protein-coupled receptors in explicit lipid and water suggests that these methods can be trusted for studying the mechanism of activation and the design of subtype-specific agonists and antagonists.

The role of extra-membranous inter-helical loops in helix–helix interactions

The effect of a short loop connecting two transmembrane alpha-helices was studied using molecular dynamics simulations. Helices F and G from bacteriorhodopsin and two corresponding polyalanine helices were embedded in octane and POPC membranes in a transmembrane configuration both with and without the inter-helical loop. The results indicate that the membrane environment and the sequence of the loop are more influential on the dynamics and structure of the motif than the presence of a loop as such, at least for the time-scales investigated. The four residues in the FG loop are stabilized by four hydrogen bonds. These hydrogen bonds are not present in the polyalanine loop, causing it to be more flexible than the FG loop. This effect was observed independently of the protein environment, stressing the importance of the sequence. The structural analysis indicates that the loop has weak stabilizing properties in all environments. The stabilization due to the presence of the loop was strongest in a simulation of the FG fragment in a membrane-mimetic octane slab. In the simulations of the helix-loop-helix motif embedded in an explicit lipid bilayer model, the lipid bilayer interface compensates to a large extent for the absence of the loop.

Toward theoretical analysis of long-range proton transfer kinetics in biomolecular pumps.

Motivated by the long-term goal of theoretically analyzing long-range proton transfer (PT) kinetics in biomolecular pumps, researchers made a number of technical developments in the framework of quantum mechanics-molecular mechanics (QM/MM) simulations. A set of collective reaction coordinates is proposed for characterizing the progress of long-range proton transfers; unlike previous suggestions, the new coordinates can describe PT along highly nonlinear three-dimensional pathways. Calculations using a realistic model of carbonic anhydrase demonstrated that adiabatic mapping using these collective coordinates gives reliable energetics and critical geometrical parameters as compared to minimum energy path calculations, which suggests that the new coordinates can be effectively used as reaction coordinate in potential of mean force calculations for long-range PT in complex systems. In addition, the generalized solvent boundary potential was implemented in the QM/MM framework for rectangular geometries, which is useful for studying reactions in membrane systems. The resulting protocol was found to produce water structure in the interior of aquaporin consistent with previous studies including a much larger number of explicit solvent and lipid molecules. The effect of electrostatics for PT through a membrane protein was also illustrated with a simple model channel embedded in different dielectric continuum environments. The encouraging results observed so far suggest that robust theoretical analysis of long-range PT kinetics in biomolecular pumps can soon be realized in a QM/MM framework.

Folding is not required for bilayer insertion: Replica exchange simulations of an α‐helical peptide with an explicit lipid bilayer

We implement the replica exchange molecular dynamics algorithm to study the interactions of a model peptide (WALP-16) with an explicitly represented DPPC membrane bilayer. We observe the spontaneous, unbiased insertion of WALP-16 into the DPPC bilayer and its folding into an alpha-helix with a transbilayer orientation. The free energy surface suggests that the insertion of the peptide into the DPPC bilayer precedes secondary structure formation. Although the peptide has some propensity to form a partially helical structure in the interfacial region of the DPPC/water system, this state is not a productive intermediate but rather an off-pathway trap for WALP-16 insertion. Equilibrium simulations show that the observed insertion/folding pathway mirrors the potential of mean force (PMF). Calculation of the enthalpic and entropic contributions to this PMF show that the surface bound conformation of WALP-16 is significantly lower in energy than other conformations, and that the insertion of WALP-16 into the bilayer without regular secondary structure is enthalpically unfavorable by 5-10 kcal/mol/residue. The observed insertion/folding pathway disagrees with the dominant conceptual model, which is that a surface-bound helix is an obligatory intermediate for the insertion of alpha-helical peptides into lipid bilayers. In our simulations, the observed insertion/folding pathway is favored because of a large (>100 kcal/mol) increase in system entropy that occurs when the unstructured WALP-16 peptide enters the lipid bilayer interior. The insertion/folding pathway that is lowest in free energy depends sensitively on the near cancellation of large enthalpic and entropic terms. This suggests the possibility that intrinsic membrane peptides may have a diversity of insertion/folding behaviors depending on the exact system of peptide and lipid under consideration.

A generalized Born formalism for heterogeneous dielectric environments: Application to the implicit modeling of biological membranes

Reliable computer simulations of complex biological environments such as integral membrane proteins with explicit water and lipid molecules remain a challenging task. We propose a modification of the standard generalized Born theory of homogeneous solvent for modeling the heterogeneous dielectric environments such as lipid/water interfaces. Our model allows the representation of biological membranes in the form of multiple layered dielectric regions with dielectric constants that are different from the solute cavity. The proposed new formalism is shown to predict the electrostatic component of solvation free energy with a relative error of 0.17% compared to exact finite-difference solutions of the Poisson equation for a transmembrane helix test system. Molecular dynamics simulations of melittin and bacteriorhodopsin are carried out and performed over 10 ns and 7 ns of simulation time, respectively. The center of melittin along the membrane normal in these stable simulations is in excellent agreement with the relevant experimental data. Simulations of bacteriorhodopsin started from the experimental structure remained stable and in close agreement with experiment. We also examined the free energy profiles of water and amino acid side chain analogs upon membrane insertion. The results with our implicit membrane model agree well with the experimental transfer free energy data from cyclohexane to water as well as explicit solvent simulations of water and selected side chain analogs.

Agonist-Induced Conformational Changes in Thyrotropin-Releasing Hormone Receptor Type I: Disulfide Cross-Linking and Molecular Modeling Approaches

The conformational changes at the cytoplasmic ends of transmembrane helices 5 and 6 (TMH5 and TMH6) of thyrotropin-releasing hormone (TRH) receptor type I (TRH-R1) during activation were analyzed by cysteine-scanning mutagenesis followed by disulfide cross-linking and molecular modeling. Sixteen double cysteine mutants were constructed by substitution of one residue at the cytoplasmic end of TMH5 and the other at that of TMH6. The cross-linking experiments indicate that four mutants, Q263C/G212C, Q263C/Y211C, T265C/G212C, and T265C/Y211C, exhibited disulfide bond formation that was sensitive to TRH occupancy. We refined our previous TRH-R1 models by embedding them into a hydrated explicit lipid bilayer. Molecular dynamics simulations of the models, as well as in silico double cysteine models, generated trajectories that were in agreement with experimental results. Our findings suggest that TRH binding induces a separation of the cytoplasmic ends of TMH5 and TMH6 and a rotation of TMH6. These changes likely increase the surface accessible area at the juxtamembrane region of intracellular loop 3 that could promote interactions between G proteins and key residues within the receptor.

Molecular modeling of the three-dimensional structure of dopamine 3 (D3) subtype receptor: discovery of novel and potent D3 ligands through a hybrid pharmacophore- and structure-based database searching approach.

The dopamine 3 (D3) subtype receptor has been implicated in several neurological conditions, and potent and selective D3 ligands may have therapeutic potential for the treatment of drug addiction, Parkinson's disease, and schizophrenia. In this paper, we report computational homology modeling of the D3 receptor based upon the high-resolution X-ray structure of rhodopsin, extensive structural refinement in the presence of explicit lipid bilayer and water environment, and validation of the refined D3 structural models using experimental data. We further describe the development, validation, and application of a hybrid computational screening approach for the discovery of several classes of novel and potent D3 ligands. This computational approach employs stepwise pharmacophore and structure-based searching of a large three-dimensional chemical database for the identification of potential D3 ligands. The obtained hits are then subjected to structural novelty screening, and the most promising compounds are tested in a D3 binding assay. Using this approach we identified four compounds with K(i) values better than 100 nM and eight compounds with K(i) values better than 1 microM out of 20 compounds selected for testing in the D3 receptor binding assay. Our results suggest that the D3 structural models obtained from this study may be useful for the discovery and design of novel and potent D3 ligands. Furthermore, the employed hybrid approach may be more effective for lead discovery from a large chemical database than either pharmacophore-based or structure-based database screening alone.

Structure of gramicidin a in a lipid bilayer environment determined using molecular dynamics simulations and solid-state NMR data.

Two different high-resolution structures recently have been proposed for the membrane-spanning gramicidin A channel: one based on solid-state NMR experiments in oriented phospholipid bilayers (Ketchem, R. R.; Roux, B.; Cross, T. A. Structure 1997, 5, 1655-1669; Protein Data Bank, PDB:1MAG); and one based on two-dimensional NMR in detergent micelles (Townsley, L. E.; Tucker, W. A.; Sham, S.; Hinton, J. F. Biochemistry 2001, 40, 11676-11686; PDB:1JNO). Despite overall agreement, the two structures differ in peptide backbone pitch and the orientation of several side chains; in particular that of the Trp at position 9. Given the importance of the peptide backbone and Trp side chains for ion permeation, we undertook an investigation of the two structures using molecular dynamics simulation with an explicit lipid bilayer membrane, similar to the system used for the solid-state NMR experiments. Based on 0.1 micros of simulation, both backbone structures converge to a structure with 6.25 residues per turn, in agreement with X-ray scattering, and broad agreement with SS backbone NMR observables. The side chain of Trp 9 is mobile, more so than Trp 11, 13, and 15, and undergoes spontaneous transitions between the orientations in 1JNO and 1MAG. Based on empirical fitting to the NMR results, and umbrella sampling calculations, we conclude that Trp 9 spends 80% of the time in the 1JNO orientation and 20% in the 1MAG orientation. These results underscore the utility of molecular dynamics simulations in the analysis and interpretation of structural information from solid-state NMR.

Computational studies of the gramicidin channel.

Ion channels are highly specific membrane-spanning protein structures which serve to facilitate the passage of selected ions across the lipid barrier. In the past decade, molecular dynamics simulations based on atomic models and realistic microscopic interactions with explicit solvent and membrane lipids have been used to gain insight into the function of these complex systems. These calculations have considerably expanded our view of ion permeation at the microscopic level. This Account will mainly focus on computational studies of the gramicidin A channel, one of the simplest and best characterized molecular pore.