Organization of Carotenoid Aggregates in Membranes Studied Selectively using Resonance Raman Optical Activity.

In living organisms, carotenoids are incorporated in biomembranes, remarkably modulating their mechanical characteristics, fluidity, and permeability. Significant resonance enhancement of Raman optical activity (ROA) signals of carotenoid chiral aggregates makes resonance ROA (RROA), a highly selective tool to study exclusively carotenoid assemblies in model membranes. Hence, RROA is combined with electronic circular dichroism (ECD), dynamic light scattering (DLS), molecular dynamics, and quantum-chemical calculations to shed new light on the carotenoid aggregation in dipalmitoylphosphatidylcholine (DPPC) liposomes. Using representative members of the carotenoid family: apolar α-carotene and more polar fucoxanthin and zeaxanthin, the authors demonstrate that the stability of carotenoid aggregates is directly linked with their orientation in membranes and the monomer structures inside the assemblies. In particular, polyene chain distortion of α-carotene molecules is an important feature of J-aggregates that show increased orientational freedom and stability inside liposomes compared to H-assemblies of more polar xanthophylls. In light of these results, RROA emerges as a new tool to study active compounds and drugs embedded in membranes.

Gold nanoshells with magnetic cores and a urea-based receptor for SERS sensing of fluoride anions: experimental and computational study.

The study demonstrates that a combination of plasmonic nanostructures and artificial receptors can be applied for sensing small molecular species. Gold nanoshells containing magnetic cores are used as the SERS-active substrates, which opens the way for the development of multimodal contrast agents with applicability extended to sensing or for the separation of analytes by magnetic solid-phase extraction. Disubstituted ureas forming hydrogen-bonded complexes with certain anions can be employed as molecular sensors. In this case study, gold nanoshells with silica-coated Mn-Zn ferrite cores were prepared by a multistep procedure. The nanoshells were co-functionalized with an N-(4-mercaptophenyl)-N'-(4-nitrophenyl)urea sensor synthesized directly on the gold surface, and with 4-nitrothiophenol, which is adopted as an internal standard. SERS measurements were carried out with acetonitrile solutions of tetrabutylammonium fluoride (Bu4NF) over a concentration range of 10-10-10-1 mol L-1. The spectral response of the sensor is dependent on the fluoride concentration in the range of 10-5-10-1 mol L-1. To investigate further the SERS mechanism, a model sensor, N-(4-bromophenyl)-N'-(4-nitrophenyl)urea, was synthesized and used in Raman spectroscopy with solutions of Bu4NF, up to a molar ratio of 1 : 20. The spectra and the interactions between the sensors and fluoride anions were also studied by extensive DFT computations.

Molecular dynamics and Raman optical activity spectra reveal nucleotide conformation ratios in solution.

Nucleotide conformational flexibility affects their biological functions. Although the spectroscopy of Raman optical activity (ROA) is well suited to structural analyses in aqueous solutions, the link between the spectral shape and the nucleotide geometry is not fully understood. We recorded the Raman and ROA spectra of model nucleotides (rAMP, rGMP, rCMP, and dTMP) and interpreted them on the basis of molecular dynamics (MD) combined with density functional theory (DFT). The relation between the sugar puckering, base conformation and spectral intensities is discussed. Hydrogen bonds between the sugar's C3' hydroxyl and the phosphate groups were found to be important for the sugar puckering. The simulated spectra correlated well with the experimental data and provided an understanding of the dependence of the spectral shapes on conformational dynamics. Most of the strongest spectral bands could be assigned to vibrational molecular motions. Decomposition of the experimental spectra into calculated subspectra based on arbitrary maps of free energies provided experimental conformer populations, which could be used to verify and improve the MD predictions. The analyses indicate some flaws of common MD force fields, such as being unable to describe the fine conformer distribution. Also the accuracy of conformer populations obtained from the spectroscopic data depends on the simulations, improvement of which is desirable for gaining a more detailed insight in the future. Improvement of the spectroscopic and computational methodology for nucleotides also provides opportunities for its application to larger nucleic acids.

Combination of Resonance and Non-Resonance Chiral Raman Scattering in a Cobalt(III) Complex.

Resonance Raman optical activity (RROA) spectra with high sensitivity reveal details on molecular structure, chirality, and excited electronic properties. Despite the difficulty of the measurements, the recorded data for the Co(III) complex with S,S-N,N-ethylenediaminedisuccinic acid are of exceptional quality and, coupled with the theory, spectacularly document the molecular behavior in resonance. This includes a huge enhancement of the chiral scattering, contribution of the antisymmetric polarizabilities to the signal, and the Herzberg-Teller effect significantly shaping the spectra. The chiral component is by about one order of magnitude bigger than for an analogous aluminum complex. The band assignment and intensity profile were confirmed by simulations based on density functional and vibronic theories. The resonance was attributed to the S0 →S3 transition, with the strongest signal enhancement of Raman and ROA spectral bands below about 800 cm-1 . For higher wavenumbers, other excited electronic states contribute to the scattering in a less resonant way. RROA spectroscopy thus appears as a unique tool to study the structure and electronic states of absorbing molecules in analytical chemistry, biology, and material science.

Low-frequency Raman optical activity provides insight into the structure of chiral liquids.

Vibrational frequencies of modes involving intermolecular motions in liquids are relatively small, in the Raman scattering close to the excitation frequency, and the bands may merge into a diverging uninterpretable signal. Raman optical activity (ROA) spectral shapes in this region, however, are structured more and may better reflect the nature of the studied systems. To understand the origin of the signal and its relation to the molecules, ROA spectra of six chiral neat liquids are recorded and analyzed on the basis of molecular dynamics and density functional theory computations. The theory of Raman scattering of liquids is discussed and adapted for modeling based on clusters and periodic boundary conditions. A plain cluster approach is compared to a crystal-like model. The results show that the low-frequency optical activity can be reliably modeled and related to the structure. However, momentary arrangement of molecules leads to large variations of optical activity, and a relatively large number of geometries need to be averaged for accurate simulations. The intermolecular modes are intertwined with intramolecular ones and start to dominate as the frequency goes down. The low-frequency ROA signal thus reflects the chemical composition and coupled with the modeling it provides a welcome means to study the structure and interactions of chiral liquids.

Polymorphism of Amyloid Fibrils Induced by Catalytic Seeding: A Vibrational Circular Dichroism Study.

Amyloidal protein fibrils occur in many biological events, but their formation and structural variability are understood rather poorly. We systematically explore fibril polymorphism for polyglutamic acid (PGA), insulin and hen egg white lysozyme. The fibrils were grown in the presence of "seeds", that is fibrils of the same or different protein. The seeds in concentrations higher than about 5 % of the total protein amount fully determined the structure of the final fibrils. Fibril structure was monitored by vibrational circular dichroism (VCD) spectroscopy and other techniques. The VCD shapes significantly differ for different fibril samples. Infrared (IR) and VCD spectra of PGA were also simulated using density functional theory (DFT) and a periodic model. The simulation provides excellent basis for data interpretation and reveals that the spectral shapes and signs depend both on fibril length and twist. The understanding of fibril formation and interactions may facilitate medical treatment of protein misfolding diseases in the future.

Intense chiral signal from α-helical poly-L-alanine observed in low-frequency Raman optical activity.

Raman optical activity (ROA) spectral features reliably indicate the structure of peptides and proteins, but the signal is often weak. However, we observed significantly enhanced low-frequency bands for α-helical poly-L-alanine (PLA) in solution. The biggest ROA signal at ∼100 cm-1 is about 10 times stronger than higher-frequency bands described previously, which facilitates the detection. The low-frequency bands of PLA were compared to those of α-helical proteins. For PLA, density functional simulations well reproduced the experimental spectra and revealed that about 12 alanine residues within two turns of the α-helix generate the strong ROA band. Averaging based on molecular dynamics (MD) provided an even more realistic spectrum compared to the static model. The low-frequency bands could be largely related to a collective motion of the α-helical backbone, partially modulated by the solvent. Helical and intermolecular vibrational coordinates have been introduced and the helical unwinding modes were assigned to the strongest ROA signal at 101-128 cm-1. Further analysis indicated that the helically arranged amide and methyl groups are important for the strong chiral signal of PLA, while the local chiral centers CαH contribute in a minor way only. The strong low-frequency ROA can thus provide precious information about the motions of the peptide backbone and facilitate future protein studies.

α-Synuclein conformations followed by vibrational optical activity. Simulation and understanding of the spectra.

α-Synuclein is a neuronal protein which adopts multiple conformations. These can be conveniently studied by the spectroscopy of vibrational optical activity (VOA). However, the interpretation of VOA spectra based on quantum-chemical simulations is difficult. To overcome the hampering of the computations by the protein size, we used the Cartesian tensor transfer technique to investigate links between the spectral shapes and protein structure. Vibrational circular dichroism (VCD) and Raman optical activity (ROA) spectra of α-synuclein in disordered, α-helical and ß-sheet (fibril) forms were measured and analyzed on the basis of molecular dynamics and density functional theory computations. For the disordered and α-helical conformers, a high fidelity of the simulated spectra with a reasonable computational cost was achieved. Most experimental spectral features could be assigned to the structure. So far unreported ROA marker bands of the secondary structure were found for the lower-frequency and CH stretching vibrations. Fibril VCD spectra were simulated with a rigid periodic model of the geometry and the results are consistent with previous studies based on cryogenic electron microscopy. The fibrils also give a specific ROA signal, but unlike VCD it is currently not fully explicable by the simulations. In connection with the computational modeling the VOA spectroscopy thus appears as an extremely useful tool for monitoring α-synuclein and other proteins in solutions.

Two Spectroscopies in One: Interference of Circular Dichroism and Raman Optical Activity.

Previously, we and other laboratories have reported an unusual and strong Raman optical activity (ROA) induced in solvents by chiral dyes. Various theories of the phenomenon appeared, but they were not capable of explaining fully the observed ROA band signs and intensities. In this work, an analysis based both on the light scattering theory and dedicated experiments provides a more complete understanding. For example, double-cell magnetic circular dichroism and magnetic ROA experiments with copper-porphyrin complex show that the induced chirality is observed without any contact of the solvents with the complex. The results thus indicate that a combination of electronic circular dichroism (ECD) with the polarized Raman scattering is responsible for the effect. The degree of circularity of solvent vibrational bands is a principal molecular property participating in the event. The insight and the possibility to predict the chirality transfer promise future applications in spectroscopy, chemical analysis and polarized imaging.

Transfer and Amplification of Chirality Within the "Ring of Fire" Observed in Resonance Raman Optical Activity Experiments.

We report extremely strong chirality transfer from a chiral nickel complex to solvent molecules detected as Raman optical activity (ROA). Electronic energies of the complex were in resonance with the excitation-laser light. The phenomenon was observed for a wide range of achiral and chiral solvents. For chiral 2-butanol, the induced ROA was even stronger than the natural one. The observations were related to so-called quantum (molecular) plasmons that enable a strong chiral Rayleigh scattering of the resonating complex. According to a model presented here, the maximal induced ROA intensity occurs at a certain distance from the solute, in a three-dimensional "ring of fire", even after rotational averaging. Most experimental ROA signs and relative intensities could be reproduced. The effect might significantly increase the potential of ROA spectroscopy in bioimaging and sensitive detection of chiral molecules.

Insight into vibrational circular dichroism of proteins by density functional modeling.

Vibrational circular dichroism (VCD) spectroscopy is an excellent method to determine the secondary structure of proteins in solution. Comparison of experimental spectra with quantum-chemical simulations represents a convenient and objective way to extract information on the structure. This has been difficult for such large molecules where approximate theoretical models have to be used. In the present study we applied the Cartesian-coordinate based tensor transfer (CCT) making it possible to extend the density functional theory (DFT) and model spectral intensities of large globular proteins nearly at quantum-chemical precision. Indeed, comparison with experiment provided a better understanding of the dependence of VCD spectral shapes on the geometry, their sensitivity to fine structural details and interactions with the environment. On a model set of globular proteins the simulated spectra correlated well with experimental data and revealed which structural information can (and cannot) be obtained from this kind of spectroscopy. Although the VCD technique has been regarded as being rather insensitive to side-chain variations, we found that the spectra of human and hen lysozyme differing by a few amino acids only are quite distinct. This has been explained by long-distance coupling of the amide vibrations. Likewise, the modeling reproduced some spectral changes caused by protein deuteration even when the protein structure was conserved.

Structure of supramolecular astaxanthin aggregates revealed by molecular dynamics and electronic circular dichroism spectroscopy.

Biomolecular aggregation is omnipresent in nature and important for metabolic processes or in medical treatment; however, the phenomenon is rather difficult to predict or understand on the basis of computational models. Recently, we found that electronic circular dichroism (ECD) spectroscopy and closely related resonance Raman optical activity (RROA) are extremely sensitive to the aggregation mechanism and structure of the astaxanthin dye. In the present study, molecular dynamics (MD) and quantum chemical (QC) computations (ZIndo/S, TDDFT) are used to link the aggregate structure with ECD spectral shapes. Realistic absorption and ECD intensities were obtained and the simulations reproduced many trends observed experimentally, such as the prevalent sign pattern and dependence of the aggregate structure on the solvent type. The computationally cheaper ZIndo/S method provided results very similar to those obtained by TDDFT. In the future, the accuracy of the combined MD/QC methodology of spectra interpretation should be improved to provide more detailed information on astaxanthin aggregates and similar macromolecular systems.

Establishing the link between fibril formation and Raman optical activity spectra of insulin.

Folding of proteins into insoluble amyloidal fibrils is implicated in a number of biological processes. Optical spectroscopy represents a convenient tool to monitor such structural variations. Recently, characteristic changes in Raman optical activity (ROA) spectra of insulin during a pre-fibrillar stage were reported but not supported by a theoretical model. In the present study, molecular dynamics and the density functional theory are used to simulate the spectra and understand the connection between the structure, and ROA and Raman spectral intensities. Theoretical results are consistent with the observations and only confirm exceptional ROA sensitivity to the protein tertiary structure. Surprisingly, this sensitivity reflects local conformational changes in the peptide main and side chains, rather than a direct through-space interaction of the protein components. Side chains providing strong ROA signals, such as tyrosine, can additionally report on local conformational features. Theoretical modeling helps in explaining the observed spectral changes and is likely to enable future applications of ROA spectroscopy in protein structural studies.

Detection of Sugars via Chirality Induced in Europium(III) Compounds.

Detection and resolution of simple monosaccharides are difficult tasks because their structure is quite similar. The present study shows that circularly polarized luminescence (CPL) induced in europium complexes provides very specific spectral patterns for fructose, mannose, glucose, and galactose. Differences were also observed between bare Eu(3+) ion and its complexes, when interacting with these sugars. The CPL spectra were measured on a Raman optical activity (ROA) spectrometer, which ensured high fluorescence intensity owing to the strong 532 nm laser excitation. The induced fluorescence was recorded in the same spectrum as the vibrational Raman bands. On the basis of the ligand field theory, most fluorescence spectral peaks could be assigned to f-shell europium transitions. Additional information on the interaction of the lanthanide with the sugar component was provided by measurement of time-dependent fluorescence, as formation of different complexes led to variations in fluorescence decay times. In nuclear magnetic resonance (NMR), the paramagnetic metal ion interacting with the sugars caused specific changes in (13)C chemical shifts. The spectroscopic data and molecular dynamics modeling showed that the interaction between the monosaccharides and Eu ion is rather weak due to the competition of the OH sugar groups with water molecules. However, multiple binding modes are possible, which contributes to the complexity and specificity of the spectra. The induced chirality and fluorescence spectra thus appear to be convenient means for monosaccharide detection and identification.

Chiral sensing of amino acids and proteins chelating with Eu(III) complexes by Raman optical activity spectroscopy.

Chiroptical spectroscopy of lanthanides sensitively reflects their environment and finds various applications including probing protein structures. However, the measurement is often hampered by instrumental detection limits. In the present study circularly polarized luminescence (CPL) of a europium complex induced by amino acids is monitored by Raman optical activity (ROA) spectroscopy, which enables us to detect weak CPL bands invisible to conventional CPL spectrometers. In detail, the spectroscopic response to the protonation state could be studied, e.g. histidine at pH = 2 showed an opposite sign of the strongest CPL band in contrast to that at pH = 7. The spectra were interpreted qualitatively on the basis of the ligand-field theory and related to CPL induced by an external magnetic field. Free energy profiles obtained by molecular dynamic simulations for differently charged alanine and histidine forms are in qualitative agreement with the spectroscopic data. The sensitivity and specificity of the detection promise future applications in probing peptide and protein side chains, chemical imaging and medical diagnosis. This potential is observed for human milk and hen egg-white lysozymes; these proteins have a similar structure, but very different induced CPL spectra.

Correction: Chiral sensing of amino acids and proteins chelating with Eu(III) complexes by Raman optical activity spectroscopy.

Correction for 'Chiral sensing of amino acids and proteins chelating with Eu(III) complexes by Raman optical activity spectroscopy' by Tao Wu et al., Phys. Chem. Chem. Phys., 2016, DOI: 10.1039/c6cp03968e.

Simulation of Raman optical activity of multi-component monosaccharide samples.

Determination of the saccharide structure in solution is a laborious process that can be significantly enhanced by optical spectroscopies. Raman optical activity (ROA) spectra are particularly sensitive to the chirality and conformation. However, the interpretation of them is largely dependent on computational tools providing a limited precision only. To understand the limitations and the link between spectral shapes and the structure, in the present study we measured and interpreted using a combination of molecular dynamics (MD) and density functional theory (DFT) Raman and ROA spectra of glucose and mannose solutions. Factors important for analyses of mixtures of conformers, anomers, and different monosaccharides are discussed as well. The accuracy of the simulations was found to be strongly dependent on the quality of the hydration model; the dielectric continuum solvent model provided lower accuracy than averaging of many solvent-solute clusters. This was due to different conformer weighting rather than direct involvement of water molecules in scattering recorded as ROA. However, the cluster-based simulations also failed to correctly reproduce the ratios of principal monosaccharide forms. The best results were obtained by a combined MD/DFT simulation, with the ratio of α- and ß-anomers and the -CH2OH group rotamers determined experimentally by NMR. Then a decomposition of experimental spectra into calculated subspectra provided realistic results even for the glucose and mannose mixtures. Raman spectra decomposition provided a better overall accuracy (∼5%) than ROA (∼10%). The combination of vibrational spectroscopy with theoretical simulations represents a powerful tool for analysing the saccharide structure. Conversely, the ROA and Raman data can be used to verify the quality of MD force fields and other parameters of computational modeling.

Circular Dichroism is Sensitive to Monovalent Cation Binding in Monensin Complexes.

Monensin is a natural antibiotic that exhibits high affinity to certain metal ions. In order to explore its potential in coordination chemistry, circular dichroism (CD) spectra of monensic acid A (MonH) and its derivatives containing monovalent cations (Li(+) , Na(+) , K(+) , Rb(+) , Ag(+) , and Et4 N(+) ) in methanolic solutions were measured and compared to computational models. Whereas the conventional CD spectroscopy allowed recording of the transitions down to 192 nm, synchrotron radiation circular dichroism (SRCD) revealed other bands in the 178-192 nm wavelength range. CD signs and intensities significantly varied in the studied compounds, in spite of their similar crystal structure. Computational modeling based on the Density Functional Theory (DFT) and continuum solvent model suggests that the solid state monensin structure is largely conserved in the solutions as well. Time-dependent Density Functional Theory (TDDFT) simulations did not allow band-to-band comparison with experimental spectra due to their limited precision, but indicated that the spectral changes were caused by a combination of minor conformational changes upon the monovalent cation binding and a direct involvement of the metal electrons in monensin electronic transitions. Both the experiment and simulations thus show that the CD spectra of monensin complexes are very sensitive to the captured ions and can be used for their discrimination. Chirality 28:420-428, 2016. © 2016 Wiley Periodicals, Inc.

Determination of absolute configuration in chiral solvents with nuclear magnetic resonance. A combined molecular dynamics/quantum chemical study.

Nuclear magnetic resonance (NMR) spectroscopy is omnipresent in chemical analysis. However, chirality of a molecule can only be detected indirectly by NMR, e.g., by monitoring its interaction with another chiral object. In the present study, we investigate the spectroscopic behavior of chiral molecules placed into a chiral solvent. In this case, the solvent-solute interaction is much weaker, but the application range of such NMR analysis is wider than for a specific chemical shift agent. Two alcohols and an amine were used as model systems, and differences in NMR chemical shifts dependent on the solute-solvent chirality combination were experimentally detected. Combined quantum mechanic/molecular mechanic (QM/MM) computations were applied to reveal the underlying solute-solvent interactions. NMR shielding was calculated using the density functional theory (DFT). While the experimental observations could not be reproduced quantitatively, the modeling provided a qualitative agreement and detailed insight into the essence of solvent-solute chiral interactions. The potentials of mean force (PMF) obtained using molecular dynamics (MD) and the weighted histogram analysis method (WHAM) indicate that the chiral interaction brings about differences in conformer ratios, which are to a large extent responsible for the NMR shifts. The MD results also predicted slight changes in the solvent structure, including the radial distribution function (RDF), to depend on the solvent/solute chirality combination. Apart from the conformer distribution, an effective average solvent electrostatic field was tested as another major factor contributing to the chiral NMR effect. The possibility to simulate spectral effects of chiral solvents from the first-principles opens up the way to NMR spectroscopic determination of the absolute configuration for a larger scale of compounds, including those not forming specific complexes.

Molecular dynamics with helical periodic boundary conditions.

Helical symmetry is often encountered in nature and thus also in molecular dynamics (MD) simulations. In many cases, an approximation based on infinite helical periodicity can save a significant amount of computer time. However, standard simulations with the usual periodic boundary conditions (PBC) are not easily compatible with it. In the present study, we propose and investigate an algorithm comprising infinitely propagated helicity, which is compatible with commonly used MD software. The helical twist is introduced as a parametric geometry constraint, and the translational PBC are modified to allow for the helical symmetry via a transitional solvent volume. The algorithm including a parallel code was implemented within the Tinker software. The viability of the helical periodic boundary conditions (HPBC) was verified in test simulations including α-helical and polyproline II like peptide structures. For an insulin-based model, the HPBC dynamics made it possible to simulate a fibrillar structure, otherwise not stable within PBC.