Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecu- lar simulation program. It has been developed over the last three decades with a primary focus on molecules of bio- logical interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estima- tors, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numer- ous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

/pdf/charmm-the-biomolecular-simulation-program-4kexem0xkx.pdf

CHARMM: the biomolecular simulation program.

저작근 근전도에 관한 정상교합자와 ii급 부정교합자의 비교 연구

Voltage-gated sodium (NaV) channels initiate electrical signalling in excitable cells and are the molecular targets for drugs and disease mutations, but the structural basis for their voltage-dependent activation, ion selectivity and drug block is unknown. Here we report the crystal structure of a voltage-gated Na+ channel from Arcobacter butzleri (NavAb) captured in a closed-pore conformation with four activated voltage sensors at 2.7 A resolution. The arginine gating charges make multiple hydrophilic interactions within the voltage sensor, including unanticipated hydrogen bonds to the protein backbone. Comparisons to previous open-pore potassium channel structures indicate that the voltage-sensor domains and the S4–S5 linkers dilate the central pore by pivoting together around a hinge at the base of the pore module. The NavAb selectivity filter is short, ∼4.6 A wide, and water filled, with four acidic side chains surrounding the narrowest part of the ion conduction pathway. This unique structure presents a high-field-strength anionic coordination site, which confers Na+ selectivity through partial dehydration via direct interaction with glutamate side chains. Fenestrations in the sides of the pore module are unexpectedly penetrated by fatty acyl chains that extend into the central cavity, and these portals are large enough for the entry of small, hydrophobic pore-blocking drugs. This structure provides the template for understanding electrical signalling in excitable cells and the actions of drugs used for pain, epilepsy and cardiac arrhythmia at the atomic level. The X-ray crystal structure of a voltage-gated sodium channel from Arcobacter butzleri has been determined, with the channel in the closed-pore conformation. Channels of this type initiate electrical signalling in excitable cells and are the molecular targets for many drugs, but the structural basis for their voltage-dependent activation and ion selectivity is not known. The selectivity filter in this sodium channel is found to be quite short, compared with those in open-pore potassium channels, and the voltage-sensor domains and linkers between segments S4 and S5 seem to dilate the central pore by pivoting together.

The crystal structure of a voltage-gated sodium channel

The goal of this review is to establish a broad and rigorous theoretical framework to describe ion permeation through biological channels. This framework is developed in the context of atomic models on the basis of the statistical mechanical projection-operator formalism of Mori and Zwanzig. The review is divided into two main parts. The first part introduces the fundamental concepts needed to construct a hierarchy of dynamical models at different level of approximation. In particular, the potential of mean force (PMF) as a configuration-dependent free energy is introduced, and its significance concerning equilibrium and non-equilibrium phenomena is discussed. In addition, fundamental aspects of membrane electrostatics, with a particular emphasis on the influence of the transmembrane potential, as well as important computational techniques for extracting essential information from all-atom molecular dynamics (MD) simulations are described and discussed. The first part of the review provides a theoretical formalism to 'translate' the information from the atomic structure into the familiar language of phenomenological models of ion permeation. The second part is aimed at reviewing and contrasting results obtained in recent computational studies of three very different channels: the gramicidin A (gA) channel, which is a narrow one-ion pore (at moderate concentration), the KcsA channel from Streptomyces lividans, which is a narrow multi-ion pore, and the outer membrane matrix porin F (OmpF) from Escherichia coli, which is a trimer of three beta-barrel subunits each forming wide aqueous multi-ion pores. Comparison with experiments demonstrates that current computational models are approaching semi-quantitative accuracy and are able to provide significant insight into the microscopic mechanisms of ion conduction and selectivity. We conclude that all-atom MD with explicit water molecules can represent important structural features of complex biological channels accurately, including such features as the location of ion-binding sites along the permeation pathway. We finally discuss the broader issue of the validity of ion permeation models and an outlook to the future.

/pdf/theoretical-and-computational-models-of-biological-ion-3gr016lg5p.pdf

Theoretical and computational models of biological ion channels

The free energy governing K+ conduction through gramicidin A channels is characterized by using over 0.1 μs of all-atom molecular dynamics simulations with explicit solvent and membrane. The results provide encouraging agreement with experiments and insights into the permeation mechanism. The free energy surface of K+, as a function of both axial and radial coordinates, is calculated. Correcting for simulation artifacts due to periodicity and the lack of hydrocarbon polarizability, the calculated single-channel conductance for K+ ions is 0.8 pS, closer to experiment than any previous calculation. In addition, the estimated single ion dissociation constants are within the range of experimental determinations. The relatively small free energy barrier to ion translocation arises from a balance of large opposing contributions from protein, single-file water, bulk electrolyte, and membrane. Mean force decomposition reveals a remarkable ability of the single-file water molecules to stabilize K+ by –40 kcal/mol, roughly half the bulk solvation free energy. The importance of the single-file water confirms the conjecture of Mackay et al. [Mackay, D. H. J., Berens, P. H., Wilson, K. R. & Hagler, A. T. (1984) Biophys. J. 46, 229–248]. Ion association with the channel involves gradual dehydration from approximately six to seven water molecules in the first shell, to just two inside the narrow pore. Ion permeation is influenced by the orientation of the single-file water column, which can present a barrier to conduction and give rise to long-range coupling of ions on either side of the pore. Small changes in the potential function, including contributions from electronic polarization, are likely to be sufficient to obtain quantitative agreement with experiments.

https://www.pnas.org/content/pnas/101/1/117.full.pdf

Energetics of ion conduction through the gramicidin channel

Biological membranes consist of bilayer arrangements of lipids forming a hydrophobic core that presents a physical barrier to all polar and charged molecules. This long-held notion has recently been challenged by biological translocon-based experiments that report small apparent free energies to insert charged side chains near the center of a transmembrane (TM) helix. We have carried out fully atomistic simulations to provide the free-energy profile for moving a TM helix containing a protonated Arg side chain across a lipid bilayer. Our results reveal the fundamental thermodynamics governing the stability of charged side chains in membranes and the microscopic interactions involved. Despite local membrane deformations, where large amounts of water and lipid head groups are pulled into the bilayer to interact with Arg, the free-energy barrier is 17 kcal/mol. We provide a rationale for the differences in our microscopic free energies and cell biological experiments using free-energy calculations that indicate that a protonated Arg at the central residue of a TM helix of the Leader peptidase might reside close to the interface and not at the membrane center. Our findings have implications for the gating mechanisms of voltage-gated ion channels, suggesting that movements of protonated Arg residues through the membrane will be prohibited.

/pdf/on-the-thermodynamic-stability-of-a-charged-arginine-side-xxszqtvb11.pdf

On the thermodynamic stability of a charged arginine side chain in a transmembrane helix

The basic amino acids lysine (Lys) and arginine (Arg) play important roles in membrane protein activity, the sensing of membrane voltages, and the actions of antimicrobial, toxin, and cell-penetrating peptides. These roles are thought to stem from the strong interactions and disruptive influences of these amino acids on lipid membranes. In this study, we employ fully atomistic molecular dynamics simulations to observe, quantify, and compare the interactions of Lys and Arg with saturated phosphatidylcholine membranes of different thickness. We make use of both charged (methylammonium and methylguanidinium) and neutral (methylamine and methylguanidine) analogue molecules, as well as Lys and Arg side chains on transmembrane helix models. We find that the free energy barrier experienced by a charged Lys crossing the membrane is strikingly similar to that of a charged Arg (to within 2 kcal/mol), despite the two having different chemistries, H-bonding capability, and hydration free energies that differ by ∼10 k...

The Different Interactions of Lysine and Arginine Side Chains with Lipid Membranes

In a recent article published by Bastug and Kuyucak [J. Chem. Phys.126, 105103 (2007)] investigated the microscopic factors affecting double ion occupancy in the gramicidin channel. The analysis relied largely on the one-dimensional potential of mean force of ions along the axis of the channel (the so-called free energy profile of the ion along the channel axis), as well as on the calculation of the equilibrium association constant of the ions in the channel binding sites. It is the purpose of this communication to clarify this issue.

Toby W. Allen

Papers

Theoretical and computational models of biological ion channels

Energetics of ion conduction through the gramicidin channel

On the thermodynamic stability of a charged arginine side chain in a transmembrane helix

The Different Interactions of Lysine and Arginine Side Chains with Lipid Membranes

Comment on "free energy simulations of single and double ion occupancy in gramicidin A"