Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Optimization of Alchemical Pathways Using Extended Thermodynamic Integration.

ABSTRACT: Thermodynamic integration (TI) is a commonly used method to determine free-energy differences. One of its disadvantages is that many intermediate λ-states need to be sampled in order to be able to integrate accurately over ⟨∂H/∂λ⟩. Here, we use the recently introduced extended TI to study alternative parameterizations of H(λ) and its influence on the smoothness of the ⟨∂H/∂λ⟩ curves as well as the efficiency of the simulations. We find that the extended TI approach can be used to select curves of low curvature. An optimal parameterization is suggested for the calculation of hydration free energies. For calculations of relative binding free energies, we show that optimized parameterizations of the Hamiltonian in the unbound state also effectively lower the curvature in the bound state of the ligand.

SUBMITTER: de Ruiter A

PROVIDER: S-EPMC7872317 | biostudies-literature | 2021 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Optimization of Alchemical Pathways Using Extended Thermodynamic Integration.

de Ruiter Anita A Petrov Drazen D Oostenbrink Chris C

Journal of chemical theory and computation 20201222 1

Thermodynamic integration (TI) is a commonly used method to determine free-energy differences. One of its disadvantages is that many intermediate λ-states need to be sampled in order to be able to integrate accurately over ⟨∂H/∂λ⟩. Here, we use the recently introduced extended TI to study alternative parameterizations of H(λ) and its influence on the smoothness of the ⟨∂H/∂λ⟩ curves as well as the efficiency of the simulations. We find that the extended TI approach can be us ...[more]

PMID: 33351609

Similar Datasets

Optimization of Protein-Ligand Electrostatic Interactions Using an Alchemical Free-Energy Method.

Project description:We present an explicit solvent alchemical free-energy method for optimizing the partial charges of a ligand to maximize the binding affinity with a receptor. This methodology can be applied to known ligand-protein complexes to determine an optimized set of ligand partial atomic changes. Three protein-ligand complexes have been optimized in this work: FXa, P38, and the androgen receptor. The sets of optimized charges can be used to identify design principles for chemical changes to the ligands which improve the binding affinity for all three systems. In this work, beneficial chemical mutations are generated from these principles and the resulting molecules tested using free-energy perturbation calculations. We show that three quarters of our chemical changes are predicted to improve the binding affinity, with an average improvement for the beneficial mutations of approximately 1 kcal/mol. In the cases where experimental data are available, the agreement between prediction and experiment is also good. The results demonstrate that charge optimization in explicit solvent is a useful tool for predicting beneficial chemical changes such as pyridinations, fluorinations, and oxygen to sulfur mutations.

| S-EPMC7007198 | biostudies-literature

Quantum algorithm for alchemical optimization in material design.

Project description:The development of tailored materials for specific applications is an active field of research in chemistry, material science and drug discovery. The number of possible molecules obtainable from a set of atomic species grow exponentially with the size of the system, limiting the efficiency of classical sampling algorithms. On the other hand, quantum computers can provide an efficient solution to the sampling of the chemical compound space for the optimization of a given molecular property. In this work, we propose a quantum algorithm for addressing the material design problem with a favourable scaling. The core of this approach is the representation of the space of candidate structures as a linear superposition of all possible atomic compositions. The corresponding 'alchemical' Hamiltonian drives the optimization in both the atomic and electronic spaces leading to the selection of the best fitting molecule, which optimizes a given property of the system, e.g., the interaction with an external potential as in drug design. The quantum advantage resides in the efficient calculation of the electronic structure properties together with the sampling of the exponentially large chemical compound space. We demonstrate both in simulations and with IBM Quantum hardware the efficiency of our scheme and highlight the results in a few test cases. This preliminary study can serve as a basis for the development of further material design quantum algorithms for near-term quantum computers.

| S-EPMC8179438 | biostudies-literature

RNA secondary structure prediction using deep learning with thermodynamic integration.

Project description:Accurate predictions of RNA secondary structures can help uncover the roles of functional non-coding RNAs. Although machine learning-based models have achieved high performance in terms of prediction accuracy, overfitting is a common risk for such highly parameterized models. Here we show that overfitting can be minimized when RNA folding scores learnt using a deep neural network are integrated together with Turner's nearest-neighbor free energy parameters. Training the model with thermodynamic regularization ensures that folding scores and the calculated free energy are as close as possible. In computational experiments designed for newly discovered non-coding RNAs, our algorithm (MXfold2) achieves the most robust and accurate predictions of RNA secondary structures without sacrificing computational efficiency compared to several other algorithms. The results suggest that integrating thermodynamic information could help improve the robustness of deep learning-based predictions of RNA secondary structure.

| S-EPMC7878809 | biostudies-literature

Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis.

Project description:Replica-exchange molecular dynamics (REMD) has been proven to efficiently improve the convergence of free-energy perturbation (FEP) calculations involving considerable reorganization of their surrounding. We previously introduced the FEP/(λ,H)-REMD algorithm for ligand binding, in which replicas along the alchemical thermodynamic coupling axis λ were expanded as a series of Hamiltonian boosted replicas along a second axis to form a two-dimensional replica-exchange exchange map [Jiang, W.; Roux, B., J. Chem. Theory Comput. 2010, 6 (9), 2559-2565]. Aiming to achieve a similar performance at a lower computational cost, we propose here a modified version of this algorithm in which only the end-states along the alchemical axis are augmented by boosted replicas. The reduced FEP/(λ,H)-REMD method with one-dimensional unbiased alchemical thermodynamic coupling axis λ is implemented on the basis of generic multiple copy algorithm (MCA) module of the biomolecular simulation program NAMD. The flexible MCA framework of NAMD enables a user to design customized replica-exchange patterns through Tcl scripting in the context of a highly parallelized simulation program without touching the source code. Two Hamiltonian tempering boosting scheme were examined with the new algorithm: a first one based on potential energy rescaling of a preidentified "solute" and a second one via the introduction of flattening torsional free-energy barriers. As two illustrative examples with reliable experiment data, the absolute binding free energies of p-xylene and n-butylbenzene to the nonpolar cavity of the L99A mutant of T4 lysozyme were calculated. The tests demonstrate that the new protocol efficiently enhances the sampling of torsional motions for backbone and side chains around the binding pocket and accelerates the convergence of the free-energy computations.

| S-EPMC6339808 | biostudies-literature

Rapid and accurate structure-based therapeutic peptide design using GPU accelerated thermodynamic integration.

Project description:Peptide-based therapeutics are an alternative to small molecule drugs as they offer superior specificity, lower toxicity, and easy synthesis. Here we present an approach that leverages the dramatic performance increase afforded by the recent arrival of GPU accelerated thermodynamic integration (TI). GPU TI facilitates very fast, highly accurate binding affinity optimization of peptides against therapeutic targets. We benchmarked TI predictions using published peptide binding optimization studies. Prediction of mutations involving charged side-chains was found to be less accurate than for non-charged, and use of a more complex 3-step TI protocol was found to boost accuracy in these cases. Using the 3-step protocol for non-charged side-chains either had no effect or was detrimental. We use the benchmarked pipeline to optimize a peptide binding to our recently discovered cancer target: EME1. TI calculations predict beneficial mutations using both canonical and non-canonical amino acids. We validate these predictions using fluorescence polarization and confirm that binding affinity is increased. We further demonstrate that this increase translates to a significant reduction in pancreatic cancer cell viability.

| S-EPMC8091910 | biostudies-literature

Thermodynamic analysis of biodegradation pathways.

Project description:Microorganisms provide a wealth of biodegradative potential in the reduction and elimination of xenobiotic compounds in the environment. One useful metric to evaluate potential biodegradation pathways is thermodynamic feasibility. However, experimental data for the thermodynamic properties of xenobiotics is scarce. The present work uses a group contribution method to study the thermodynamic properties of the University of Minnesota Biocatalysis/Biodegradation Database. The Gibbs free energies of formation and reaction are estimated for 914 compounds (81%) and 902 reactions (75%), respectively, in the database. The reactions are classified based on the minimum and maximum Gibbs free energy values, which accounts for uncertainty in the free energy estimates and a feasible concentration range relevant to biodegradation. Using the free energy estimates, the cumulative free energy change of 89 biodegradation pathways (51%) in the database could be estimated. A comparison of the likelihood of the biotransformation rules in the Pathway Prediction System and their thermodynamic feasibility was then carried out. This analysis revealed that when evaluating the feasibility of biodegradation pathways, it is important to consider the thermodynamic topology of the reactions in the context of the complete pathway. Group contribution is shown to be a viable tool for estimating, a priori, the thermodynamic feasibility and the relative likelihood of alternative biodegradation reactions. This work offers a useful tool to a broad range of researchers interested in estimating the feasibility of the reactions in existing or novel biodegradation pathways.

| S-EPMC4073797 | biostudies-literature

Protein thermostability calculations using alchemical free energy simulations.

Project description:Thermal stability of proteins is crucial for both biotechnological and therapeutic applications. Rational protein engineering therefore frequently aims at increasing thermal stability by introducing stabilizing mutations. The accurate prediction of the thermodynamic consequences caused by mutations, however, is highly challenging as thermal stability changes are caused by alterations in the free energy of folding. Growing computational power, however, increasingly allows us to use alchemical free energy simulations, such as free energy perturbation or thermodynamic integration, to calculate free energy differences with relatively high accuracy. In this article, we present an automated protocol for setting up alchemical free energy calculations for mutations of naturally occurring amino acids (except for proline) that allows an unprecedented, automated screening of large mutant libraries. To validate the developed protocol, we calculated thermodynamic stability differences for 109 mutations in the microbial Ribonuclease Barnase. The obtained quantitative agreement with experimental data illustrates the potential of the approach in protein engineering and design.

| S-EPMC2872215 | biostudies-literature

Computing Relative Free Energies of Solvation using Single Reference Thermodynamic Integration Augmented with Hamiltonian Replica Exchange.

Project description:This paper introduces an efficient single-topology variant of Thermodynamic Integration (TI) for computing relative transformation free energies in a series of molecules with respect to a single reference state. The presented TI variant that we refer to as Single-Reference TI (SR-TI) combines well-established molecular simulation methodologies into a practical computational tool. Augmented with Hamiltonian Replica Exchange (HREX), the SR-TI variant can deliver enhanced sampling in select degrees of freedom. The utility of the SR-TI variant is demonstrated in calculations of relative solvation free energies for a series of benzene derivatives with increasing complexity. Noteworthy, the SR-TI variant with the HREX option provides converged results in a challenging case of an amide molecule with a high (13-15 kcal/mol) barrier for internal cis/trans interconversion using simulation times of only 1 to 4 ns.

| S-EPMC2998072 | biostudies-literature

Adsorption of phenol using adsorbent derived from Saccharum officinarum biomass: optimization, isotherms, kinetics, and thermodynamic study.

Project description:The present research shows the application of Taguchi's design of experiment approach to optimize the process parameters for the removal of phenol onto surface of Saccharum officinarum biomass activated carbon (SBAC) from an aqueous solution to maximize adsorption capacity of SBAC. The effect of adsorption parameters viz. adsorbent dose (m), temperature (T), initial concentration (C0) and mixing time (t) on response characteristics i.e., adsorption capacity (qt) has been studied at three levels by using L9 orthogonal array (OA) which further analyzed by variance analysis (ANOVA) for adsorption data and signal/noise (S/N) ratio data by using 'larger the better' characteristics. Using ANOVA, the optimum parameters are found to be m = 2 g/L, C0 = 150 mg/L, T = 313 K and t = 90 min, resulting in a maximum adsorption capacity of 64.59 mg/g. Adopting ANOVA, the percentage contribution of each process parameter in descending order of sequence is adsorbent dose 59.97% > initial phenol concentration 31.70% > contact time 4.28% > temperature 4.04%. The phenol adsorption onto SBAC was best fitted with the pseudo-second-order kinetic model and follows the Radke-Prausnitz isotherm model. Thermodynamic parameters suggested a spontaneous, exothermic nature and the adsorption process approaches physisorption followed by chemisorption. Hence the application of Taguchi orthogonal array design is a cost-effective and time-efficient approach for carrying out experiments and optimizing procedures for adsorption of phenol and improve the adsorption capacity of SBAC.

| S-EPMC10603077 | biostudies-literature

Application of referenced thermodynamic integration to Bayesian model selection.

Project description:Evaluating normalising constants is important across a range of topics in statistical learning, notably Bayesian model selection. However, in many realistic problems this involves the integration of analytically intractable, high-dimensional distributions, and therefore requires the use of stochastic methods such as thermodynamic integration (TI). In this paper we apply a simple but under-appreciated variation of the TI method, here referred to as referenced TI, which computes a single model's normalising constant in an efficient way by using a judiciously chosen reference density. The advantages of the approach and theoretical considerations are set out, along with pedagogical 1 and 2D examples. The approach is shown to be useful in practice when applied to a real problem -to perform model selection for a semi-mechanistic hierarchical Bayesian model of COVID-19 transmission in South Korea involving the integration of a 200D density.

| S-EPMC10424863 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data