1
|
Odinokov A, Son WJ, Yakubovich A, Park JY, Jung Y. Ab Initio Prediction of Vapor Pressure for Diverse Atomic Layer Deposition Precursors. J Chem Theory Comput 2024; 20:6144-6151. [PMID: 38991151 DOI: 10.1021/acs.jctc.3c01416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
Understanding the saturated vapor pressure (Pvap) is vital for evaluating atomic layer deposition (ALD) precursors, as it directly influences the ALD temperature window and, by extension, the processability of compounds. The early estimation of vapor pressure ranges is crucial during the initial stages of novel precursor design, reducing the reliance on empirical synthesis or experimentation. However, predicting vapor pressure through computer simulations is often impeded by the scarcity of suitable empirical force fields for molecular dynamics simulations. This challenge is further compounded by the diverse chemical substances and the introduction of new elements into modern ALD processes, necessitating robust force fields that can accommodate metals, organics, and halides. In response, this study introduces a novel approach utilizing a quantum mechanically derived force field for the prediction of vapor pressure across a wide spectrum of potential ALD precursors. This approach enables the creation of system-specific force fields through parametrization based on ab initio calculations for a single molecule. We develop a comprehensive workflow to simulate both liquid and gaseous equilibrium phases, allowing the calculation of vapor pressure across a wide temperature range. Our methodology has been validated with a diverse set of ALD precursors, demonstrating its robustness in predicting Pvap at specified temperatures. The approach yields a Pearson's correlation coefficient (R2) greater than 0.9 on a logarithmic scale and a root-mean-squared deviation in self-solvation-free energies as low as 1.3 kcal mol-1. This innovative workflow, which does not require any prior experimental data, marks a significant advancement in the computer-aided design of novel ALD precursors, paving the way for accelerating developments in technology.
Collapse
Affiliation(s)
- Alexey Odinokov
- Samsung Advanced Institute of Technology, Samsung Electronics, 130 Samsung-ro, Suwon 16678, Korea
| | - Won-Joon Son
- CSE Team, Innovation Center, Samsung Electronics, Hwaseong 18448, Korea
| | - Alexander Yakubovich
- Samsung Advanced Institute of Technology, Samsung Electronics, 130 Samsung-ro, Suwon 16678, Korea
| | - Ji Young Park
- CSE Team, Innovation Center, Samsung Electronics, Hwaseong 18448, Korea
| | - Yongsik Jung
- Samsung Advanced Institute of Technology, Samsung Electronics, 130 Samsung-ro, Suwon 16678, Korea
| |
Collapse
|
2
|
Cabral VÁ, Govoni B, Verli H. Unravelling carbohydrate binding module 21 (CBM21) dynamics of interaction with amylose. Carbohydr Polym 2024; 330:121792. [PMID: 38368081 DOI: 10.1016/j.carbpol.2024.121792] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Revised: 12/21/2023] [Accepted: 01/05/2024] [Indexed: 02/19/2024]
Abstract
The carbohydrate binding module 21 (CBM21) from Rhizopus oryzae is a dual-site CBM proposed to disrupt polysaccharide structures. Additionally, it serves as a purification tag in industry. CBM21 crystal structure features a Glc residue in an unusual 1S3 conformation, whose relevance for the CBM mechanism of action is unclear. In this context, we seek to contribute for the understanding of CBM21 mechanism of action by: i) investigating the role of the 1S3 conformation on carbohydrate recognition, and ii) characterize the protein-carbohydrate binding dynamics using molecular dynamics and metadynamics simulations at MM and QM/MM levels. Results indicate the 1S3 Glc conformation is unlikely to occur under biological conditions, being originated from the crystallographic environment. CBM21 binding to small ligands appears transient and unstable, while protein dimerization and polysaccharide chain size influence complex stability. In interactions with amylose, CBM21 exhibits a repeated unbinding followed by re-binding, while simultaneously alternating between binding sites I and II. These results suggest that CBM21 acts through transient interactions, directing carbohydrates to the catalytic center rather than forming strong and long-lasting bonds with carbohydrates. Accordingly, we expect such atomistic depiction of CBM21 mechanism could aid in CBM design targeting biotechnological applications.
Collapse
Affiliation(s)
- Vinicius Ávila Cabral
- Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul, Av. Bento Gonçalves, 9500, Porto Alegre 91500-970, RS, Brazil
| | - Bruna Govoni
- Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul, Av. Bento Gonçalves, 9500, Porto Alegre 91500-970, RS, Brazil
| | - Hugo Verli
- Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul, Av. Bento Gonçalves, 9500, Porto Alegre 91500-970, RS, Brazil.
| |
Collapse
|
3
|
Roccatano D. A molecular dynamics simulation study of glycine/serine octapeptides labeled with 2,3-diazabicyclo[2.2.2]oct-2-ene fluorophore. J Chem Phys 2024; 160:145101. [PMID: 38587229 DOI: 10.1063/5.0190073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2023] [Accepted: 03/19/2024] [Indexed: 04/09/2024] Open
Abstract
The compound 2,3-diazabicyclo[2.2.2]oct-2-ene (DBO) is a versatile fluorophore widely used in Förster resonance energy transfer (FRET) spectroscopy studies due to its remarkable sensitivity, enabling precise donor-acceptor distance measurements, even for short peptides. Integrating time-resolved and FRET spectroscopies with molecular dynamics simulations provides a robust approach to unravel the structure and dynamics of biopolymers in a solution. This study investigates the structural behavior of three octapeptide variants: Trp-(Gly-Ser)3-Dbo, Trp-(GlyGly)3-Dbo, and Trp-(SerSer)3-Dbo, where Dbo represents the DBO-containing modified aspartic acid, using molecular dynamics simulations. Glycine- and serine-rich amino acid fragments, common in flexible protein regions, play essential roles in functional properties. Results show excellent agreement between end-to-end distances, orientational factors from simulations, and the available experimental and theoretical data, validating the reliability of the GROMOS force field model. The end-to-end distribution, modeled using three Gaussian distributions, reveals a complex shape, confirmed by cluster analysis highlighting a limited number of significant conformations dominating the peptide landscape. All peptides predominantly adopt a disordered state in the solvent, yet exhibit a compact shape, aligning with the model of disordered polypeptide chains in poor solvents. Conformations show marginal dependence on chain composition, with Ser-only chains exhibiting slightly more elongation. This study enhances our understanding of peptide behavior, providing valuable insights into their structural dynamics in solution.
Collapse
Affiliation(s)
- Danilo Roccatano
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln LN6 7TS, United Kingdom
| |
Collapse
|
4
|
Fontana C, de Meirelles JL, Verli H. Theoretical models of staurosporine and analogs uncover detailed structural information in biological solution. J Mol Graph Model 2024; 126:108653. [PMID: 37922640 DOI: 10.1016/j.jmgm.2023.108653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 10/12/2023] [Accepted: 10/12/2023] [Indexed: 11/07/2023]
Abstract
Staurosporine and its analogs (STA-analogs) are indolocarbazoles (ICZs) compounds able to inhibit kinase proteins in a non-specific way, while present antimicrobial and cytostatic properties. The knowledge of molecular features associated to the complexation, including the ligand shape in solution and thermodynamics of complexation, is substantial to the development of new bioactive ICZs with improved therapeutic properties. In this context, the empirical approach of GROMOS force field is able to accurately reproduce condensed phase physicochemical properties of molecular systems after parameterization. Hence, through parameterization under GROMOS force field and molecular simulations, we assessed STA-analogs dynamics in aqueous solution, as well as its interaction with water to probe conformational and structural features involved in complexation to therapeutic targets. The coexistence of multiple conformers observed in simulations, and confirmed by metadynamics calculations, expanding the conformational space knowledge of these ligands with potential implications in understanding the ligand conformational selection during complexation. Also, changes in availability to H-bonding concerning the different substituents and water can reflect on effects at complexation free energy due to variation at the desolvation energetic costs. Based on these results, we expect the obtained structural data provide systemic framework for rational chemical modification of STA-analogs.
Collapse
Affiliation(s)
- Crisciele Fontana
- Universidade Federal do Rio Grande do Sul, Centro de Biotecnologia, Av. Bento Gonçalves, 9500 (Caixa Postal 15005), Porto Alegre, CEP 91501-970, RS, Brazil
| | - João Luiz de Meirelles
- Universidade Federal do Rio Grande do Sul, Centro de Biotecnologia, Av. Bento Gonçalves, 9500 (Caixa Postal 15005), Porto Alegre, CEP 91501-970, RS, Brazil
| | - Hugo Verli
- Universidade Federal do Rio Grande do Sul, Centro de Biotecnologia, Av. Bento Gonçalves, 9500 (Caixa Postal 15005), Porto Alegre, CEP 91501-970, RS, Brazil.
| |
Collapse
|
5
|
Gonçalves YMH, Horta BAC. gmak: A Parameter-Space Mapping Strategy for Force-Field Calibration. J Chem Theory Comput 2023; 19:605-618. [PMID: 36634285 DOI: 10.1021/acs.jctc.2c00955] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
In the context of classical molecular simulations, the accuracy of a force field is highly influenced by the values of the relevant simulation parameters. In this work, a parameter-space mapping (PSM) workflow is proposed to aid in the calibration of force-field parameters, based mainly on the following features: (i) regular-grid discretization of the search space; (ii) partial sampling of the search-space grid; (iii) training of surrogate models to predict the estimates of the target properties for nonsampled parameter sets; (iv) post hoc interpretation of the results in terms of multiobjective optimization concepts; (v) attenuation of statistical errors achieved via empiric extension of the duration of the simulations; (vi) iterative search-space translation according to a user-defined scalar objective function that measures the accuracy of the force field (e.g., the weighted root-mean-square deviation of the target properties relative to the reference data). This combination of features results in a hybrid of a single- and a multiobjective optimization strategy, allowing for the approximate determination of both a local minimum of the chosen objective function and its neighboring Pareto efficient points. The PSM workflow is implemented in the extensible Python program gmak, which is made available in the Git repository at http://github.com/mssm-labmmol/gmak. Using this implementation, the PSM workflow was tested in a proof-of-concept fashion in the recalibration of the Lennard-Jones parameters of the 3-point Optimal Point Charge (OPC3) water model for compatibility with the GROMOS treatment of nonbonded interactions. The recalibrated model reproduces typical pure-liquid properties with an accuracy similar to the original OPC3 model and represents a significant improvement relative to the Simple Point Charge (SPC) model, which is the official recommendation for simulations using GROMOS force fields.
Collapse
Affiliation(s)
- Yan M H Gonçalves
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
- Peers Consulting & Technology, Av. Ibirapuera, 1753-18° andar, Moema, São Paulo, São Paulo 04029-90, Brazil
| | - Bruno A C Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
- Peers Consulting & Technology, Av. Ibirapuera, 1753-18° andar, Moema, São Paulo, São Paulo 04029-90, Brazil
- Laboratory of Applied Intelligence, University of Vale do Itajaí, Itajaí, Santa Catarina 88302-901, Brazil
| |
Collapse
|
6
|
Oliveira MP, Hünenberger PH. Force fields optimized against experimental data for large compound families using CombiFF: Validation considering non-target properties and polyfunctional compounds. J Mol Graph Model 2023; 118:108312. [PMID: 36252318 DOI: 10.1016/j.jmgm.2022.108312] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Revised: 08/09/2022] [Accepted: 08/16/2022] [Indexed: 10/14/2022]
Abstract
The CombiFF scheme is a workflow for the automated calibration of force-field parameters against condensed-phase experimental data considering simultaneously entire classes of organic molecules. The main steps of this scheme are: (i) selection of a molecule family; (ii) enumeration of all isomers; (iii) query for experimental data; (iv) automatic construction of the molecular topologies; (v) iterative refinement of the force-field parameters considering the entire family. In two recent articles, CombiFF was applied to the design of GROMOS-compatible united-atom force fields for the saturated acyclic haloalkanes and for saturated acyclic compounds involving eight common chemical functional groups of oxygen and nitrogen. This calibration and the subsequent initial validation involved two limitations: (i) the experimental data considered was restricted to values for the pure-liquid density ρliq and the vaporization enthalpy ΔHvap of the compounds; (ii) beyond monofunctional compounds, the training set only involved homo-polyhaloalkanes (possibly mixing halogen types) in the first study, and homo-polyfunctional compounds of the considered oxygen or nitrogen functional groups (no mixing of different group types) in the second one. The goal of this article is to further test the accuracy of CombiFF-generated force fields by extending the validation to: (i) nine additional properties that were not used as optimization targets (pure-liquid thermodynamic, dielectric and transport properties, as well as solvation properties); (ii) hetero-polyfunctional molecules that were not included in the calibration and initial validation sets. The results for the nine additional properties show good agreement with experiment, except for the shear viscosity and the dielectric permittivity. There, larger discrepancies are observed, likely due to the united-atom representation adopted for the aliphatic groups and to the implicit treatment of electronic polarization effects. The results for the hetero-polyfunctional molecules also show reasonable agreement with experiment in terms of the monitored properties.
Collapse
Affiliation(s)
- Marina P Oliveira
- 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
|
7
|
Valente ÉC, Polêto MD, de Oliveira TV, Soares LDS, dos Reis Coimbra JS, Guimarães AP, de Oliveira EB. Effects of the Cations Li+, Na+, K+, Mg2+, or Ca2+ on Physicochemical Properties of Xanthan Gum in Aqueous Medium – A view from Computational Molecular Dynamics Calculations. FOOD BIOPHYS 2022. [DOI: 10.1007/s11483-022-09773-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
8
|
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
|
9
|
Capability of MXene 2D material as an amoxicillin, ampicillin, and cloxacillin adsorbent in wastewater. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118545] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
|
10
|
Gonçalves YMH, Kashefolgheta S, Oliveira MP, Hünenberger PH, Horta BAC. Simultaneous parametrization of torsional and third-neighbor interaction terms in force-field development: The LLS-SC algorithm. J Comput Chem 2022; 43:644-653. [PMID: 35133016 DOI: 10.1002/jcc.26819] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2021] [Revised: 01/11/2022] [Accepted: 01/13/2022] [Indexed: 01/10/2023]
Abstract
The calibration of torsional interaction terms by fitting relative gas-phase conformational energies against their quantum-mechanical values is a common procedure in force-field development. However, much less attention has been paid to the optimization of third-neighbor nonbonded interaction parameters, despite their strong coupling with the torsions. This article introduces an algorithm termed LLS-SC, aimed at simultaneously parametrizing torsional and third-neighbor interaction terms based on relative conformational energies. It relies on a self-consistent (SC) procedure where each iteration involves a linear least-squares (LLS) regression followed by a geometry optimization of the reference structures. As a proof-of-principle, this method is applied to obtain torsional and third-neighbor interaction parameters for aliphatic chains in the context of the GROMOS 53A6 united-atom force field. The optimized parameter set is compared to the original one, which has been fitted manually against thermodynamic properties for small linear alkanes. The LLS-SC implementation is freely available under http://github.com/mssm-labmmol/profiler.
Collapse
Affiliation(s)
- Yan M H Gonçalves
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | | | | | | | - Bruno A C Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| |
Collapse
|
11
|
Hajibabaei A, Ha M, Pourasad S, Kim J, Kim KS. Machine Learning of First-Principles Force-Fields for Alkane and Polyene Hydrocarbons. J Phys Chem A 2021; 125:9414-9420. [PMID: 34657427 DOI: 10.1021/acs.jpca.1c05819] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Machine learning (ML) interatomic potentials (ML-IAPs) are generated for alkane and polyene hydrocarbons using on-the-fly adaptive sampling and a sparse Gaussian process regression (SGPR) algorithm. The ML model is generated based on the PBE+D3 level of density functional theory (DFT) with molecular dynamics (MD) for small alkane and polyene molecules. Intermolecular interactions are also trained with clusters and condensed phases of small molecules. It shows excellent transferability to long alkanes and closely describes the ab inito potential energy surface for polyenes. Simulation of liquid ethane also shows reasonable agreement with experimental reports. This is a promising initiative toward a universal ab initio quality force-field for hydrocarbons and organic molecules.
Collapse
Affiliation(s)
- Amir Hajibabaei
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Korea
| | - Miran Ha
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Korea
| | - Saeed Pourasad
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Korea
| | - Junu Kim
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Korea
| | - Kwang S Kim
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Korea
| |
Collapse
|
12
|
P. Oliveira M, Hünenberger PH. Systematic optimization of a fragment-based force field against experimental pure-liquid properties considering large compound families: application to oxygen and nitrogen compounds. Phys Chem Chem Phys 2021; 23:17774-17793. [PMID: 34350931 PMCID: PMC8386690 DOI: 10.1039/d1cp02001c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 06/30/2021] [Indexed: 12/04/2022]
Abstract
The CombiFF approach is a workflow for the automated refinement of force-field parameters against experimental condensed-phase data, considering entire classes of organic molecules constructed using a fragment library via combinatorial isomer enumeration. One peculiarity of this approach is that it relies on an electronegativity-equalization scheme to account for induction effects within molecules, with values of the atomic hardness and electronegativity as electrostatic parameters, rather than the partial charges themselves. In a previous article [M. P. Oliveira, M. Andrey, S. R. Rieder, L. Kern, D. F. Hahn, S. Riniker, B. A. C. Horta and P. H. Hünenberger, J. Chem. Theory. Comput. 2020, 16, 7525], CombiFF was introduced and applied to calibrate a GROMOS-compatible united-atom force field for the saturated acyclic (halo-)alkane family. Here, this scheme is employed for the construction of a corresponding force field for saturated acyclic compounds encompassing eight common chemical functional groups involving oxygen and/or nitrogen atoms, namely: ether, aldehyde, ketone, ester, alcohol, carboxylic acid, amine, and amide. Monofunctional as well as homo-polyfunctional compounds are considered. A total of 1712 experimental liquid densities ρliq and vaporization enthalpies ΔHvap concerning 1175 molecules are used for the calibration (339 molecules) and validation (836 molecules) of the 102 non-bonded interaction parameters of the force field. Using initial parameter values based on the GROMOS 2016H66 parameter set, convergence is reached after five iterations. Given access to one processor per simulated system, this operation only requires a few days of wall-clock computing time. After optimization, the root-mean-square deviations from experiment are 29.9 (22.4) kg m-3 for ρliq and 4.1 (5.5) kJ mol-1 for ΔHvap for the calibration (validation) set. Thus, a very good level of agreement with experiment is achieved in terms of these two properties, although the errors are inhomogeneously distributed across the different chemical functional groups.
Collapse
Affiliation(s)
- Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 5503
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 5503
| |
Collapse
|
13
|
Kashefolgheta S, Wang S, Acree WE, Hünenberger PH. Evaluation of nine condensed-phase force fields of the GROMOS, CHARMM, OPLS, AMBER, and OpenFF families against experimental cross-solvation free energies. Phys Chem Chem Phys 2021; 23:13055-13074. [PMID: 34105547 PMCID: PMC8207520 DOI: 10.1039/d1cp00215e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Accepted: 04/28/2021] [Indexed: 12/02/2022]
Abstract
Experimental solvation free energies are nowadays commonly included as target properties in the validation of condensed-phase force fields, sometimes even in their calibration. In a previous article [Kashefolgheta et al., J. Chem. Theory. Comput., 2020, 16, 7556-7580], we showed how the involved comparison between experimental and simulation results could be made more systematic by considering a full matrix of cross-solvation free energies . For a set of N molecules that are all in the liquid state under ambient conditions, such a matrix encompasses N×N entries for considering each of the N molecules either as solute (A) or as solvent (B). In the quoted study, a cross-solvation matrix of 25 × 25 experimental value was introduced, considering 25 small molecules representative for alkanes, chloroalkanes, ethers, ketones, esters, alcohols, amines, and amides. This experimental data was used to compare the relative accuracies of four popular condensed-phase force fields, namely GROMOS-2016H66, OPLS-AA, AMBER-GAFF, and CHARMM-CGenFF. In the present work, the comparison is extended to five additional force fields, namely GROMOS-54A7, GROMOS-ATB, OPLS-LBCC, AMBER-GAFF2, and OpenFF. Considering these nine force fields, the correlation coefficients between experimental values and simulation results range from 0.76 to 0.88, the root-mean-square errors (RMSEs) from 2.9 to 4.8 kJ mol-1, and average errors (AVEEs) from -1.5 to +1.0 kJ mol-1. In terms of RMSEs, GROMOS-2016H66 and OPLS-AA present the best accuracy (2.9 kJ mol-1), followed by OPLS-LBCC, AMBER-GAFF2, AMBER-GAFF, and OpenFF (3.3 to 3.6 kJ mol-1), and then by GROMOS-54A7, CHARM-CGenFF, and GROMOS-ATB (4.0 to 4.8 kJ mol-1). These differences are statistically significant but not very pronounced, and are distributed rather heterogeneously over the set of compounds within the different force fields.
Collapse
Affiliation(s)
- Sadra Kashefolgheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| | - Shuzhe Wang
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| | - William E. Acree
- Department of Chemistry, University of North Texas1155 Union Circle Drive #305070DentonTexas 76203USA
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| |
Collapse
|
14
|
Ramos MC, Horta BAC. Drug-Loading Capacity of PAMAM Dendrimers Encapsulating Quercetin Molecules: A Molecular Dynamics Study with the 2016H66 Force Field. J Chem Inf Model 2021; 61:987-1000. [PMID: 33502188 DOI: 10.1021/acs.jcim.0c00960] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The complexation of quercetin molecules with poly(amidoamine) (PAMAM) dendrimers of generation 0-3 was studied by molecular dynamics simulations. Three main points were addressed: (i) the effect of starting from different initial structures; (ii) the performance of the 2016H66 force field (recently validated in the context of dendrimer simulations) in predicting the experimental drug(quercetin)-loading capacity of PAMAM dendrimers; and (iii) the stability of quercetin-PAMAM complexes and their interactions. Initial structures generated by different restraint protocols led to faster convergence compared to initial structures generated by randomly placing the drug molecules in the simulation box. The simulations yielded meta-stable complexes where the loading numbers have converged to average values and were compared to experimentally obtained values. Once the first meta-stable state was reached, the drug-dendrimer complexes did not deviate significantly throughout the simulation. They were characterized in terms of structural properties, such as the radius of gyration and radial distribution functions. The results suggest that quercetin molecules interact mostly with the internal dendrimer monomers rather than to their surface.
Collapse
Affiliation(s)
- Mayk C Ramos
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Bruno A C Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| |
Collapse
|
15
|
van Gunsteren WF, Daura X, Fuchs PFJ, Hansen N, Horta BAC, Hünenberger PH, Mark AE, Pechlaner M, Riniker S, Oostenbrink C. On the Effect of the Various Assumptions and Approximations used in Molecular Simulations on the Properties of Bio-Molecular Systems: Overview and Perspective on Issues. Chemphyschem 2020; 22:264-282. [PMID: 33377305 DOI: 10.1002/cphc.202000968] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Indexed: 12/14/2022]
Abstract
Computer simulations of molecular systems enable structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic or supra-molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio-molecular systems are considered. What the developers of bio-molecular force fields and simulation software can do to facilitate and broaden research involving bio-molecular simulations is also discussed.
Collapse
Affiliation(s)
- Wilfred F van Gunsteren
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Xavier Daura
- Institute of Biotechnology and Biomedicine, Universitat Autonoma de Barcelona (UAB), 08193, Barcelona, Spain.,Catalan Institution for Research and Advanced Studies (ICREA), 08010, Barcelona, Spain
| | - Patrick F J Fuchs
- Sorbonne Université, Ecole Normale Supérieure, PSL Research University, CNRS, Laboratoire des Biomolécules (LBM), F-75005, Paris, France.,Université de Paris, UFR Sciences du Vivant, F-75013, Paris, France
| | - Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569, Stuttgart, Germany
| | - Bruno A C Horta
- Instituto de Química, Universidade Federal de Rio de Janeiro, Rio de Janeiro, 21941-909, Brazil
| | - Philippe H Hünenberger
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Alan E Mark
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, QLD, 4072, Australia
| | - Maria Pechlaner
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Chris Oostenbrink
- Institute of Molecular Modelling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| |
Collapse
|
16
|
Diem M, Oostenbrink C. The effect of different cutoff schemes in molecular simulations of proteins. J Comput Chem 2020; 41:2740-2749. [PMID: 33026106 PMCID: PMC7756334 DOI: 10.1002/jcc.26426] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 09/06/2020] [Accepted: 09/09/2020] [Indexed: 01/02/2023]
Abstract
Molecular simulations of nanoscale systems invariably involve assumptions and approximations to describe the electrostatic interactions, which are long-ranged in nature. One approach is the use of cutoff schemes with a reaction-field contribution to account for the medium outside the cutoff scheme. Recent reports show that macroscopic properties may depend on the exact choice of cutoff schemes in modern day simulations. In this work, a systematic analysis of the effects of different cutoff schemes was performed using a set of 52 proteins. We find no statistically significant differences between using a twin-range or a single-range cutoff scheme. Applying the cutoff based on charge groups or based on atomic positions, does lead to significant differences, which is traced to the cutoff noise for energies and forces. While group-based cutoff schemes show increased cutoff noise in the potential energy, applying an atomistic cutoff leads to artificial structure in the solvent at the cutoff distance. Carefully setting the temperature control, or using an atomistic cutoff for the solute and a group-based cutoff for the solvent significantly reduces the effects of the cutoff noise, without introducing structure in the solvent. This study aims to deepen the understanding of the implications different cutoffs have on molecular dynamics simulations.
Collapse
Affiliation(s)
- Matthias Diem
- Institute for Molecular Modeling and SimulationUniversity of Natural Resources and Life SciencesViennaAustria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and SimulationUniversity of Natural Resources and Life SciencesViennaAustria
| |
Collapse
|
17
|
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
|
18
|
Kashefolgheta S, Oliveira MP, Rieder SR, Horta BAC, Acree WE, Hünenberger PH. Evaluating Classical Force Fields against Experimental Cross-Solvation Free Energies. J Chem Theory Comput 2020; 16:7556-7580. [DOI: 10.1021/acs.jctc.0c00688] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Sadra Kashefolgheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Marina P. Oliveira
- 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
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - William E. Acree
- Department of Chemistry, University of North Texas, 1155 Union Circle Drive #305070, Denton, Texas 76203, United States
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
19
|
Diem M, Oostenbrink C. The Effect of Using a Twin-Range Cutoff Scheme for Nonbonded Interactions: Implications for Force-Field Parametrization? J Chem Theory Comput 2020; 16:5985-5990. [PMID: 32813524 PMCID: PMC7558304 DOI: 10.1021/acs.jctc.0c00509] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2020] [Indexed: 02/01/2023]
Abstract
Recently, concerns have been voiced regarding the validity of the GROMOS force fields, being parametrized using a twin-range cutoff scheme, in which longer ranged nonbonded forces and energies are updated less frequently than shorter ranged ones. Here we demonstrate that the influence of such a scheme on the thermodynamic, structural, and dynamic properties used in the parametrization of the GROMOS force fields is minor. We find root-mean-square differences of maximally 0.5 kJ/mol for the solvation free energy and heat of vaporization and of maximally 0.4% for the density. Slightly larger differences are observed when switching from a group-based to an atom-based cutoff scheme. In cases where the twin-range cutoff scheme does result in minor differences compared to a single-range cutoff these are well within the deviation from the experimentally measured values.
Collapse
Affiliation(s)
- Matthias Diem
- Institute for Molecular Modeling
and Simulation, University of Natural Resources
and Life Sciences, 1190 Vienna, Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling
and Simulation, University of Natural Resources
and Life Sciences, 1190 Vienna, Austria
| |
Collapse
|
20
|
Pannuzzo M, Horta BAC, La Rosa C, Decuzzi P. Predicting the Miscibility and Rigidity of Poly(lactic- co-glycolic acid)/Polyethylene Glycol Blends via Molecular Dynamics Simulations. Macromolecules 2020; 53:3643-3654. [PMID: 32831403 PMCID: PMC7428138 DOI: 10.1021/acs.macromol.0c00110] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 04/07/2020] [Indexed: 11/30/2022]
Abstract
![]()
The
addition of polyethylene glycol (PEG) chains to poly(lactic-co-glycolic acid) (PLGA) matrices is extensively used to
modulate the biodegradation, drug loading and release, mechanical
properties, and chemical stability of the original system. Multiple
parameters, including the molecular weight, relative concentration,
polarity, and solubility, affect the physicochemical properties of
the polymer blend. Here, molecular dynamics simulations with the united-atom
2016H66 force field are used to model the behavior of PLGA and PEG chains and thus predict the overall
physicochemical features of the resulting blend. First, the model
accuracy is validated against fundamental properties of pure PLGA
and PEG samples. In agreement with previous experimental and theoretical
observations, the PLGA solubility results to be higher in acetonitrile
than in water, with Flory parameters νACN = 0.63
± 0.01 and νW = 0.21 ± 0.02, and the Young’s
modulus of PLGA and PEG equal to Y = 2.0 ± 0.43
and 0.32 ± 0.34 GPa, respectively. Next, four PEG/PLGA blending
regimes are identified by varying the relative concentrations and
molecular weights of the individual polymers. The computational results
demonstrate that at low PEG concentrations (<8% w/w), homogeneous
blends are generated for both low and high PEG molecular weights.
In contrast, at comparable PEG and PLGA concentrations (∼50%
w/w), short PEG chains are only partially miscible whereas long PEG
chains segregate within the PLGA matrix. This behavior has been confirmed
experimentally via differential scanning calorimetry and is in agreement
with previous observations. Finally, the computed Young’s modulus
of PLGA/PEG blends is observed to decrease with the PEG content returning
the lowest values for the partial and fully segregated regimens (Y ≈ 1.3 GPa). This work proposes a computational
scheme for predicting the physicochemical properties of PLGA/PEG blends
paving the way toward the rational design of polymer mixtures for
biomedical applications.
Collapse
Affiliation(s)
- Martina Pannuzzo
- Laboratory of Nanotechnology for Precision Medicine, Fondazione Istituto Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| | - Bruno A C Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Carmelo La Rosa
- Department of Chemical Sciences, University of Catania, Viale A. Doria 6, Catania 95125, Italy
| | - Paolo Decuzzi
- Laboratory of Nanotechnology for Precision Medicine, Fondazione Istituto Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| |
Collapse
|
21
|
De Oliveira TV, Guimarães AP, Bressan GC, Maia ER, Coimbra JSDR, Polêto MD, De Oliveira EB. Structural and molecular bases of angiotensin-converting enzyme inhibition by bovine casein-derived peptides: an in silico molecular dynamics approach. J Biomol Struct Dyn 2020; 39:1386-1403. [PMID: 32066337 DOI: 10.1080/07391102.2020.1730243] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The angiotensin-converting enzyme (ACE) plays a key role in blood pressure regulation process, and its inhibition is one of the main drug targets for the treatment of hypertension. Though various peptides from milk proteins are well-known for their ACE-inhibitory capacity, research devoted to understand the molecular bases of such property remain scarce, specifically for such peptides. Therefore, in this work, computational molecular docking and molecular dynamics calculations were performed to enlighten the intermolecular interactions involved in ACE inhibition by six different casein-derived peptides (FFVAPFPEVFGK, FALPQYLK, ALNEINQFYQK, YLGYLEQLLR, HQGLPQEVLNENLLR and NAVPITPTLNR). Two top ranked docking poses for each peptide (one with N- and the other C-terminal peptide extremity oriented towards the ACE active site) were selected for dynamic simulations (50 ns; GROMOS53A6 force field), and the results were correlated to in vitro ACE inhibition capacity. Two molecular features appeared to be essential for peptides to present high ACE inhibition capacity in vitro: i) to interact with the S1 active site residues (Ala354, Glu384, and Tyr523) by hydrogen bonds; ii) to interact with Zn2+ coordinated residues (His383, His387, and Glu411) by short-lenght hydrogen bonds, as observed in the cases of ALNEINQFYQK (IACE = 80.7%), NAVPITPTLNR (IACE = 80.7%), and FALPQYLK (IACE = 79.0%). Regardless of the temporal stability of these strong interactions, they promoted some disruption of Zn2+ tetrahedral coordination during the molecular dynamics trajectories, and were pointed as the main reason for the greatest ACE inhibition by these peptides. On the other hand, peptides with intermediate inhibition capacity (50% < IACE < 45%) interacted mainly by weaker interactions (e.g.: electrostatic and hydrophobic) with the Zn2+ coordinated residues, and were not able to change significantly its tetrahedral coordination structure. These findings may: i) assist the discrimination in silico of "good" and "bad" ACE-inhibitory peptides from other food sources, and/or ii) aid in designing de novo new molecules with ACE-inhibitory capacity. Communicated by Ramaswamy Sarma.
Collapse
Affiliation(s)
| | - Ana Paula Guimarães
- Departamento de Química (DEQ), Universidade Federal de Viçosa (UFV), Viçosa, MG, Brazil
| | - Gustavo Costa Bressan
- Departamento de Bioquímica e Biologia Molecular (DBB), Universidade Federal de Viçosa (UFV), Viçosa, MG, Brazil
| | - Elaine Rose Maia
- Laboratório de Estudos Estruturais Moleculares (LEEM), Instituto de Química, Universidade de Brasília, Brasília, DF, Brazil
| | | | - Marcelo Depólo Polêto
- Departamento de Biologia Geral (DBG), Universidade Federal de Viçosa (UFV), Viçosa, MG, Brazil
| | | |
Collapse
|
22
|
Arantes PR, Pedebos C, Polêto MD, Pol-Fachin L, Verli H. The Lazy Life of Lipid-Linked Oligosaccharides in All Life Domains. J Chem Inf Model 2020; 60:631-643. [PMID: 31769974 DOI: 10.1021/acs.jcim.9b00904] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Lipid-linked oligosaccharides (LLOs) play an important role in the N-glycosylation pathway as the donor substrate of oligosaccharyltransferases (OSTs), which are responsible for the en bloc transfer of glycan chains onto a nascent polypeptide. The lipid component of LLO in both eukarya and archaea consists of a dolichol, and an undecaprenol in prokarya, whereas the number of isoprene units may change between species. Given the potential relevance of LLOs and their related enzymes to diverse biotechnological applications, obtaining reliable LLO models from distinct domains of life could support further studies on complex formation and their processing by OSTs, as well as protein engineering on such systems. In this work, molecular modeling techniques, such as quantum mechanics calculations, molecular dynamics simulations, and metadynamics were employed to study eukaryotic (Glc3-Man9-GlcNAc2-PP-Dolichol), bacterial (Glc1-GalNAc5-Bac1-PP-Undecaprenol), and archaeal (Glc1-Man1-Gal1-Man1-Glc1-Gal1-Glc1-P-Dolichol) LLOs in membrane bilayers. Microsecond molecular dynamics simulations and metadynamics calculations of LLOs revealed that glycan chains are more prone to interact with the membrane lipid head groups, while the PP linkages are positioned at the lipid phosphate head groups level. The dynamics of isoprenoid chains embedded within the bilayer are described, and membrane dynamics and related properties are also investigated. Overall, there are similarities regarding the structure and dynamics of the eukaryotic, the bacterial, and the archaeal LLOs in bilayers, which can support the comprehension of their association with OSTs. These data may support future studies on the transferring mechanism of the oligosaccharide chain to an acceptor protein.
Collapse
Affiliation(s)
- Pablo R Arantes
- Centro de Biotecnologia , Universidade Federal do Rio Grande do Sul , Porto Alegre , RS 91509-900 , Brazil.,Laboratory of Genetic Toxicology , Federal University of Health Sciences of Porto Alegre - UFCSPA, Sarmento Leite, 245, Lab.714 , Porto Alegre , RS 90050-170 , Brazil
| | - Conrado Pedebos
- Centro de Biotecnologia , Universidade Federal do Rio Grande do Sul , Porto Alegre , RS 91509-900 , Brazil.,School of Chemistry , University of Southampton , Southampton , SO17 1BJ , U.K
| | - Marcelo D Polêto
- Centro de Biotecnologia , Universidade Federal do Rio Grande do Sul , Porto Alegre , RS 91509-900 , Brazil.,Departamento de Biologia Geral , Universidade Federal de Viçosa , Viçosa , MG 36570-000 , Brazil
| | | | - Hugo Verli
- Centro de Biotecnologia , Universidade Federal do Rio Grande do Sul , Porto Alegre , RS 91509-900 , Brazil
| |
Collapse
|
23
|
Kubincová A, Riniker S, Hünenberger PH. Reaction-field electrostatics in molecular dynamics simulations: development of a conservative scheme compatible with an atomic cutoff. Phys Chem Chem Phys 2020; 22:26419-26437. [DOI: 10.1039/d0cp03835k] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Shifting and switching schemes are developed to enable strict energy conservation in molecular dynamics simulations relying on reaction-field electrostatic (as well as Lennard-Jones) interactions with an atom-based cutoff truncation.
Collapse
Affiliation(s)
| | - Sereina Riniker
- Laboratory of Physical Chemistry
- ETH Zurich
- 8093 Zurich
- Switzerland
| | | |
Collapse
|
24
|
dos Santos TJ, Abreu CR, Horta BA, Tavares FW. Self-diffusion coefficients of methane/n-hexane mixtures at high pressures: An evaluation of the finite-size effect and a comparison of force fields. J Supercrit Fluids 2020. [DOI: 10.1016/j.supflu.2019.104639] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
25
|
Ramos MC, Horta VAC, Horta BAC. Molecular Dynamics Simulations of PAMAM and PPI Dendrimers Using the GROMOS-Compatible 2016H66 Force Field. J Chem Inf Model 2019; 59:1444-1457. [DOI: 10.1021/acs.jcim.8b00911] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Mayk C. Ramos
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, 21941-909, Brazil
| | - Vitor A. C. Horta
- Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Juiz de Fora, 36036-900, Brazil
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, 21941-909, Brazil
| |
Collapse
|