scholarly journals Oxidation of Isodrimeninol with PCC Yields Drimane Derivatives with Activity against Candida Yeast by Inhibition of Lanosterol 14-Alpha Demethylase

Biomolecules ◽  
2020 ◽  
Vol 10 (8) ◽  
pp. 1101
Author(s):  
Victor Marin ◽  
Andres Iturra ◽  
Andres Opazo ◽  
Bernd Schmidt ◽  
Matthias Heydenreich ◽  
...  
Keyword(s):  
Free Energy ◽  
Antifungal Activity ◽  
Candida Species ◽  
Binding Free Energy ◽  
Free Energy Calculations ◽  
Antifungal Drugs ◽  
Energy Estimates ◽  
Stable Complex ◽  
Pyridinium Chlorochromate ◽  
2D Nmr Spectroscopy

Candida species cause an opportunistic yeast infection called Candidiasis, which is responsible for more than 50,000 deaths every year around the world. Effective treatments against candidiasis caused by non-albicans Candida species such as C. glabrata, C. parapsilosis, C. aureus, and C. krusei are limited due to severe resistance to conventional antifungal drugs. Natural drimane sesquiterpenoids have shown promising antifungal properties against Candida yeast and have emerged as valuable candidates for developing new candidiasis therapies. In this work, we isolated isodrimeninol (C1) from barks of Drimys winteri and used it as starting material for the hemi-synthesis of four sesquiterpenoids by oxidation with pyridinium chlorochromate (PCC). The structure of the products (C2, C3, C4, and C5) was elucidated by 1D and 2D NMR spectroscopy resulting in C4 being a novel compound. Antifungal activity assays against C. albicans, C. glabrata, and C. krusei revealed that C4 exhibited an increased activity (IC50 of 75 μg/mL) compared to C1 (IC50 of 125 μg/mL) in all yeast strains. The antifungal activity of C1 and C4 was rationalized in terms of their capability to inhibit lanosterol 14-alpha demethylase using molecular docking, molecular dynamics simulations, and MM/GBSA binding free energy calculations. In silico analysis revealed that C1 and C4 bind to the outermost region of the catalytic site of 14-alpha demethylase and block the entrance of lanosterol (LAN) to the catalytic pocket. Binding free energy estimates suggested that C4 forms a more stable complex with the enzyme than C1, in agreement with the experimental evidence. Based on this new approach it is possible to design new drimane-type sesquiterpenoids for the control of Candida species as inhibitors of 14-alpha demethylase.

scholarly journals The SAMPL6 SAMPLing challenge: Assessing the reliability and efficiency of binding free energy calculations

2019 ◽  
Author(s):  
Andrea Rizzi ◽  
Travis Jensen ◽  
David R. Slochower ◽  
Matteo Aldeghi ◽  
Vytautas Gapsys ◽  
...  
Keyword(s):  
Free Energy ◽  
Convergence Rates ◽  
Binding Free Energy ◽  
Potential Of Mean Force ◽  
Free Energy Calculations ◽  
Energy Estimates ◽  
Initial Population ◽  
Energy Calculations ◽  
Partial Charges ◽  
Expanded Ensemble

AbstractApproaches for computing small molecule binding free energies based on molecular simulations are now regularly being employed by academic and industry practitioners to study receptor-ligand systems and prioritize the synthesis of small molecules for ligand design. Given the variety of methods and implementations available, it is natural to ask how the convergence rates and final predictions of these methods compare. In this study, we describe the concept and results for the SAMPL6 SAMPLing challenge, the first challenge from the SAMPL series focusing on the assessment of convergence properties and reproducibility of binding free energy methodologies. We provided parameter files, partial charges, and multiple initial geometries for two octa-acid (OA) and one cucurbit[8]uril (CB8) host-guest systems. Participants submitted binding free energy predictions as a function of the number of force and energy evaluations for seven different alchemical and physical-pathway (i.e., potential of mean force and weighted ensemble of trajectories) methodologies implemented with the GROMACS, AMBER, NAMD, or OpenMM simulation engines. To rank the methods, we developed an efficiency statistic based on bias and variance of the free energy estimates. For the two small OA binders, the free energy estimates computed with alchemical and potential of mean force approaches show relatively similar variance and bias as a function of the number of energy/force evaluations, with the attach-pull-release (APR), GROMACS expanded ensemble, and NAMD double decoupling submissions obtaining the greatest efficiency. The differences between the methods increase when analyzing the CB8-quinine system, where both the guest size and correlation times for system dynamics are greater. For this system, nonequilibrium switching (GROMACS/NS-DS/SB) obtained the overall highest efficiency. Surprisingly, the results suggest that specifying force field parameters and partial charges is insufficient to generally ensure reproducibility, and we observe differences between seemingly converged predictions ranging approximately from 0.3 to 1.0 kcal/mol, even with almost identical simulations parameters and system setup (e.g., Lennard-Jones cutoff, ionic composition). Further work will be required to completely identify the exact source of these discrepancies. Among the conclusions emerging from the data, we found that Hamiltonian replica exchange—while displaying very small variance—can be affected by a slowly-decaying bias that depends on the initial population of the replicas, that bidirectional estimators are significantly more efficient than unidirectional estimators for nonequilibrium free energy calculations for systems considered, and that the Berendsen barostat introduces non-negligible artifacts in expanded ensemble simulations.

Large-Scale Assessment of Binding Free Energy Calculations in Active Drug Discovery Projects

2020 ◽  
Author(s):  
Christina Schindler ◽  
Hannah Baumann ◽  
Andreas Blum ◽  
Dietrich Böse ◽  
Hans-Peter Buchstaller ◽  
...  
Keyword(s):  
Free Energy ◽  
Drug Discovery ◽  
Large Scale ◽  
Active Drug ◽  
Binding Free Energy ◽  
Free Energy Calculations ◽  
Energy Calculation ◽  
Large Scale Assessment ◽  
New Public ◽  
Binding Free Energy Calculations

Here we present an evaluation of the binding affinity prediction accuracy of the free energy calculation method FEP+ on internal active drug discovery projects and on a large new public benchmark set.<br>

SAMPL7 TrimerTrip Host-Guest Binding Poses and Binding Affinities from Spherical-Coordinates-Biased Simulations

2020 ◽  
Author(s):  
Zhaoxi Sun
Keyword(s):  
Free Energy ◽  
Computational Modelling ◽  
Free Energy Calculations ◽  
Energy Estimates ◽  
Binding Affinities ◽  
Spherical Coordinates ◽  
Collective Variables ◽  
Guest Binding ◽  
Starting Point ◽  
Sampl Challenge

Host-guest binding remains a major challenge in modern computational modelling. The newest 7<sup>th</sup> statistical assessment of the modeling of proteins and ligands (SAMPL) challenge contains a new series of host-guest systems. The TrimerTrip host binds to 16 structurally diverse guests. Previously, we have successfully employed the spherical coordinates as the collective variables coupled with the enhanced sampling technique metadynamics to enhance the sampling of the binding/unbinding event, search for possible binding poses and predict the binding affinities in all three host-guest binding cases of the 6<sup>th</sup> SAMPL challenge. In this work, we employed the same protocol to investigate the TrimerTrip host in the SAMPL7 challenge. As no binding pose is provided by the SAMPL7 host, our simulations initiate from randomly selected configurations and are proceeded long enough to obtain converged free energy estimates and search for possible binding poses. The predicted binding affinities are in good agreement with the experimental reference, and the obtained binding poses serve as a nice starting point for end-point or alchemical free energy calculations.

Automated, Accurate, and Scalable Relative Protein-Ligand Binding Free Energy Calculations using Lambda Dynamics

2020 ◽  
Author(s):  
E. Prabhu Raman ◽  
Thomas J. Paul ◽  
Ryan L. Hayes ◽  
Charles L. Brooks III
Keyword(s):  
Free Energy ◽  
Ligand Binding ◽  
Binding Affinity ◽  
Binding Free Energy ◽  
Computational Cost ◽  
Combinatorial Libraries ◽  
Free Energy Calculations ◽  
Lead Optimization ◽  
Efficient Estimation ◽  
Lead Compound

<p>Accurate predictions of changes to protein-ligand binding affinity in response to chemical modifications are of utility in small molecule lead optimization. Relative free energy perturbation (FEP) approaches are one of the most widely utilized for this goal, but involve significant computational cost, thus limiting their application to small sets of compounds. Lambda dynamics, also rigorously based on the principles of statistical mechanics, provides a more efficient alternative. In this paper, we describe the development of a workflow to setup, execute, and analyze Multi-Site Lambda Dynamics (MSLD) calculations run on GPUs with CHARMm implemented in BIOVIA Discovery Studio and Pipeline Pilot. The workflow establishes a framework for setting up simulation systems for exploratory screening of modifications to a lead compound, enabling the calculation of relative binding affinities of combinatorial libraries. To validate the workflow, a diverse dataset of congeneric ligands for seven proteins with experimental binding affinity data is examined. A protocol to automatically tailor fit biasing potentials iteratively to flatten the free energy landscape of any MSLD system is developed that enhances sampling and allows for efficient estimation of free energy differences. The protocol is first validated on a large number of ligand subsets that model diverse substituents, which shows accurate and reliable performance. The scalability of the workflow is also tested to screen more than a hundred ligands modeled in a single system, which also resulted in accurate predictions. With a cumulative sampling time of 150ns or less, the method results in average unsigned errors of under 1 kcal/mol in most cases for both small and large combinatorial libraries. For the multi-site systems examined, the method is estimated to be more than an order of magnitude more efficient than contemporary FEP applications. The results thus demonstrate the utility of the presented MSLD workflow to efficiently screen combinatorial libraries and explore chemical space around a lead compound, and thus are of utility in lead optimization.</p>

scholarly journals Automation of absolute protein-ligand binding free energy calculations for docking refinement and compound evaluation

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Germano Heinzelmann ◽  
Michael K. Gilson
Keyword(s):  
Free Energy ◽  
Early Stage ◽  
Binding Free Energy ◽  
Low Cost ◽  
Free Energy Calculations ◽  
Free Energies ◽  
Energy Calculations ◽  
Drug Candidates ◽  
Binding Free Energy Calculations ◽  
Binding Free Energies

AbstractAbsolute binding free energy calculations with explicit solvent molecular simulations can provide estimates of protein-ligand affinities, and thus reduce the time and costs needed to find new drug candidates. However, these calculations can be complex to implement and perform. Here, we introduce the software BAT.py, a Python tool that invokes the AMBER simulation package to automate the calculation of binding free energies for a protein with a series of ligands. The software supports the attach-pull-release (APR) and double decoupling (DD) binding free energy methods, as well as the simultaneous decoupling-recoupling (SDR) method, a variant of double decoupling that avoids numerical artifacts associated with charged ligands. We report encouraging initial test applications of this software both to re-rank docked poses and to estimate overall binding free energies. We also show that it is practical to carry out these calculations cheaply by using graphical processing units in common machines that can be built for this purpose. The combination of automation and low cost positions this procedure to be applied in a relatively high-throughput mode and thus stands to enable new applications in early-stage drug discovery.

scholarly journals Molecular modeling study on the drug resistance mechanism of NS5B polymerase to TMC647055

2016 ◽  
Vol 94 (2) ◽  
pp. 147-158 ◽  
Author(s):  
Huiqun Wang ◽  
Wei Cui ◽  
Chenchen Guo ◽  
Bo-Zhen Chen ◽  
Mingjuan Ji
Keyword(s):  
Drug Resistance ◽  
Free Energy ◽  
Rational Design ◽  
Binding Free Energy ◽  
Resistance Mechanism ◽  
Free Energy Calculations ◽  
Side Chain ◽  
Free Energy Decomposition ◽  
Hcv Ns5b ◽  
Ns5b Polymerase

NS5B polymerase plays an important role in viral replication machinery. TMC647055 (TMC) is a novel and potent non-nucleoside inhibitor of the HCV NS5B polymerase. However, mutations that result in drug resistance to TMC have been reported. In this study, we used molecular dynamics (MD) simulations, binding free energy calculations, and free energy decomposition to investigate the drug resistance mechanism of HCV to TMC resulting from L392I, P495T, P495S, and P495L mutations in NS5B polymerase. From the calculated results we determined that the decrease in the binding affinity between TMC and NS5BL392I polymerase is mainly caused by the extra methyl group at the CB atom of Ile. The polarity of the side-chain of residue 495 has no distinct influence on residue 495 binding with TMC, whereas the smaller size of the side-chain of residue 495 causes a substantial decrease in the van der Walls interaction between TMC and residue 495. Moreover, the longer length of the side-chain of residue 495 has a significant effect on the electrostatic interaction between TMC and Arg-503. Finally, we performed the same calculations and detailed analysis on other 3 mutations (L392V, P495V, and P495I). The results further confirmed our conclusions. The computational results not only reveal the drug resistance mechanism between TMC647055 and NS5B polymerase, but also provide valuable information for the rational design of more potent non-nucleoside inhibitors targeting HCV NS5B polymerase.

Assessing the ligand selectivity of sphingosine kinases using molecular dynamics and MM-PBSA binding free energy calculations

2016 ◽  
Vol 12 (4) ◽  
pp. 1174-1182 ◽  
Author(s):  
Liang Fang ◽  
Xiaojian Wang ◽  
Meiyang Xi ◽  
Tianqi Liu ◽  
Dali Yin
Keyword(s):  
Molecular Dynamics ◽  
Free Energy ◽  
Model Building ◽  
Binding Free Energy ◽  
Homology Model ◽  
Free Energy Calculations ◽  
Dynamics Simulation ◽  
Ligand Selectivity ◽  
Sphingosine Kinases ◽  
Binding Free Energy Calculations

Three residues of SK1 were identified important for selective SK1 inhibitory activity via SK2 homology model building, molecular dynamics simulation, and MM-PBSA studies.

scholarly journals Alchemical Absolute Protein-Ligand Binding Free Energies for Drug Design

2021 ◽  
Author(s):  
Yuriy Khalak ◽  
Gary Tresdern ◽  
Matteo Aldeghi ◽  
Hannah Magdalena Baumann ◽  
David L. Mobley ◽  
...  
Keyword(s):  
Free Energy ◽  
Drug Discovery ◽  
Drug Design ◽  
Ligand Binding ◽  
Binding Free Energy ◽  
Free Energy Calculations ◽  
Free Energies ◽  
Energy Calculations ◽  
Binding Free Energy Calculations ◽  
Binding Free Energies

The recent advances in relative protein-ligand binding free energy calculations have shown the value of alchemical methods in drug discovery. Accurately assessing absolute binding free energies, although highly desired, remains...

Sign in / Sign up

Export Citation Format

Share Document