Vibrational Optical Activity of Cysteine in Aqueous Solution: A Comparison of Theoretical and Experimental Spectra

Vibrational circular dichroism spectra of proline in water at different pH values

Recording VCD spectra of aqueous solution poses a particular challenge as water is a strong infrared absorber. Likewise, the computational analysis of VCD spectra by means of DFT-based spectral calculations requires the consideration of explicit solvent molecules, thus posing an even greater challenge. Several studies suggested that by modeling the solvent environment with a few water molecules in a micro-solvation approach would be sufficient to describe experimental spectra. For example, using proline at different pH values, we herein show that a change in the relative spatial orientation of a single water molecule in five-fold solvated structures strongly affects the computed VCD spectral signatures and that Boltzmann-weighted spectra do not correctly reproduce the experiment. We thus explored an approach based on molecular dynamics and subsequent DFT-calculations, in which we considered 30 water molecules (about 1.5 solvation shells). Once again, it was found that the Boltzmann-weighted spectra obtained on the basis of several hundred structures did not correctly reproduce experimental signatures, and a simple averaging scheme resulted in well-matching spectra with comparable bandwidths. The rationale behind the procedure was that sampling the configurational space of the solvent molecules is as equally important as the conformational sampling of the solute. For conformationally more flexible molecules, it is assumed that a much larger set of structures will have to be computed in order to properly sample the conformational space of both solute and solvent.

Effective fully polarizable QM/MM approaches to compute Raman and Raman Optical Activity spectra in aqueous solution

Raman and Raman Optical Activity (ROA) signals are amply affected by solvent effects, especially in the presence of strongly solute-solvent interactions such as Hydrogen Bonding (HB). In this work, we extend the fully atomistic polarizable Quantum Mechanics/Molecular Mechanics approach, based on the Fluctuating Charges and Fluctuating Dipoles force field to the calculation of Raman and ROA spectra. Such an approach is able to accurately describe specific HB interactions, by also accounting for anisotropic contributions due to the inclusion of fluctuating dipoles. To highlight the potentiality of the novel approach, Raman and ROA spectra of L-Serine and L-Cysteine dissolved in aqueous solution are computed and compared both with alternative theoretical approaches and experimental measurements.

UV-vis and ECD spectroelectrochemistry of atropisomeric naphthalenediimide derivative

The UV-vis and ECD spectroelectrochemistry (SEC) of a chiral binaphthalenylamine derivative of the N-butyl naphthalenediimide (NDIB-NH₂) enantiomers were applied to measure UV-vis and ECD spectra of NDIB-NH₂ radicals and dianion formed in the reduction and oxidation processes observed in cyclic voltammetry (CV). The CV curves and EPR spectroelectrochemistry enabled us to establish conditions at which a radical-anion [NDIB-NH₂]̇^.-, a dianion [NDIB-NH₂]^2-, and a radical-cation [NDIB-NH₂]̇^.+ are formed. The DFT restricted open-shell CAM-B3LYP-D3/def2TZVP/PCM calculations demonstrated that in the radical-anion [NDIB-NH₂]̇^.-, spin is spread over the NDI system while in the radical-cation [NDIB-NH₂]̇⁺ it is spread over the aminonaphthalene moiety. The UV-vis spectra of radical-anion and dianion show the most significant changes in the 400-800 nm range. In that range, the ECD spectra varied with the change of electrode potential more than the UV-vis did and enabled the identification of a new ECD band of [NDIB-NH₂]̇^.- at ca. 400 nm hidden in the background in the UV spectra at -1000 mV. A broad structured ECD pattern with a maximum at ca. 530 nm was observed for [NDIB-NH₂]̇^.- (-1000 mV), while a single smooth ECD band of [NDIB-NH₂]^2- was located at 520 nm (-1750 mV). For the first time, an isosbestic point (455 nm) was found in ECD spectroelectrochemical measurements for the radical-cation [NDIB-NH₂]̇^.+ in equilibrium with the NDIB-NH₂ neutral form. The TD-DFT CAM-B3LYP-D3/6-31G** calculations combined with the hybrid (explicit combined with implicit) solvation model fairly well reproduced the UV-vis and ECD SEC of neutral and redox forms of NDIB-NH₂ but the ECD spectrum of [NDIB-NH₂]̇^.+ above 390 nm.

Chirality Discrimination at Binary Organic|Water Interfaces Monitored by Interfacial Tension Measurements with Preliminary Comparison with Molecular Dynamics Simulations

Chirality discrimination at a binary toluene (organic)/water(aqueous) interface between R- or S-Tol-BINAP (2,2'-Bis(di-p-tolylphosphino)-1,1'-binaphthyl) molecules and the water-soluble serine chiral specie is examined for the first time, using a combination of interfacial tension measurements and molecular dynamic simulations. Experimental interfacial measurements exhibit a clear chirality-controlled difference when a homochiral versus a heterochiral enantiomeric pairs are introduced at the interfaces. The related molecular dynamics simulations support the experimental results and provide further molecular insight of intermolecular interactions at the interfaces. The results indicate that interfacial tension measurements can capture the preferential interactions which exist between different pairs of enantiomers at the binary interfaces, opening up a new way for probing chirality discrimination at liquid-liquid interfaces.

Dissecting Bonding Interactions in Cysteine Dimers

Neutral (n) and zwitterionic (z) forms of cysteine monomers are combined in this work to extensively explore the potential energy surfaces for the formation of cysteine dimers in aqueous environments represented by a continuum. A simulated annealing search followed by optimization and characterization of the candidate structures afforded a total of 746 structurally different dimers held together via 80 different types of intermolecular contacts in 2894 individual non-covalent interactions as concluded from Natural Bond Orbitals (NBO), Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interactions (NCI) analyses. This large pool of interaction possibilities includes the traditional primary hydrogen bonds and salt bridges which actually dictate the structures of the dimers, as well as the less common secondary hydrogen bonds, exotic X⋯Y (X = C, N, O, S) contacts, and H⋯H dihydrogen bonds. These interactions are not homogeneous but have rather complex distributions of strengths, interfragment distances and overall stabilities. Judging by their Gibbs bonding energies, most of the structures located here are suitable for experimental detection at room conditions.

Chiral chromatography and surface chirality of carbon nanoparticles

Chiral carbon nanoparticles (CNPs) represent a rapidly evolving area of research for optical and biomedical technologies. Similar to small molecules, applications of CNPs as well as fundamental relationships between their optical activity and structural asymmetry would greatly benefit from their enantioselective separations by chromatography. However, this technique remains in its infancy for chiral carbon and other nanoparticles. The possibility of effective separations using high performance liquid chromatography (HPLC) with chiral stationary phases remains an open question whose answer can also shed light on the components of multiscale chirality of the nanoparticles. Herein, we report a detailed methodology of HPLC for successful separation of chiral CNPs and establish a path for its future optimization. A mobile phase of water/acetonitrile was able to achieve chiral separation of CNPs derived from L- and D-cysteine denoted as L-CNPs and D-CNPs. Molecular dynamics simulations show that the teicoplanin-based stationary phase has a higher affinity for L-CNPs than for D-CNPs, in agreement with experiments. The experimental and computational findings jointly indicate that chiral centers of chiral CNPs are present at their surface, which is essential for the multiple applications of these chiral nanostructures and equally essential for interactions with biomolecules and circularly polarized photons.

Revisiting conformations of methyl lactate in water and methanol

The recently developed efficient protocols to implicit [Grimme et al., J. Phys. Chem. A 125, 4039-4054 (2021)] and explicit quantum mechanical modeling of non-rigid molecules in solution [Katsyuba et al., J. Phys. Chem. B 124, 6664-6670 (2020)] are applied to methyl lactate (ML). Building upon this work, a new combination scheme is proposed to incorporate solvation effects for the computation of infrared (IR) absorption spectra. Herein, Boltzmann populations calculated for implicitly solvated single conformers are used to weight the IR spectra of explicitly solvated clusters with a size of typically ten solvent molecules, i.e., accounting for the first solvation shell. It is found that in water and methanol, the most abundant conformers of ML are structurally modified relative to the gas phase, where the major form is ML1, in which the syn conformation of the -OH moiety is stabilized by a OH⋯O=C intramolecular hydrogen bond (HB). In solution, this syn conformation transforms to the gauche form because the intramolecular HB is disrupted by explicit water molecules that form intermolecular HBs with the hydroxyl and carbonyl groups. Similar changes induced by the gas-solution transition are observed for the minor conformers, ML2 and/or ML3, characterized by OH⋯OCH₃ intramolecular HB in the gas phase. The relative abundance of ML1 is shown to decrease from ∼96% in gas to ∼51% in water and ∼92% in methanol. The solvent strongly influences frequencies, IR intensities, and normal modes, resulting in qualitatively different spectra compared to the gas phase. Some liquid-state conformational markers in the fingerprint region of IR spectra are revealed.

A Raman, SERS and UV-circular dichroism spectroscopic study of N-acetyl-l-cysteine in aqueous solutions

The aim of this work is to evaluate the vibrational and structural properties of N-acetyl-l-cysteine (NAC), and its molecular structure and electronic properties in relation to the action of thiol and amine groups at different pH.

Drop coating deposition Raman spectroscopy of proteinogenic amino acids compared with their solution and crystalline state

The Raman spectra of 20 proteinogenic amino acids were recorded in the solution, glass phase (as drop coating deposition Raman (DCDR) samples) and crystalline forms in the wide spectral range of 200-3200cm^-1. The most apparent spectral differences between the Raman spectra of the crystalline forms, glass phases and aqueous solutions of amino acids were briefly discussed and described in the frame of published works. The possible density dependencies of spectral bands were noted. In some cases, a strong influence of the sample density, as well as of the organization of the water envelope, was observed. The most apparent changes were observed for Ser and Thr. Nevertheless, for the majority of amino acids, the DCDR sample form is an intermediate between the solution and crystalline forms. In contrast, aromatic amino acids have only a small sensitivity to the form of the sample. Our reference set of Raman spectra is useful for revealing discrepancies between the SERS and solid/solution spectra of amino acids. We also found that some previously published Raman spectra of polycrystalline samples resemble glassy state rather than crystalline spectra. Therefore, this reference set of spectra will find application in every branch of Raman spectroscopy where the spectra of biomolecules are collected from coatings.

Solvation of a chiral carboxylic acid: effects of hydrogen bonding on the IR and VCD spectra of α-methoxyphenylacetic acid

Taking the title compound as a representative example of a chiral carboxylic acid, we investigate how its VCD spectral pattern is affected by hydrogen bonding of different solvents.

The Clusters-in-a-Liquid Approach for Solvation: New Insights from the Conformer Specific Gas Phase Spectroscopy and Vibrational Optical Activity Spectroscopy

Vibrational optical activity spectroscopies, namely vibrational circular dichroism (VCD) and Raman optical activity (ROA), have been emerged in the past decade as powerful spectroscopic tools for stereochemical information of a wide range of chiral compounds in solution directly. More recently, their applications in unveiling solvent effects, especially those associated with water solvent, have been explored. In this review article, we first select a few examples to demonstrate the unique sensitivity of VCD spectral signatures to both bulk solvent effects and explicit hydrogen-bonding interactions in solution. Second, we discuss the induced solvent chirality, or chiral transfer, VCD spectral features observed in the water bending band region in detail. From these chirality transfer spectral data, the related conformer specific gas phase spectroscopic studies of small chiral hydration clusters, and the associated matrix isolation VCD experiments of hydrogen-bonded complexes in cold rare gas matrices, a general picture of solvation in aqueous solution emerges. In such an aqueous solution, some small chiral hydration clusters, rather than the chiral solutes themselves, are the dominant species and are the ones that contribute mainly to the experimentally observed VCD features. We then review a series of VCD studies of amino acids and their derivatives in aqueous solution under different pHs to emphasize the importance of the inclusion of the bulk solvent effects. These experimental data and the associated theoretical analyses are the foundation for the proposed "clusters-in-a-liquid" approach to account for solvent effects effectively. We present several approaches to identify and build such representative chiral hydration clusters. Recent studies which applied molecular dynamics simulations and the subsequent snapshot averaging approach to generate the ROA, VCD, electronic CD, and optical rotatory dispersion spectra are also reviewed. Challenges associated with the molecular dynamics snapshot approach are discussed and the successes of the seemingly random "ad hoc explicit solvation" reported before are also explained. To further test and improve the "clusters-in-a-liquid" model in practice, future work in terms of conformer specific gas phase spectroscopy of sequential solvation of a chiral solute, matrix isolation VCD measurements of small chiral hydration clusters, and more sophisticated models for the bulk solvent effects would be highly valuable.

Conformational dynamics of carbohydrates: Raman optical activity of D-glucuronic acid and N-acetyl-D-glucosamine using a combined molecular dynamics and quantum chemical approach

As two biologically and medically relevant monosaccharides, the constituents of hyaluronic acid, d-glucuronic acid and N-acetyl-d-glucosamine, constitute perfect test cases for the development of carbohydrate-specific structural methods. These two molecules have been analysed by Raman optical activity (ROA), a spectroscopic technique exhibiting exquisite sensitivity to stereochemistry. We show that it is possible to support the experiment with a simulation approach combining density functional theory (DFT) and molecular dynamics (MD), both using explicit solvation. Thus, we have gained new insight into the crucial hydration effects that contribute to the conformational dynamics of carbohydrates and managed to characterize in detail the poorly understood vibrational nature of this class of biomolecules. Experimental and calculated ROA spectra of these two molecules are reported and excellent agreement has been found. More specifically, comparison has been made with the more commonly used gas phase and implicitly solvated calculation approaches, which offer poor or zero modelling of solvent interactions. The calculated spectra have been used to resolve the structural origins of the observed bands, a current challenge in the study of carbohydrates due to a lack of definitive vibrational assignments. We report and analyse major features in the fingerprint region of the ROA spectra, with recurrent structural and spectral features between the two monosaccharides observed.

The Raman optical activity of β-D-xylose: where experiment and computation meet

Besides its applications in bioenergy and biosynthesis, β-d-xylose is a very simple monosaccharide that exhibits relatively high rigidity. As such, it provides the best basis to study the impact of different solvation shell radii on the computation of its Raman optical activity (ROA) spectrum. Indeed, this chiroptical spectroscopic technique provides exquisite sensitivity to stereochemistry, and benefits much from theoretical support for interpretation. Our simulation approach combines density functional theory (DFT) and molecular dynamics (MD) in order to efficiently account for the crucial hydration effects in the simulation of carbohydrates and their spectroscopic response predictions. Excellent agreement between the simulated spectrum and the experiment was obtained with a solvation radius of 10 Å. Vibrational bands have been resolved from the computed ROA data, and compared with previous results on different monosaccharides in order to identify specific structure-spectrum relationships and to investigate the effect of the solvation environment on the conformational dynamics of small sugars. From the comparison with ROA analytical results, a shortcoming of the classical force field used for the MD simulations has been identified and overcome, again highlighting the complementary role of experiment and theory in the structural characterisation of complex biomolecules. Indeed, due to unphysical puckering, a spurious ring conformation initially led to erroneous conformer ratios, which are used as weights for the averaging of the spectral average, and only by removing this contribution was near perfect comparison between theory and experiment achieved.

A configurational and conformational study of (−)-Oseltamivir using a multi-chiroptical approach

Four chiroptical methods, i.e. electronic circular dichroism (ECD), optical rotatory dispersion (ORD), vibrational circular dichroism (VCD), and Raman optical activity (ROA) were employed to discover a set of the most probable conformations of (−)-Oseltamivir in solution.

A consistent force field parameter set for zwitterionic amino acid residues

Isolated amino acids play an important role in biochemistry and are therefore an interesting object of study. Atomistic molecular dynamics (MD) simulations can provide a high-resolution picture of the dynamic features of these species, especially in their biological environment. Unfortunately, most standard force field packages lack libraries for isolated amino acids in their zwitterionic form. Although several studies have used ad-hoc parameterizations for single amino acids, a consistent force-field parameter set for these molecules is still missing. Here, we present such a parameter library derived from the widely used parm99SB set from the AMBER program package. The parameter derivation for all 20 proteinogenic amino acids transparently followed established procedures with histidine treated in three different protonation states. All amino acids were subjected to MD simulations in four different forms for comparison: zwitterionic, N-teminally capped with acetyl, C-terminally capped with N-methyl, and capped at both termini. Simulation results show similarities between the different forms. Five zwitterionic amino acids-arginine, glutamate, glycine, phenylalanine, leucine-were simulated in a protein environment. Proteins and ligands generally retained their initial structure. The new parameter set will thus facilitate future atomistic simulations of these species.

Deducing the molecular properties of zwitterionic, protonated, deprotonated, and double-deprotonated forms of L-cysteine from vibrational spectroscopy (IR, Raman, VCD) and quantum chemical calculations

The behavior of L-cysteine (C3H7NO2S, (2R)-2-amino-3-sulfanylpropanoic acid) in water at different pH values was analyzed both experimentally and theoretically. The behavior was studied at pH values of 5.21 (at this pH, L-cysteine is a zwitterionic species), 1.00 (protonated species), 8.84 (monodeprotonated species), and 13.00 (dideprotonated species). We carried out a vibrational study using nonchiroptical (IR-Raman) and chiroptical (VCD) techniques complemented by quantum chemical calculations. We adopted a dual strategy, as follows. (i) The hybrid density functionals B3LYP and M062X and the ab initio MP2 method were employed, with the same 6-311++G (d,p) basis set, in order to characterize the relative energies and structures of an extensive set of conformers of L-cysteine. The presence of water was included by utilizing the IEF-PCM implicit solvation model. (ii) The vibrational analysis was made using a chirality-sensitive using a chirality-sensitive technique (VCD) and chirality-insensitive techniques (IR, including MIR and FIR, and Raman), especially in aqueous solution. The results obtained theoretically and experimentally were compared in order to deduce the most stable structures at each pH. Moreover, for the first time, the monodeprotonated anion of L-cysteine was detected in aqueous solution by means of IR, Raman and vibrational circular dichroism (VCD). Finally, analysis of the low-frequency region using the IR and Raman techniques was shown to be a very important way to understanding the conformational preference of the zwitterionic species.

Where does the Raman optical activity of [Rh(en)3]3+ come from? Insight from a combined experimental and theoretical approach

Raman optical activity (ROA) spectra are measured and calculated for Δ- and Λ-tris-(ethylenediamine)rhodium(iii) chloride in aqueous solution.