Spectroscopic detection of DNA quadruplexes by vibrational circular dichroism

Guanosine hydrogels in focus: A comprehensive analysis through mid-infrared spectroscopy

Guanosine nucleosides and nucleotides have the peculiar ability to self-assemble in water to form supramolecular complex architectures from G-quartets to G-quadruplexes. G-quadruplexes exhibit in turn a large liquid crystalline lyotropic polymorphism, but they eventually cross-link or entangle to form a densely connected 3D network (a molecular hydrogel), able to entrap very large amount of water (up to the 99% v/v). This high water content of the hydrogels enables tunable softness, deformability, self-healing, and quasi-liquid properties, making them ideal candidates for different biotechnological and biomedical applications. In order to fully exploit their possible applications, Attenuated Total Reflection-Fourier Transform InfraRed (ATR-FTIR) spectroscopy was used to unravel the vibrational characteristics of supramolecular guanosine structures. First, the characteristic vibrations of the known quadruplex structure of guanosine 5'-monophosphate, potassium salt (GMP/K), were investigated: the identified peaks reflected both the chemical composition of the sample and the formation of quartets, octamers, and quadruplexes. Second, the role of K+ and Na+ cations in promoting the quadruplex formation was assessed: infrared spectra confirmed that both cations induce the formation of G-quadruplexes and that GMP/K is more stable in the G-quadruplex organization. Finally, ATR-FTIR spectroscopy was used to investigate binary mixtures of guanosine (Gua) and GMP/K or GMP/Na, both systems forming G-hydrogels. The same G-quadruplex-based structure was found in both mixtures, but the proportion of Gua and GMP affected some features, like sugar puckering, guanine vibrations, and base stacking, reflecting the known side-to-side aggregation and bundle formation occurring in these binary systems.

What are the minimal folding seeds in proteins? Experimental and theoretical assessment of secondary structure propensities of small peptide fragments

Certain peptide sequences, some of them as short as amino acid triplets, are significantly overpopulated in specific secondary structure motifs in folded protein structures. For example, 74% of the EAM triplet is found in α-helices, and only 3% occurs in the extended parts of proteins (typically β-sheets). In contrast, other triplets (such as VIV and IYI) appear almost exclusively in extended parts (79% and 69%, respectively). In order to determine whether such preferences are structurally encoded in a particular peptide fragment or appear only at the level of a complex protein structure, NMR, VCD, and ECD experiments were carried out on selected tripeptides: EAM (denoted as pro-'α-helical' in proteins), KAM(α), ALA(α), DIC(α), EKF(α), IYI(pro-β-sheet or more generally, pro-extended), and VIV(β), and the reference α-helical CATWEAMEKCK undecapeptide. The experimental data were in very good agreement with extensive quantum mechanical conformational sampling. Altogether, we clearly showed that the pro-helical vs. pro-extended propensities start to emerge already at the level of tripeptides and can be fully developed at longer sequences. We postulate that certain short peptide sequences can be considered minimal "folding seeds". Admittedly, the inherent secondary structure propensity can be overruled by the large intramolecular interaction energies within the folded and compact protein structures. Still, the correlation of experimental and computational data presented herein suggests that the secondary structure propensity should be considered as one of the key factors that may lead to understanding the underlying physico-chemical principles of protein structure and folding from the first principles.

Light rare earth elements stabilize G-quadruplex structure in variants of human telomeric sequences

Recently, light rare earth elements (LREEs) are gaining importance in modern-day technologies. Thus, the entry of LREEs into biochemical pathways cannot be ignored, which might affect the conformation of biomacromolecules. Herein, for the first time, we discover the G-quadruplex formation in the human telomeric variants in presence of micromolar concentrations of LREEs. Thermal melting show that the LREE-induced unimolecular G-quadruplex structure. Isothermal titration calorimetry, UV-vis, and CD spectroscopy results suggest the binding stoichiometry of lanthanide ions to telomeric variants is 2:1. The data confirms that the LREE ions coordinate between adjacent G-quartets. The excess LREE ions are most likely binding to quadruplex loops. The CD spectra revealed that the LREE-induced quadruplex in human telomere and its variant have antiparallel orientation. The binding equilibria of LREEs have been studied both in the presence and absence of competing metal cations. Addition of LREEs to the Na+ or K+-induced G-quadruplexes led to conformational change, which may be ascribed to the displacement of K+ or Na+ ions by LREE ions and formation of a more compact LREE-induced G-quadruplex structure in human telomeric variant. Moreover, the thymine in the central loop of the human telomeric sequence stabilizes LREE induced G-quadruplex.

Molecular Vibrations in Chiral Europium Complexes Revealed by Near-Infrared Raman Optical Activity

Raman optical activity (ROA) is commonly measured with green light (532 nm) excitation. At this wavelength, however, Raman scattering of europium complexes is masked by circularly polarized luminescence (CPL). This can be avoided using near-infrared (near-IR, 785 nm) laser excitation, as demonstrated here by Raman and ROA spectra of three chiral europium complexes derived from camphor. Since luminescence is strongly suppressed, many vibrational bands can be detected. They carry a wealth of structural information about the ligand and the metal core, and can be interpreted based on density functional theory (DFT) simulations of the spectra. For example, jointly with ROA experimental data, the simulations make it possible to determine absolute configuration of chiral lanthanide compounds in solution.

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.

A Raman optical activity spectrometer can sensitively detect lanthanide circularly polarized luminescence

Recently, many studies have appeared in which the Raman optical activity (ROA) instrument was found to be convenient for measuring circularly polarized luminescence (CPL). Typically, weak lanthanide luminescence including circular polarization could be detected. The new detection scheme is referred to as ROA-CPL spectroscopy. It is particularly useful when also the vibrational (ROA) itself is detectable as the molecule structure can be examined more reliably. In this review, development of this chiroptical approach and its applications in structural studies of biomolecules are summarized.

Chiral Sum Frequency Generation Spectroscopy Detects Double-Helix DNA at Interfaces

Many DNA-based technologies involve the immobilization of DNA and therefore require a fundamental understanding of the DNA structure-function relationship at interfaces. We present three immobilization methods compatible with chiral sum frequency generation (SFG) spectroscopy at interfaces. They are the "anchor" method for covalently attaching DNA on a glass surface, the "island" method for dropcasting DNA on solid substrates, and the "buoy" method using a hydrocarbon moiety for localizing DNA at the air-water interface. Although SFG was previously used to probe DNA, the chiral and achiral SFG responses of single-stranded and double-stranded DNA have not been compared systemically. Using the three immobilization methods, we obtain the achiral and chiral C-H stretching spectra. The results introduce four potential applications of chiral SFG. First, chiral SFG gives null response from single-stranded DNA but prominent signals from double-stranded DNA, providing a simple binary readout for label-free detection of DNA hybridization. Second, with heterodyne detection, chiral SFG gives an opposite-signed spectral response useful for distinguishing native (D-) right-handed double helix from non-native (L-) left-handed double helix. Third, chiral SFG captures the aromatic C-H stretching modes of nucleobases that emerge upon hybridization, revealing the power of chiral SFG to probe highly localized molecular structures within DNA. Finally, chiral SFG is sensitive to macroscopic chirality but not local chiral centers and thus can detect not only canonical antiparallel double helix but also other DNA secondary structures, such as a poly-adenine parallel double helix. Our work benchmarks the SFG responses of DNA immobilized by the three distinct methods, building a basis for new chiral SFG applications to solve fundamental and biotechnological problems.

Aptamers targeting amyloidogenic proteins and their emerging role in neurodegenerative diseases

Aptamers are oligonucleotides selected from large pools of random sequences based on their affinity for bioactive molecules and are used in similar ways to antibodies. Aptamers provide several advantages over antibodies, including their small size, facile, large-scale chemical synthesis, high stability, and low immunogenicity. Amyloidogenic proteins, whose aggregation is relevant to neurodegenerative diseases, such as Alzheimer's, Parkinson's, and prion diseases, are among the most challenging targets for aptamer development due to their conformational instability and heterogeneity, the same characteristics that make drug development against amyloidogenic proteins difficult. Recently, chemical tethering of aptagens (equivalent to antigens) and advances in high-throughput sequencing-based analysis have been used to overcome some of these challenges. In addition, internalization technologies using fusion to cellular receptors and extracellular vesicles have facilitated central nervous system (CNS) aptamer delivery. In view of the development of these techniques and resources, here we review antiamyloid aptamers, highlighting preclinical application to CNS therapy.

Noble Metallic Pyramidal Substrate for Surface-Enhanced Raman Scattering Detection of Plasmid DNA Based on Template Stripping Method

In this paper, a new method for manufacturing flexible and repeatable sensors made of silicon solar cells is reported. The method involves depositing the noble metal film directly onto the Si template and stripping out the substrate with a pyramid morphology by using an adhesive polymer. In order to evaluate the enhancement ability of the substrate, Rhodamine 6G (R6G) were used as surface-enhanced Raman scattering (SERS) probe molecules, and the results showed a high sensitivity and stability. The limit of detection was down to 10−12 M for R6G. The finite-difference time domain (FDTD) was used to reflect the distribution of the electromagnetic field, and the electric field was greatly enhanced on the surface of the inverted pyramidal substrate, especially in pits. The mechanism of Raman enhancement of two types of pyramidal SERS substrate, before and after stripping of the noble metal film, is discussed. By detecting low concentrations of plasmid DNA, the identification of seven characteristic peaks was successfully realized using a noble metallic pyramidal substrate.

Variable temperature FTIR spectra of polycrystalline purine nucleobases and estimating strengths of individual hydrogen bonds

In the first part of this work, we report the FTIR spectra of pure NH and isotopically substituted ND (10-15% D and 80-90% D) polycrystalline hypoxanthine, xanthine, adenine and guanine recorded in the 400-4000 cm^-1 range, as a function of temperature (10-300 K). We provide assignments of the stretching and out-of-plane bending amine (NH₂) and imine (NH) bands to the distinct H-bonds present in the crystal, based on the temperature sensitivity and isotopic exchange behavior. Empirical correlations between spectral and thermodynamic or structural parameters enabled us to estimate the energies and lengths of H-bonds in the studied nucleobase crystals and to correlate them with literature data. The empirical H-bonding energies are compared with H-bonding and stacking energies computed for hypoxanthine. In the second part, strategies for using the empirical correlations together with information extracted from quantum mechanical data (in particular from the Bader's quantum theory of atoms in molecules, QTAIM) for the evaluation of hydrogen bonding properties are discussed, and their advantages and drawbacks pointed out. The justification for a cooperative use of quantum-mechanical calculations with empirical spectra-energy correlations is discussed.

How important are the intermolecular hydrogen bonding interactions in methanol solvent for interpreting the chiroptical properties?

Two crispine A analogs and tetrahydrofuro[2,3-b]furan-3,3a(6aH)-diol, endowed with hydroxyl groups that can participate in intramolecular hydrogen bonding, have been synthesized and experimental vibrational circular dichroism (VCD) spectra and optical rotatory dispersion (ORD) data have been measured in CD₃OD/CH₃OH solvents. The absolute configurations (ACs) of these compounds have been determined using their synthetic schemes, supplemented wherever possible with X-ray diffraction data. The ACs are also analyzed with quantum chemical (QC) calculations of VCD and ORD utilizing implicit solvation as well as explicit solvation models, with the later employing classical molecular dynamics (MD) simulations. It is found that VCD calculations with implicit solvation model are adequate for determining the ACs, despite propensity of studied compounds for intermolecular hydrogen bonding between solute and solvent molecules. This observation is important because time-consuming MD simulations may not be necessary in the type of situations studied here. Additionally, it is found that the QC predicted VCD spectra provided enough diastereomer discrimination for determining the correct AC of studied compounds independently. The same observation did not apply to ORD.

Proximity Ligation Assay Detection of Protein-DNA Interactions-Is There a Link between Heme Oxygenase-1 and G-quadruplexes?

G-quadruplexes (G4) are stacked nucleic acid structures that are stabilized by heme. In cells, they affect DNA replication and gene transcription. They are unwound by several helicases but the composition of the repair complex and its heme sensitivity are unclear. We found that the accumulation of G-quadruplexes is affected by heme oxygenase-1 (Hmox1) expression, but in a cell-type-specific manner: hematopoietic stem cells (HSCs) from Hmox1^−/− mice have upregulated expressions of G4-unwinding helicases (e.g., Brip1, Pif1) and show weaker staining for G-quadruplexes, whereas Hmox1-deficient murine induced pluripotent stem cells (iPSCs), despite the upregulation of helicases, have more G-quadruplexes, especially after exposure to exogenous heme. Using iPSCs expressing only nuclear or only cytoplasmic forms of Hmox1, we found that nuclear localization promotes G4 removal. We demonstrated that the proximity ligation assay (PLA) can detect cellular co-localization of G-quadruplexes with helicases, as well as with HMOX1, suggesting the potential role of HMOX1 in G4 modifications. However, this colocalization does not mean a direct interaction was detectable using the immunoprecipitation assay. Therefore, we concluded that HMOX1 influences G4 accumulation, but rather as one of the proteins regulating the heme availability, not as a rate-limiting factor. It is noteworthy that cellular G4–protein colocalizations can be quantitatively analyzed using PLA, even in rare cells.

Infrared Spectroscopy Reveals the Preferred Motif Size and Local Disorder in Parallel Stranded DNA G-Quadruplexes

Infrared spectroscopy detects the formation of G-quadruplexes in guanine-rich nucleic acid sequences through shifts in the guanine C=O stretch mode. Here, we use ultrafast 2D infrared (IR) spectroscopy and isotope substitution to show that these shifts arise from vibrational delocalization among stacked G-quartets. This provides a direct measure of the sizes of locally ordered motifs in heterogeneous samples with substantial disordered regions. We find that parallel-stranded, potassium-bound DNA G-quadruplexes are limited to five consecutive G-quartets and 3-4 consecutive layers are preferred for longer polyguanine tracts. The resulting potassium-dependent G-quadruplex assembly landscape reflects the polyguanine tract lengths found in genomes, the ionic conditions prevalent in healthy mammalian cells, and the onset of structural disorder in disease states. Our study describes spectral markers that can be used to probe other G-quadruplex structures and provides insight into the fundamental limits of their formation in biological and artificial systems.

Deconvolution of conformational exchange from Raman spectra of aqueous RNA nucleosides

Ribonucleic acids (RNAs) are key to the central dogma of molecular biology. While Raman spectroscopy holds great potential for studying RNA conformational dynamics, current computational Raman prediction and assignment methods are limited in terms of system size and inclusion of conformational exchange. Here, a framework is presented that predicts Raman spectra using mixtures of sub-spectra corresponding to major conformers calculated using classical and ab initio molecular dynamics. Experimental optimization allowed purines and pyrimidines to be characterized as predominantly syn and anti, respectively, and ribose into exchange between equivalent south and north populations. These measurements are in excellent agreement with Raman spectroscopy of ribonucleosides, and previous experimental and computational results. This framework provides a measure of ribonucleoside solution populations and conformational exchange in RNA subunits. It complements other experimental techniques and could be extended to other molecules, such as proteins and carbohydrates, enabling biological insights and providing a new analytical tool.

An RNA aptamer with potent affinity for a toxic dimer of amyloid β42 has potential utility for histochemical studies of Alzheimer's disease

Oligomers of β-amyloid 42 (Aβ42), rather than fibrils, drive the pathogenesis of Alzheimer's disease (AD). In particular, toxic oligomeric species called protofibrils (PFs) have attracted significant attention. Herein, we report RNA aptamers with higher affinity toward PFs derived from a toxic Aβ42 dimer than toward fibrils produced from WT Aβ42 or from a toxic, conformationally constrained Aβ42 variant, E22P-Aβ42. We obtained these RNA aptamers by using the preincubated dimer model of E22P-Aβ42, which dimerized via a linker located at Val-40, as the target of in vitro selection. This dimer formed PFs during incubation. Several physicochemical characteristics of an identified aptamer, E22P-AbD43, suggested that preferential affinity of this aptamer toward PFs is due to its higher affinity for the toxic dimer unit (K_D = 20 ± 6.0 nm) of Aβ42 than for less-toxic Aβ40 aggregates. Comparison of CD data from the full-length and random regions of E22P-AbD43 suggested that the preferential binding of E22P-AbD43 toward the dimer might be related to the formation of a G-quadruplex structure. E22P-AbD43 significantly inhibited the nucleation phase of the dimer and its associated neurotoxicity in SH-SY5Y human neuroblastoma cells. Of note, E22P-AbD43 also significantly protected against the neurotoxicity of WT Aβ42 and E22P-Aβ42. Furthermore, in an AD mouse model, E22P-AbD43 preferentially recognized diffuse aggregates, which likely originated from PFs or higher-order oligomers with curvilinear structures, compared with senile plaques formed from fibrils. We conclude that the E22P-AbD43 aptamer is a promising research and diagnostic tool for further studies of AD etiology.

Single-stranded DNA structural diversity: TAGGGT from monomers to dimers to tetramer formation

RATIONALE

DNA quadruplex structures have emerged as novel drug targets due to their role in preventing abnormal gene transcription and maintaining telomere stability. Trapped Ion Mobility Spectrometry-Mass Spectrometry (TIMS-MS), combined with theoretical modeling, is a powerful tool for studying the kinetic intermediates of DNA complexes formed in solution and interrogated in the gas phase after desolvation.

METHODS

A TAGGGT ssDNA sequence was purchased and studied in 10 mM ammonium acetate using nanospray electrospray ionization (nESI)-TIMS-MS in positive and negative ion mode. Collisional cross section (CCS) profiles were measured using internal calibration (Tune Mix). Theoretical structures were proposed based on molecular dynamics, charge location and geometry optimization for the most intense IMS bands based on the number of TAGGGT units, adduct form and charge states.

RESULTS

A distribution of monomeric, dimeric and tetrameric TAGGGT structures were formed in solution and separated in the gas phase based on their mobility and m/z value (e.g., [M + 2H]⁺² , [2M + 3H]⁺³ , [M - 2H]^-2 , [2M - 3H]^-3 , [4M + 4H]⁺⁴ , [4M + 3H + NH₄ ]⁺⁴ , [4M + 2H + 2NH₄ ]⁺⁴ and [4M + H + 3NH₄ ]⁺⁴ ). The high mobility resolution of the TIMS-MS analyzer permitted the observation of multiple CCS bands per molecular ion form. Comparison with theoretical candidate structures suggests that monomeric TAGGGT species are stabilized by A-T and G⁺ -G interactions, with the size of the conformer influenced by the proton location. In the case of the TAGGGT quadruplex, the protonated species displayed a broad CCS distribution, while six discrete conformers were stabilized by the presence of ammonium ions (n = 1-3).

CONCLUSIONS

This is the first observation of multiple conformations of TAGGGT complexes (n = 1, 2 and 4) in 10 mM ammonium acetate. Candidate structures with intramolecular interactions of the form of G⁺ -G and traditional A-T base pairing agreed with the experimental trends. Our results demonstrate the structural diversity of TAGGGT monomers, dimers and tetramers in the gas phase beyond the previously reported solution structure, using 10 mM ammonium acetate to replicate biological conditions.

Influence of the hydrophobic domain on the self-assembly and hydrogen bonding of hydroxy-amphiphiles

The amphiphiles 1-octadecanol (octadecyl (stearyl) alcohol, ODA) and 1,2-dioleoylglycerol (DOG) were studied by IR spectroscopy and X-ray diffraction combined with multiscale theoretical modeling. The computations allowed us to rationalize the experimental findings and deduce the supramolecular structure of the formed assemblies while providing a fairly detailed insight into their hydrogen-bonding patterns. IR spectra revealed that the amphiphilic assemblies dramatically differ in structural order and hydrogen-bond strength, both being high in ODA and low in DOG. On the other hand, both compounds demonstrated common features, namely a splitting of the IR bands arising from O-H stretching vibrations (νOH) as well as complete hydrophobicity. However, the observed phenomena have different origins in the two amphiphiles. While the νOH split in ODA occurs due to a vibrational coupling along the string of inter-layer O-HO hydrogen bonds, in DOG it arises from different types of hydrogen bonds (intra- and intermolecular). The hydrophobicity of ODA stems from the very tight O-HO hydrogen bonding network connecting the opposite monolayers in a densely packed tilted crystalline phase (L_c'), whereas in DOG it occurs because the polar sites are locked inside reverted micellar-like assemblies. ODA and DOG illustrate that, in the assemblies of amphiphilic hydroxyl compounds, hydrogen bonds can be formed in a wide structural latitude, which is primarily governed by the chemical nature of apolar chains. Such a wide structural variability of OH-involving hydrogen bonds can be essential for the biological functioning of relevant molecules, such as glycolipids, acylglycerols, and, potentially, glycoproteins and carbohydrates.

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.

A Vibrational Circular Dichroism Microsampling Accessory: Mapping Enhanced Vibrational Circular Dichroism in Amyloid Fibril Films

We report the first vibrational circular dichroism (VCD) measurement of spatial heterogeneity in a sample using infrared (IR) microsampling. Vibrational circular dichroism spectra are typically measured using a standard IR cell with an IR beam diameter of 10 mm or greater making it impossible to investigate the spatial heterogeneity of a solid film sample. We have constructed a VCD sampling assembly with either 3 mm or 1 mm spatial resolution. An XY-translation stage was used to measure spectra at different spatial locations producing IR and VCD maps of the sample. In addition, a rotating sample stage was employed using a dual photoelastic modulator (PEM) setup to suppress artifacts due to linear birefringence in solid-phase or film samples. Infrared and VCD mapping of an insulin fibril film has been carried out at both 3 and 1 mm spatial resolution, and lysozyme films were mapped at 1 mm resolution. The IR spectra of different spots vary in intensity due primarily to sample thickness. The changes in the VCD intensity across the map largely correlate to corresponding changes in the IR map. Closer inspection of the insulin map revealed changes in the relative intensities of the VCD spectra not present in the parent IR spectra, which indicated differences in the degree of supramolecular chirality of the fibrils in the various spatial regions. For lysozyme films, in addition to different degrees of supramolecular chirality, reversal of the net fibril chirality was observed. The large signal-to-noise ratio observed at 1 mm resolution implies the feasibility of further increasing the spatial resolution by one or two orders of magnitude for protein fibril film samples.

Silver colloids as plasmonic substrates for direct label-free surface-enhanced Raman scattering analysis of DNA

Unraveling the role played by the surface chemistry of silver colloids in the direct SERS analysis of DNA.

Stereochemical analysis of glycerophospholipids by vibrational circular dichroism

The stereochemistry of glycerophospholipids (GPLs) has been of interest for its roles in the evolution of life and in their biological activity. However, because of their structural complexity, no convenient method to determine their configuration has been reported. In this work, through the first systematic application of vibrational circular dichroism (VCD) spectroscopy to various diacylated GPLs, we have revealed that their chirality can be assigned by the sign of a VCD exciton couplet generated by the interaction of two carbonyl groups. This paper also presents spectroscopic evidence for the stereochemistry of GPLs isolated from bacteria, eukaryotes, and mitochondria.

Vibrational Properties of the Phosphate Group Investigated by Molecular Dynamics and Density Functional Theory

The phosphate group (PO2(-)) is an important building block occurring in many components of living matter including nucleic acids. It provides distinct features in vibrational spectra and is useful as a local probe of NA conformation and interactions with the environment. For this purpose, it is desirable to explore in detail various factors influencing spectral shapes of characteristic phosphate vibrations. In the present study, effects of the solvent and conformational averaging are analyzed for simple model molecules, dimethylphosphate, ethylmethylphosphate, and ethylmethylthiophosphate. Infrared absorption (IR) and Raman spectra were measured and calculated using a combination of molecular dynamics (MD) and density functional theory (DFT). To fully understand the link between the structure and the spectra, the solvent has to be explicitly included in the computational modeling. The results indicate that vibrational properties of the phosphate moiety are very sensitive to its conformation and interactions with the aqueous environment indeed. Polarizable continuum solvent models without explicit water molecules provided significantly worse agreement with the experiment. The combined MD/DFT approach captures well spectral characteristics for the model systems and constitutes the most reliable basis for exploration of phosphate vibrational properties in biomolecular structural studies.

Determination of absolute configuration and conformation of a cyclic dipeptide by NMR and chiral spectroscopic methods

Increasing precision of contemporary computational methods makes spectroscopies such as vibrational (VCD) and electronic (ECD) circular dichroism attractive for determination of absolute configurations (AC) of organic compounds. This is, however, difficult for polar, flexible molecules with multiple chiral centers. Typically, a combination of several methods provides the best picture of molecular behavior. As a test case, all possible stereoisomers with known AC (RS, SR, SS, and RR) of the cyclic dipeptide cyclo(Arg-Trp) (CAT) were synthesized, and the performances of the ECD, infrared (IR), VCD, Raman, Raman optical activity (ROA), and nuclear magnetic resonance (NMR) techniques for AC determination were investigated. The spectra were interpreted with the aid of density functional theory (DFT) calculations. Folded geometries stabilized by van der Waals and electrostatic interactions between the diketopiperazine (DKP) ring and the indole group are predicted to be preferred for CAT, with more pronounced folding due to Arg-Trp stacking in the case of SS/RR-CAT. The RS/SR isomers prefer a twist-boat puckering of the DKP ring, which is relatively independent of the orientation of the side chains. Calculated conformer-averaged VCD and ECD spectra explain most of the experimentally observed bands and allow for AC determination of the tryptophan side-chain, whereas the stereochemical configuration of the arginine side-chain is visible only in VCD. NMR studies provide characteristic long-range (2)J(C,H) and (3)J(C,H) coupling constants, and nuclear Overhauser effect (NOE) correlations, which in combination with either ECD or VCD also allow for complete AC determination of CAT.

Detection of mercury-TpT dinucleotide binding by Raman spectra: a computational study

The Hg(2+) ion stabilizes the thymine-thymine mismatched base pair and provides new ways of creating various DNA structures. Recently, such T-Hg-T binding was detected by the Raman spectroscopy. In this work, detailed differences in vibrational frequencies and Raman intensity patterns in the free TpT dinucleotide and its metal-mediated complex (TpT·Hg)(2) are interpreted on the basis of quantum chemical modeling. The computations verified specific marker Raman bands indicating the effect of mercury binding to DNA. Although the B3LYP functional well-describes the Raman frequencies, a dispersion correction had to be added for all atoms including mercury to obtain realistic geometry of the (TpT·Hg)(2) dimer. Only then, the DFT complex structure agreed with those obtained with the wave function-based MP2 method. The aqueous solvent modeled as a polarizable continuum had a minor effect on the dispersion interaction, but it stabilized conformations of the sugar and phosphate parts. A generalized definition of internal coordinate force field was introduced to monitor covalent bond mechanical strengthening and weakening upon the Hg(2+) binding. Induced vibrational frequency shifts were rationalized in terms of changes in electronic structure. The simulations thus also provided reliable insight into the complex structure and stability.

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.