The all-atom optimized potentials for liquid simulations (OPLS-AA) force field is a popular force field for simulating biomolecules. However, the current OPLS parameters for hydrocarbons developed using short alkanes cannot reproduce the liquid properties of long alkanes in molecular dynamics simulations. Therefore, the extension of OPLS-AA to (phospho)lipid molecules required for the study of biological membranes was hampered in the past. Here, we optimized the OPLS-AA force field for both short and long hydrocarbons. Following the framework of the OPLS-AA parametrization, we refined the torsional parameters for hydrocarbons by fitting to the gas-phase ab initio energy profiles calculated at the accurate MP2/aug-cc-pVTZ theory level. Additionally, the depth of the Lennard-Jones potential for methylene hydrogen atoms was adjusted to reproduce the densities and the heats of vaporization of alkanes and alkenes of different lengths. Optimization of partial charges finally allowed to reproduce the gel-to-liquid-phase transition temperature for pentadecane and solvation free energies. It is shown that the optimized parameter set (L-OPLS) yields improved hydrocarbon diffusion coefficients, viscosities, and gauche-trans ratios. Moreover, its applicability for lipid bilayer simulations is shown for a GMO bilayer in its liquid-crystalline phase.

In all species tested, except humans, biological differences between vitamins D2 and D3 are accepted as fact. To test the presumption of equivalence in humans, we compared the ability of equal molar quantities of vitamin D2 or D3 to increase serum 25-hydroxyvitamin D [25(OH)D], the measure of vitamin D nutrition. Subjects took 260 nmol (approximately 4000 IU) vitamin D2 (n=17) or vitamin D3 (n=55) daily for 14 d. 25(OH)D was assayed with a method that detects both the vitamin D2 and D3 forms. With vitamin D3, mean (+/-SD) serum 25(OH)D increased from 41.3+/-17.7 nmol/L before to 64.6+/-17.2 nmol/L after treatment. With vitamin D2, the 25(OH)D concentration went from 43.7+/-17.7 nmol/L before to 57.4+/-13.0 nmol/L after. The increase in 25(OH)D with vitamin D3 was 23.3+/-15.7 nmol/L, or 1.7 times the increase obtained with vitamin D2 (13.7+/-11.4 nmol/L; P=0.03). There was an inverse relation between the increase in 25(OH)D and the initial 25(OH)D concentration. The lowest 2 tertiles for basal 25(OH)D showed larger increases in 25(OH)D: 30.6 and 25.5 nmol/L, respectively, for the first and second tertiles. In the highest tertile [25(OH)D >49 nmol/L] the mean increase in 25(OH)D was 13.3 nmol/L (P < 0.03 for comparison with each lower tertile). Although the 1.7-times greater efficacy for vitamin D3 shown here may seem small, it is more than what others have shown for 25(OH)D increases when comparing 2-fold differences in vitamin D3 dose. The assumption that vitamins D2 and D3 have equal nutritional value is probably wrong and should be reconsidered.

Phospholipid force fields are of ample importance for the simulation of artificial bilayers, membranes, and also for the simulation of integral membrane proteins. Here, we compare the two most applied atomic force fields for phospholipids, the all-atom CHARMM27 and the united atom Berger force field, with a newly developed all-atom generalized AMBER force field (GAFF) for dioleoylphosphatidylcholine molecules. Only the latter displays the experimentally observed difference in the order of the C2 atom between the two acyl chains. The interfacial water dynamics is smoothly increased between the lipid carbonyl region and the bulk water phase for all force fields; however, the water order and with it the electrostatic potential across the bilayer showed distinct differences between the force fields. Both Berger and GAFF underestimate the lipid self-diffusion. GAFF offers a consistent force field for the atomic scale simulation of biomembranes.

By means of molecular dynamics simulations with an all-atom force field, we investigated the affinities of alkali cations and halide anions for the dioleoylphosphatidylcholine lipid membrane in aqueous salt solutions. In addition, changes in phospholipid lateral diffusion and in headgroup mobility upon adding NaCl were observed using fluorescence spectroscopy. The simulations revealed that sodium is attracted to the headgroup region with its concentration being maximal in the vicinity of the phosphate groups. Potassium and cesium, however, do not preferentially adsorb to the membrane. Similarly, halide anions do not exhibit a strong affinity for the lipid headgroups but merely compensate for the positive charge of the sodium countercations. Nevertheless, larger halides such as bromide and iodide penetrate deeper into the headgroup region toward the boundary with the hydrophobic alkyl chain, this effect being likely underestimated within the present nonpolarizable force field. Addition of alkali halide salts modifies physical properties of the bilayer including the electronic density profiles, the electrostatic potential, and the area per lipid headgroup.

We developed a SVM-trained, quantitative matrix-based method for the prediction of MHC class I binding peptides, in which the features of the scoring matrix are energy terms retrieved from molecular dynamics simulations. At the same time we used the equilibrated structures obtained from the same simulations in a simple and efficient docking procedure. Our method consists of two steps: First, we predict potential binders from sequence data alone and second, we construct protein-peptide complexes for the predicted binders. So far, we tested our approach on the HLA-A0201 allele. We constructed two prediction models, using local, position-dependent (DynaPred(POS)) and global, position-independent (DynaPred) features. The former model outperformed the two sequence-based methods used in our evaluation; the latter shows a much higher generalizability towards other alleles than the position-dependent models. The constructed peptide structures can be refined within seconds to structures with an average backbone RMSD of 1.53 A from the corresponding experimental structures.

Protein-ligand docking is an essential step in modern drug discovery process. The challenge here is to accurately predict and e±ciently optimize the position and orientation of ligands in the binding pocket of a target protein. In this paper, we present a new method called PSOVina which combined the particle swarm optimization (PSO) algorithm with the e±cient BroydenFletcher-Goldfarb-Shannon (BFGS) local search method adopted in AutoDock Vina to tackle the conformational search problem in docking. Using a diverse data set of 201 protein-ligand complexes from the PDBbind database and a full set of ligands and decoys for four representative targets from the directory of useful decoys (DUD) virtual screening data set, we assessed the docking performance of PSOVina in comparison to the original Vina program. Our results showed that PSOVina achieves a remarkable execution time reduction of 51-60% without compromising the prediction accuracies in the docking and virtual screening experiments. This improvement in time e±ciency makes PSOVina a better choice of a docking tool in large-scale protein-ligand docking applications. Our work lays the foundation for the future development of swarm-based algorithms in molecular docking programs. PSOVina is freely available to noncommercial users at http://cbbio.cis.umac.mo.

Trans-soluble N-ethylmaleimide-sensitive factor attachment protein (SNAP) receptor (SNARE) complexes formed between the SNARE motifs of synaptobrevin II, SNAP-25, and syntaxin play an essential role in Ca 2ϩ -regulated exocytosis. Apart from the well studied interactions of the SNARE domains, little is known about the functional relevance of other evolutionarily conserved structures in the SNARE proteins. Here, we show that substitution of two highly conserved tryptophan residues within the juxtamembrane domain (JMD) of the vesicular SNARE Synaptobrevin II (SybII) profoundly impairs priming of granules in mouse chromaffin cells without altering catecholamine release from single vesicles. Using molecular dynamic simulations of membrane-embedded SybII, we show that Trp residues of the JMD influence the electrostatic surface potential by controlling the position of neighboring lysine and arginine residues at the membrane-water interface. Our observations indicate a decisive role of the tryptophan moiety of SybII in keeping the vesicles in the release-ready state and support a model wherein tryptophan-mediated protein-lipid interactions assist in bridging the apposing membranes before fusion.

No abstract