RESUMEN
Atomic partial charges are crucial parameters in molecular dynamics simulation, dictating the electrostatic contributions to intermolecular energies and thereby the potential energy landscape. Traditionally, the assignment of partial charges has relied on surrogates of ab initio semiempirical quantum chemical methods such as AM1-BCC and is expensive for large systems or large numbers of molecules. We propose a hybrid physical/graph neural network-based approximation to the widely popular AM1-BCC charge model that is orders of magnitude faster while maintaining accuracy comparable to differences in AM1-BCC implementations. Our hybrid approach couples a graph neural network to a streamlined charge equilibration approach in order to predict molecule-specific atomic electronegativity and hardness parameters, followed by analytical determination of optimal charge-equilibrated parameters that preserve total molecular charge. This hybrid approach scales linearly with the number of atoms, enabling for the first time the use of fully consistent charge models for small molecules and biopolymers for the construction of next-generation self-consistent biomolecular force fields. Implemented in the free and open source package EspalomaCharge, this approach provides drop-in replacements for both AmberTools antechamber and the Open Force Field Toolkit charging workflows, in addition to stand-alone charge generation interfaces. Source code is available at https://github.com/choderalab/espaloma-charge.
RESUMEN
A methodology that combines alchemical free energy calculations (FEP) with machine learning (ML) has been developed to compute accurate absolute hydration free energies. The hybrid FEP/ML methodology was trained on a subset of the FreeSolv database and retrospectively shown to outperform most submissions from the SAMPL4 competition. Compared to pure machine-learning approaches, FEP/ML yields more precise estimates of free energies of hydration and requires a fraction of the training set size to outperform standalone FEP calculations. The ML-derived correction terms are further shown to be transferable to a range of related FEP simulation protocols. The approach may be used to inexpensively improve the accuracy of FEP calculations and to flag molecules which will benefit the most from bespoke force field parametrization efforts.
Asunto(s)
Aprendizaje Automático , Simulación por Computador , Entropía , Estudios Retrospectivos , TermodinámicaRESUMEN
The Zika virus (ZIKV), discovered in Africa in 1947, swiftly spread across continents, causing significant concern due to its recent association with microcephaly in newborns and Guillain-Barré syndrome in adults. Despite a decrease in prevalence, the potential for a resurgence remains, necessitating urgent therapeutic interventions. Like other flaviviruses, ZIKV presents promising drug targets within its replication machinery, notably the NS3 helicase (NS3Hel) protein, which plays critical roles in viral replication. However, a lack of structural information impedes the development of specific inhibitors targeting NS3Hel. Here we applied high-throughput crystallographic fragment screening on ZIKV NS3Hel, which yielded structures that reveal 3D binding poses of 46 fragments at multiple sites of the protein, including 11 unique fragments in the RNA-cleft site. These fragment structures provide templates for direct design of hit compounds and should thus assist the development of novel direct-acting antivirals against ZIKV and related flaviviruses, thus opening a promising avenue for combating future outbreaks.
RESUMEN
Alchemical free energy calculations are a useful tool for predicting free energy differences associated with the transfer of molecules from one environment to another. The hallmark of these methods is the use of "bridging" potential energy functions representing alchemical intermediate states that cannot exist as real chemical species. The data collected from these bridging alchemical thermodynamic states allows the efficient computation of transfer free energies (or differences in transfer free energies) with orders of magnitude less simulation time than simulating the transfer process directly. While these methods are highly flexible, care must be taken in avoiding common pitfalls to ensure that computed free energy differences can be robust and reproducible for the chosen force field, and that appropriate corrections are included to permit direct comparison with experimental data. In this paper, we review current best practices for several popular application domains of alchemical free energy calculations performed with equilibrium simulations, in particular relative and absolute small molecule binding free energy calculations to biomolecular targets.