Modelling the Encapsulation of 3-Hydroxyflavone with Cyclodextrin and Octa Acid and Comparing Their Differences.

The 3-hydroxyflavone (3-HF) is one of the common fluorescence probes. It has two distinct fluorescence bands: normal form and tautomer form. However, 3-hydroxyflavone has poor performance in water because of hydrogen bonding perturbation. The utilization of supramolecular chemistry would improve the fluorescence performance of 3-hydroxyflavone in water. In this paper, it reviews supramolecular chemistry of 3-hydroxyflavone with cyclodextrin and octa acid. Past research has found that the addition of ß-cyclodextrin to 3-hydroxyflavone in water would slightly improve the fluorescence intensity of the tautomer form. When adding γ-cyclodextrin to 3-hydroxyflavone in water, the green fluorescence intensity would be enhanced. Finally, the addition of octa acid creates a dry environment for the 3-hydroxyflavone, and it would only have a tautomer form. The ONIOM calculation shows the ways of self-assembly of ß- and γ-cyclodextrin. It can explain the difference in ratio between the tautomer form and normal form after understanding the interaction.

Spectroscopic Insights of an Emissive Complex between 4'-N,N-Diethylaminoflavonol in Octa-Acid Deep-Cavity Cavitand and Rhodamine 6G.

Excited-state chemistry relies on the communication between molecules, making it a crucial aspect of the field. One important question that arises is whether intermolecular communication and its rate can be modified when a molecule is confined. To explore the interaction in such systems, we investigated the ground and excited states of 4'-N,N-diethylaminoflavonol (DEA3HF) in an octa acid-based (OA) confined medium and in ethanolic solution, both in the presence of Rhodamine 6G (R6G). Despite the observed spectral overlap between the flavonol emission and the R6G absorption, as well as the fluorescence quenching of the flavonol in the presence of R6G, the almost constant fluorescence lifetime at different amounts of R6G discards the presence of FRET in the studied systems. Steady-state and time-resolved fluorescence indicate the formation of an emissive complex between the proton transfer dye encapsulated within water-soluble supramolecular host octa acid (DEA3HF@(OA)2) and R6G. A similar result was observed between DEA3HF:R6G in ethanolic solution. The respective Stern-Volmer plots corroborate with these observations, suggesting a static quenching mechanism for both systems.

The SAMPL6 SAMPLing challenge: assessing the reliability and efficiency of binding free energy calculations.

Approaches 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.

Host-guest interactions between octa acid and cations/nucleobases.

The nature of host-guest interaction in between octa acid cavitand (OA) and some representative cationic guests (Li+ , Na+ , K+ , Be+2 , Mg+2 , Ca+2 , Li3 O+ , Na3 O+ , K3 O+ ) as well as heterocyclic moieties like [adenine (A), guanine (G), cytosine (C), thymine (T), uracil (U), and tetrathiafulvalene (TTF)] has been examined with the aid of density functional theory (DFT)-based computations. Thermochemical results indicate that all the guests bind with OA in a thermodynamically favorable fashion at 298.15 K temperature and one atmospheric pressure. OA exhibits high selectivity in binding the lighter cations/metal cluster cations as compared to the heavier congeners along each given series. Moreover, OA exhibits enhanced affinity as well as selectivity in binding A/G/TTF molecules as compared to C/T/U. Noncovalent interaction and energy decomposition analyses reveal that in addition to the van der Waals interaction, significant contribution from electrostatic as well as orbital interactions dictate the outcome in all the host-guest complexes. Time dependent DFT calculations have been carried out to assess the role of the guests in tuning the electronic properties as well as absorption spectrum of OA. © 2017 Wiley Periodicals, Inc.

Overview of the SAMPL6 host-guest binding affinity prediction challenge.

Accurately predicting the binding affinities of small organic molecules to biological macromolecules can greatly accelerate drug discovery by reducing the number of compounds that must be synthesized to realize desired potency and selectivity goals. Unfortunately, the process of assessing the accuracy of current computational approaches to affinity prediction against binding data to biological macromolecules is frustrated by several challenges, such as slow conformational dynamics, multiple titratable groups, and the lack of high-quality blinded datasets. Over the last several SAMPL blind challenge exercises, host-guest systems have emerged as a practical and effective way to circumvent these challenges in assessing the predictive performance of current-generation quantitative modeling tools, while still providing systems capable of possessing tight binding affinities. Here, we present an overview of the SAMPL6 host-guest binding affinity prediction challenge, which featured three supramolecular hosts: octa-acid (OA), the closely related tetra-endo-methyl-octa-acid (TEMOA), and cucurbit[8]uril (CB8), along with 21 small organic guest molecules. A total of 119 entries were received from ten participating groups employing a variety of methods that spanned from electronic structure and movable type calculations in implicit solvent to alchemical and potential of mean force strategies using empirical force fields with explicit solvent models. While empirical models tended to obtain better performance than first-principle methods, it was not possible to identify a single approach that consistently provided superior results across all host-guest systems and statistical metrics. Moreover, the accuracy of the methodologies generally displayed a substantial dependence on the system considered, emphasizing the need for host diversity in blind evaluations. Several entries exploited previous experimental measurements of similar host-guest systems in an effort to improve their physical-based predictions via some manner of rudimentary machine learning; while this strategy succeeded in reducing systematic errors, it did not correspond to an improvement in statistical correlation. Comparison to previous rounds of the host-guest binding free energy challenge highlights an overall improvement in the correlation obtained by the affinity predictions for OA and TEMOA systems, but a surprising lack of improvement regarding root mean square error over the past several challenge rounds. The data suggests that further refinement of force field parameters, as well as improved treatment of chemical effects (e.g., buffer salt conditions, protonation states), may be required to further enhance predictive accuracy.

On the fly estimation of host-guest binding free energies using the movable type method: participation in the SAMPL5 blind challenge.

We review our performance in the SAMPL5 challenge for predicting host-guest binding affinities using the movable type (MT) method. The challenge included three hosts, acyclic Cucurbit[2]uril and two octa-acids with and without methylation at the entrance to their binding cavities. Each host was associated with 6-10 guest molecules. The MT method extrapolates local energy landscapes around particular molecular states and estimates the free energy by Monte Carlo integration over these landscapes. Two blind submissions pairing MT with variants of the KECSA potential function yielded mean unsigned errors of 1.26 and 1.53 kcal/mol for the non-methylated octa-acid, 2.83 and 3.06 kcal/mol for the methylated octa-acid, and 2.77 and 3.36 kcal/mol for Cucurbit[2]uril host. While our results are in reasonable agreement with experiment, we focused on particular cases in which our estimates gave incorrect results, particularly with regard to association between the octa-acids and an adamantane derivative. Working on the hypothesis that differential solvation effects play a role in effecting computed binding affinities for the parent octa-acid and the methylated octa-acid and that the ligands bind inside the pockets (rather than on the surface) we devised a new solvent accessible surface area term to better quantify solvation energy contributions in MT based studies. To further explore this issue a, molecular dynamics potential of mean force (PMF) study indicates that, as found by our docking calculations, the stable binding mode for this ligand is inside (rather than surface bound) the octa-acid cavity whether the entrance is methylated or not. The PMF studies also obtained the correct order for the methylation-induced change in binding affinities and associated the difference, to a large extent to differential solvation effects. Overall, the SAMPL5 challenge yielded in improvements our solvation modeling and also demonstrated the need for thorough validation of input data integrity prior to any computational analysis.

Converging ligand-binding free energies obtained with free-energy perturbations at the quantum mechanical level.

In this article, the convergence of quantum mechanical (QM) free-energy simulations based on molecular dynamics simulations at the molecular mechanics (MM) level has been investigated. We have estimated relative free energies for the binding of nine cyclic carboxylate ligands to the octa-acid deep-cavity host, including the host, the ligand, and all water molecules within 4.5 Å of the ligand in the QM calculations (158-224 atoms). We use single-step exponential averaging (ssEA) and the non-Boltzmann Bennett acceptance ratio (NBB) methods to estimate QM/MM free energy with the semi-empirical PM6-DH2X method, both based on interaction energies. We show that ssEA with cumulant expansion gives a better convergence and uses half as many QM calculations as NBB, although the two methods give consistent results. With 720,000 QM calculations per transformation, QM/MM free-energy estimates with a precision of 1 kJ/mol can be obtained for all eight relative energies with ssEA, showing that this approach can be used to calculate converged QM/MM binding free energies for realistic systems and large QM partitions. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.