Search papers, labs, and topics across Lattice.
This study introduces a novel stepwise construction approach to modeling supramolecular polymers using sequential replica exchange with solute tempering (REST), addressing the inefficiencies of traditional global REST methods. By incrementally adding monomers and applying REST to optimize their binding positions, the method significantly enhances computational efficiency and improves replica traversal during simulations. The results demonstrate that this approach can effectively generate polymer structures with characteristic hydrogen-bonding patterns at a reduced computational cost compared to conventional methods.
Sequential REST not only accelerates the modeling of supramolecular polymers but also enhances efficiency by optimizing monomer binding positions one at a time.
Predicting detailed atomistic structures of self-assembling systems remains a challenge for all-atom molecular dynamics simulations. Replica exchange with solute tempering (REST) has been used to study those systems by accelerating all monomers in a global and uniform manner. While such a global approach can in principle predict any morphology of the system, it has computational drawbacks such as inefficient replica traversal due to order-disorder transitions and the growing number of replicas with system size. To address these issues, here we propose an alternative, stepwise construction approach to modeling supramolecular polymers under the assumption of one-dimensional polymerization. Specifically, we generate polymer structures by adding new monomers one by one to the system and applying REST to the new monomers to find their optimal binding positions based on an energy-based scoring function. The monomer addition and enhanced sampling are repeated sequentially until a polymer of desired length is obtained. We test the above procedure using a model supramolecular polymer in explicit solvent, and show that it can generate a polymer structure with characteristic H-bonding patterns at reduced computational costs, while also improving the efficiency of replica traversal significantly. We thus expect that the sequential REST will be useful for modeling supramolecular polymers, particularly for cases where global REST simulations are too demanding computationally.
