Supplementary data are available at Bioinformatics online.

CHARMM27 is a widespread and popular force field for biomolecular simulation, and several recent algorithms such as implicit solvent models have been developed specifically for it. We have here implemented the CHARMM force field and all necessary extended functional forms in the GROMACS molecular simulation package, to make CHARMM-specific features available and to test them in combination with techniques for extended time steps, to make all major force fields available for comparison studies in GROMACS, and to test various solvent model optimizations, in particular the effect of Lennard-Jones interactions on hydrogens. The implementation has full support both for CHARMM-specific features such as multiple potentials over the same dihedral angle and the grid-based energy correction map on the ϕ, ψ protein backbone dihedrals, as well as all GROMACS features such as virtual hydrogen interaction sites that enable 5 fs time steps. The medium-to-long time effects of both the correction maps and virtual sites have been tested by performing a series of 100 ns simulations using different models for water representation, including comparisons between CHARMM and traditional TIP3P. Including the correction maps improves sampling of near native-state conformations in our systems, and to some extent it is even able to refine distorted protein conformations. Finally, we show that this accuracy is largely maintained with a new implicit solvent implementation that works with virtual interaction sites, which enables performance in excess of 250 ns/day for a 900-atom protein on a quad-core desktop computer.

Voltage-gated ion channels open and close in response to changes in membrane potential, thereby enabling electrical signaling in excitable cells. The voltage sensitivity is conferred through four voltage-sensor domains (VSDs) where positively charged residues in the fourth transmembrane segment (S4) sense the potential. While an open state is known from the Kv1.2/2.1 X-ray structure, the conformational changes underlying voltage sensing have not been resolved. We present 20 additional interactions in one open and four different closed conformations based on metal-ion bridges between all four segments of the VSD in the voltage-gated Shaker K channel. A subset of the experimental constraints was used to generate Rosetta models of the conformations that were subjected to molecular simulation and tested against the remaining constraints. This achieves a detailed model of intermediate conformations during VSD gating. The results provide molecular insight into the transition, suggesting that S4 slides at least 12 Å along its axis to open the channel with a 3 10 helix region present that moves in sequence in S4 in order to occupy the same position in space opposite F290 from open through the three first closed states. V oltage-gated ion channels are critical for biological signaling, and they are able to regulate ion flux on a millisecond time scale. To sense changes in membrane voltage, each ion channel is equipped with four voltage-sensor domains (VSDs) connected to a central ion-conducting pore domain. The fourth transmembrane segment (S4) of each VSD carries several positively charged amino-acid residues responsible for VSD gating (1). At least three elementary charges per VSD must traverse outwards through the membrane electric field to open a channel that corresponds to a considerable displacement of the S4 helix ( Fig. 1A) (2-4). The positive charges in S4 make salt bridges with negative countercharges on their move through the VSD (4-8). It has even been proposed that the VSD undergoes a conformational alteration following the opening, when the channel relaxes to an inactivated, that is closed, state (9). In addition to conferring voltage dependence to ion channels, VSDs also regulate enzymes (10), act as voltage-gated proton channels (11,12), are susceptible to disease-causing mutations (13,14), and serve as targets for drugs and toxins (1,(15)(16)(17)(18). Therefore, it is of crucial interest to understand the details underlying voltage sensing by VSDs.Few, if any, segments of membrane proteins have received more attention than the S4 helix of voltage sensors. In addition to their paramount biological importance, they can help us understand fundamental biophysical problems such as why some membrane protein segments can be hydrophilic (19), how charges effectively move through a membrane, or how a potential triggers structural changes on a microsecond time scale. These questions are inherently linked to transient conformations and contacts that can be difficult to capture in a single structure. Although active...

We describe how membrane proteins diffuse laterally in the membrane plane together with the lipids surrounding them. We find a number of intriguing phenomena. The lateral displacements of the protein and the lipids are strongly correlated, as the protein and the neighboring lipids form a dynamical protein-lipid complex, consisting of approximately 50-100 lipids. The diffusion of the lipids in the complex is much slower compared to the rest of the lipids. We also find a strong directional correlation between the movements of the protein and the lipids in its vicinity. The results imply that in crowded membrane environments there are no "free" lipids, as they are all influenced by the protein structure and dynamics. Our results indicate that, in studies of cell membranes, protein and lipid dynamics have to be considered together.

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins—not to mention numerous applications in drug design. Here, we present a full 1 µs atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120° rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation (∼35°) of the extracellular end of all S4 segments is present also in a reference 0.5 µs simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 310 helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4–lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5–1 µs). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations.

The activation of voltage-gated ion channels is controlled by the S4 helix, with arginines every third residue. The x-ray structures are believed to reflect an open-inactivated state, and models propose combinations of translation, rotation, and tilt to reach the resting state. Recently, experiments and simulations have independently observed occurrence of 3(10)-helix in S4. This suggests S4 might make a transition from α- to 3(10)-helix in the gating process. Here, we show 3(10)-helix structure between Q1 and R3 in the S4 segment of a voltage sensor appears to facilitate the early stage of the motion toward a down state. We use multiple microsecond-steered molecular simulations to calculate the work required for translating S4 both as α-helix and transformed to 3(10)-helix. The barrier appears to be caused by salt-bridge reformation simultaneous to R4 passing the F233 hydrophobic lock, and it is almost a factor-two lower with 3(10)-helix. The latter facilitates translation because R2/R3 line up to face E183/E226, which reduces the requirement to rotate S4. This is also reflected in a lower root mean-square deviation distortion of the rest of the voltage sensor. This supports the 3(10) hypothesis, and could explain some of the differences between the open-inactivated- versus activated-states.

SUMMARYFor the purpose of developing a national system for outbreak surveillance, local outbreak signals were compared in three sources of syndromic data – telephone triage of acute gastroenteritis, web queries about symptoms of gastrointestinal illness, and over-the-counter (OTC) pharmacy sales of antidiarrhoeal medication. The data sources were compared against nine known waterborne and foodborne outbreaks in Sweden in 2007–2011. Outbreak signals were identified for the four largest outbreaks in the telephone triage data and the two largest outbreaks in the data on OTC sales of antidiarrhoeal medication. No signals could be identified in the data on web queries. The signal magnitude for the fourth largest outbreak indicated a tenfold larger outbreak than officially reported, supporting the use of telephone triage data for situational awareness. For the two largest outbreaks, telephone triage data on adult diarrhoea provided outbreak signals at an early stage, weeks and months in advance, respectively, potentially serving the purpose of early event detection. In conclusion, telephone triage data provided the most promising source for surveillance of point-source outbreaks.

Molecular-dynamics simulations were carried out to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistent with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive. The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be nonconducting. The conductance of the hexamer was estimated to be 115 ± 34 pS and 74 ± 20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis, we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K⁺ and Cl⁻ with their first solvation shells intact. The free energy barrier encountered by K⁺ is only 2.2 kcal/mol whereas Cl⁻ encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics.