In this paper, we investigate the interconversion processes of the major flame retardant - 1,2,5,6,9,10-hexabromocyclododecane (HBCD) - by the means of statistical thermodynamics based on classical force-fields. Three ideas will be presented. First, the application of classical hybrid Monte-Carlo simulations for quantum mechanical processes will be justified. Second, the problem of insufficient convergence properties of hybrid Monte-Carlo methods for the generation of low temperature canonical ensembles will be solved by an interpolation approach. Furthermore, it will be shown how free energy differences can be used for a rate matrix computation. The results of our numerical simulations will be compared to experimental results.
A novel approach for the simulation of host.guest systems by systematically scanning the host moleculefs orientations within the guest cavity is presented along with a thermodynamic strategy for determining preferential binding modes and corresponding optimal interaction energies between host and guest molecules. By way of example, the elution order of hexabromocyclododecane stereoisomers from high performance liquid chromatography separation on a permethylated β-cyclcodextrin stationary phase has been computed using classical molecular dynamics simulations with the explicit solvents water and acetonitrile. Comparison of estimated with experimental separation data reveals remarkable squared coefficients of correlation with R2=0.87 and a very high correlation RLOO2 . 0:72 using the leaveone-out cross-validation method and water as solvent. In particular, the approach presented shapes up as very robust in terms of the evaluated time range under consideration, reflecting well thermodynamic equilibria. These and further observations correlating with experimental results suggest the suitability of the underlying force fields and our multi-mode approach for the estimation of relative binding affinities for host–guest systems with unknown binding modes.
Ergopeptides, like ergocornine and a-ergocryptine, exist in an S- and in an R-configuration. Kinetic experiments imply that certain configurations are preferred depending on the solvent. The experimental methods are explained in this article. Furthermore, computational methods are used to understand this configurational preference. Standard quantum chemical methods can predict the favored configurations by using minimum energy calculations on the potential energy landscape. However, the explicit role of the solvent is not revealed by this type of methods. In order to better understand its influence, classical mechanical molecular simulations are applied. It appears from our research that 'folding' the ergopeptide molecules into an intermediate state (between the S- and the R-configuration) is mechanically hindered for the preferred configurations.