Solvation and the secondary structure of a proline-containing dipeptide: insights from VCD spectroscopy

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.

Synthesis and VCD Spectroscopic Characterization of a Series of Azacryptands from a Chiral Valine-Based Derivative of Tris(2-aminoethyl)amine (TREN)

Utilizing experimental and computational vibrational circular dichroism (VCD) spectroscopy, we explored the conformational preferences of a series of chiral C3 -symmetric octaazacryptands with tris(2-aminoethyl)-amine head groups derived from valine. While the spectra of the smallest azacryptand with p-phenyl linkers and its elongated derivative with p-biphenyls linker were found to match well with the computed spectra, the computed conformational preferences of the m-biphenyl-based azacryptand did not seem to reflect the conformations dominating in chloroform solution. A detailed analysis revealed that structural changes resulting in a collapsed cage structure gave a notably better match with the experiment. It could subsequently be concluded from the VCD analysis, that the octaazacryptands prefer a collapsed structure, which is not predicted by density functional theory (DFT) calculations as the global minimum structures. These findings are expected to have consequences also for future studies on inclusion complexes of such azacryptands.

Modelling solute-solvent interactions in VCD spectra analysis with the micro-solvation approach

Vibrational circular dichroism (VCD) spectroscopy has become an important part of the (stereo-)chemists' toolbox as a reliable method for the determination of absolute configurations. Being the chiroptical version of infrared spectroscopy, it has also been recognized as being very sensitive to conformational changes and intermolecular interactions. This sensitivity originates from the fact that the VCD spectra of individual conformers are often more different than their IR spectra, so that changes in conformational distributions or band positions and intensities become more pronounced. What is an advantage for studies focussing on intermolecular interactions can, however, quickly turn into a major obstacle during AC determinations: solute-solvent interactions can have a strong influence on spectral signatures and they must be accurately treated when simulating VCD and IR spectra. In this perspective, we showcase selected examples which exhibit particularly pronounced solvent effects. It is demonstrated that it is typically sufficient to model solute-solvent interactions by placing single solvent molecules near hydrogen bonding sites of the solute and subsequently use the optimized structures for spectra simulations. This micro-solvation approach works reasonably well for medium-sized, not too conformationally flexible molecules. We thus also discuss its limitations and outline the next steps that method development needs to take in order to further improve the workflows for VCD spectra predictions.

Diverging conformations guide dipeptide self-assembly into crystals or hydrogels

The prediction of dipeptide assembly into crystals or gels is challenging. This work reveals the diverging conformational landscape that guides self-organization towards different outcomes. In silico and experimental data enabled deciphering of the electronic circular dichroism (ECD) spectra of self-assembling dipeptides to reveal folded or extended conformers as key players.

Impact of conformation and intramolecular interactions on vibrational circular dichroism spectra identified with machine learning

Vibrational Circular Dichroism (VCD) spectra often differ strongly from one conformer to another, even within the same absolute configuration of a molecule. Simulated molecular VCD spectra typically require expensive quantum chemical calculations for all conformers to generate a Boltzmann averaged total spectrum. This paper reports whether machine learning (ML) can partly replace these quantum chemical calculations by capturing the intricate connection between a conformer geometry and its VCD spectrum. Three hypotheses concerning the added value of ML are tested. First, it is shown that for a single stereoisomer, ML can predict the VCD spectrum of a conformer from solely the conformer geometry. Second, it is found that the ML approach results in important time savings. Third, the ML model produced is unfortunately hardly transferable from one stereoisomer to another.

Can the absolute configuration of cyclic peptides be determined with vibrational circular dichroism?

Cyclic peptides show a wide range of biological activities, among others as antibacterial agents. These peptides are often large and flexible with multiple chiral centers. The determination of the stereochemistry of molecules with multiple chiral centers is a challenging and important task in drug development. Chiroptical spectroscopies such as vibrational circular dichroism (VCD) can distinguish between different stereoisomers. The absolute configuration (AC) of a stereoisomer can be determined by comparing its experimental spectra to computed spectra of stereoisomers with known AC. In this way, the AC of rigid molecules with up to seven chiral centers has been assigned (Bogaerts et al., Phys. Chem. Chem. Phys., 2020, 22, 18014). The question arises whether this is possible with more conformationally flexible molecules such as cyclic peptides. We here investigate to what extent the AC of cyclic peptides can be determined with VCD. More specifically, we investigate the maximum number of chiral centers a cyclic peptide can have in order to be able to unambiguously assign the AC with VCD. We present experimental and computed IR and VCD spectra for a series of eight tetrapeptides and hexapeptides with two, three and four chiral centers. We use our recently developed computational protocol with a conformational search based on sampling with meta-dynamics. We use visual inspection to compare the computed spectra of different stereoisomers with an experimental spectrum of the corresponding cyclic peptide with known AC. We find that the AC of the investigated cyclic peptides with two chiral centers can be unambiguously assigned with VCD. This is however not possible for all of the cyclic peptides with three chiral centers and for none of those with four chiral centers. At best, one can limit the number of possible stereoisomers in those cases. Our work shows that other techniques are needed to assign the AC of cyclic peptides with three or more chiral centers. Our study also constitutes a warning that the spectra of all stereoisomers should be computed before attempting to match to an experimental spectrum, to avoid an accidental erroneous match.

Non-covalent interactions reveal the protein chain δ conformation in a flexible single-residue model

The δ conformation is a local secondary structure in proteins that implicates a π_amide N-H⋯N interaction between a backbone N atom and the NH of the following residue. Small-molecule models thereof have been limited so far to rigid proline-type compounds. We show here that in derivatives of a cyclic amino acid with a sulphur atom in the γ-position, specific side-chain/backbone N-H⋯S interactions stabilize the δ conformation sufficiently to allow it to compete with classical C5 and C7 H-bonded conformers.

Exhaustive Mapping of the Conformational Space of Natural Dipeptides by the DFT-D3//COSMO-RS Method

We extensively mapped energy landscapes and conformations of 22 (including three His protonation states) proteinogenic α-amino acids in trans configuration and the corresponding 484 (22²) dipeptides. To mimic the environment in a protein chain, the N- and C-termini of the studied systems were capped with acetyl and N-methylamide groups, respectively. We systematically varied the main chain dihedral angles (ϕ, ψ) by 40° steps and all side chain angles by 90° or 120° steps. We optimized the molecular geometries with the GFN2-xTB semiempirical (SQM) method and performed single point density functional theory calculations at the BP86-D3/DGauss-DZVP//COSMO-RS level in water, 1-octanol, N,N-dimethylformamide, and n-hexane. For each restrained (nonequilibrium) structure, we also calculated energy gradients (in water) and natural atomic charges. The exhaustive and unprecedented QM-based sampling enabled us to construct Ramachandran plots of quantum mechanical (QM(BP86-D3)//COSMO-RS) energies calculated on SQM structures, for all 506 (484 dipeptides and 22 amino acids) studied systems. We showed how the character of an amino acid side chain influences the conformational space of single amino acids and dipeptides. With clustering techniques, we were able to identify unique minima of amino acids and dipeptides (i.e., minima on the GFN2-xTB potential energy surfaces) and analyze the distribution of their BP86-D3//COSMO-RS conformational energies in all four solvents. We also derived an empirical formula for the number of unique minima based on the overall number of rotatable bonds within each peptide. The final peptide conformer data set (PeptideCs) comprises over 400 million structures, all of them annotated with QM(BP86-D3)//COSMO-RS energies. Thanks to its completeness and unbiased nature, the PeptideCs can serve, inter alia, as a data set for the validation of new methods for predicting the energy landscapes of protein structures. This data set may also prove to be useful in the development and reparameterization of biomolecular force fields. The data set is deposited at Figshare (10.25452/figshare.plus.19607172) and can be accessed using a simple web interface at http://peptidecs.uochb.cas.cz.

A Computational Protocol for Vibrational Circular Dichroism Spectra of Cyclic Oligopeptides

Cyclic peptides are a promising class of compounds for next-generation antibiotics as they may provide new ways of limiting antibiotic resistance development. Although their cyclic structure will introduce some rigidity, their conformational space is large and they usually have multiple chiral centers that give rise to a wide range of possible stereoisomers. Chiroptical spectroscopies such as vibrational circular dichroism (VCD) are used to assign stereochemistry and discriminate enantiomers of chiral molecules, often in combination with electronic structure methods. The reliable determination of the absolute configuration of cyclic peptides will require robust computational methods than can identify all significant conformers and their relative population and reliably assign their stereochemistry from their chiroptical spectra by comparison with ab initio calculated spectra. We here present a computational protocol for the accurate calculation of the VCD spectra of a series of flexible cyclic oligopeptides. The protocol builds on the Conformer-Rotamer Ensemble Sampling Tool (CREST) developed by Grimme and co-workers ( Phys. Chem. Chem. Phys. 2020, 22, 7169-7192 and J. Chem. Theory. Comput. 2019, 15, 2847-2862) in combination with postoptimizations using B3LYP and moderately sized basis sets. Our recommended computational protocol for the computation of VCD spectra of cyclic oligopeptides consists of three steps: (1) conformational sampling with CREST and tight-binding density functional theory (xTB); (2) energy ranking based on single-point energy calculations as well as geometry optimization and VCD calculations of conformers that are within 2.5 kcal/mol of the most stable conformer using B3LYP/6-31+G*/CPCM; and (3) VCD spectra generation based on Boltzmann weighting with Gibbs free energies. Our protocol provides a feasible basis for generating VCD spectra also for larger cyclic peptides of biological/pharmaceutical interest and can thus be used to investigate promising compounds for next-generation antibiotics.

Insights in the vibrational optical activity spectra of the antibiotic vancomycin in DMSO

Vibrational circular dichroism (VCD) and Raman optical activity (ROA) are two spectroscopic techniques that are sensitive towards the conformational behaviour of molecules, and are often complementary herein. In this work we pursue the determination of the conformational ensemble of the antibiotic glycopeptide vancomycin in DMSO through comparison of experimental and computational spectra, both for VCD and ROA. ROA is found to be highly suitable for the task, identifying an ensemble that strongly resembles the NMR conformation. In the case of VCD, however, a too high sensitivity of the intensities with respect to minor conformational changes hampers a reliable conformational analysis. Whence attempting to improve the match between the VCD experiment and calculations by any means - e.g., by inducing minor conformational changes or including solvent effects in the calculations - we show that there is the risk of going down the rabbit hole. In conclusion, this work contributes to the broader understanding of where, when and how VCD and ROA can be deployed as techniques for conformational analysis.

Solvation of the Boc-Val-Phe-nPr peptide characterized by VCD spectroscopy and DFT calculations

The conformational preferences of peptides are strongly determined by hydrogen bonding interactions. Intermolecular solute-solvent interactions compete with intramolecular interactions, which typically stabilize the secondary structure of the peptide. The analysis of vibrational circular dichroism (VCD) spectra can give insights into solvation-induced changes in the conformational distribution of small peptides. Here we describe the VCD spectroscopic characterization of the model peptide Boc-Val-Phe-nPr in chloroform as representative for a weakly interacting solvent and dimethyl sulfoxide (DMSO-d₆) as a strongly hydrogen bonding solvent. We show that the conformational preferences of the peptide in chloroform are well-described by the computationally predicted distribution of the isolated molecule assuming only implicit solvation effects through a continuum solvation model. In order to simulate the spectra recorded in DMSO-d₆, solvation was accounted for explicitly by computed microsolvated structures containing one to three solvent molecules. A good match of the computed spectra with the experimental data is obtained by this method. Comparing the conformational distributions in deuterated chloroform-d₁ and DMSO-d₆, structures with intramolecular hydrogen bonds such as the (δ,δ)-conformer family contribute to the conformational distribution only when there is no strong interaction with the solvent. This is in contrast to the results for the related Boc-Pro-Phe-nPr studied before, for which the intramolecular interaction was found to persist in DMSO-d₆. Furthermore, we discuss the influence of hydrogen bonding to different numbers of solvent molecules on the spectral signatures and show that the structure of the peptide in DMSO-d₆ is best described as a mixture of twofold-solvated (δ,β)- and threefold-solvated (β,β)-conformers.

Conformational preferences induced by cyclization in orbitides: a vibrational CD study

Orbitides are bioactive head-to-tail natural cyclic peptides from plant species. Their bioactivity is intrinsically related to the main conformations adopted in solution, whose correct characterization represents an important bottleneck for...

Anion-binding of a chiral tris(2-aminoethyl)amine-based tripodal thiourea: A spectroscopic and computational study

Thioureas are well-known structural motifs in supramolecular anion recognition. Their conformational preferences are typically characterized by detailed NMR spectroscopy and crystallography, which are often complemented with computational results from geometry...

VCD spectroscopy distinguishes the enamine and iminium ion of a 1,1’-binaphthyl azepine

We present a VCD spectroscopic characterization of a chiral 1,1’-binaphthyl azepine catalyst and show that the VCD spectra of an in-situ generated enamine and an ex-situ prepared iminium ion are...

Induced VCD and conformational chirality in host-guest complexes of a chiral ammonium salt with crown ethers

The hydrogen bonded complexes of the chiral ammonium salt α-methylbenzyl ammonium chloride (MBA-H⁺Cl^-) and the achiral crown ethers 18c6 and 15c5 serve as model systems to investigate the effect of host-guest complex formation on the conformational preferences of the macrocycles. We demonstrate that the intermolecular interactions result in new VCD signatures, that can be assigned to vibrational modes of the crown ethers. Based on a detailed conformational analysis, we investigate the origin of these signatures and discuss induced VCD (iVCD) and conformational chirality as possible sources of VCD intensity. The macrocycle in the MBA-H⁺/18c6 complex prefers either an achiral D_3d-symmetric conformation, which gives rise to iVCD, or chiral conformations, that feature individual contributions to the VCD spectrum. For the MBA-H⁺/15c5 complex, the contributions of the macrocycle to the VCD signatures are less pronounced and found to arise solely from conformational chirality. Therefore, analysis of the VCD signatures confirms that the small chiral guest molecule is able to affect the conformational preferences of a macrocyclic host. The study thus demonstrates the suitability of VCD spectroscopy for the characterization of analogous supramolecular host-guest complexes.

The important role of non-covalent interactions for the vibrational circular dichroism of lactic acid in aqueous solution

We present a computational study of vibrational circular dichroism (VCD) in solutions of (S)-lactic acid, relying on ab initio molecular dynamics (AIMD) and full solvation with bulk water. We discuss the effect of the hydrogen bond network on the aggregation behaviour of the acid: while aggregates of the solute represent conditions encountered in a weakly interacting solvent, the presence of water drastically interferes with the clusters - more strongly than originally anticipated. For both scenarios we computed the VCD spectra by means of nuclear velocity perturbation theory (NVPT). The comparison with experimental data allows us to establish a VCD-structure relationship that includes the solvent network around the chiral solute. We suggest that fundamental modes with strong polarisation such as the carbonyl stretching vibration can borrow VCD from the chirally restructured solvent cage, which extends the common explanatory models of VCD generation in aqueous solution.

Paving the way to conformationally unravel complex glycopeptide antibiotics by means of Raman optical activity

It is crucial for fundamental physical chemistry techniques to find their application in tackling real-world challenges. Hitherto, Raman optical activity (ROA) spectroscopy is one of the examples where a promising future within the pharmaceutical sector is foreseen, but has not yet been established. Namely, the technique is believed to be able to contribute in investigating the conformational behaviour of drug candidates. We, herein, strive towards the alignment of the ROA analysis outcome and the pharmaceutical expectations by proposing a fresh strategy that ensures a more complete, reliable, and transferable ROA study. The strategy consists of the treatment of the conformational space by means of a principal component analysis (PCA) and a clustering algorithm, succeeded by a thorough ROA spectral analysis and a novel way of estimating the contributions of the different chemical fragments to the total ROA spectral intensities. Here, vancomycin, an antibiotic glycopeptide, has been treated; it is the first antibiotic glycopeptide studied by means of ROA and is a challenging compound in ROA terms. By applying our approach we discover that ROA is capable of independently identifying the correct conformation of vancomycin in aqueous solution. In addition, we have a clear idea of what ROA can and cannot tell us regarding glycopeptides. Finally, the glycopeptide class turns out to be a spectroscopically curious case, as its spectral responses are unlike the typical ROA spectral responses of peptides and carbohydrates. This preludes future ROA studies of this intriguing molecular class.

VCD spectroscopy reveals that a water molecule determines the conformation of azithromycin in solution

We report the IR and VCD spectra of azithromycin, a macrolide antibiotic with a total of 18 stereogenic centers. The computational analysis of the spectra reveals that a single water molecule has to be considered in the conformational search. Its key role is the stabilization of an extended hydrogen bonding network and an otherwise unstable conformation that determines the VCD spectral signatures.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.