Coarse-grained ͑CG͒ models provide a computationally efficient method for rapidly investigating the long time-and length-scale processes that play a critical role in many important biological and soft matter processes. Recently, Izvekov and Voth introduced a new multiscale coarse-graining ͑MS-CG͒ method ͓J. Phys. Chem. B 109, 2469 ͑2005͒; J. Chem. Phys. 123, 134105 ͑2005͔͒ for determining the effective interactions between CG sites using information from simulations of atomically detailed models. The present work develops a formal statistical mechanical framework for the MS-CG method and demonstrates that the variational principle underlying the method may, in principle, be employed to determine the many-body potential of mean force ͑PMF͒ that governs the equilibrium distribution of positions of the CG sites for the MS-CG models. A CG model that employs such a PMF as a "potential energy function" will generate an equilibrium probability distribution of CG sites that is consistent with the atomically detailed model from which the PMF is derived. Consequently, the MS-CG method provides a formal multiscale bridge rigorously connecting the equilibrium ensembles generated with atomistic and CG models. The variational principle also suggests a class of practical algorithms for calculating approximations to this many-body PMF that are optimal. These algorithms use computer simulation data from the atomically detailed model. Finally, important generalizations of the MS-CG method are introduced for treating systems with rigid intramolecular constraints and for developing CG models whose equilibrium momentum distribution is consistent with that of an atomically detailed model.

Coarse-grained (CG) models of biomolecules have recently attracted considerable interest because they enable the simulation of complex biological systems on length-scales and timescales that are inaccessible for atomistic molecular dynamics simulation. A CG model is defined by a map that transforms an atomically detailed configuration into a CG configuration. For CG models of relatively small biomolecules or in cases that the CG and all-atom models have similar resolution, the construction of this map is relatively straightforward and can be guided by chemical intuition. However, it is more challenging to construct a CG map when large and complex domains of biomolecules have to be represented by relatively few CG sites. This work introduces a new and systematic methodology called essential dynamics coarse-graining (ED-CG). This approach constructs a CG map of the primary sequence at a chosen resolution for an arbitrarily complex biomolecule. In particular, the resulting ED-CG method variationally determines the CG sites that reflect the essential dynamics characterized by principal component analysis of an atomistic molecular dynamics trajectory. Numerical calculations illustrate this approach for the HIV-1 CA protein dimer and ATP-bound G-actin. Importantly, since the CG sites are constructed from the primary sequence of the biomolecule, the resulting ED-CG model may be better suited to appropriately explore protein conformational space than those from other CG methods at the same degree of resolution.

No abstract

Epigenetics contributes to the pathogenesis of immune-mediated diseases like rheumatoid arthritis (RA). Here we show the first comprehensive epigenomic characterization of RA fibroblast-like synoviocytes (FLS), including histone modifications (H3K27ac, H3K4me1, H3K4me3, H3K36me3, H3K27me3, and H3K9me3), open chromatin, RNA expression and whole-genome DNA methylation. To address complex multidimensional relationship and reveal epigenetic regulation of RA, we perform integrative analyses using a novel unbiased method to identify genomic regions with similar profiles. Epigenomically similar regions exist in RA cells and are associated with active enhancers and promoters and specific transcription factor binding motifs. Differentially marked genes are enriched for immunological and unexpected pathways, with “Huntington’s Disease Signaling” identified as particularly prominent. We validate the relevance of this pathway to RA by showing that Huntingtin-interacting protein-1 regulates FLS invasion into matrix. This work establishes a high-resolution epigenomic landscape of RA and demonstrates the potential for integrative analyses to identify unanticipated therapeutic targets.