Amber free energy tools: Interoperable software for free energy simulations using generalized quantum mechanical/molecular mechanical and machine learning potentials.

We report the development and testing of new integrated cyberinfrastructure for performing free energy simulations with generalized hybrid quantum mechanical/molecular mechanical (QM/MM) and machine learning potentials (MLPs) in Amber. The Sander molecular dynamics program has been extended to leverage fast, density-functional tight-binding models implemented in the DFTB+ and xTB packages, and an interface to the DeePMD-kit software enables the use of MLPs. The software is integrated through application program interfaces that circumvent the need to perform "system calls" and enable the incorporation of long-range Ewald electrostatics into the external software's self-consistent field procedure. The infrastructure provides access to QM/MM models that may serve as the foundation for QM/MM-ΔMLP potentials, which supplement the semiempirical QM/MM model with a MLP correction trained to reproduce ab initio QM/MM energies and forces. Efficient optimization of minimum free energy pathways is enabled through a new surface-accelerated finite-temperature string method implemented in the FE-ToolKit package. Furthermore, we interfaced Sander with the i-PI software by implementing the socket communication protocol used in the i-PI client-server model. The new interface with i-PI allows for the treatment of nuclear quantum effects with semiempirical QM/MM-ΔMLP models. The modular interoperable software is demonstrated on proton transfer reactions in guanine-thymine mispairs in a B-form deoxyribonucleic acid helix. The current work represents a considerable advance in the development of modular software for performing free energy simulations of chemical reactions that are important in a wide range of applications.

Software Infrastructure for Next-Generation QM/MM-ΔMLP Force Fields.

We present software infrastructure for the design and testing of new quantum mechanical/molecular mechanical and machine-learning potential (QM/MM-ΔMLP) force fields for a wide range of applications. The software integrates Amber's molecular dynamics simulation capabilities with fast, approximate quantum models in the xtb package and machine-learning potential corrections in DeePMD-kit. The xtb package implements the recently developed density-functional tight-binding QM models with multipolar electrostatics and density-dependent dispersion (GFN2-xTB), and the interface with Amber enables their use in periodic boundary QM/MM simulations with linear-scaling QM/MM particle-mesh Ewald electrostatics. The accuracy of the semiempirical models is enhanced by including machine-learning correction potentials (ΔMLPs) enabled through an interface with the DeePMD-kit software. The goal of this paper is to present and validate the implementation of this software infrastructure in molecular dynamics and free energy simulations. The utility of the new infrastructure is demonstrated in proof-of-concept example applications. The software elements presented here are open source and freely available. Their interface provides a powerful enabling technology for the design of new QM/MM-ΔMLP models for studying a wide range of problems, including biomolecular reactivity and protein-ligand binding.

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials.

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

Surface-Accelerated String Method for Locating Minimum Free Energy Paths.

We present a surface-accelerated string method (SASM) to efficiently optimize low-dimensional reaction pathways from the sampling performed with expensive quantum mechanical/molecular mechanical (QM/MM) Hamiltonians. The SASM accelerates the convergence of the path using the aggregate sampling obtained from the current and previous string iterations, whereas approaches like the string method in collective variables (SMCV) or the modified string method in collective variables (MSMCV) update the path only from the sampling obtained from the current iteration. Furthermore, the SASM decouples the number of images used to perform sampling from the number of synthetic images used to represent the path. The path is optimized on the current best estimate of the free energy surface obtained from all available sampling, and the proposed set of new simulations is not restricted to being located along the optimized path. Instead, the umbrella potential placement is chosen to extend the range of the free energy surface and improve the quality of the free energy estimates near the path. In this manner, the SASM is shown to improve the exploration for a minimum free energy pathway in regions where the free energy surface is relatively flat. Furthermore, it improves the quality of the free energy profile when the string is discretized with too few images. We compare the SASM, SMCV, and MSMCV using 3 QM/MM applications: a ribozyme methyltransferase reaction using 2 reaction coordinates, the 2'-O-transphosphorylation reaction of Hammerhead ribozyme using 3 reaction coordinates, and a tautomeric reaction in B-DNA using 5 reaction coordinates. We show that SASM converges the paths using roughly 3 times less sampling than the SMCV and MSMCV methods. All three algorithms have been implemented in the FE-ToolKit package made freely available.

Improving rice nitrogen-use efficiency by modulating a novel monouniquitination machinery for optimal root plasticity response to nitrogen.

Plant nitrogen (N)-use efficiency (NUE) is largely determined by the ability of root to take up external N sources, whose availability and distribution in turn trigger the modification of root system architecture (RSA) for N foraging. Therefore, improving N-responsive reshaping of RSA for optimal N absorption is a major target for developing crops with high NUE. In this study, we identified RNR10 (REGULATOR OF N-RESPONSIVE RSA ON CHROMOSOME 10) as the causal gene that underlies the significantly different root developmental plasticity in response to changes in N level exhibited by the indica (Xian) and japonica (Geng) subspecies of rice. RNR10 encodes an F-box protein that interacts with a negative regulator of auxin biosynthesis, DNR1 (DULL NITROGEN RESPONSE1). Interestingly, RNR10 monoubiquitinates DNR1 and inhibits its degradation, thus antagonizing auxin accumulation, which results in reduced root responsivity to N and nitrate (NO3-) uptake. Therefore, modulating the RNR10-DNR1-auxin module provides a novel strategy for coordinating a desirable RSA and enhanced N acquisition for future sustainable agriculture.

Modern semiempirical electronic structure methods and machine learning potentials for drug discovery: Conformers, tautomers, and protonation states.

Modern semiempirical electronic structure methods have considerable promise in drug discovery as universal "force fields" that can reliably model biological and drug-like molecules, including alternative tautomers and protonation states. Herein, we compare the performance of several neglect of diatomic differential overlap-based semiempirical (MNDO/d, AM1, PM6, PM6-D3H4X, PM7, and ODM2), density-functional tight-binding based (DFTB3, DFTB/ChIMES, GFN1-xTB, and GFN2-xTB) models with pure machine learning potentials (ANI-1x and ANI-2x) and hybrid quantum mechanical/machine learning potentials (AIQM1 and QDπ) for a wide range of data computed at a consistent ωB97X/6-31G* level of theory (as in the ANI-1x database). This data includes conformational energies, intermolecular interactions, tautomers, and protonation states. Additional comparisons are made to a set of natural and synthetic nucleic acids from the artificially expanded genetic information system that has important implications for the design of new biotechnology and therapeutics. Finally, we examine the acid/base chemistry relevant for RNA cleavage reactions catalyzed by small nucleolytic ribozymes, DNAzymes, and ribonucleases. Overall, the hybrid quantum mechanical/machine learning potentials appear to be the most robust for these datasets, and the recently developed QDπ model performs exceptionally well, having especially high accuracy for tautomers and protonation states relevant to drug discovery.

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery.

We report QDπ-v1.0 for modeling the internal energy of drug molecules containing H, C, N, and O atoms. The QDπ model is in the form of a quantum mechanical/machine learning potential correction (QM/Δ-MLP) that uses a fast third-order self-consistent density-functional tight-binding (DFTB3/3OB) model that is corrected to a quantitatively high-level of accuracy through a deep-learning potential (DeepPot-SE). The model has the advantage that it is able to properly treat electrostatic interactions and handle changes in charge/protonation states. The model is trained against reference data computed at the ωB97X/6-31G* level (as in the ANI-1x data set) and compared to several other approximate semiempirical and machine learning potentials (ANI-1x, ANI-2x, DFTB3, MNDO/d, AM1, PM6, GFN1-xTB, and GFN2-xTB). The QDπ model is demonstrated to be accurate for a wide range of intra- and intermolecular interactions (despite its intended use as an internal energy model) and has shown to perform exceptionally well for relative protonation/deprotonation energies and tautomers. An example application to model reactions involved in RNA strand cleavage catalyzed by protein and nucleic acid enzymes illustrates QDπ has average errors less than 0.5 kcal/mol, whereas the other models compared have errors over an order of magnitude greater. Taken together, this makes QDπ highly attractive as a potential force field model for drug discovery.

Growth-Promoting Effects of Zhenqi Granules on Finishing Pigs.

Developing nonantibiotic livestock growth promoters attracts intensive interest in the post-antibiotic era. In this study, we investigated the growth-promoting efficacy of Zhenqi granules (ZQ) in pigs and further explored the possible mechanisms by transcriptomics analysis. Weaned piglets (52 days old with an average body weight of 17.92 kg) were fed with diets supplemented with different doses of ZQ (0 g/kg, 1 g/kg, and 2 g/kg) for 30 days and continued observations for an additional 32 days after removing ZQ from the diets. Compared with the control group, the average daily gain, carcass weight, average back fat thickness, and fat meat percentage of the group supplemented with 1 g/kg of ZQ showed a significant increase, and the feed/gain ratio was lower. The group supplemented with 2 g/kg of ZQ also showed a significant increase in average daily gain and average backfat thickness. A transcriptomics analysis revealed that the supplementation of ZQ at 1 g/kg upregulated the expression of genes related to collagen biosynthesis and lipid biosynthesis in skeletal muscle and liver. This effect was primarily through upregulating the mRNA levels of structural proteins and lipid-related enzymes. This study demonstrates the growth-promoting efficacy of ZQ and provides some insights of the mechanism of growth promotion.

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery.

Predicting protein-ligand binding affinities and the associated thermodynamics of biomolecular recognition is a primary objective of structure-based drug design. Alchemical free energy simulations offer a highly accurate and computationally efficient route to achieving this goal. While the AMBER molecular dynamics package has successfully been used for alchemical free energy simulations in academic research groups for decades, widespread impact in industrial drug discovery settings has been minimal because of the previous limitations within the AMBER alchemical code, coupled with challenges in system setup and postprocessing workflows. Through a close academia-industry collaboration we have addressed many of the previous limitations with an aim to improve accuracy, efficiency, and robustness of alchemical binding free energy simulations in industrial drug discovery applications. Here, we highlight some of the recent advances in AMBER20 with a focus on alchemical binding free energy (BFE) calculations, which are less computationally intensive than alternative binding free energy methods where full binding/unbinding paths are explored. In addition to scientific and technical advances in AMBER20, we also describe the essential practical aspects associated with running relative alchemical BFE calculations, along with recommendations for best practices, highlighting the importance not only of the alchemical simulation code but also the auxiliary functionalities and expertise required to obtain accurate and reliable results. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative BFE simulations in AMBER20, with a focus on real-world drug discovery applications.

Improved Alchemical Free Energy Calculations with Optimized Smoothstep Softcore Potentials.

Progress in the development of GPU-accelerated free energy simulation software has enabled practical applications on complex biological systems and fueled efforts to develop more accurate and robust predictive methods. In particular, this work re-examines concerted (a.k.a., one-step or unified) alchemical transformations commonly used in the prediction of hydration and relative binding free energies (RBFEs). We first classify several known challenges in these calculations into three categories: endpoint catastrophes, particle collapse, and large gradient-jumps. While endpoint catastrophes have long been addressed using softcore potentials, the remaining two problems occur much more sporadically and can result in either numerical instability (i.e., complete failure of a simulation) or inconsistent estimation (i.e., stochastic convergence to an incorrect result). The particle collapse problem stems from an imbalance in short-range electrostatic and repulsive interactions and can, in principle, be solved by appropriately balancing the respective softcore parameters. However, the large gradient-jump problem itself arises from the sensitivity of the free energy to large values of the softcore parameters, as might be used in trying to solve the particle collapse issue. Often, no satisfactory compromise exists with the existing softcore potential form. As a framework for solving these problems, we developed a new family of smoothstep softcore (SSC) potentials motivated by an analysis of the derivatives along the alchemical path. The smoothstep polynomials generalize the monomial functions that are used in most implementations and provide an additional path-dependent smoothing parameter. The effectiveness of this approach is demonstrated on simple yet pathological cases that illustrate the three problems outlined. With appropriate parameter selection, we find that a second-order SSC(2) potential does at least as well as the conventional approach and provides vast improvement in terms of consistency across all cases. Last, we compare the concerted SSC(2) approach against the gold-standard stepwise (a.k.a., decoupled or multistep) scheme over a large set of RBFE calculations as might be encountered in drug discovery.

Validation of Free Energy Methods in AMBER.

Herein we provide high-precision validation tests of the latest GPU-accelerated free energy code in AMBER. We demonstrate that consistent free energy results are obtained in both the gas phase and in solution. We first show, in the context of thermodynamic integration (TI), that the results are invariant with respect to "split" (e.g., stepwise decharge-vdW-recharge) versus "unified" protocols. This brought to light a subtle inconsistency in previous versions of AMBER that was traced to the improper treatment of 1-4 vdW and electrostatic interactions involving atoms across the softcore boundary. We illustrate that under the assumption that the ensembles produced by different legs of the alchemical transformation between molecules A and B in the gas phase and aqueous phase are very small, the inconsistency in the relative hydration free energy ΔΔGhydr[A → B] = ΔGaq[A → B] - ΔGgas[A → B] is minimal. However, for general cases where the ensembles are shown to be substantially different, as expected in ligand-protein binding applications, these errors can be large. Finally, we demonstrate that results for relative hydration free energy simulations are independent of TI or multistate Bennett's acceptance ratio (MBAR) analysis, invariant to the specific choice of the softcore region, and in agreement with results derived from absolute hydration free energy values.

Effective flocculation of target microalgae with self-flocculating microalgae induced by pH decrease.

A flocculation method was developed to harvest target microalgae with self-flocculating microalgae induced by decreasing pH to just below isoelectric point. The flocculation efficiencies of target microalgae were much higher than those flocculated only via pH decrease. The mechanism could be that negatively charged self-flocculating microalgal cells became positively charged during pH decrease, subsequently attracted negatively charged target microalgae cells to form flocs and settled down due to gravity. Microalgal biomass concentration and released polysaccharide (RPS) from target microalgae influenced flocculation efficiencies, while multivalent metal ions in growth medium could not. Furthermore, neutralizing pH and then supplementing nutrients allowed flocculated medium to be recycled for cultivation. Finally, Spearman's Rank Correlation Coefficients (Rs) between flocculation efficiency and key factors were also investigated. These results suggest that this method is effective, simple to operate and allows the reuse of flocculated medium, thereby contributing to the economic production from microalgae to biodiesel.

Freshwater microalgae harvested via flocculation induced by pH decrease.

BACKGROUND: Recent studies have demonstrated that microalga has been widely regarded as one of the most promising raw materials of biofuels. However, lack of an economical, efficient and convenient method to harvest microalgae is a bottleneck to boost their full-scale application. Many methods of harvesting microalgae, including mechanical, electrical, biological and chemical based, have been studied to overcome this hurdle. RESULTS: A new flocculation method induced by decreasing pH value of growth medium was developed for harvesting freshwater microalgae. The flocculation efficiencies were as high as 90% for Chlorococcum nivale, Chlorococcum ellipsoideum and Scenedesmus sp. with high biomass concentrations (>1g/L). The optimum flocculation efficiency was achieved at pH 4.0. The flocculation mechanism could be that the carboxylate ions of organic matters adhering on microalgal cells accepted protons when pH decreases and the negative charges were neutralized, resulting in disruption of the dispersing stability of cells and subsequent flocculation of cells. A linear correlation between biomass concentration and acid dosage was observed. Furthermore, viability of flocculated cells was determined by Evans Blue assay and few cells were found to be damaged with pH decrease. After neutralizing pH and adding nutrients to the flocculated medium, microalgae were proved to maintain a similar growth yield in the flocculated medium comparing with that in the fresh medium. The recycling of medium could contribute to the economical production from algae to biodiesel. CONCLUSIONS: The study provided an economical, efficient and convenient method to harvest fresh microalgae. Advantages include capability of treating high cell biomass concentrations (>1g/L), excellent flocculation efficiencies (≥ 90%), operational simplicity, low cost and recycling of medium. It has shown the potential to overcome the hurdle of harvesting microalgae to promote full-scale application to biofuels from microalgae.