Concurrent and Adaptive Extreme Scale Binding Free Energy Calculations

The efficacy of drug treatments depends on how tightly small molecules bind to their target proteins. The rapid and accurate quantification of the strength of these interactions (as measured by binding affinity) is a grand challenge of computational chemistry, surmounting which could revolutionize drug design and provide the platform for patient-specific medicine. Recent evidence suggests that molecular dynamics (MD) can achieve useful predictive accuracy (< 1 kcal/mol). For this predictive accuracy to impact clinical decision making, binding free energy computational campaigns must provide results rapidly and without loss of accuracy. This demands advances in algorithms, scalable software systems, and efficient utilization of supercomputing resources. We introduce a framework called HTBAC, designed to support accurate and scalable drug binding affinity calculations, while marshaling large simulation campaigns. We show that HTBAC supports the specification and execution of free-energy protocols at scale. This paper makes three main contributions: (1) shows the importance of adaptive execution for ensemble-based free energy protocols to improve binding affinity accuracy; (2) presents and characterizes HTBAC -- a software system that enables the scalable and adaptive execution of binding affinity protocols at scale; and (3) for a widely used free-energy protocol (TIES), shows improvements in the accuracy of simulations for a fixed amount of resource, or reduced resource consumption for a fixed accuracy as a consequence of adaptive execution

Comparison of Equilibrium and Nonequilibrium Approaches for Relative Binding Free Energy Predictions

Alchemical relative binding free energy calculations have recently found important applications in drug optimization. A series of congeneric compounds are generated from a preidentified lead compound, and their relative binding affinities to a protein are assessed in order to optimize candidate drugs. While methods based on equilibrium thermodynamics have been extensively studied, an approach based on nonequilibrium methods has recently been reported together with claims of its superiority. However, these claims pay insufficient attention to the basis and reliability of both methods. Here we report a comparative study of the two approaches across a large data set, comprising more than 500 ligand transformations spanning in excess of 300 ligands binding to a set of 14 diverse protein targets. Ensemble methods are essential to quantify the uncertainty in these calculations, not only for the reasons already established in the equilibrium approach but also to ensure that the nonequilibrium calculations reside within their domain of validity. If and only if ensemble methods are applied, we find that the nonequilibrium method can achieve accuracy and precision comparable to those of the equilibrium approach. Compared to the equilibrium method, the nonequilibrium approach can reduce computational costs but introduces higher computational complexity and longer wall clock times. There are, however, cases where the standard length of a nonequilibrium transition is not sufficient, necessitating a complete rerun of the entire set of transitions. This significantly increases the computational cost and proves to be highly inconvenient during large-scale applications. Our findings provide a key set of recommendations that should be adopted for the reliable implementation of nonequilibrium approaches to relative binding free energy calculations in ligand-protein systems

Development and analysis of Tinker-OpenMM as a GPU-based free energy perturbation engine

The utilization of computational technologies for the lead optimization process is one of the biggest challenges in the computational chemistry field. In this dissertation, I describe the addition of GPU-based absolute and relative free energy calculation methods using polarizable force field AMOEBA to Tinker-OpenMM. I then proceed to test the capabilities of this platform by studying the binding free energy and binding structures of derivatives of the MELK inhibitor IN17. Also, I present the implementation of virial-based pressure control to the Tinker-OpenMM platform that is needed for performing isobaric simulations.Cellular and Molecular Biolog

Effect of Oxidative Damage on the Stability and Dimerization of Superoxide Dismutase 1

During their life cycle, proteins are subject to different modifications involving reactive oxygen species. Such oxidative damage to proteins may lead to the formation of insoluble aggregates and cytotoxicity and is associated with age-related disorders including neurodegenerative diseases, cancer, and diabetes. Superoxide dismutase 1 (SOD1), a key antioxidant enzyme in human cells, is particularly susceptible to such modifications. Moreover, this homodimeric metalloenzyme has been directly linked to both familial and sporadic amyotrophic lateral sclerosis (ALS), a devastating, late-onset motor neuronal disease, with more than 150 ALS-related mutations in the SOD1 gene. Importantly, oxidatively damaged SOD1 aggregates have been observed in both familial and sporadic forms of the disease. However, the molecular mechanisms as well as potential implications of oxidative stress in SOD1-induced cytotoxicity remain elusive. In this study, we examine the effects of oxidative modification on SOD1 monomer and homodimer stability, the key molecular properties related to SOD1 aggregation. We use molecular dynamics simulations in combination with thermodynamic integration to study microscopic-level site-specific effects of oxidative "mutations" at the dimer interface, including lysine, arginine, proline and threonine carbonylation, and cysteine oxidation. Our results show that oxidative damage of even single residues at the interface may drastically destabilize the SOD1 homodimer, with several modifications exhibiting a comparable effect to that of the most drastic ALS-causing mutations known. Additionally, we show that the SOD1 monomer stability decreases upon oxidative stress, which may lead to partial local unfolding and consequently to increased aggregation propensity. Importantly, these results suggest that oxidative stress may play a key role in development of ALS, with the mutations in the SOD1 gene being an additional factor

Automated free-energy calculation from atomistic simulations

We devise automated workflows for the calculation of Helmholtz and Gibbs free energies and their temperature and pressure dependence and provide the corresponding computational tools. We employ nonequilibrium thermodynamics for evaluating the free energy of solid and liquid phases at a given temperature and reversible scaling for computing free energies over a wide range of temperatures, including the direct integration of P−T coexistence lines. By changing the chemistry and the interatomic potential, alchemical and upscaling free energy calculations are possible. Several examples illustrate the accuracy and efficiency of our implementation

TIES 20: Relative Binding Free Energy with a Flexible Superimposition Algorithm and Partial Ring Morphing

The TIES (Thermodynamic Integration with Enhanced Sampling) protocol is a formally exact alchemical approach in computational chemistry to the calculation of relative binding free energies. The validity of TIES relies on the correctness of matching atoms across compared pairs of ligands, laying the foundation for the transformation along an alchemical pathway. We implement a flexible topology superimposition algorithm which uses an exhaustive joint-traversal for computing the largest common component(s). The algorithm is employed to enable matching and morphing of partial rings in the TIES protocol along with a validation study using 55 transformations and five different proteins from our previous work. We find that TIES 20 with the RESP charge system, using the new superimposition algorithm, reproduces the previous results with mean unsigned error of 0.75 kcal/mol with respect to the experimental data. Enabling the morphing of partial rings decreases the size of the alchemical region in the dual-topology transformations resulting in a significant improvement in the prediction precision. We find that increasing the ensemble size from 5 to 20 replicas per λ window only has a minimal impact on the accuracy. However, the non-normal nature of the relative free energy distributions underscores the importance of ensemble simulation. We further compare the results with the AM1-BCC charge system and show that it improves agreement with the experimental data by slightly over 10%. This improvement is partly due to AM1-BCC affecting only the charges of the atoms local to the mutation, which translates to even fewer morphed atoms, consequently reducing issues with sampling and therefore ensemble averaging. TIES 20, in conjunction with the enablement of ring morphing, reduces the size of the alchemical region and significantly improves the precision of the predicted free energies