New parameter sets of the GROMOS biomolecular force field, 54A7 and 54B7, are introduced. These parameter sets summarise some previously published force field modifications: The 53A6 helical propensities are corrected through new u/w torsional angle terms and a modification of the N-H, C=O repulsion, a new atom type for a charged -CH 3 in the choline moiety is added, the Na ? and Cl -ions are modified to reproduce the free energy of hydration, and additional improper torsional angle types for free energy calculations involving a chirality change are introduced. The new helical propensity modification is tested using the benchmark proteins hen egg-white lysozyme, fox1 RNA binding domain, chorismate mutase and the GCN4-p1 peptide. The stability of the proteins is improved in comparison with the 53A6 force field, and good agreement with a range of primary experimental data is obtained.

A large variety of antimicrobial peptides have been shown to act, at least in vitro, by poration of the lipid membrane. The nanometre size of these pores, however, complicates their structural characterization by experimental techniques. Here we use molecular dynamics simulations, to study the interaction of a specific class of antimicrobial peptides, melittin, with a dipalmitoylphosphatidylcholine bilayer in atomic detail. We show that transmembrane pores spontaneously form above a critical peptide to lipid ratio. The lipid molecules bend inwards to form a toroidally shaped pore but with only one or two peptides lining the pore. This is in strong contrast to the traditional models of toroidal pores in which the peptides are assumed to adopt a transmembrane orientation. We find that peptide aggregation, either prior or after binding to the membrane surface, is a prerequisite to pore formation. The presence of a stable helical secondary structure of the peptide, however is not. Furthermore, results obtained with modified peptides point to the importance of electrostatic interactions in the poration process. Removing the charges of the basic amino-acid residues of melittin prevents pore formation. It was also found that in the absence of counter ions pores not only form more rapidly but lead to membrane rupture. The rupture process occurs via a novel recursive poration pathway, which we coin the Droste mechanism.

Molecular dynamics simulations of the magainin MG-H2 peptide interacting with a model phospholipid membrane have been used to investigate the mechanism by which antimicrobial peptides act. Multiple copies of the peptide were randomly placed in solution close to the membrane. The peptide readily bound to the membrane, and above a certain concentration, the peptide was observed to cooperatively induce the formation of a nanometer-sized, toroidally shaped pore in the bilayer. In sharp contrast with the commonly accepted model of a toroidal pore, only one peptide was typically found near the center of the pore. The remaining peptides lay close to the edge of the pore, maintaining a predominantly parallel orientation with respect to the membrane.

A new force field for the simulation of dipalmitoylphosphatidylcholine (DPPC) in the liquid-crystalline, fluid phase at zero surface tension is presented. The structure of the bilayer with the area per lipid (0.629 nm(2); experiment 0.629-0.64 nm(2)), the volume per lipid (1.226 nm(3); experiment 1.229-1.232 nm(3)), and the ordering of the palmitoyl chains (order parameters) are all in very good agreement with experiment. Experimental electron density profiles are well reproduced in particular with regard to the penetration of water into the bilayer. The force field was further validated by simulating the spontaneous assembly of DPPC into a bilayer in water. Notably, the timescale on which membrane sealing was observed using this model appears closer to the timescales for membrane resealing suggested by electroporation experiments than previous simulations using existing models.

Methods to compute free energy differences between different states of a molecular system are reviewed with the aim of identifying their basic ingredients and their utility when applied in practice to biomolecular systems. A free energy calculation is comprised of three basic components: (i) a suitable model or Hamiltonian, (ii) a sampling protocol with which one can generate a representative ensemble of molecular configurations, and (iii) an estimator of the free energy difference itself. Alternative sampling protocols can be distinguished according to whether one or more states are to be sampled. In cases where only a single state is considered, six alternative techniques could be distinguished: (i) changing the dynamics, (ii) deforming the energy surface, (iii) extending the dimensionality, (iv) perturbing the forces, (v) reducing the number of degrees of freedom, and (vi) multi-copy approaches. In cases where multiple states are to be sampled, the three primary techniques are staging, importance sampling, and adiabatic decoupling. Estimators of the free energy can be classified as global methods that either count the number of times a given state is sampled or use energy differences. Or, they can be classified as local methods that either make use of the force or are based on transition probabilities. Finally, this overview of the available techniques and how they can be best used in a practical context is aimed at helping the reader choose the most appropriate combination of approaches for the biomolecular system, Hamiltonian and free energy difference of interest.

Signaling from JAK (Janus kinase) protein kinases to STAT (signal transducers and activators of transcription) transcription factors is key to many aspects of biology and medicine, yet the mechanism by which cytokine receptors initiate signaling is enigmatic. We present a complete mechanistic model for activation of receptor-bound JAK2, based on an archetypal cytokine receptor, the growth hormone receptor. For this, we used fluorescence resonance energy transfer to monitor positioning of the JAK2 binding motif in the receptor dimer, substitution of the receptor extracellular domains with Jun zippers to control the position of its transmembrane (TM) helices, atomistic modeling of TM helix movements, and docking of the crystal structures of the JAK2 kinase and its inhibitory pseudokinase domain with an opposing kinase-pseudokinase domain pair. Activation of the receptor dimer induced a separation of its JAK2 binding motifs, driven by a ligand-induced transition from a parallel TM helix pair to a left-handed crossover arrangement. This separation leads to removal of the pseudokinase domain from the kinase domain of the partner JAK2 and pairing of the two kinase domains, facilitating trans-activation. This model may well generalize to other class I cytokine receptors.

Nanosecond molecular dynamics simulations of bovine pancreatic trypsin inhibitor and lysozyme in water are analyzed in terms of backbone atomic Physikalische Chemie positional fluctuations and dynamical cross-correlations. It is found that although the molecular systems are stable, B-factors calculated over a time Zü rich, Switzerland period as long as 500 ps are not representative for the motions within the proteins. This is especially true for the most mobile residues. On a nanosecond time-scale, the B-factors calculated from the simulations of the proteins in solution are considerably larger than those obtained by structure refinement of the proteins in crystals, based on X-ray data. The time evolution of the atomic fluctuations shows that for large portions of the proteins under study, atomic positional fluctuations are not yet converged after a nanosecond. Cross-correlations do not converge faster than the fluctuations themselves. Most display very erratic behavior if the sampling covers less than about 200 ps. It is also shown that inclusion of mobile atoms into the procedure used to remove rigid-body motion from the simulation can lead to spurious correlations between the motions of the atoms at the surface of the protein.

Molecular dynamics simulations of fully hydrated pure bilayers of four widely studied phospholipids, 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) using a recent revision of the GROMOS96 force field are reported. It is shown that the force field reproduces the structure and the hydration of bilayers formed by each of the four lipids with high accuracy. Specifically, the solvation and the orientation of the dipole of the phosphocholine headgroup and of the ester carbonyls show that the structure of the primary hydration shell in the simulations closely matches experimental findings. This work highlights the need to reproduce a broad range of properties beyond the area per lipid, which is poorly defined experimentally, and to consider the effect of system size and sampling times well beyond those commonly used.