More Like this

1. Analysis of the conformational properties of amine ligands at the gold/water interface with QM, MM and QM/MM simulations

   https://doi.org/10.1039/C7CP06709G  

   Liang, Dongyue ; Hong, Jiewei ; Fang, Dong ; Bennett, Joseph W. ; Mason, Sara E. ; Hamers, Robert J. ; Cui, Qiang ( January 2018 , Physical Chemistry Chemical Physics)

   We describe a strategy of integrating quantum mechanical (QM), hybrid quantum mechanical/molecular mechanical (QM/MM) and MM simulations to analyze the physical properties of a solid/water interface. This protocol involves using a correlated ab initio (CCSD(T)) method to first calibrate Density Functional Theory (DFT) as the QM approach, which is then used in QM/MM simulations to compute relevant free energy quantities at the solid/water interface using a mean-field approximation of Yang et al. that decouples QM and MM thermal fluctuations; gas-phase QM/MM and periodic DFT calculations are used to determine the proper QM size in the QM/MM simulations. Finally, the QM/MM free energy results are compared with those obtained from MM simulations to directly calibrate the force field model for the solid/water interface. This protocol is illustrated by examining the orientations of an alkyl amine ligand at the gold/water interface, since the ligand conformation is expected to impact the chemical properties ( e.g. , charge) of the solid surface. DFT/MM and MM simulations using the INTERFACE force field lead to consistent results, suggesting that the effective gold/ligand interactions can be adequately described by a van der Waals model, while electrostatic and induction effects are largely quenched by solvation. The observed differences among periodic DFT, QM/MM and MM simulations, nevertheless, suggest that explicitly including electronic polarization and potentially charge transfer in the MM model can be important to the quantitative accuracy. The strategy of integrating multiple computational methods to cross-validate each other for complex interfaces is applicable to many problems that involve both inorganic/metallic and organic/biomolecular components, such as functionalized nanoparticles. 

   more » « less
2. Refining amino acid hydrophobicity for dynamics simulation of membrane proteins

   https://doi.org/10.7717/peerj.4230  

   Hills, Jr ( January 2018 , PeerJ)

   Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems. 

   more » « less
3. Coarse-Grained Force Field Calibration Based on Multi-Objective Bayesian Optimization to Simulate Water Diffusion in Poly-ɛ-caprolactone

   https://doi.org/10.1021/acs.jpca.0c01939  

   Sestito, Jesse M. ; Thatcher, Mary L ; Shu, Leshi ; Harris, Tequila A. ; Wang, Yan ( June 2020 , The Journal of Physical Chemistry A)

   Molecular dynamics at the atomistic scale is increasingly being used to predict material properties and speed up the materials design and development process. However, the accuracy of molecular dynamics predictions is sensitively dependent on the force fields. In the traditional force field calibration process, a specific property, predicted by the model, is compared with the experimental observation and the force field parameters are adjusted to minimize the difference. This leads to the issue that the calibrated force fields are not generic and robust enough to predict different properties. Here, a new calibration method based on multi-objective Bayesian optimization is developed to speed up the development of molecular dynamics force fields that are capable of predicting multiple properties accurately. This is achieved by reducing the number of simulation runs to generate the Pareto front with an efficient sequential sampling strategy. The methodology is demonstrated by generating a new coarse-grained force field for polycaprolactone, where the force field can predict mechanical properties and water diffusion in the polymer. 

   more » « less
4. A multiscale polymerization framework towards network structure and fracture of double-network hydrogels

   https://doi.org/10.1038/s41524-021-00509-5  

   Zhang, Mingzhen ; Zhang, Dong ; Chen, Hong ; Zhang, Yanxian ; Liu, Yonglan ; Ren, Baiping ; Zheng, Jie ( March 2021 , npj Computational Materials)

   Double-network (DN) hydrogels, consisting of two contrasting and interpenetrating polymer networks, are considered as perhaps the toughest soft-wet materials. Current knowledge of DN gels from synthesis methods to toughening mechanisms almost exclusively comes from chemically-linked DN hydrogels by experiments. Molecular modeling and simulations of inhomogeneous DN structure in hydrogels have proved to be extremely challenging. Herein, we developed a new multiscale simulation platform to computationally investigate the early fracture of physically-chemically linked agar/polyacrylamide (agar/PAM) DN hydrogels at a long timescale. A “random walk reactive polymerization” (RWRP) was developed to mimic a radical polymerization process, which enables to construct a physically-chemically linked agar/PAM DN hydrogel from monomers, while conventional and steered MD simulations were conducted to examine the structural-dependent energy dissipation and fracture behaviors at the relax and deformation states. Collective simulation results revealed that energy dissipation of agar/PAM hydrogels was attributed to a combination of the pulling out of agar chains from the DNs, the disruption of massive hydrogen bonds between and within DN structures, and the strong association of water molecules with both networks, thus explaining a different mechanical enhancement of agar/PAM hydrogels. This computational work provided atomic details of network structure, dynamics, solvation, and interactions of a hybrid DN hydrogel, and a different structural-dependent energy dissipation mode and fracture behavior of a hybrid DN hydrogel, which help to design tough hydrogels with new network structures and efficient energy dissipation modes. Additionally, the RWRP algorithm can be generally applied to construct the radical polymerization-produced hydrogels, elastomers, and polymers.

    

   more » « less
5. The UA?CG Workflow: High Performance Molecular Dynamics of Coarse-Grained Polymers

   https://doi.org/10.1109/PDP.2016.127  

   Ozog, David ; Malony, Allen D. ; Guenza, Marina ( February 2016 , Parallel, Distributed, and Network-Based Processing (DPD), 2016 24th Euromicro International Conference on)

   Our analytically based technique for coarse-graining (CG) polymer simulations dramatically improves spatial and temporal scaling while preserving thermodynamic quantities and bulk properties. The purpose of CG codes is to run more efficient molecular dynamics simulations, yet the research field generally lacks thorough analysis of how such codes scale with respect to full-atom representations. This paper conducts an in-depth performance study of highly realistic polymer melts on modern supercomputing systems. We also present a workflow that integrates our analytical solution for calculating CG forces with new high-performance techniques for mapping back and forth between the atomistic and CG descriptions in LAMMPS. The workflow benefits from the performance of CG, while maintaining full-atom accuracy. Our results show speedups up to 12x faster than atomistic simulations. 

   more » « less