Dataset Information

Enhanced Computational Sampling of Perylene and Perylothiophene Packing with Rigid-Body Models.

ABSTRACT: Molecular simulations have the potential to advance the understanding of how the structure of organic materials can be engineered through the choice of chemical components but are limited by computational costs. The computational costs can be significantly lowered through the use of modeling approximations that capture the relevant features of a system, while lowering algorithmic complexity or by decreasing the degrees of freedom that must be integrated. Such methods include coarse-graining techniques, approximating long-range electrostatics with short-range potentials, and the use of rigid bodies to replace flexible bonded constraints between atoms. To understand whether and to what degree these techniques can be leveraged to enhance the understanding of planar organic molecules, we investigate the morphologies predicted by molecular dynamic simulations using simplified molecular models of perylene and perylothiophene. Approximately, 10 000 wall-clock hours of graphics processing unit-accelerated simulations are performed using both rigid and flexible models to test their efficiency and predictive capability with the two chemistries. We characterize the 1191 resulting morphologies using simulated X-ray diffraction and cluster analysis to distinguish structural transitions, summarized by four phase diagrams. We find that the morphologies generated by the rigid model of perylene and perylothiophene match with those generated by the flexible model. We find that ordered, hexagonally packed columnar phases are thermodynamically favored over a wide range of densities and temperatures for both molecules, in qualitative agreement with experiments. Furthermore, we find the rigid model to be more computationally efficient for both molecules, providing more samples per second and shorter times to equilibrium. Owing to the structural accuracy and improved computational efficiency of modeling polyaromatic groups as rigid bodies, we recommend this modeling choice for enhancing the sampling in polyaromatic molecular simulations.

SUBMITTER: Miller ED

PROVIDER: S-EPMC6640971 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:We present an alchemical enhanced sampling (ACES) method implemented in the GPU-accelerated AMBER free energy MD engine. The methods hinges on the creation of an "enhanced sampling state" by reducing or eliminating selected potential energy terms and interactions that lead to kinetic traps and conformational barriers while maintaining those terms that curtail the need to otherwise sample large volumes of phase space. For example, the enhanced sampling state might involve transforming regions of a ligand and/or protein side chain into a noninteracting "dummy state" with internal electrostatics and torsion angle terms turned off. The enhanced sampling state is connected to a real-state end point through a Hamiltonian replica exchange (HREMD) framework that is facilitated by newly developed alchemical transformation pathways and smoothstep softcore potentials. This creates a counterdiffusion of real and enhanced-sampling states along the HREMD network. The effect of a differential response of the environment to the real and enhanced-sampling states is minimized by leveraging the dual topology framework in AMBER to construct a counterbalancing HREMD network in the opposite alchemical direction with the same (or similar) real and enhanced sampling states at inverted end points. The method has been demonstrated in a series of test cases of increasing complexity where traditional MD, and in several cases alternative REST2-like enhanced sampling methods, are shown to fail. The hydration free energy for acetic acid was shown to be independent of the starting conformation, and the values for four additional edge case molecules from the FreeSolv database were shown to have a significantly closer agreement with experiment using ACES. The method was further able to handle different rotamer states in a Cdk2 ligand identified as fractionally occupied in crystal structures. Finally, ACES was applied to T4-lysozyme and demonstrated that the side chain distribution of V111χ1 could be reliably reproduced for the apo state, bound to p-xylene, and in p-xylene→ benzene transformations. In these cases, the ACES method is shown to be highly robust and superior to a REST2-like enhanced sampling implementation alone.

Dataset Information

Enhanced Computational Sampling of Perylene and Perylothiophene Packing with Rigid-Body Models.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets