1
|
Oliveira MP, Gonçalves YMH, Ol Gheta SK, Rieder SR, Horta BAC, Hünenberger PH. Comparison of the United- and All-Atom Representations of (Halo)alkanes Based on Two Condensed-Phase Force Fields Optimized against the Same Experimental Data Set. J Chem Theory Comput 2022; 18:6757-6778. [PMID: 36190354 DOI: 10.1021/acs.jctc.2c00524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The level of accuracy that can be achieved by a force field is influenced by choices made in the interaction-function representation and in the relevant simulation parameters. These choices, referred to here as functional-form variants (FFVs), include for example the model resolution, the charge-derivation procedure, the van der Waals combination rules, the cutoff distance, and the treatment of the long-range interactions. Ideally, assessing the effect of a given FFV on the intrinsic accuracy of the force-field representation requires that only the specific FFV is changed and that this change is performed at an optimal level of parametrization, a requirement that may prove extremely challenging to achieve in practice. Here, we present a first attempt at such a comparison for one specific FFV, namely the choice of a united-atom (UA) versus an all-atom (AA) resolution in a force field for saturated acyclic (halo)alkanes. Two force-field versions (UA vs AA) are optimized in an automated way using the CombiFF approach against 961 experimental values for the pure-liquid densities ρliq and vaporization enthalpies ΔHvap of 591 compounds. For the AA force field, the torsional and third-neighbor Lennard-Jones parameters are also refined based on quantum-mechanical rotational-energy profiles. The comparison between the UA and AA resolutions is also extended to properties that have not been included as parameterization targets, namely the surface-tension coefficient γ, the isothermal compressibility κT, the isobaric thermal-expansion coefficient αP, the isobaric heat capacity cP, the static relative dielectric permittivity ϵ, the self-diffusion coefficient D, the shear viscosity η, the hydration free energy ΔGwat, and the free energy of solvation ΔGche in cyclohexane. For the target properties ρliq and ΔHvap, the UA and AA resolutions reach very similar levels of accuracy after optimization. For the nine other properties, the AA representation leads to more accurate results in terms of η; comparably accurate results in terms of γ, κT, αP, ϵ, D, and ΔGche; and less accurate results in terms of cP and ΔGwat. This work also represents a first step toward the calibration of a GROMOS-compatible force field at the AA resolution.
Collapse
Affiliation(s)
- Marina P Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Yan M H Gonçalves
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - S Kashef Ol Gheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A C Horta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Philippe H Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
2
|
Mavra A, Petrou CC, Vlasiou MC. Ligand and Structure-Based Virtual Screening in Combination, to Evaluate Small Organic Molecules as Inhibitors for the XIAP Anti-Apoptotic Protein: The Xanthohumol Hypothesis. Molecules 2022; 27:molecules27154825. [PMID: 35956774 PMCID: PMC9369490 DOI: 10.3390/molecules27154825] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Revised: 07/25/2022] [Accepted: 07/26/2022] [Indexed: 11/30/2022] Open
Abstract
Herein, we propose two chalcone molecules, (E)-1-(4-methoxyphenyl)-3-(p-tolyl) prop-2-en-1-one and (E)-3-(4-hydroxyphenyl)-1-(2,4,6-trihydroxyphenyl) prop-2-en-1-one, based on the anticancer bioactive molecule Xanthohumol, which are suitable for further in vitro and in vivo studies. Their ability to create stable complexes with the antiapoptotic X-linked IAP (XIAP) protein makes them promising anticancer agents. The calculations were based on ligand-based and structure-based virtual screening combined with the pharmacophore build. Additionally, the structures passed Lipinski’s rule for drug use, and their reactivity was confirmed using density functional theory studies. ADMET studies were also performed to reveal the pharmacokinetic potential of the compounds. The candidates were chosen from 10,639,400 compounds, and the docking protocols were evaluated using molecular dynamics simulations.
Collapse
|
3
|
Partial charges for molecular-mechanical models of heterocyclic compounds: pyridine nitrogen. J Mol Model 2021; 27:118. [PMID: 33818696 DOI: 10.1007/s00894-021-04731-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Accepted: 03/14/2021] [Indexed: 10/21/2022]
Abstract
This paper deals with the effect of introducing an additional interaction site onto molecular-mechanical models of nitrogen-containing heterocyclic compounds. The introduction of only one additional site next to nitrogen atoms is shown to result in significant improvement of the quality of the models along with negligible slowdown of calculation speed. Concretely, it was proposed to introduce the site inside the aromatic ring at a distance of 0.4 Å from the nitrogen atom center. All the parametrization can be completely automated. The proposed force field allows predicting heats of evaporation of liquids of the compound under investigation with an accuracy of 1 kcal/mol.
Collapse
|
4
|
Oliveira MP, Andrey M, Rieder SR, Kern L, Hahn DF, Riniker S, Horta BAC, Hünenberger PH. Systematic Optimization of a Fragment-Based Force Field against Experimental Pure-Liquid Properties Considering Large Compound Families: Application to Saturated Haloalkanes. J Chem Theory Comput 2020; 16:7525-7555. [DOI: 10.1021/acs.jctc.0c00683] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Maurice Andrey
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R. Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Leyla Kern
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - David F. Hahn
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
5
|
Finney AR, Innocenti Malini R, Freeman CL, Harding JH. Amino Acid and Oligopeptide Effects on Calcium Carbonate Solutions. CRYSTAL GROWTH & DESIGN 2020; 20:3077-3092. [PMID: 32581657 PMCID: PMC7304842 DOI: 10.1021/acs.cgd.9b01693] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2019] [Revised: 04/08/2020] [Indexed: 05/04/2023]
Abstract
Biological organisms display sophisticated control of nucleation and crystallization of minerals. In order to mimic living systems, deciphering the mechanisms by which organic molecules control the formation of mineral phases from solution is a key step. We have used computer simulations to investigate the effects of the amino acids arginine, aspartic acid, and glycine on species that form in solutions of calcium carbonate (CaCO3) at lower and higher levels of supersaturation. This provides net positive, negative, and neutral additives. In addition, we have prepared simulations containing hexapeptides of the amino acids to consider the effect of additive size on the solution species. We find that additives have limited impact on the formation of extended, liquid-like CaCO3 networks in supersaturated solutions. Additives control the amount of (bi)carbonate in solution, but more importantly, they are able to stabilize these networks on the time scales of the simulations. This is achieved by coordinating the networks and assembled additive clusters in solutions. The association leads to subtle changes in the coordination of CaCO3 and reduced mobility of the cations. We find that the number of solute association sites and the size and topology of the additives are more important than their net charge. Our results help to understand why polymer additives are so effective at stabilizing dense liquid CaCO3 phases.
Collapse
Affiliation(s)
- Aaron R. Finney
- Department
of Materials Science and Engineering, University
of Sheffield, Sheffield S1 3JD, United Kingdom
- Department
of Chemical Engineering, University College
London, London WC1E 6BT, United Kingdom
- E-mail:
| | - Riccardo Innocenti Malini
- Laboratory
for Biomimetic Membranes and Textiles, EMPA,
Swiss Federal Laboratories for Materials Science and Technology, St. Gallen 9014, Switzerland
| | - Colin L. Freeman
- Department
of Materials Science and Engineering, University
of Sheffield, Sheffield S1 3JD, United Kingdom
| | - John H. Harding
- Department
of Materials Science and Engineering, University
of Sheffield, Sheffield S1 3JD, United Kingdom
| |
Collapse
|
6
|
Işık M, Bergazin TD, Fox T, Rizzi A, Chodera JD, Mobley DL. Assessing the accuracy of octanol-water partition coefficient predictions in the SAMPL6 Part II log P Challenge. J Comput Aided Mol Des 2020; 34:335-370. [PMID: 32107702 PMCID: PMC7138020 DOI: 10.1007/s10822-020-00295-0] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 01/24/2020] [Indexed: 12/12/2022]
Abstract
The SAMPL Challenges aim to focus the biomolecular and physical modeling community on issues that limit the accuracy of predictive modeling of protein-ligand binding for rational drug design. In the SAMPL5 log D Challenge, designed to benchmark the accuracy of methods for predicting drug-like small molecule transfer free energies from aqueous to nonpolar phases, participants found it difficult to make accurate predictions due to the complexity of protonation state issues. In the SAMPL6 log P Challenge, we asked participants to make blind predictions of the octanol-water partition coefficients of neutral species of 11 compounds and assessed how well these methods performed absent the complication of protonation state effects. This challenge builds on the SAMPL6 p[Formula: see text] Challenge, which asked participants to predict p[Formula: see text] values of a superset of the compounds considered in this log P challenge. Blind prediction sets of 91 prediction methods were collected from 27 research groups, spanning a variety of quantum mechanics (QM) or molecular mechanics (MM)-based physical methods, knowledge-based empirical methods, and mixed approaches. There was a 50% increase in the number of participating groups and a 20% increase in the number of submissions compared to the SAMPL5 log D Challenge. Overall, the accuracy of octanol-water log P predictions in SAMPL6 Challenge was higher than cyclohexane-water log D predictions in SAMPL5, likely because modeling only the neutral species was necessary for log P and several categories of method benefited from the vast amounts of experimental octanol-water log P data. There were many highly accurate methods: 10 diverse methods achieved RMSE less than 0.5 log P units. These included QM-based methods, empirical methods, and mixed methods with physical modeling supported with empirical corrections. A comparison of physical modeling methods showed that QM-based methods outperformed MM-based methods. The average RMSE of the most accurate five MM-based, QM-based, empirical, and mixed approach methods based on RMSE were 0.92 ± 0.13, 0.48 ± 0.06, 0.47 ± 0.05, and 0.50 ± 0.06, respectively.
Collapse
Affiliation(s)
- Mehtap Işık
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA.
- Tri-Institutional PhD Program in Chemical Biology, Weill Cornell Graduate School of Medical Sciences, Cornell University, New York, NY, 10065, USA.
| | | | - Thomas Fox
- Computational Chemistry, Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co KG, 88397, Biberach, Germany
| | - Andrea Rizzi
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
- Tri-Institutional Training Program in Computational Biology and Medicine, New York, NY, 10065, USA
| | - John D Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
| | - David L Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine, CA, 92697, USA
- Department of Chemistry, University of California, Irvine, CA, 92697, USA
| |
Collapse
|
7
|
Magnitude of attraction in CF4-CH4 interactions: Are CF4-CH4 interactions weaker than average of CF4-CF4 and CH4-CH4 interactions? J Fluor Chem 2020. [DOI: 10.1016/j.jfluchem.2020.109468] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
8
|
Nikitin A. Non-zero Lennard-Jones parameters for the Toukan–Rahman water model: more accurate calculations of the solvation free energy of organic substances. J Comput Aided Mol Des 2019; 34:437-441. [DOI: 10.1007/s10822-019-00256-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2019] [Accepted: 11/19/2019] [Indexed: 11/29/2022]
|
9
|
The role of hydration effects in 5-fluorouridine binding to SOD1: insight from a new 3D-RISM-KH based protocol for including structural water in docking simulations. J Comput Aided Mol Des 2019; 33:913-926. [PMID: 31686367 DOI: 10.1007/s10822-019-00239-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2019] [Accepted: 10/17/2019] [Indexed: 12/13/2022]
Abstract
Misfolded Cu/Zn superoxide dismutase enzyme (SOD1) shows prion-like propagation in neuronal cells leading to neurotoxic aggregates that are implicated in amyotrophic lateral sclerosis (ALS). Tryptophan-32 (W32) in SOD1 is part of a potential site for templated conversion of wild type SOD1. This W32 binding site is located on a convex, solvent exposed surface of the SOD1 suggesting that hydration effects can play an important role in ligand recognition and binding. A recent X-ray crystal structure has revealed that 5-Fluorouridine (5-FUrd) binds at the W32 binding site and can act as a pharmacophore scaffold for the development of anti-ALS drugs. In this study, a new protocol is developed to account for structural (non-displaceable) water molecules in docking simulations and successfully applied to predict the correct docked conformation binding modes of 5-FUrd at the W32 binding site. The docked configuration is within 0.58 Å (RMSD) of the observed configuration. The docking protocol involved calculating a hydration structure around SOD1 using molecular theory of solvation (3D-RISM-KH, 3D-Reference Interaction Site Model-Kovalenko-Hirata) whereby, non-displaceable water molecules are identified for docking simulations. This protocol was also used to analyze the hydrated structure of the W32 binding site and to explain the role of solvation in ligand recognition and binding to SOD1. Structural water molecules mediate hydrogen bonds between 5-FUrd and the receptor, and create an environment favoring optimal placement of 5-FUrd in the W32 binding site.
Collapse
|
10
|
Bhadra P, Siu SWI. Refined Empirical Force Field to Model Protein-Self-Assembled Monolayer Interactions Based on AMBER14 and GAFF. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:9622-9633. [PMID: 31246036 DOI: 10.1021/acs.langmuir.9b01367] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Understanding protein interaction with material surfaces is important for the development of nanotechnological devices. The structures and dynamics of proteins can be studied via molecular dynamics (MD) if the protein-surface interactions can be accurately modeled. To answer this question, we computed the adsorption free energies of peptides (representing eleven different amino acids) on a hydrophobic self-assembled monolayer (CH3-SAM) and compared them to the benchmark experimental data set. Our result revealed that existing biomolecular force fields, GAFF and AMBER ff14sb, cannot reproduce the experimental peptide adsorption free energies by Wei and Latour (Langmuir, 2009, 25, 5637-5646). To obtain the improved force fields, we systematically tuned the Lennard-Jones parameters of selected amino acid sidechains and the functional group of SAM with repeated metadynamics and umbrella sampling simulations. The final parameter set has yielded a significant improvement in the free energy values with R = 0.83 and MSE = 0.65 kcal/mol. We applied the refined force field to predict the initial adsorption orientation of lysozyme on CH3-SAM. Two major orientations-face-down and face-up-were predicted. Our analysis on the protein structure, solvent accessible surface area, and binding of native ligand NAG3 suggested that lysozyme in the face-up orientation can remain active after initial adsorption. However, because of its weaker affinity (ΔΔG = 7.86 kcal/mol) for the ligand, the bioactivity of the protein is expected to reduce. Our work facilitates the use of MD for the study of protein-SAM systems. The refined force field compatible with GROMACS is available at https://cbbio.cis.um.edu.mo/software/SAMFF .
Collapse
Affiliation(s)
- Pratiti Bhadra
- Department of Computer and Information Science , University of Macau , Avenida da Universidade , Taipa , Macau
| | - Shirley W I Siu
- Department of Computer and Information Science , University of Macau , Avenida da Universidade , Taipa , Macau
| |
Collapse
|
11
|
Nikitin A, Del Frate G. Development of Nonbonded Models for Metal Cations Using the Electronic Continuum Correction. J Comput Chem 2019; 40:2464-2472. [PMID: 31301182 DOI: 10.1002/jcc.26021] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Revised: 05/14/2019] [Accepted: 06/06/2019] [Indexed: 12/24/2022]
Abstract
The parametrization of classical nonbonded models of metal ions has been widely addressed in the recent years. Despite the continuous development of novel and more physically inspired functional forms, the 12-6 Lennard-Jones plus Coulomb potential is still the most adopted force field in molecular dynamics (MD) codes, owing to its simple form and easy implementation. However, due to the integer formal charge, unpolarizable force fields of ions may suffer from overestimated interatomic electrostatic interactions, leading to nonphysical clustering or repulsion between such full charges. The electronic continuum correction (ECC) can fix this problem through a simple inclusion of solvent polarization effects via ionic charge rescaling. In this work, the development of novel nonbonded models for mono, divalent, and highly charged metal ions is presented. For each metal species, the ionic charge has been scaled, according to the ECC. Lennard-Jones parameters have been optimized using experimental structural and thermodynamic properties as target quantities. Performances of the proposed models are discussed and compared with the literature data, while transferability attitudes among different and well-known water models are evaluated. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Alexei Nikitin
- Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Moscow, 119991, Russian Federation.,Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Gianluca Del Frate
- IMT School for Advanced Studies Lucca, Piazza S. Francesco 19, I-55100, Lucca, Italy
| |
Collapse
|
12
|
Gonçalves YMH, Senac C, Fuchs PFJ, Hünenberger PH, Horta BAC. Influence of the Treatment of Nonbonded Interactions on the Thermodynamic and Transport Properties of Pure Liquids Calculated Using the 2016H66 Force Field. J Chem Theory Comput 2019; 15:1806-1826. [DOI: 10.1021/acs.jctc.8b00425] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Yan M. H. Gonçalves
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Caroline Senac
- Sorbonne
Université,
CNRS, INSERM, Laboratoire d’Imagerie Biomédicale (LIB), F-75006 Paris, France
| | - Patrick F. J. Fuchs
- Sorbonne Université,
École normale supérieure, PSL University, CNRS, Laboratoire
des biomolécules, LBM, 75005 Paris, France
- Université
Paris Diderot, 75013 Paris, France
| | | | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| |
Collapse
|
13
|
Hagler AT. Force field development phase II: Relaxation of physics-based criteria… or inclusion of more rigorous physics into the representation of molecular energetics. J Comput Aided Mol Des 2018; 33:205-264. [DOI: 10.1007/s10822-018-0134-x] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2018] [Accepted: 07/18/2018] [Indexed: 01/04/2023]
|
14
|
Del Frate G, Nikitin A. Including Electronic Screening in Classical Force Field of Zinc Ion for Biomolecular Simulations. ChemistrySelect 2018. [DOI: 10.1002/slct.201802864] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Gianluca Del Frate
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa Italy
- Present address: IMT School for Advanced Studies Lucca Piazza S.Francesco 19 Lucca 55100 Italy
| | - Alexei Nikitin
- Engelhardt Institute of Molecular BiologyRussian Academy of Sciences Moscow 119991 Russia
- Scuola Normale Superiore Piazza dei Cavalieri 7 6126 Pisa Italy
| |
Collapse
|
15
|
Ermilova I, Lyubartsev AP. Extension of the Slipids Force Field to Polyunsaturated Lipids. J Phys Chem B 2016; 120:12826-12842. [DOI: 10.1021/acs.jpcb.6b05422] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Inna Ermilova
- Department of Materials and
Environmental Chemistry, Stockholm University, SE 106 91 Stockholm, Sweden
| | - Alexander P. Lyubartsev
- Department of Materials and
Environmental Chemistry, Stockholm University, SE 106 91 Stockholm, Sweden
| |
Collapse
|
16
|
Horta BAC, Merz PT, Fuchs PFJ, Dolenc J, Riniker S, Hünenberger PH. A GROMOS-Compatible Force Field for Small Organic Molecules in the Condensed Phase: The 2016H66 Parameter Set. J Chem Theory Comput 2016; 12:3825-50. [DOI: 10.1021/acs.jctc.6b00187] [Citation(s) in RCA: 88] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bruno A. C. Horta
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Pascal T. Merz
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Patrick F. J. Fuchs
- Institut Jacques Monod, UMR 7592 CNRS, Université Paris-Diderot, Sorbonne Paris Cité, F-75205 Paris, France
| | - Jozica Dolenc
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Chemistry,
Biology and Pharmacy Information Center, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | | |
Collapse
|
17
|
The good, the bad and the user in soft matter simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2529-2538. [PMID: 26862882 DOI: 10.1016/j.bbamem.2016.02.004] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/16/2015] [Revised: 02/01/2016] [Accepted: 02/02/2016] [Indexed: 11/21/2022]
Abstract
Molecular dynamics (MD) simulations have become popular in materials science, biochemistry, biophysics and several other fields. Improvements in computational resources, in quality of force field parameters and algorithms have yielded significant improvements in performance and reliability. On the other hand, no method of research is error free. In this review, we discuss a few examples of errors and artifacts due to various sources and discuss how to avoid them. Besides bringing attention to artifacts and proper practices in simulations, we also aim to provide the reader with a starting point to explore these issues further. In particular, we hope that the discussion encourages researchers to check software, parameters, protocols and, most importantly, their own practices in order to minimize the possibility of errors. The focus here is on practical issues. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.
Collapse
|