Abstract:Precise regulation of the peptide self-assembly into ordered nanostructures with intriguing properties has attracted intense attention. However, predicting peptide assembly at atomic resolution is a challenge due to both the structural flexibility of peptides and the associated huge computational costs. A machine learning-guided adaptive parametrization method was proposed for developing a mixed atomic and coarse-grained (CG) model through a multiobjective optimization strategy. Our model incorporates the united-atom (UA) model for diphenylalanine (P) and the polarizable electrostatic-variable coarse-grained (VaCG) model for aqueous ionic liquid [BMIM]+[BF4]- solution. In this mixed model, the coupling van der Waals (vdW) interaction is addressed by introducing virtual sites (VS) in the UA model to interact with solvent CG beads. The coupling parameters, including the electrostatic parameter and vdW parameters, are automatically optimized through ML-guided adaptive parametrization. The performance of this model was tested by some microstructural properties, e.g., the average number of P-P intermolecular hydrogen bonds (HBs) and radius distribution functions (RDFs) between P and different fragments of IL, in comparison with all-atom (AA) simulations. The computational cost is significantly reduced using such a parametrization scheme, which could search tens of thousands of force-field parameter sets, while needing only a small fraction of them to be assessed with molecular dynamics (MD) simulations. We used such a mixed resolution model to investigate the self-assembly in IL-water mixtures with variants of IL concentration (X). The long-range-ordered fibril structure is formed in a pure water system (X = 0). With an increase of IL concentrations, the formation of an ordered self-assembly nanostructure is prohibited, instead forming branched fibril at X = 2 mol % or amorphous aggregates when X > 10 mol %, resulting from the interplay between π-stacking and HB interactions between P and IL. The qualitative agreement between the simulated structures and the observed morphologies in experiments indicates the applicability of ML-guided parametrization strategy in the study of complex systems, such as polymers, lipid bilayers, and polysaccharides.

Parameterization of Pace Force Field for Membrane Environment and Simulation of Helical Peptides and Helix-Helix Association

Pace Force Field for Protein Simulations. 2. Folding Simulations of Peptides

Pace Force Field for Protein Simulations. 1. Full Parameterization of Version 1 and Verification

Further optimization of a hybrid united-atom and coarse-grained force field for folding simulations: Improved backbone hydration and interactions between charged side chains.

Development of a Hybrid-Resolution Force Field for Peptide Self-Assembly Simulations: Optimizing Peptide-Peptide and Peptide-Solvent Interactions

A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

Residue-Speci fi c Force Field Based on the Protein Coil Library

CHARMM-GUI PACE CG Builder for Solution, Micelle, and Bilayer Coarse-Grained Simulations

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Machine Learning-Guided Adaptive Parametrization for Coupling Terms in a Mixed United-Atom/Coarse-Grained Model for Diphenylalanine Self-Assembly in Aqueous Ionic Liquids

Folding and insertion thermodynamics of the transmembrane WALP peptide

Universal Implementation of a Residue-Specific Force Field Based on CMAP Potentials and Free Energy Decomposition

Benchmarking polarizable and non-polarizable force fields for Ca2+-peptides against a comprehensive QM dataset

A Test of AMBER Force Fields in Predicting the Secondary Structure of Α-Helical and Β-Hairpin Peptides

Hybrid QM/MM simulation of the hydration phenomena of dipalmitoylphosphatidylcholine headgroup.

A Top‐down and Bottom‐up Combined Strategy for Parameterization of Coarse‐grained Force Fields for Phospholipids

Phosphorylation Modification Force Field FB18CMAP Improving Conformation Sampling of Phosphoproteins

Balancing Force Field Protein–Lipid Interactions To Capture Transmembrane Helix–Helix Association