Abstract
ABSTRACTThe calculation of Gibbs free energy of solvation ΔG is commonly applied to predict the solubility of solutes in different solvents. The transition path between the thermodynamic end states of full and non-existent solute–solvent interactions is described by the spacing of intermediate states and the soft-core potential employed to scale the interactions. The choice of both, the soft-core potential and the distribution of intermediate states, has a large impact on the simulated ΔG. In this work, we determine the free energy of hydration for 10 neutral amino acid side chain analogues to analyse the impact of different soft-core potentials on the consistency of results from two different free energy methods, i.e. thermodynamic integration and free energy perturbation. Additionally, we propose a new method to improve the alchemical pathway by adaptive spacing of intermediate states. Our technique is easy to implement and not limited to any simulation package or free energy technique in particular, making it easily adaptable to a variety of workflows. Exploiting the ability to systematically optimise alchemical pathways, we introduce a scheme to define a suitable number of windows based on prescribed configurational space overlap. With this, excessive sampling is avoided and computational costs can be reduced.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
