Dataset Information

Optimization of Slipids Force Field Parameters Describing Headgroups of Phospholipids.

ABSTRACT: The molecular mechanics force field Slipids developed in a series of works by Jämbeck and Lyubartsev (J. Phys. Chem. B 2012, 116, 3164-3179; J. Chem. Theory Comput. 2012, 8, 2938-2948) generally provides a good description of various lipid bilayer systems. However, it was also found that order parameters of C-H bonds in the glycerol moiety of the phosphatidylcholine headgroup deviate significantly from NMR results. In this work, the dihedral force field parameters have been reparameterized in order to improve the agreement with experiment. For this purpose, we have computed energies for a large amount of lipid headgroup conformations using density functional theory on the B3P86/cc-pvqz level and optimized dihedral angle parameters simultaneously to provide the best fit to the quantum chemical energies. The new parameter set was validated for three lipid bilayer systems against a number of experimental properties including order parameters, area per lipid, scattering form factors, bilayer thickness, area compressibility and lateral diffusion coefficients. In addition, the order parameter dependence on cholesterol content in the POPC bilayer was investigated. It is shown that the new force field significantly improves agreement with the experimental order parameters for the lipid headgroup while keeping good agreement with other experimentally measured properties.

SUBMITTER: Grote F

PROVIDER: S-EPMC7584372 | biostudies-literature | 2020 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Optimization of Slipids Force Field Parameters Describing Headgroups of Phospholipids.

Grote Fredrik F Lyubartsev Alexander P AP

The journal of physical chemistry. B 20200925 40

The molecular mechanics force field Slipids developed in a series of works by Jämbeck and Lyubartsev (J. Phys. Chem. B 2012, 116, 3164-3179; J. Chem. Theory Comput. 2012, 8, 2938-2948) generally provides a good description of various lipid bilayer systems. However, it was also found that order parameters of C-H bonds in the glycerol moiety of the phosphatidylcholine headgroup deviate significantly from NMR results. In this work, the dihedral force field ...[more]

PMID: 32926619

Similar Datasets

Project description:Molecular dynamics (MD) simulations based on atomic models play an important role in the drug-discovery process to screen molecules, estimate binding free energies, and optimize lead compounds in chemical space. Accurate computations of thermodynamic and kinetic properties using MD simulations are highly dependent on the accuracy of the underlying atomic force field. In this context, going beyond the nonpolarizable fixed-charge model by accounting explicitly for induced polarization is highly desirable. The CHARMM polarizable force field based on classical Drude oscillators, in which an auxiliary charged particle is attached via a harmonic spring to its parent nucleus, offers both a computationally convenient and rigorous framework to model explicitly induced electronic polarization in MD simulations. For any molecule of interest, electrostatic partial charges, atomic polarizabilities, and Thole shielding factors, as well as bonded parameters can either be determined from ab initio calculations or ascribed from the knowledge-based library of the CHARMM Generalized force field. While this approach is fairly reliable in general, it is well understood that the overall accuracy of the models with respect to thermodynamic properties such as bulk density, enthalpies, and solvation free energies is particularly sensitive to the nonbonded Lennard-Jones (LJ) parameters. In the present study, we systematically refined the set of LJ parameters for the atom types available in the Drude force field to best match the experimental thermodynamic properties for 416 small druglike organic molecules. To further test the transferability of the optimized parameters, the hydration free energy of 372 molecules was computed. The calculations resulted in a small average error of 0.46 kcal/mol and a Pearson R of 0.9, representing a significant improvement over the additive GAFF force field in our previous study, where an average error of ∼2 kcal/mol was obtained. Such an improvement is consistent with the ability of the polarizable Drude model to more accurately model interactions in different environments. The effort provides a roadmap for the global optimization of force field parameters using experimental data. It is hoped that the present effort will further the application of the Drude polarizable force field in molecular simulations including drug design and discovery.

Project description:Accurate representations of van der Waals dispersion-repulsion interactions play an important role in high-quality molecular dynamics simulations. Training the force field parameters used in the Lennard Jones (LJ) potential typically used to represent these interactions is challenging, generally requiring adjustment based on simulations of macroscopic physical properties. The large computational expense of these simulations, especially when many parameters must be trained simultaneously, limits the size of training data set and number of optimization steps that can be taken, often requiring modelers to perform optimizations within a local parameter region. To allow for more global LJ parameter optimization against large training sets, we introduce a multi-fidelity optimization technique which uses Gaussian process surrogate modeling to build inexpensive models of physical properties as a function of LJ parameters. This approach allows for fast evaluation of approximate objective functions, greatly accelerating searches over parameter space and enabling the use of optimization algorithms capable of searching more globally. In this study, we use an iterative framework which performs global optimization with differential evolution at the surrogate level, followed by validation at the simulation level and surrogate refinement. Using this technique on two previously studied training sets, containing up to 195 physical property targets, we refit a subset of the LJ parameters for the OpenFF 1.0.0 (Parsley) force field. We demonstrate that this multi-fidelity technique can find improved parameter sets compared to a purely simulation-based optimization by searching more broadly and escaping local minima. Additionally, this technique often finds significantly different parameter minima that have comparably accurate performance. In most cases, these parameter sets are transferable to other similar molecules in a test set. Our multi-fidelity technique provides a platform for rapid, more global optimization of molecular models against physical properties, as well as a number of opportunities for further refinement of the technique.

Project description:Optimization of exposure parameters (target, filter, and kVp) in digital mammography necessitates maximization of the image signal-to-noise ratio (SNR), while simultaneously minimizing patient dose. The goal of this study is to compare, for each of the major commercially available full field digital mammography (FFDM) systems, the impact of the selection of technique factors on image SNR and radiation dose for a range of breast thickness and tissue types. This phantom study is an update of a previous investigation and includes measurements on recent versions of two of the FFDM systems discussed in that article, as well as on three FFDM systems not available at that time. The five commercial FFDM systems tested, the Senographe 2000D from GE Healthcare, the Mammomat Novation DR from Siemens, the Selenia from Hologic, the Fischer Senoscan, and Fuji's 5000MA used with a Lorad M-IV mammography unit, are located at five different university test sites. Performance was assessed using all available x-ray target and filter combinations and nine different phantom types (three compressed thicknesses and three tissue composition types). Each phantom type was also imaged using the automatic exposure control (AEC) of each system to identify the exposure parameters used under automated image acquisition. The figure of merit (FOM) used to compare technique factors is the ratio of the square of the image SNR to the mean glandular dose. The results show that, for a given target/filter combination, in general FOM is a slowly changing function of kVp, with stronger dependence on the choice of target/filter combination. In all cases the FOM was a decreasing function of kVp at the top of the available range of kVp settings, indicating that higher tube voltages would produce no further performance improvement. For a given phantom type, the exposure parameter set resulting in the highest FOM value was system specific, depending on both the set of available target/filter combinations, and on the receptor type. In most cases, the AECs of the FFDM systems successfully identified exposure parameters resulting in FOM values near the maximum ones, however, there were several examples where AEC performance could be improved.

Project description:Computational prediction of noncovalent binding free energies with methods based on molecular mechanical force fields has become increasingly routine in drug discovery projects, where they promise to speed the discovery of small molecule ligands to bind targeted proteins with high affinity. Because the reliability of free energy methods still has significant room for improvement, new force fields, or modifications of existing ones, are regularly introduced with the aim of improving the accuracy of molecular simulations. However, comparatively little work has been done to systematically assess how well force fields perform, particularly in relation to the calculation of binding affinities. Hardware advances have made these calculations feasible, but comprehensive force field assessments for protein-ligand sized systems still remain costly. Here, we turn to cyclodextrin host-guest systems, which feature many hallmarks of protein-ligand binding interactions but are generally much more tractable due to their small size. We present absolute binding free energy and enthalpy calculations, using the attach-pull-release (APR) approach, on a set of 43 cyclodextrin-guest pairs for which experimental ITC data are available. The test set comprises both α- and β-cyclodextrin hosts binding a series of small organic guests, each with one of three functional groups: ammonium, alcohol, or carboxylate. Four water models are considered (TIP3P, TIP4Pew, SPC/E, and OPC), along with two partial charge assignment procedures (RESP and AM1-BCC) and two cyclodextrin host force fields. The results suggest a complex set of considerations when choosing a force field for biomolecular simulations. For example, some force field combinations clearly outperform others at the binding enthalpy calculations but not for the binding free energy. Additionally, a force field combination which we expected to be the worst performer gave the most accurate binding free energies - but the least accurate binding enthalpies. The results have implications for the development of improved force fields, and we propose this test set, and potential future elaborations of it, as a powerful validation suite to evaluate new force fields and help guide future force field development.

Dataset Information

Optimization of Slipids Force Field Parameters Describing Headgroups of Phospholipids.

Publications

Optimization of Slipids Force Field Parameters Describing Headgroups of Phospholipids.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets