These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
Preprints are manuscripts made publicly available before they have been submitted for formal peer review and publication. They might contain new research findings or data. Preprints can be a draft or final version of an author's research but must not have been accepted for publication at the time of submission.
Recent improvements to free energy perturbation (FEP) calculations, especiallyFEP+, established their utility for pharmaceutical lead optimization. However, to dateFEP has typically been helpful only when (1) high-quality X-ray data is available and(2) the target protein does not undergo significant conformational changes. Also, alack of systematic studies on determining an adequate sampling time is often one ofthe primary limitations of FEP calculations. Herein, we propose a modified versionof the FEP/REST (i.e., replica exchange with solute tempering) sampling protocol,based on systematic studies on several targets by probing a large number of permutations with different sampling schemes. Improved FEP+ binding affinity predictions for regular flexible-loop (F-loop) motions and considerable structural changes can be obtained by extending the pre-REST sampling time from 0.24 ns to 5 ns/λand 2×10 ns/λ, respectively. We obtained much more precise ∆∆G calculations of the individual perturbations, including the sign of the transformations and less error. We extended the REST simulations from 5 ns to 8 ns to achieve reasonable free energy convergence.Implementing REST to the entire ligand as opposed to solely the perturbed region, and also some important flexible protein residues (pREST region) in ligand binding domain (LBD) , also considerably improved the FEP+ results in most of the studied cases. Preliminary molecular dynamics (MD) runs were useful for establishing the correct binding mode of the compounds and thus precise alignment for FEP+.