Icosahedral packing of polymer-tethered nanospheres and stabilization of the gyroid phase

We present results of molecular simulations that predict the phases formed by the self-assembly of model nanospheres functionalized with a single polymer "tether", including double gyroid, perforated lamella and crystalline bilayer phases. We show that microphase separation of the immiscible tethers and nanospheres causes confinement of the nanoparticles, which promotes local icosahedral packing that stabilizes the gyroid and perforated lamella phases. We present a new metric for determining the local arrangement of particles based on spherical harmonic "fingerprints", which we use to quantify the extent of icosahedral ordering.Comment: 8 pages, 4 figure

CHARMM: The biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological 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 estimators, 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 numerous 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. © 2009 Wiley Periodicals, Inc.J Comput Chem, 2009.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/63074/1/21287_ftp.pd

Generalized solvent boundary potential for computer simulations

Copyright 2001 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in The Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1336570.A general approach has been developed to allow accurate simulations of a small region part of a large macromolecular system while incorporating the influence of the remaining distant atoms with an effective boundary potential. The method is called the Generalized Solvent Boundary Potential (GSBP). By representing the surrounding solvent as a continuum dielectric, both the solvent-shielded static field from the distant atoms of the macromolecule and the reaction field from the dielectricsolvent acting on the atoms in the region of interest are included. The static field is calculated once, using the finite-difference Poisson–Boltzmann (PB) equation, and the result is stored on a discrete grid for efficient simulations. The solventreaction field is developed using a basis-set expansion whose coefficients correspond to generalized electrostatic multipoles. A matrix representing the reaction field Green’s function between those generalized multipoles is calculated only once using the PB equation and stored for efficient simulations. In the present work, the formalism is applied to both spherical and orthorhombic simulation regions for which orthonormal basis-sets exist based on spherical harmonics or cartesian Legendre polynomials. The GSBP method is also tested and illustrated with simple model systems and two detailed atomic systems: the active site region of aspartyl-tRNA synthetase (spherical region) and the interior of the KcsA potassium channel (orthorhombic region). Comparison with numerical finite-difference PB calculations shows that GSBP can accurately describe all long-range electrostatic interactions and remain computationally inexpensive

Evaluating the Effects of Cutoffs and Treatment of Long-range Electrostatics in Protein Folding Simulations

The use of molecular dynamics simulations to provide atomic-level descriptions of biological processes tends to be computationally demanding, and a number of approximations are thus commonly employed to improve computational efficiency. In the past, the effect of these approximations on macromolecular structure and stability has been evaluated mostly through quantitative studies of small-molecule systems or qualitative observations of short-timescale simulations of biological macromolecules. Here we present a quantitative evaluation of two commonly employed approximations, using a test system that has been the subject of a number of previous protein folding studies–the villin headpiece. In particular, we examined the effect of (i) the use of a cutoff-based force-shifting technique rather than an Ewald summation for the treatment of electrostatic interactions, and (ii) the length of the cutoff used to determine how many pairwise interactions are included in the calculation of both electrostatic and van der Waals forces. Our results show that the free energy of folding is relatively insensitive to the choice of cutoff beyond 9 Å, and to whether an Ewald method is used to account for long-range electrostatic interactions. In contrast, we find that the structural properties of the unfolded state depend more strongly on the two approximations examined here

Protein structural variation in computational models and crystallographic data

Normal mode analysis offers an efficient way of modeling the conformational flexibility of protein structures. Simple models defined by contact topology, known as elastic network models, have been used to model a variety of systems, but the validation is typically limited to individual modes for a single protein. We use anisotropic displacement parameters from crystallography to test the quality of prediction of both the magnitude and directionality of conformational variance. Normal modes from four simple elastic network model potentials and from the CHARMM forcefield are calculated for a data set of 83 diverse, ultrahigh resolution crystal structures. While all five potentials provide good predictions of the magnitude of flexibility, the methods that consider all atoms have a clear edge at prediction of directionality, and the CHARMM potential produces the best agreement. The low-frequency modes from different potentials are similar, but those computed from the CHARMM potential show the greatest difference from the elastic network models. This was illustrated by computing the dynamic correlation matrices from different potentials for a PDZ domain structure. Comparison of normal mode results with anisotropic temperature factors opens the possibility of using ultrahigh resolution crystallographic data as a quantitative measure of molecular flexibility. The comprehensive evaluation demonstrates the costs and benefits of using normal mode potentials of varying complexity. Comparison of the dynamic correlation matrices suggests that a combination of topological and chemical potentials may help identify residues in which chemical forces make large contributions to intramolecular coupling.Comment: 17 pages, 4 figure