Redox potential replica exchange molecular dynamics at constant pH in AMBER: Implementation and validation

Abstract: Redox processes are important in chemistry, with applications in biomedicine, chemical analysis, among others. As many redox experiments are also performed at a fixed value of pH, having an efficient computational method to support experimental measures at both constant redox potential and pH is very important. Such computational techniques have the potential to validate experimental observations performed under these conditions and to provide additional information unachievable experimentally such as an atomic level description of macroscopic measures. We present the implementation of discrete redox and protonation states methods for constant redox potential Molecular Dynamics (CEMD), for coupled constant pH and constant redox potential MD (C(pH,E)MD), and for Replica Exchange MD along the redox potential dimension (E-REMD) on the AMBER software package. Validation results are presented for a small system that contains a single heme group: N-acetylmicroperoxidase-8 (NAcMP8) axially connected to a histidine peptide. The methods implemented allow one to make standard redox potential (Eo) predictions with the same easiness and accuracy as pKa predictions using the constant pH molecular dynamics and pH-REMD methods currently available on AMBER. In our simulations, we can correctly describe, in agreement also with theoretical predictions, the following behaviors: when a redox-active group is reduced, the pKa of a near pH-active group increases because it becomes easier for a proton to be attached; equivalently, when a pH-active group is protonated, the Eo of an adjacent redox-active group rises. Furthermore, our results also show that E-REMD is able to achieve faster statistical convergence than CEMD or C(pH,E)MD. Moreover, computational benchmarks using our methodologies show high-performance of GPU accelerated calculations in comparison to conventional CPU calculations.