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

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.