Dataset Information

LIPID11: a modular framework for lipid simulations using amber.

ABSTRACT: Accurate simulation of complex lipid bilayers has long been a goal in condensed phase molecular dynamics (MD). Structure and function of membrane-bound proteins are highly dependent on the lipid bilayer environment and are challenging to study through experimental methods. Within Amber, there has been limited focus on lipid simulations, although some success has been seen with the use of the General Amber Force Field (GAFF). However, to date there are no dedicated Amber lipid force fields. In this paper we describe a new charge derivation strategy for lipids consistent with the Amber RESP approach and a new atom and residue naming and type convention. In the first instance, we have combined this approach with GAFF parameters. The result is LIPID11, a flexible, modular framework for the simulation of lipids that is fully compatible with the existing Amber force fields. The charge derivation procedure, capping strategy, and nomenclature for LIPID11, along with preliminary simulation results and a discussion of the planned long-term parameter development are presented here. Our findings suggest that LIPID11 is a modular framework feasible for phospholipids and a flexible starting point for the development of a comprehensive, Amber-compatible lipid force field.

SUBMITTER: Skjevik AA

PROVIDER: S-EPMC3771533 | biostudies-literature | 2012 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

LIPID11: a modular framework for lipid simulations using amber.

Skjevik Åge A ÅA Madej Benjamin D BD Walker Ross C RC Teigen Knut K

The journal of physical chemistry. B 20120904 36

Accurate simulation of complex lipid bilayers has long been a goal in condensed phase molecular dynamics (MD). Structure and function of membrane-bound proteins are highly dependent on the lipid bilayer environment and are challenging to study through experimental methods. Within Amber, there has been limited focus on lipid simulations, although some success has been seen with the use of the General Amber Force Field (GAFF). However, to date there are no dedicated Amber lipid force fields. In th ...[more]

PMID: 22916730

Similar Datasets

Project description:Computational cardiac models have been extensively used to study different cardiac biomechanics; specifically, finite-element analysis has been one of the tools used to study the internal stresses and strains in the cardiac wall during the cardiac cycle. Cubic-Hermite finite element meshes have been used for simulating cardiac biomechanics due to their convergence characteristics and their ability to capture smooth geometries compactly-fewer elements are needed to build the cardiac geometry-compared to linear tetrahedral meshes. Such meshes have previously been used only with simple ventricular geometries with non-physiological boundary conditions due to challenges associated with creating cubic-Hermite meshes of the complex heart geometry. However, it is critical to accurately capture the different geometric characteristics of the heart and apply physiologically equivalent boundary conditions to replicate the in vivo heart motion. In this work, we created a four-chamber cardiac model utilizing cubic-Hermite elements and simulated a full cardiac cycle by coupling the 3D finite element model with a lumped circulation model. The myocardial fiber-orientations were interpolated within the mesh using the Log-Euclidean method to overcome the singularity associated with interpolation of orthogonal matrices. Physiologically equivalent rigid body constraints were applied to the nodes along the valve plane and the accuracy of the resulting simulations were validated using open source clinical data. We then simulated a complete cardiac cycle of a healthy heart and a heart with acute myocardial infarction. We compared the pumping functionality of the heart for both cases by calculating the ventricular work. We observed a 20% reduction in acute work done by the heart immediately after myocardial infarction. The myocardial wall displacements obtained from the four-chamber model are comparable to actual patient data, without requiring complicated non-physiological boundary conditions usually required in truncated ventricular heart models.

Project description:We designed a strategy for microarray analysis that is based on the identification of transcriptional modules formed by genes coordinately expressed in multiple disease data sets. Mapping changes in gene expression at the module level generated disease-specific transcriptional fingerprints that provide a stable framework for the visualization and functional interpretation of microarray data. The first step of the module-construction process analyzes expression patterns of transcripts across samples for individual diseases: sets of coordinately expressed transcripts were identified with an unsupervised clustering algorithm; in this case, the GeneSpring Version 7.1 (Agilent) implementation of the K-Means algorithm (k = 30). All transcripts detected in at least one sample were used as input; no screening for differential expression was performed. The second step of the module-construction process analyzed the “clustering behavior” of transcripts across diseases, taking into account the possibility that genes may cocluster in some diseases but not in others. Also, in our example, the transcripts that clustered together across all eight diseases were grouped to form a set of modules (round 1 of selection), and the stringency of the analysis was then decreased gradually to identify transcripts that belong to a similar K-means cluster in only a subset of diseases (round 2: seven out of eight diseases; round 3: six out of eight diseases). It is important to note that the module-selection process is “data-driven” and does not involve manual selection of genes by the investigator. We implemented the module-construction strategy described above, using as input a total of 239 peripheral-blood mononuclear cell (PBMC) samples obtained from individuals with one of the following conditions: systemic juvenile idiopathic arthritis (n = 47), systemic lupus erythematosus (n = 40), type I diabetes (n = 20), metastatic melanoma (n = 39), acute infections (Escherichia coli [n = 22], Staphylococcus aureus [n = 18], Influenza A [n = 16]) or liver-transplant recipients undergoing immunosuppressive therapy (n = 37). Transcriptional profiles were generated with Affymetrix U133A and U133B GeneChips (> 44,000 probe sets). A total of 4742 transcripts, distributed among 28 sets, were selected after running of the module-construction algorithm described above. Each module is assigned a unique identifier indicating the round and order of selection (i.e., M3.1 is the first module identified in the third round of selection). The stringency of this algorithm was tested statistically by implementation of the same module-construction procedure after randomization of the original data set. This process was repeated 200 times, without a single module identified. Therefore, the analysis of gene-cluster membership across multiple diseases provided a stringent means to identify PBMC transcriptional modules.

Dataset Information

LIPID11: a modular framework for lipid simulations using amber.

Publications

LIPID11: a modular framework for lipid simulations using amber.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets