Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases.

Polarizable force fields are an essential component for the chemically accurate modeling of complex molecular systems with a significant degree of fluxionality, beyond harmonic or perturbative approximations. In this contribution we examine the performance of such an approach for the vibrational spectroscopy of the alanine amino acid, in the gas and condensed phases, from the Fourier transform of appropriate time correlation functions generated along molecular dynamics (MD) trajectories. While the infrared (IR) spectrum only requires the electric dipole moment, the vibrational circular dichroism (VCD) spectrum further requires knowledge of the magnetic dipole moment, for which we provide relevant expressions to be used with polarizable force fields. The AMOEBA force field was employed here to model alanine in the neutral and zwitterionic isolated forms, solvated by water or nitrogen, and as a crystal. Within this framework, comparison of the electric and magnetic dipole moments to those obtained with nuclear velocity perturbation theory based on density-functional theory for the same MD trajectories are found to agree well with one another. The statistical convergence of the IR and VCD spectra is examined and found to be more demanding in the latter case. Comparisons with experimental frequencies are also provided for the condensed phases.

Influence of the environment on the infrared spectrum of alanine: An effective mode analysis.

The vibrational spectrum of the alanine amino acid was computationally determined in the infrared range 1000-2000 cm-1, under various environments encompassing the gas, hydrated, and crystalline phases, by means of classical molecular dynamics trajectories, carried out with the Atomic Multipole Optimized Energetics for Biomolecular Simulation polarizable force field. An effective mode analysis was performed, in which the spectra are optimally decomposed into different absorption bands arising from well-defined internal modes. In the gas phase, this analysis allows us to unravel the significant differences between the spectra obtained for the neutral and zwitterionic forms of alanine. In condensed phases, the method provides invaluable insight into the molecular origins of the vibrational bands and further shows that peaks with similar positions can be traced to rather different molecular motions.

Control of Crystallinity and Stereocomplexation of Synthetic Carbohydrate Polymers from d- and l-Xylose.

Manipulating the stereochemistry of polymers is a powerful method to alter their physical properties. Despite the chirality of monosaccharides, reports on the impact of stereochemistry in natural polysaccharides and synthetic carbohydrate polymers remain absent. Herein, we report the cocrystallisation of regio- and stereoregular polyethers derived from d- and l-xylose, leading to enhanced thermal properties compared to the enantiopure polymers. To the best of our knowledge, this is the first example of a stereocomplex between carbohydrate polymers of opposite chirality. In contrast, atactic polymers obtained from a racemic mixture of monomers are amorphous. We also show that the polymer hydroxyl groups are amenable to post-polymerisation functionalization. These strategies afford a family of carbohydrate polyethers, the physical and chemical properties of which can both be controlled, and which opens new possibilities for polysaccharide mimics in biomedical applications or as advanced materials.

Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches.

Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra. Several approaches have been proposed to include them into computational models for calculating VCD signals, in particular those resting on the "cluster-in-a-liquid" model. Here we examine the capabilities of this ansatz on the example of flexible (1S,2S)-trans-1-amino-2-indanol solvated in dimethyl sulfoxide (DMSO). We compare cluster sets obtained from static calculations with results from explicit molecular dynamics (MD) trajectories based on either force field (FF) or first-principles (FP) methods. While the FFMD approach provides a broader sampling of configurational space, FPMD and time-correlation functions of dipole moments account for anharmonicity and entropy effects in the VCD calculation. They provide a means to evaluate the immediate effect of the solvent on the spectrum. This survey singles out several challenges associated with the use of clusters to describe solvation effects in systems showing shallow potential energy surfaces and non-covalent interactions. Static structures of clusters involving a limited number of solvent molecules satisfactorily capture the main effects of solvation in the bulk limit on the VCD spectra, if these structures are correctly weighted. The importance of taking into consideration their fluxionality, i.e. different solvent conformations sharing a same hydrogen bond pattern, and the limitations of small clusters for describing the solvent dynamics are discussed.