Cavity Closure Dynamics of Peracetylated β-Cyclodextrins in Supercritical Carbon Dioxide

Structural properties of peracetylated β-cyclodextrin in supercritical carbon dioxide were investigated by means of molecular dynamics simulations. The study indicated a strong reduction of the cavity accessibility to guest molecules, compared to native β-cyclodextrin in water. Indeed, the cavity is self-closed during the largest part of the simulation, which agrees well with suggestions made on the basis on high-pressure NMR experiments. Self-closure happens because one glucose unit undergoes a main conformational change (from chair to skew) that brings one of the acetyl groups in the wide rim of the cyclodextrin to the cavity interior. This arrangement turns out to be quite favorable, persisting for several nanoseconds. In addition to the wide rim self-closure, a narrow rim self-closure may also occur, though it is less likely and exhibits short duration (<1 ns). Therefore, the number of solvent molecules reaching the cavity interior is much smaller than that found in the case of native β-cyclodextrin in water after correction to account for different molar densities. These findings support the weak tendency of the macromolecule to form host–guest complexes in this nonconventional medium, as reported by some experiments. Finally, Lewis acid/base interactions between the acetyl carbonyl groups and the solvent CO₂ molecules were analyzed through ab initio calculations that revealed the existence of a quite favorable four-member ring structure not yet reported. The ensemble of these results can contribute to establish general thermodynamic principles controlling the formation of inclusion complexes in supercritical CO₂, where the hydrophilicity/hydrophobicity balance is not applicable

Taste for Chiral Guests: Investigating the Stereoselective Binding of Peptides to β‑Cyclodextrins

Obtaining compounds of diastereomeric purity is extremely important in the field of biological and pharmaceutical industry, where amino acids and peptides are widely employed. In this work, we theoretically investigate the possibility of chiral separation of peptides by β-cyclodextrins (β-CDs), providing a description of the associated interaction mechanisms by means of molecular dynamics (MD) simulations. The formation of host/guest complexes by including a model peptide in the macrocycle cavity is analyzed and discussed. We consider the terminally blocked phenylalanine dipeptide (Ace-Phe-Nme), in the l- and d-configurations, to be involved in the host/guest recognition process. The CD–peptide free energies of binding for the two enantiomers are evaluated through a combined approach that assumes: (1) extracting a set of independent molecular structures from the MD simulation, (2) evaluating the interaction energies for the host/guest complexes by hybrid quantum mechanics/molecular mechanics (QM/MM) calculations carried out on each structure, for which we also compute, (3) the solvation energies through the Poisson–Boltzmann surface area method. We find that chiral discrimination by the CD macrocycle is of the order of 1 kcal/mol, which is comparable to experimental data for similar systems. According to our results, the Ace-(d)­Phe-Nme isomer leads to a more stable complex with a β-CD compared to the Ace-(l)­Phe-Nme isomer. Nevertheless, we show that the chiral selectivity of β-CDs may strongly depend on the secondary structure of larger peptides. Although the free energy differences are relatively small, the predicted selectivities can be rationalized in terms of host/guest hydrogen bonds and hydration effects. Indeed, the two enantiomers display different interaction modes with the cyclodextrin macrocavity and different mobility within the cavity. This finding suggests a new interpretation for the interactions that play a key role in chiral recognition, which may be exploited to design more efficient and selective chiral separations of peptides

Hydration Effect on Amide I Infrared Bands in Water: An Interpretation Based on an Interaction Energy Decomposition Scheme

The sensitivity of some infrared bands to the local environment can be exploited to shed light on the structure and the dynamics of biological systems. In particular, the amide I band, which is specifically related to vibrations within the peptide bonds, can give information on the ternary structure of proteins, and can be used as a probe of energy transfer. In this work, we propose a model to quantitatively interpret the frequency shift on the amide I band of a model peptide induced by the formation of hydrogen bonds in the first solvation shell. This method allows us to analyze to what extent the electrostatic interaction, electronic polarization and charge transfer affect the position of the amide I band. The impact of the anharmoniticy of the pontential energy surface on the hydration induced shift is elucidated as well

Vibrational Energy Relaxation of the Amide I Mode of <i>N</i>‑Methylacetamide in D₂O Studied through Born–Oppenheimer Molecular Dynamics

The vibrational relaxation of the amide I mode of deuterated <i>N</i>-methylacetamide in D₂O solution is studied through nonequilibrium simulations using the semiempirical Born–Oppenheimer molecular dynamics (SEBOMD) approach to describe the whole solute–solvent system. Relaxation pathways and lifetimes are determined using the instantaneous normal mode (INM) analysis. The relaxation of the amide I mode is characterized by three different time scales; most of the excess energy (80%) is redistributed through intramolecular vibrational energy redistribution processes, with a smaller contribution (20%) of intermolecular energy flowing into the solvent. The amide II mode is found to contribute modestly (7%) to the relaxation mechanism. The amide I mode and the total vibrational energy decay curves obtained using SEBOMD and INM are in satisfactory agreement with the experimental measurements