Study of the electrostatics treatment in molecular dynamics simulations.

This article considers the treatment of long-range interactions in molecular dynamics simulations. We investigate the effects of using different cutoff distances, constant versus distance-dependent dielectric, and different smoothing methods. In contrast to findings of earlier studies, we find that increasing the cutoff over 8 A does not significantly improve the accuracy (Arnold and Ornstein, Proteins 1994;18:19-33), and using a distance-dependent dielectric instead of a constant dielectric also does not improve accuracy (Guenot and Kollman, Protein Sci 1992;1:1185-1205). This might depend on differences in simulation protocols or force fields, or both, because we use the CHARMM22 force field with stochastic boundary conditions, whereas earlier studies used other protocols and energy functions. We also note that the stability of the simulations is highly dependent on the starting structure, showing that accurate molecular simulations not only depend on a realistic simulation protocol but also on correct initial conditions.

Architecture of beta-barrel membrane proteins: analysis of trimeric porins.

We have analyzed the known three-dimensional structures of trimeric porins from bacterial outer membranes. The distribution of surface-exposed residues in a direction perpendicular to the membrane is similar to that in helical membrane proteins, with aliphatic residues concentrated in the central 20 A of the bilayer. Outside these residues is a layer of aromatic residues, followed by polar and charged residues. Residues in the trimer interface are more conserved than residues not in the interface. By comparing the interface and noninterface residues, an interface preference scale has been derived that may be used as a basis for predicting interface surfaces in monomer models.