Interactions of Monovalent and Divalent Cations at Palmitoyl-Oleoyl-Phosphatidylcholine Interface.

Li+ is a biologically active and medically important cation. Experiments show that Li+ modulates some phospholipid bilayer properties in a manner similar to divalent cations, rather than other monovalent cations. We previously performed a comparative simulation study of the interaction of several monovalent cations with palmitoyl-oleoyl-phosphatidylcholine bilayers and reported that Li+ exhibited the highest association with lipids and formed a unique tetrahedral coordinated structure with lipid head groups. Here we extend these studies to two biologically important divalent cations, Mg2+ and Ca2+, and observe that, just like monovalent cations, Mg2+ and Ca2+ reduce bilayer areas and increase chain order. Bilayer area changes induced by cations are strongly correlated with the amount of charge inside the headgroup region; however, Mg2+ and Li+ are clear outliers. At the same time though, Mg2+ adsorption in the bilayer is the smallest among all cations, which is in contrast to Li+ that binds strongly to lipids. In fact, in contrast to all other cations, Mg2+ remains fully hydrated in the lipid headgroup region. However, Li+ and Mg2+ share high overlap between their inner-shell coordination topologies. This suggests that Li+ can structurally replace Mg2+, which is bound to other biomolecules with up to fourfold coordination, provided such replacement is energetically feasible. We compute structural topologies and compare them quantitatively using a new weighted-graphs-based method. Finally, we find that the specificity of cation interaction with lipid head groups exhibit consistent trend with the solvation shell energetics of ions in lipid headgroup and bulk water regions.

Towards a Unified Understanding of Lithium Action in Basic Biology and its Significance for Applied Biology.

Lithium has literally been everywhere forever, since it is one of the three elements created in the Big Bang. Lithium concentration in rocks, soil, and fresh water is highly variable from place to place, and has varied widely in specific regions over evolutionary and geologic time. The biological effects of lithium are many and varied. Based on experiments in which animals are deprived of lithium, lithium is an essential nutrient. At the other extreme, at lithium ingestion sufficient to raise blood concentration significantly over 1 mM/, lithium is acutely toxic. There is no consensus regarding optimum levels of lithium intake for populations or individuals-with the single exception that lithium is a generally accepted first-line therapy for bipolar disorder, and specific dosage guidelines for sufferers of that condition are generally agreed on. Epidemiological evidence correlating various markers of social dysfunction and disease vs. lithium level in drinking water suggest benefits of moderately elevated lithium compared to average levels of lithium intake. In contrast to other biologically significant ions, lithium is unusual in not having its concentration in fluids of multicellular animals closely regulated. For hydrogen ions, sodium ions, potassium ions, calcium ions, chloride ions, and magnesium ions, blood and extracellular fluid concentrations are closely and necessarily regulated by systems of highly selective channels, and primary and secondary active transporters. Lithium, while having strong biological activity, is tolerated over body fluid concentrations ranging over many orders of magnitude. The lack of biological regulation of lithium appears due to lack of lithium-specific binding sites and selectivity filters. Rather lithium exerts its myriad physiological and biochemical effects by competing for macromolecular sites that are relatively specific for other cations, most especially for sodium and magnesium. This review will consider what is known about the nature of this competition and suggest using and extending this knowledge towards the goal of a unified understanding of lithium in biology and the application of that understanding in medicine and nutrition.

Effects of Lithium and Other Monovalent Ions on Palmitoyl Oleoyl Phosphatidylcholine Bilayer.

Interactions of monovalent salts with lipid membranes are explored with molecular dynamics (MD) simulations. The simulations included the monovalent ions Na+ and K+, for their importance in physiology, Li+ for its small size and importance in several medical conditions including bipolar disorder, and Rb+ for its large size. All simulations included Cl- as counterions. One bilayer was simulated without salt as a control. Palmitoyl oleoyl phosphatidylcholine (POPC) bilayers experienced reductions in area per lipid with the addition of salt; the smaller the ion the smaller the area, with the exception of Li+. Li+ exhibited unique binding affinities between phosphates and sn-2 carbonyls that lowered the order of the top part of sn-2 chain, which increased the area per lipid, compared to other ionic simulations. Further, we observe that monovalent salts alter bilayer properties through structural changes and not so much through the changes in surface potential.

Impact of sterol tilt on membrane bending rigidity in cholesterol and 7DHC-containing DMPC membranes.

Cholesterol is so essential to the proper function of mammalian cell membranes that even strikingly small inborn errors in cholesterol synthesis can be devastating. Here we combine molecular dynamics simulations with small angle x-ray diffraction experiments to compare mixed sterol/DMPC membranes over a wide range of sterol compositions for two types of sterols: cholesterol and its immediate metabolic precursor 7DHC, that differs from cholesterol by one double bond. We find that while most membrane properties are only slightly affected by the replacement of one sterol by the other, the tilt degree of freedom, as gauged by the tilt modulus, is significantly larger for cholesterol than for 7DHC over a large range of concentrations. In silico mutations of one sterol into the other further support these findings. Moreover, bending rigidities calculated from simulations and estimated in experiments show that cholesterol stiffens membranes to a larger extent than 7DHC. We discuss the possible mechanistic link between sterol tilt and the way it impacts the membrane mechanical properties, and comment on how this link may shed light on the way replacement of cholesterol by 7DHC leads to disease.

Improving predicted protein loop structure ranking using a Pareto-optimality consensus method.

BACKGROUND: Accurate protein loop structure models are important to understand functions of many proteins. Identifying the native or near-native models by distinguishing them from the misfolded ones is a critical step in protein loop structure prediction. RESULTS: We have developed a Pareto Optimal Consensus (POC) method, which is a consensus model ranking approach to integrate multiple knowledge- or physics-based scoring functions. The procedure of identifying the models of best quality in a model set includes: 1) identifying the models at the Pareto optimal front with respect to a set of scoring functions, and 2) ranking them based on the fuzzy dominance relationship to the rest of the models. We apply the POC method to a large number of decoy sets for loops of 4- to 12-residue in length using a functional space composed of several carefully-selected scoring functions: Rosetta, DOPE, DDFIRE, OPLS-AA, and a triplet backbone dihedral potential developed in our lab. Our computational results show that the sets of Pareto-optimal decoys, which are typically composed of approximately 20% or less of the overall decoys in a set, have a good coverage of the best or near-best decoys in more than 99% of the loop targets. Compared to the individual scoring function yielding best selection accuracy in the decoy sets, the POC method yields 23%, 37%, and 64% less false positives in distinguishing the native conformation, indentifying a near-native model (RMSD < 0.5A from the native) as top-ranked, and selecting at least one near-native model in the top-5-ranked models, respectively. Similar effectiveness of the POC method is also found in the decoy sets from membrane protein loops. Furthermore, the POC method outperforms the other popularly-used consensus strategies in model ranking, such as rank-by-number, rank-by-rank, rank-by-vote, and regression-based methods. CONCLUSIONS: By integrating multiple knowledge- and physics-based scoring functions based on Pareto optimality and fuzzy dominance, the POC method is effective in distinguishing the best loop models from the other ones within a loop model set.

An improved united atom force field for simulation of mixed lipid bilayers.

We introduce a new force field (43A1-S3) for simulation of membranes by the Gromacs simulation package. Construction of the force fields is by standard methods of electronic structure computations for bond parameters and charge distribution and specific volumes and heats of vaporization for small-molecule components of the larger lipid molecules for van der Waals parameters. Some parameters from the earlier 43A1 force field are found to be correct in the context of these calculations, while others are modified. The validity of the force fields is demonstrated by correct replication of X-ray form factors and NMR order parameters over a wide range of membrane compositions in semi-isotropic NTP 1 atm simulations. 43-A1-S3 compares favorably with other force fields used in conjunction with the Gromacs simulation package with respect to the breadth of phenomena that it accurately reproduces.

Automated optimization of water-water interaction parameters for a coarse-grained model.

We have developed an automated parameter optimization software framework (ParOpt) that implements the Nelder-Mead simplex algorithm and applied it to a coarse-grained polarizable water model. The model employs a tabulated, modified Morse potential with decoupled short- and long-range interactions incorporating four water molecules per interaction site. Polarizability is introduced by the addition of a harmonic angle term defined among three charged points within each bead. The target function for parameter optimization was based on the experimental density, surface tension, electric field permittivity, and diffusion coefficient. The model was validated by comparison of statistical quantities with experimental observation. We found very good performance of the optimization procedure and good agreement of the model with experiment.

The effect of curcumin on the stability of Aß dimers.

Aß oligomers are potential targets for the diagnosis and therapy of Alzheimer's disease (AD). On the other hand, the molecule curcumin has been shown to possess significant therapeutic potential in many areas. In this paper, we use all-atom explicit solvent molecular dynamics simulations to study the effect of curcumin on the stability of Aß amyloid protein oligomers. We observed that curcumin decreases the ß-sheet secondary structural content within the Aß oligomers without reducing the contacts between the monomers. The breaking of the ß-sheet is found to be preceded by a deformation of the ß-sheet structure due to hydrophobic interaction from the nearby curcumin. Furthermore, the π-stacking interaction between curcumin (keto ring and enol ring) and the aromatic residues of Aß, which exists throughout the simulations, has also contributed to the diminishing of the ß-sheet structure. Our analysis of the underwrapped amide-carbonyl hydrogen bonds reveals several stable dehydrons of the oligomer, especially the dehydron pair 34L and 41I, which curcumin tends to hover over. We have examined the paths of curcumin on the Aß proteins and determined the common routes where curcumin lingers as it traverses around the Aß. In consequence, our study has provided a detailed interaction picture between curcumin and the Aß oligomers.

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study.

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.

Amyloid ß peptides aggregation in a mixed membrane bilayer: a molecular dynamics study.

The aggregation of amyloid ß peptides resulting in neurotoxic oligomers is an important but yet mysterious process in Alzheimer's disease development. Molecular dynamics simulations were performed to investigate the self-assembly of three full-length amyloid peptides in the zwitterionic dipalmitoylphosphatidylcholine and cholesterol mixed lipid bilayer. During the 1000 ns simulation, the residues 1-27 were found to interact preferentially with the lipid-aqueous interface region, while residues 28-42 show an inclination to remain inside the bilayer hydrophobic tail region. The interaction between peptides and lipids has facilitated the association of Aß peptides. However, the interaction between cholesterol and peptides is inversely correlated with the extent of the peptide-peptide interactions. Our simulation has uncovered the formation of a short segment of parallel ß-sheet between two peptide chains. In another chain, the N- and C-termini came close to each other. All the structural transitions indicate that our simulation has caught a glimpse of the complicated peptide oligomerization process. The full understanding of the underlying mechanism still requires further experimental and theoretical studies.

Molecular dynamic simulation study of cholesterol and conjugated double bonds in lipid bilayers.

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is believed to have anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be associated with anti-obesity effects. In this paper we extend earlier molecular dynamics (MD) simulations of pure CLA-phosphatidylcholine bilayers to investigate the comparative effects of cholesterol on bilayers composed of the two respective isomers. Simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed in which, for each isomer, the simulated bilayers contained 10% and 30% cholesterol (Chol). From MD trajectories we calculate and compare structural properties of the bilayers, including areas per molecule, thickness of bilayers, tilt angle of cholesterols, order parameter profiles, and one and two-dimensional radial distribution function (RDF), as functions of Chol concentration. While the structural effect of cholesterol is approximately the same for both isomers, we find differences at an atomistic level in order parameter profiles and in two-dimensional radial distribution functions.

A Coarse-Grained Model Based on Morse Potential for Water and n-Alkanes.

In order to extend the time and distance scales of molecular dynamics simulations, it is essential to create accurate coarse-grained force fields, in which each particle contains several atoms. Coarse-grained force fields that utilize the Lennard-Jones potential form for pairwise nonbonded interactions have been shown to suffer from serious inaccuracy, notably with respect to describing the behavior of water. In this paper, we describe a coarse-grained force field for water, in which each particle contains four water molecules, based on the Morse potential form. By molecular dynamics simulations, we show that our force field closely replicates important water properties. We also describe a Morse potential force field for alkanes and a simulation method for alkanes in which individual particles may have variable size, providing flexibility in constructing complex molecules comprised partly or solely of alkane groups. We find that, in addition to being more accurate, the Morse potential also provides the ability to take larger time steps than the Lennard-Jones, because the short distance repulsion potential profile is less steep. We suggest that the Morse potential form should be considered as an alternative for the Lennard-Jones form for coarse-grained molecular dynamics simulations.

SIMULATION OF ION CONDUCTION IN α-HEMOLYSIN NANOPORES WITH COVALENTLY ATTACHED ß-CYCLODEXTRIN BASED ON BOLTZMANN TRANSPORT MONTE CARLO MODEL.

Ion channels, as natures' solution to regulating biological environments, are particularly interesting to device engineers seeking to understand how natural molecular systems realize device-like functions, such as stochastic sensing of organic analytes. What's more, attaching molecular adaptors in desired orientations inside genetically engineered ion channels, enhances the system functionality as a biosensor. In general, a hierarchy of simulation methodologies is needed to study different aspects of a biological system like ion channels. Biology Monte Carlo (BioMOCA), a three-dimensional coarse-grained particle ion channel simulator, offers a powerful and general approach to study ion channel permeation. BioMOCA is based on the Boltzmann Transport Monte Carlo (BTMC) and Particle-Particle-Particle-Mesh (P(3)M) methodologies developed at the University of Illinois at Urbana-Champaign. In this paper, we have employed BioMOCA to study two engineered mutations of α-HL, namely (M113F)(6)(M113C-D8RL2)(1)-ß-CD and (M113N)(6)(T117C-D8RL3)(1)-ß-CD. The channel conductance calculated by BioMOCA is slightly higher than experimental values. Permanent charge distributions and the geometrical shape of the channels gives rise to selectivity towards anions and also an asymmetry in I-V curves, promoting a rectification largely for cations.

Cholesterol packing around lipids with saturated and unsaturated chains: a simulation study.

The fundamental role of cholesterol in the regulation of eukaryotic membrane structure is well-established. However the manner in which atomic level interactions between cholesterol and lipids, with varying degrees of chain unsaturation and polar groups, affect the overall structure and organization of the bilayer is only beginning to be understood. In this paper we describe a series of Molecular Dynamics simulations designed to provide new insights into lipid-cholesterol interactions as a function of chain unsaturation. We have run simulations of varying concentrations of cholesterol in dipalmitoyl phosphatidylcholine (DPPC), palmitoyl-oleyol phosphatidylcholine (POPC), and dioleyol phosphatidylcholine (DOPC) bilayers. Structural analysis of the simulations reveals both atomistic and systemic details of the interactions and are presented here. In particular, we find that the minimum partial molecular area of cholesterol occurs in POPC-Chol mixtures implying the most favorable packing. Physically, this appears to be related to the fact that the two faces of the cholesterol molecule are different from each other and that the steric cross section of cholesterol molecules drops sharply near the small chain tails.

Cholesterol surrogates: a comparison of cholesterol and 16:0 ceramide in POPC bilayers.

Experimental evidence indicates that, under some circumstances, "surrogate" molecules may play the same role as cholesterol in ordering membrane lipids. The simplest molecule in this class is Ceramide. In this article, we describe atomic-level molecular dynamics simulations designed to shed light on this phenomenon. We run simulations of hydrated phosphoryl-oleoyl phosphatidylcholine (POPC) bilayers containing cholesterol, and containing ceramide, in concentrations ranging from 5% to 33%. We also perform a simulation of a pure POPC bilayer to verify the simulation force fields against experimental structural data for POPC. Our simulation data are in good agreement with experimental data for the partial molecular volumes, areas, form factors, and order parameters. These simulations suggest that ceramide and cholesterol have a very similar effect on the POPC bilayer, although ceramide is less effective in inducing order in the bilayer compared with cholesterol at the same concentrations.

The influence of amino acid protonation states on molecular dynamics simulations of the bacterial porin OmpF.

Several groups, including our own, have found molecular dynamics (MD) calculations to result in the size of the pore of an outer membrane bacterial porin, OmpF, to be reduced relative to its size in the x-ray crystal structure. At the narrowest portion of its pore, loop L3 was found to move toward the opposite face of the pore, resulting in decreasing the cross-section area by a factor of approximately 2. In an earlier work, we computed the protonation states of titratable residues for this system and obtained values different from those that had been used in previous MD simulations. Here, we show that MD simulations carried out with these recently computed protonation states accurately reproduce the cross-sectional area profile of the channel lumen in agreement with the x-ray structure. Our calculations include the investigation of the effect of assigning different protonation states to the one residue, D(127), whose protonation state could not be modeled in our earlier calculation. We found that both assumptions of charge states for D(127) reproduced the lumen size profile of the x-ray structure. We also found that the charged state of D(127) had a higher degree of hydration and it induced greater mobility of polar side chains in its vicinity, indicating that the apparent polarizability of the D(127) microenvironment is a function of the D(127) protonation state.

G3(MP2) study of the C3H6O+* isomers fragmented from 1,4-dioxane+*.

The fragmentation process of ionized 1,4-dioxane and the reactions between the C3H6O+* ions, one of the major fragments, and various reactants (including acetonitrile, formaldehyde, ethylene, and propene) have been studied experimentally with mass spectrometry. In the present work, G3(MP2) calculations were carried out to investigate these processes theoretically. In agreement with experiment, isomers CH3OCHCH2+* (1) and *CH2CH2OCH2+ (2) were found to be the C3H6O+* ions fragmented from ionized 1,4-dioxane, with 2 being the major product. The mechanisms of the formation of 1 and 2 were successfully established. In addition, the characteristic reactivities, as well as the corresponding reaction mechanisms, of both isomers were rationalized with the aid of calculations. Finally, a minor reaction between isomer 2 and propene was identified, and the presence of the product of this reaction was found to be useful in explaining the aforementioned mass spectrometric data.