Conformational preferences of two C-glycosyl analogues of Manp-(1 → 3)-Manp, were studied using a combined method of theoretical and experimental chemistry. Molecular dynamics was utilized to provide conformational behavior along C-glycosidic bonds of methyl 3-deoxy-3-C-[(α-d-mannopyranosyl)methyl]-α-d- and l-mannopyranosides. The OPLS2005 and Glycam06 force fields were used. Simulations were performed with explicit water (TIP3P) and methanol. Results were compared with a complete conformational scan at the MM4 level with the dielectric constant corresponding to methanol. In order to verify predicted conformational preferences, vicinal 3JHH NMR coupling constants were calculated by the Karplus equation on simulated potential energy surfaces (PES). A set of new parameters for the Karplus equation was also designed. Predicted 3JHH were compared with experimental data. We also used reverse methodology, in which the 3JHH coupling constants were calculated at the DFT level for each family of (ϕ, ψ)-conformers separately and then experimental values were decomposed onto calculated 3JHH couplings in order to obtain experimentally derived populations of conformers. As an alternative method of evaluation of preferred conformers, analysis of sensitive 13C chemical shifts was introduced. We were able to thoroughly discuss several fundamental issues in predictions of preferred conformers of C-saccharides, such as the solvent effect, reliability of the force field, character of empirical Karplus equation or applicability of NMR parameters in predictions of conformational preferences in general.
Clusters of a solute and a few solvent molecules obtained from molecular dynamics (MD) are a powerful tool to study solvation effects by advanced quantum chemical (QC) methods. For spectroscopic properties strongly dependent on the solvation, however, a large number of clusters are needed for a good convergence. In this work, a parallel variable selection (PVS) method is proposed that in some cases efficiently reduces the number of clusters needed for the averaging. The mass, charge, or atomic density MD distributions are used as a secondary variable to preselect the most probable cluster geometries used for averaging of solute spectral properties. When applied to nuclear magnetic resonance chemical shift of a model alcohol, the method allowed one to significantly reduce the total computational time, by a factor of 10. Even larger savings were achieved for the modeling of Raman and Raman optical activity spectra of (S)-lactamide molecule dissolved in water. The results thus suggest that the PVS method can be generally used for simulations of spectroscopic properties of solvated molecules and makes multiscale MD/QC computations more affordable.
Molecular pincers or tweezers are designed to hold and release the target molecule. Potential applications involve drug distribution in medicine, environment technologies, or microindustrial techniques. Typically, the binding is dominated by van der Waals forces. Modeling of such complexes can significantly enhance their design; yet obtaining accurate complexation energies by theory is difficult. In this study, density functional theory (DFT) computations combined with dielectric continuum solvent model are compared with the potential of mean force approach using umbrella sampling and the weighted histogram analysis method (WHAM) with molecular dynamics (MD) simulations. For DFT, functional and basis set effects are discussed. The computed results are compared to experimental data based on NMR spectroscopic measurements of five synthesized tweezers based on the Tröger's basis. Whereas the DFT computations correctly provided the observed trends in complex stability, they failed to produce realistic magnitudes of complexation energies. Typically, the binding was overestimated by DFT if compared to experiment. The simpler semiempirical PM6-DH2X scheme proposed lately yielded better magnitudes of the binding energies than DFT but not the right order. The MD-WHAM simulations provided the most realistic Gibbs binding energies, although the approximate MD force fields were not able to reproduce completely the ordering of relative stabilities of model complexes found by NMR. Yet the modeling provides interesting insight into the complex geometry and flexibility and appears as a useful tool in the tweezers' design.