Automated Preparation of Nanoscopic Structures: Graph-Based Sequence Analysis, Mismatch Detection, and pH-Consistent Protonation with Uncertainty Estimates

TL;DR

The paper introduces ASAP, a fast, robust, and modular protocol for automated analysis, pH-consistent protonation, and uncertainty estimation of nanoscopic structures, enhancing molecular simulation accuracy and efficiency.

Contribution

ASAP provides a novel, graph-based, automated pipeline for structure analysis and pH prediction with uncertainty estimates, applicable to biomolecules and other nanostructures.

Findings

Effective error assessment of input structures.

Accurate pKa prediction with uncertainty quantification.

Facilitates fast pH-state switching during simulations.

Abstract

Structure and function in nanoscale atomistic assemblies are tightly coupled, and every atom with its specific position and even every electron will have a decisive effect on the electronic structure, and hence, on the molecular properties. Molecular simulations of nanoscopic atomistic structures therefore require accurately resolved three-dimensional input structures. If extracted from experiment, these structures often suffer from severe uncertainties, of which the lack of information on hydrogen atoms is a prominent example. Hence, experimental structures require careful review and curation, which is a time-consuming and error-prone process. Here, we present a fast and robust protocol for the automated structure analysis, and pH-consistent protonation, in short, ASAP. For biomolecules as a target, the ASAP protocol integrates sequence analysis and error assessment of a given input structure. ASAP allows for pKa prediction from reference data through Gaussian process regression including uncertainty estimation and connects to system-focused atomistic modeling described in (J. Chem. Theory Comput. 16, 2020, 1646). Although focused on biomolecules, ASAP can be extended to other nanoscopic objects, because most of its design elements rely on a general graph-based foundation guaranteeing transferability. The modular character of the underlying pipeline supports different degrees of automation, which allows for (i) efficient feedback loops for human-machine interaction with a low entrance barrier and for (ii) integration into autonomous procedures such as automated force field parametrizations. This facilitates fast switching of the pH-state through on-the-fly system-focused reparametrization during a molecular simulation at virtually no extra computational cost.