抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白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.

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system's overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules must keep safe pairwise distances. Obtaining such a molecular arrangement can be considered a packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules must be greater than some specified tolerance. We have developed a code able to pack millions of atoms, grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straightforward. In addition, different atoms belonging to the same molecule may also be restricted to different spatial regions, in such a way that more ordered molecular arrangements can be built, as micelles, lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies from a few seconds to a few minutes in a personal computer. The input files are simple and currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and can be downloaded from http://www.ime.unicamp.br/~martinez/packmol/.

/pdf/packmol-a-package-for-building-initial-configurations-for-qsy9r74wrw.pdf

PACKMOL: a package for building initial configurations for molecular dynamics simulations.

CHARMM is an academic research program used widely for macromolecular mechanics and dynamics with versatile analysis and manipulation tools of atomic coordinates and dynamics trajectories. CHARMM-GUI, http://www.charmm-gui.org, has been developed to provide a web-based graphical user interface to generate various input files and molecular systems to facilitate and standardize the usage of common and advanced simulation techniques in CHARMM. The web environment provides an ideal platform to build and validate a molecular model system in an interactive fashion such that, if a problem is found through visual inspection, one can go back to the previous setup and regenerate the whole system again. In this article, we describe the currently available functional modules of CHARMM-GUI Input Generator that form a basis for the advanced simulation techniques. Future directions of the CHARMM-GUI development project are also discussed briefly together with other features in the CHARMM-GUI website, such as Archive and Movie Gallery.

https://onlinelibrary.wiley.com/doi/pdf/10.1002/jcc.20945

CHARMM-GUI: a web-based graphical user interface for CHARMM.

While the quality of the current CHARMM22/CMAP additive force field for proteins has been demonstrated in a large number of applications, limitations in the model with respect to the equilibrium between the sampling of helical and extended conformations in folding simulations have been noted. To overcome this, as well as make other improvements in the model, we present a combination of refinements that should result in enhanced accuracy in simulations of proteins. The common (non Gly, Pro) backbone CMAP potential has been refined against experimental solution NMR data for weakly structured peptides, resulting in a rebalancing of the energies of the α-helix and extended regions of the Ramachandran map, correcting the α-helical bias of CHARMM22/CMAP. The Gly and Pro CMAPs have been refitted to more accurate quantum-mechanical energy surfaces. Side-chain torsion parameters have been optimized by fitting to backbone-dependent quantum-mechanical energy surfaces, followed by additional empirical optimization targeting NMR scalar couplings for unfolded proteins. A comprehensive validation of the revised force field was then performed against data not used to guide parametrization: (i) comparison of simulations of eight proteins in their crystal environments with crystal structures; (ii) comparison with backbone scalar couplings for weakly structured peptides; (iii) comparison with NMR residual dipolar couplings and scalar couplings for both backbone and side-chains in folded proteins; (iv) equilibrium folding of mini-proteins. The results indicate that the revised CHARMM 36 parameters represent an improved model for the modeling and simulation studies of proteins, including studies of protein folding, assembly and functionally relevant conformational changes.

/pdf/optimization-of-the-additive-charmm-all-atom-protein-force-4qtc6t5kz1.pdf

Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral angles.

CHARMM is an academic research program used widely for macromolecular mechanics and dynamicswith versatile analysis and manipulation tools of atomic coordinates and dynamics trajectories. CHARMM-GUI,http://www.charmm-gui.org, has been developed to provide a web-based graphical user interface to generate variousinput ﬁles and molecular systems to facilitate and standardize the usage of common and advanced simulation techni-ques in CHARMM. The web environment provides an ideal platform to build and validate a molecular model sys-tem in an interactive fashion such that, if a problem is found through visual inspection, one can go back to the pre-vious setup and regenerate the whole system again. In this article, we describe the currently available functionalmodules of CHARMM-GUI Input Generator that form a basis for the advanced simulation techniques. Future direc-tions of the CHARMM-GUI development project are also discussed brieﬂy together with other features in theCHARMM-GUI website, such as Archive and Movie Gallery.q 2008 Wiley Periodicals, Inc. J Comput Chem 29: 1859–1865, 2008Key words: molecular dynamics; Poisson-Boltzmann; electrostatic potential; membrane; visualization; explicit sol-vent; implicit solvent; MarvinSpace

Software News and Updates CHARMM-GUI: A Web-Based Graphical User Interface for CHARMM

Molecular dynamics simulations of membrane proteins have provided deeper insights into their functions and interactions with surrounding environments at the atomic level. However, compared to solvation of globular proteins, building a realistic protein/membrane complex is still challenging and requires considerable experience with simulation software. Membrane Builder in the CHARMM-GUI website (http://www.charmm-gui.org) helps users to build such a complex system using a web browser with a graphical user interface. Through a generalized and automated building process including system size determination as well as generation of lipid bilayer, pore water, bulk water, and ions, a realistic membrane system with virtually any kinds and shapes of membrane proteins can be generated in 5 minutes to 2 hours depending on the system size. Default values that were elaborated and tested extensively are given in each step to provide reasonable options and starting points for both non-expert and expert users. The efficacy of Membrane Builder is illustrated by its applications to 12 transmembrane and 3 interfacial membrane proteins, whose fully equilibrated systems with three different types of lipid molecules (DMPC, DPPC, and POPC) and two types of system shapes (rectangular and hexagonal) are freely available on the CHARMM-GUI website. One of the most significant advantages of using the web environment is that, if a problem is found, users can go back and re-generate the whole system again before quitting the browser. Therefore, Membrane Builder provides the intuitive and easy way to build and simulate the biologically important membrane system.

/pdf/automated-builder-and-database-of-protein-membrane-complexes-5h7kcwiajj.pdf

Automated builder and database of protein/membrane complexes for molecular dynamics simulations.

CHARMM36 (C36) is the most up-to-date pairwise additive all-atom lipid force field and is able to accurately represent bilayer properties of saturated and monounsaturated lipid molecules in the natural constant particle, pressure, and temperature (NPT) ensemble. However, molecular dynamics (MD) simulations on 1-stearoyl-2-docosahexaenoyl-sn-glycerco-3-phosphocholine (SDPC) bilayers of the polyunsaturated fatty acid (PUFA) chains result in inaccuracies of the surface area per lipid (SA), deuterium order parameters (SCD), and X-ray form factors. Therefore, in this study, high-level quantum mechanical calculations are used to improve the dihedral potential of neighboring double bonds, and the corresponding force field is referred to as C36p. The SA for SDPC at 303 K increases from 63.2 ± 0.2 (C36) to 70.8 ± 0.2 (C36p) A2 and agrees favorably with X-ray diffraction results at 297 K. The resulting SCD are in excellent agreement with experimental values of both the sn-1 and sn-2 chains. Calculated NMR 13C relax...

Improving the CHARMM force field for polyunsaturated fatty acid chains.

https://kuscholarworks.ku.edu/bitstream/handle/1808/17465/Wonpil-v99_p175-183.pdf;sequence=1

Taehoon Kim

Papers

CHARMM-GUI: a web-based graphical user interface for CHARMM.

Software News and Updates CHARMM-GUI: A Web-Based Graphical User Interface for CHARMM

Automated builder and database of protein/membrane complexes for molecular dynamics simulations.

Improving the CHARMM force field for polyunsaturated fatty acid chains.

Revisiting Hydrophobic Mismatch with Free Energy Simulation Studies of Transmembrane Helix Tilt and Rotation