Molecular dynamics simulations allow a direct study of the structure and dynamics of membrane proteins and lipids. We describe the behavior of aromatic residues and lipid properties in POPE and POPC ...bilayer models with the Escherichia coli OmpF trimer, single alamethicin and Influenza M2 helices, 4-helix M2 bundles, and two alamethicin 6-helix channel models. The total simulation time is over 24 ns, of systems containing solvent, protein, and between 104 and 318 lipids. Various types of adjustment between lipids and proteins occur, depending on the size of the protein and the degree of hydrophobic mismatch between lipid and protein. Single helices cause little measurable effect on nearby lipids whereas the 4-helix bundles, 6-helix channel models, and OmpF cause a significant lowering of order parameters in nearby lipid chains, an increased difference between odd and even chain dihedrals in the magnitude of the trans dihedral fractions and dihedral transition rates, and in most cases a decreased gauche population and a decrease in bilayer thickness. An increased tilt of the lipid chains near the proteins can account for most of the observed decrease in order parameters. The orientation of tryptophans and tyrosines on the outside of the proteins is determined by packing at the protein exterior and non-specific hydrogen bonding with lipids and solvent. The tyrosines in the broad bands that delimit the hydrophobic exterior of OmpF show little change in orientation over one nanosecond. Their rings are oriented predominantly perpendicular to the bilayer plane, with the hydroxyl group pointing toward the lipid−water interface. Phenylalanines in OmpF, alamethicin, and Influenza M2 are more mobile and assume a variety of orientations.
Lipid peroxidation is an important part of the pathological pathway of membrane damage in membranes that have high levels of polyunsaturated fatty acids such as linoleic, linolenic, arachidonic, and ...docosahexaenoic acids. Neural membranes are particularly rich in polyunsaturated acids and such damage is implicated in neurological diseases, such as Alzheimer's disease. To obtain a bilayer model that represents the property of susceptibility to lipid peroxidation, we carried out molecular dynamics (MD) simulations of a bilayer of 1-palmitoyl-2-linoleyl-sn-glycero-3-phosphatidylcholine (PLPC). Parameters for the torsional potentials of the cis,cis-Δ9,12 bis-allylic region of the linoleate chain were fitted to the results of high-level ab initio calculations on model compounds. The MD simulations of the bilayer provided the structural properties of the system and show that the unsaturation induces disorder and affects the physical properties of the membrane.
The uptake of nutrients is essential for the survival of bacterial cells. Many specialized systems have evolved, such as the maltose-dependent ABC transport system that transfers oligosaccharides ...through the cytoplasmic membrane. The maltose/maltodextrin-binding protein (MBP) serves as an initial high-affinity binding component in the periplasm that delivers the bound sugar into the cognate ABC transporter MalFGK2. We have investigated the domain motions induced by the binding of the ligand maltotriose into the binding cleft using molecular dynamics simulations. We find that MBP is predominantly in the open state without ligand and in the closed state with ligand bound. Oligosaccharide binding induces a closure motion (30.0° rotation), whereas ligand removal leads to domain opening (32.6° rotation) around a well-defined hinge affecting key areas relevant for chemotaxis and transport. Our simulations suggest that a “hook-and-eye” motif is involved in the binding. A salt bridge between Glu-111 and Lys-15 forms that effectively locks the protein-ligand complex in a semiclosed conformation inhibiting any further opening and promoting complete closure. This previously unrecognized feature seems to secure the ligand in the binding site and keeps MBP in the closed conformation and suggests a role in the initial steps of substrate transport.
Molecular dynamics simulations of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal significant differences in interactions of the selectivity filter with K
+ compared ...with Na
+ ions. K
+ ions and water molecules within the filter undergo concerted single-file motion in which they translocate between adjacent sites within the filter on a nanosecond timescale. In contrast, Na
+ ions remain bound to sites within the filter and do not exhibit translocation on a nanosecond timescale. Furthermore, entry of a K
+ ion into the filter from the extracellular mouth is observed, whereas this does not occur for a Na
+ ion. Whereas K
+ ions prefer to sit within a cage of eight oxygen atoms of the filter, Na
+ ions prefer to interact with a ring of four oxygen atoms plus two water molecules. These differences in interactions in the selectivity filter may contribute to the selectivity of KcsA for K
+ ions (in addition to the differences in dehydration energy between K
+ and Na
+) and the block of KcsA by internal Na
+ ions. In our simulations the selectivity filter exhibits significant flexibility in response to changes in ion/protein interactions, with a somewhat greater distortion induced by Na
+ than by K
+ ions.
Methodological issues in molecular dynamics (MD) simulations, such as the treatment of long-range electrostatic interactions or the type of pressure coupling, have important consequences for the ...equilibrium properties observed. We report a series of long (up to 150 ns) MD simulations of dipalmitoylphosphatidylcholine (DPPC) bilayers in which the methodology of simulation is systematically varied. Comparisons of simulations with truncation schemes, Ewald summations, and modified Coulomb interactions, either by shift functions or reaction field models, to describe long-range electrostatics point out the artifacts inherent in each of these methods and above all those of straight cutoff methods. We further show that bilayer properties are less sensitive to the details of the pressure-coupling algorithm and that an increased integration time step of 5 fs can be safely used in simulations of phosphatidylcholine lipid bilayers.
Alamethicin is an antimicrobial peptide that forms stable channels with well-defined conductance levels. We have used extended molecular dynamics simulations of alamethicin bundles consisting of 4, ...5, 6, 7, and 8 helices in a palmitoyl-oleolyl-phosphatidylcholine bilayer to evaluate and analyze channel models and to link the models to the experimentally measured conductance levels. Our results suggest that four helices do not form a stable water-filled channel and might not even form a stable intermediate. The lowest measurable conductance level is likely to correspond to the pentamer. At higher aggregation numbers the bundles become less symmetrical. Water properties inside the different-sized bundles are similar. The hexamer is the most stable model with a stability comparable with simulations based on crystal structures. The simulation was extended from 4 to 20
ns or several times the mean passage time of an ion. Essential dynamics analyses were used to test the hypothesis that correlated motions of the helical bundles account for high-frequency noise observed in open channel measurements. In a 20-ns simulation of a hexameric alamethicin bundle, the main motions are those of individual helices, not of the bundle as a whole. A detailed comparison of simulations using different methods to treat long-range electrostatic interactions (a twin range cutoff, Particle Mesh Ewald, and a twin range cutoff combined with a reaction field correction) shows that water orientation inside the alamethicin channels is sensitive to the algorithms used. In all cases, water ordering due to the protein structure is strong, although the exact profile changes somewhat. Adding an extra 4-nm layer of water only changes the water ordering slightly in the case of particle mesh Ewald, suggesting that periodicity artifacts for this system are not serious.
Peptide−membrane interactions are important for understanding the binding, partitioning, and folding of membrane proteins; the activity of antimicrobial and fusion peptides; and a number of other ...processes. We describe molecular dynamics simulations (10−25 ns) of two pentapeptides Ace-WLXLL (with X = Arg or Lys side chain) (White, S. H., and Wimley, W.C. (1996) Nat. Struct. Biol. 3, 842−848) in water and three different membrane mimetic systems: (i) a water/cyclohexane interface, (ii) water-saturated octanol, and (iii) a solvated dioleoylphosphatidylcholine bilayer. A salt bridge is found between the protonated Arg or Lys side chains with the carboxyl terminus at the three interfaces. In water/cyclohexane, the salt bridge is most exposed to the water phase and least stable. In water/octanol and the lipid bilayer systems, the salt bridge once formed persists throughout the simulations. In the lipid bilayer, the salt bridge is more stable when the peptide penetrates deeper into the bilayer. In one of two peptides, a cation−π interaction between the Arg and the Trp side chains is stable in the lipid bilayer for about 15 ns before breaking. In all cases, the conformations of the peptides are restricted by their presence at the interface and can be assigned to a few major conformational clusters. Side chains facing the water phase are most mobile. In the lipid bilayer, the peptides remain in the interface area, where they overlap with the carbonyl area of the lipid bilayer and perturb the local density profile of the bilayer. The tryptophan side chain remains in the water−lipid interface, where it interacts with the lipid choline group and forms hydrogen bonds with the ester carbonyl of the lipid and with water in the interface.
The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing ...the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-
sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively “frozen” in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (Randa et al., 1999,
Biophys. J. 77:2400–2410).