Enhancing torsional sampling using fully adaptive simulated tempering.

Enhanced sampling algorithms are indispensable when working with highly disconnected multimodal distributions. An important application of these is the conformational exploration of particular internal degrees of freedom of molecular systems. However, despite the existence of many commonly used enhanced sampling algorithms to explore these internal motions, they often rely on system-dependent parameters, which negatively impact efficiency and reproducibility. Here, we present fully adaptive simulated tempering (FAST), a variation of the irreversible simulated tempering algorithm, which continuously optimizes the number, parameters, and weights of intermediate distributions to achieve maximally fast traversal over a space defined by the change in a predefined thermodynamic control variable such as temperature or an alchemical smoothing parameter. This work builds on a number of previously published methods, such as sequential Monte Carlo, and introduces a novel parameter optimization procedure that can, in principle, be used in any expanded ensemble algorithms. This method is validated by being applied on a number of different molecular systems with high torsional kinetic barriers. We also consider two different soft-core potentials during the interpolation procedure and compare their performance. We conclude that FAST is a highly efficient algorithm, which improves simulation reproducibility and can be successfully used in a variety of settings with the same initial hyperparameters.

Sire: An interoperability engine for prototyping algorithms and exchanging information between molecular simulation programs.

Sire is a Python/C++ library that is used both to prototype new algorithms and as an interoperability engine for exchanging information between molecular simulation programs. It provides a collection of file parsers and information converters that together make it easier to combine and leverage the functionality of many other programs and libraries. This empowers researchers to use sire to write a single script that can, for example, load a molecule from a PDBx/mmCIF file via Gemmi, perform SMARTS searches via RDKit, parameterize molecules using BioSimSpace, run GPU-accelerated molecular dynamics via OpenMM, and then display the resulting dynamics trajectory in a NGLView Jupyter notebook 3D molecular viewer. This functionality is built on by BioSimSpace, which uses sire's molecular information engine to interconvert with programs such as GROMACS, NAMD, Amber, and AmberTools for automated molecular parameterization and the running of molecular dynamics, metadynamics, and alchemical free energy workflows. Sire comes complete with a powerful molecular information search engine, plus trajectory loading and editing, analysis, and energy evaluation engines. This, when combined with an in-built computer algebra system, gives substantial flexibility to researchers to load, search for, edit, and combine molecular information from multiple sources and use that to drive novel algorithms by combining functionality from other programs. Sire is open source (GPL3) and is available via conda and at a free Jupyter notebook server at https://try.openbiosim.org. Sire is supported by the not-for-profit OpenBioSim community interest company.

BRIDGE: An Open Platform for Reproducible High-Throughput Free Energy Simulations.

Biomolecular Reaction and Interaction Dynamics Global Environment (BRIDGE) is an open-source web platform developed with the aim to provide an environment for the design of reliable methods to conduct reproducible biomolecular simulations. It is built on the well-known Galaxy bioinformatics platform. Through this, BRIDGE hosts computational chemistry tools on public web servers for internet use and provides machine- and operating-system-independent portability using the Docker container platform for local use. This construction improves the accessibility, shareability, and reproducibility of computational methods for molecular simulations. Here we present integrated free energy tools (or apps) to calculate absolute binding free energies (ABFEs) and relative binding free energies (RBFEs), as illustrated through use cases. We present free energy perturbation (FEP) methods contained in various software packages such as GROMACS and YANK that are independent of hardware configuration, software libraries, or operating systems when ported in the Galaxy-BRIDGE Docker container platform. By performing cyclin-dependent kinase 2 (CDK2) FEP calculations on geographically dispersed web servers (in Africa and Europe), we illustrate that large-scale computations can be performed using the exact same codes and methodology by collaborating groups through publicly shared protocols and workflows. The ease of public sharing and independent reproduction of simulations via BRIDGE makes possible an open review of methods and complete simulation protocols. This makes the discovery of inhibitors for drug targets accessible to nonexperts and the computer experiments that are used to arrive at leads verifiable by experts and reviewers. We illustrate this on ß-galactoside α-2,3-sialyltransferase I (ST3Gal-I), a breast cancer drug target, where a combination of RBFE and ABFE methods are used to compute the binding free energies of three inhibitors.

ProtoCaller: Robust Automation of Binding Free Energy Calculations.

ProtoCaller is a Python library distributed through Anaconda which automates relative protein-ligand binding free energy calculations in GROMACS. It links a number of popular specialized tools used to perform protein setup and parametrization, such as PDB2PQR, Modeller, and AmberTools. ProtoCaller supports commonly used AMBER force fields with additional cofactor parameters, and AM1-BCC is used to derive ligand charges. ProtoCaller also comes with an extensive PDB parser, an enhanced maximum common substructure algorithm providing improved ligand-ligand mapping, and a light GROMACS wrapper for running multiple molecular dynamics simulations. ProtoCaller is highly relevant to most researchers in the field of biomolecular simulation, allowing a customizable balance between automation and user intervention.

Enhancing Ligand and Protein Sampling Using Sequential Monte Carlo.

The sampling problem is one of the most widely studied topics in computational chemistry. While various methods exist for sampling along a set of reaction coordinates, many require system-dependent hyperparameters to achieve maximum efficiency. In this work, we present an alchemical variation of adaptive sequential Monte Carlo (SMC), an irreversible importance resampling method that is part of a well-studied class of methods that have been used in various applications but have been underexplored in computational biophysics. Afterward, we apply alchemical SMC on a variety of test cases, including torsional rotations of solvated ligands (butene and a terphenyl derivative), translational and rotational movements of protein-bound ligands, and protein side chain rotation coupled to the ligand degrees of freedom (T4-lysozyme, protein tyrosine phosphatase 1B, and transforming growth factor ß). We find that alchemical SMC is an efficient way to explore targeted degrees of freedom and can be applied to a variety of systems using the same hyperparameters to achieve a similar performance. Alchemical SMC is a promising tool for preparatory exploration of systems where long-timescale sampling of the entire system can be traded off against short-timescale sampling of a particular set of degrees of freedom over a population of conformers.

Sensitivity of Binding Free Energy Calculations to Initial Protein Crystal Structure.

Binding free energy calculations using alchemical free energy (AFE) methods are widely considered to be the most rigorous tool in the computational drug discovery arsenal. Despite this, the calculations suffer from accuracy, precision, and reproducibility issues. In this publication, we perform a high-throughput study of more than a thousand AFE calculations, utilizing over 220 µs of total sampling time, on three different protein systems to investigate the impact of the initial crystal structure on the resulting binding free energy values. We also consider the influence of equilibration time and discover that the initial crystal structure can have a significant effect on free energy values obtained at short timescales that can manifest itself as a free energy difference of more than 1 kcal/mol. At longer timescales, these differences are largely overtaken by important rare events, such as torsional ligand motions, typically resulting in a much higher uncertainty in the obtained values. This work emphasizes the importance of rare event sampling and long-timescale dynamics in free energy calculations even for routinely performed alchemical perturbations. We conclude that an optimal protocol should not only concentrate computational resources on achieving convergence in the alchemical coupling parameter (λ) space but also on longer simulations and multiple repeats.

Degradation Rate Observations as a Function of Drug Load in Solid-State Drug Products.

Degradation rates of solid-state drug products generally increase as the drug load decreases. A model for quantifying this effect based on surface area ratios is proposed here. This model relates the degradation rate to an estimate of the proportion of drug substance in contact with the excipient, and that the drug substance in contact with excipients degrades more quickly. Degradation data from previously published case studies and from 5 new case studies were found to be consistent with our proposed model; our model performed better than similar previously published models. It was also found that the relationship between degradation rate and drug load is largely independent of the temperature and humidity conditions, suggesting that drug load solely affects the pre-exponential factor of the Arrhenius equation and does not significantly affect the activation energy of the degradation process. A second method for calculating the proportion of the drug substance surface in contact with the excipient surface is presented in the Supplementary Material. Fundamentally, the 2 methods are very similar and provide almost identical fits to the experimental data.