Simulation of Electronic Circular Dichroism of Nucleic Acids: From the Structure to the Spectrum

A Perspective on the Simulation of Electronic Circular Dichroism and Circularly Polarized Luminescence Spectra in Chiral Solid Materials

Chiral materials have shown tremendous potential for many technological applications, such as optoelectronics, sensing, magnetism, information technology, and imaging. Characterization of these materials is mostly based on chiroptical spectroscopies, such as electronic circular dichroism (ECD) and circularly polarized luminescence (CPL). These experimental measurements would greatly benefit from theoretical simulations for interpretation of the spectra as well as predictions on new materials. While ECD and CPL simulations are well established for molecular systems, they are not for materials. In this Perspective, we describe the theoretical quantities necessary to simulate ECD and CPL spectra in oriented systems. Then, we discuss the approximate strategies currently used to perform these calculations, what computational machinery is already available to develop more general approaches, and some of the open challenges for the simulation of ECD and CPL spectra in solid materials. When methods that are as reliable and computationally efficient as those for molecules are developed, these simulations will provide invaluable insight and guidance for the rational design of optically active materials.

Experimental detection of conformational transitions between forms of DNA: problems and prospects

Under different conditions, the DNA double helix can take different geometric forms. Of the large number of its conformations, in addition to the "canonical" B form, the A, C, and Z forms are widely known, and the D, Hoogsteen, and X forms are less known. DNA locally takes the A, C, and Z forms in the cell, in complexes with proteins. We compare different methods for detecting non-canonical DNA conformations: X-ray, IR, and Raman spectroscopy, linear and circular dichroism in both the infrared and ultraviolet regions, as well as NMR (measurement of chemical shifts and their anisotropy, scalar and residual dipolar couplings and inter-proton distances from NOESY (nuclear Overhauser effect spectroscopy) data). We discuss the difficulties in applying these methods, the problems of theoretical interpretation of the experimental results, and the prospects for reliable identification of non-canonical DNA conformations.

Tidying up the conformational ensemble of a disordered peptide by computational prediction of spectroscopic fingerprints

The most advanced structure prediction methods are powerless in exploring the conformational ensemble of disordered peptides and proteins and for this reason the "protein folding problem" remains unsolved. We present a novel methodology that enables the accurate prediction of spectroscopic fingerprints (circular dichroism, infrared, Raman, and Raman optical activity), and by this allows for "tidying up" the conformational ensembles of disordered peptides and disordered regions in proteins. This concept is elaborated for and applied to a dodecapeptide, whose spectroscopic fingerprint is measured and theoretically predicted by means of enhanced-sampling molecular dynamics coupled with quantum mechanical calculations. Following this approach, we demonstrate that peptides lacking a clear propensity for ordered secondary-structure motifs are not randomly, but only conditionally disordered. This means that their conformational landscape, or phase-space, can be well represented by a basis-set of conformers including about 10 to 100 structures. The implications of this finding have profound consequences both for the interpretation of experimental electronic and vibrational spectral features of peptides in solution and for the theoretical prediction of these features using accurate and computationally expensive techniques. The here-derived methods and conclusions are expected to fundamentally impact the rationalization of so-far elusive structure-spectra relationships for disordered peptides and proteins, towards improved and versatile structure prediction methods.

Effect of the Thermal Fluctuations of the Photophysics of GC and CG DNA Steps: A Computational Dynamical Study

Here we refine and assess two computational procedures aimed to include the effect of thermal fluctuations on the electronic spectra and the ultrafast excited state dynamics of multichromophore systems, focusing on DNA duplexes. Our approach is based on a fragment diabatization procedure that, from a given Quantum Mechanical (QM) reference method, can provide the parameters (energy and coupling) of the reference diabatic states on the basis of the isolated fragments, either for a purely electronic excitonic Hamiltonian (FrDEx) or a linear vibronic coupling Hamiltonian (FrD-LVC). After having defined the most cost-effective procedure for DNA duplexes on two smaller fragments, FrDEx is used to simulate the absorption and Electronic Circular Dichroism (ECD) spectra of (GC)₅ sequences, including the coupling with the Charge Transfer (CT) states, on a number of structures extracted from classical Molecular Dynamics (MD) simulations. The computed spectra are close to the reference TD-DFT calculations and fully consistent with the experimental ones. We then couple MD simulations and FrD-LVC to simulate the interplay between local excitations and CT transitions, both intrastrand and interstrand, in GC and CG steps when included in a oligoGC or in oligoAT DNA sequence. We predict that for both sequences a substantial part of the photoexcited population on G and C decays, within 50-100 fs, to the corresponding intrastrand CT states. This transfer is more effective for GC steps that, on average, are more closely stacked than CG ones.

Organic materials repurposing, a data set for theoretical predictions of new applications for existing compounds

We present a data set of 48182 organic semiconductors, constituted of molecules that were prepared with a documented synthetic pathway and are stable in solid state. We based our search on the Cambridge Structural Database, from which we selected semiconductors with a computational funnel procedure. For each entry we provide a set of electronic properties relevant for organic materials research, and the electronic wavefunction for further calculations and/or analyses. This data set has low bias because it was not built from a set of materials designed for organic electronics, and thus it provides an excellent starting point in the search of new applications for known materials, with a great potential for novel physical insight. The data set contains molecules used as benchmarks in many fields of organic materials research, allowing to test the reliability of computational screenings for the desired application, "rediscovering" well-known molecules. This is demonstrated by a series of different applications in the field of organic materials, confirming the potential for the repurposing of known organic molecules.

Bright Frenkel Excitons in Molecular Crystals: A Survey

We computed the optical properties of a large set of molecular crystals (∼2200 structures) composed of molecules whose lowest excited states are strongly coupled and generate wide excitonic bands. Such bands are classified in terms of their dimensionality (1-, 2-, and 3-dimensional), the position of the optically allowed state in relation with the excitonic density of states, and the presence of Davydov splitting. The survey confirms that one-dimensional aggregates are rare in molecular crystals highlighting the need to go beyond the simple low-dimensional models. Furthermore, this large set of data is used to search for technologically interesting and less common properties. For instance, we considered the largest excitonic bandwidth that is achievable within known molecular crystals and identified materials with strong super-radiant states. Finally, we explored the possibility that strong excitonic coupling can be used to generate emissive states in the near-infrared region in materials formed by molecules with bright visible absorption and we could identify the maximum allowable red shift in this material class. These insights with the associated searchable database provide practical guidelines for designing materials with interesting optical properties.

Electronic Circular Dichroism of the Cas9 Protein and gRNA:Cas9 Ribonucleoprotein Complex

The Streptococcus pyogenes Cas9 protein (SpCas9), a component of CRISPR-based immune system in microbes, has become commonly utilized for genome editing. This nuclease forms a ribonucleoprotein (RNP) complex with guide RNA (gRNA) which induces Cas9 structural changes and triggers its cleavage activity. Here, electronic circular dichroism (ECD) spectroscopy was used to confirm the RNP formation and to determine its individual components. The ECD spectra had characteristic features differentiating Cas9 and gRNA, the former showed a negative/positive profile with maxima located at 221, 209 and 196 nm, while the latter revealed positive/negative/positive/negative pattern with bands observed at 266, 242, 222 and 209 nm, respectively. For the first time, the experimental ECD spectrum of the gRNA:Cas9 RNP complex is presented. It exhibits a bisignate positive/negative ECD couplet with maxima at 273 and 235 nm, and it differs significantly from individual spectrum of each RNP components. Additionally, the Cas9 protein and RNP complex retained biological activity after ECD measurements and they were able to bind and cleave DNA in vitro. Hence, we conclude that ECD spectroscopy can be considered as a quick and non-destructive method of monitoring conformational changes of the Cas9 protein as a result of Cas9 and gRNA interaction, and identification of the gRNA:Cas9 RNP complex.

Chirality of a rhodamine heterodimer linked to a DNA scaffold: an experimental and computational study

The chiroptical properties of multi-chromophoric systems are governed by the intermolecular arrangement of the monomeric units. We report on a computational and experimental study of the linear optical properties and supramolecular structure of a rhodamine heterodimer assembled on a DNA scaffold. The experimental absorption and circular dichroism (CD) profiles confirm the dimer formation. Computationally, starting from low-cost DFT/TDDFT simulations of the bare dimer we attribute the measured -/+ CD sign sequence of the S1/S2 bands to a specific chiral conformation of the heterodimer. In the monomers, as typical for rhodamine dyes, the electric transition dipole of the lowest π-π* transition is parallel to the long axis of the xanthene planes. We show that in the heterodimer the sign sequence of the two CD bands is related to the orientation of these long axes. To account explicitly for environment effects, we use molecular dynamics (MD) simulations for characterizing the supramolecular structure of the two optical isomers tethered on DNA. Average absorption and CD-profiles were modeled using ab initio TDDFT calculations at the geometries sampled along a few nanosecond MD run. The absorption profiles computed for both optical isomers are in good agreement with the experimental absorption spectrum and do not allow one to discriminate between them. The computed averaged CD profiles provide the orientation of monomers in the enantiomer that is dominant under the experimental conditions.

Studying the excited electronic states of guanine rich DNA quadruplexes by quantum mechanical methods: main achievements and perspectives

Calculations are providing more and more useful insights into the interaction between light and DNA quadruplexes.

Quantum Chemical Modeling of the Photoinduced Activity of Multichromophoric Biosystems

Multichromophoric biosystems represent a broad family with very diverse members, ranging from light-harvesting pigment-protein complexes to nucleic acids. The former are designed to capture, harvest, efficiently transport, and transform energy from sunlight for photosynthesis, while the latter should dissipate the absorbed radiation as quickly as possible to prevent photodamages and corruption of the carried genetic information. Because of the unique electronic and structural characteristics, the modeling of their photoinduced activity is a real challenge. Numerous approaches have been devised building on the theoretical development achieved for single chromophores and on model Hamiltonians that capture the essential features of the system. Still, a question remains: is a general strategy for the accurate modeling of multichromophoric systems possible? By using a quantum chemical point of view, here we review the advancements developed so far highlighting differences and similarities with the single chromophore treatment. Finally, we outline the important limitations and challenges that still need to be tackled to reach a complete and accurate picture of their photoinduced properties and dynamics.

Steady-State Linear and Non-linear Optical Spectroscopy of Organic Chromophores and Bio-macromolecules

Bio-macromolecules as DNA, lipid membranes and (poly)peptides are essential compounds at the core of biological systems. The development of techniques and methodologies for their characterization is therefore necessary and of utmost interest, even though difficulties can be experienced due to their intrinsic complex nature. Among these methods, spectroscopies, relying on optical properties are especially important to determine their macromolecular structures and behaviors, as well as the possible interactions and reactivity with external dyes—often drugs or pollutants—that can (photo)sensitize the bio-macromolecule leading to eventual chemical modifications, thus damages. In this review, we will focus on the theoretical simulation of electronic spectroscopies of bio-macromolecules, considering their secondary structure and including their interaction with different kind of (photo)sensitizers. Namely, absorption, emission and electronic circular dichroism (CD) spectra are calculated and compared with the available experimental data. Non-linear properties will be also taken into account by two-photon absorption, a highly promising technique (i) to enhance absorption in the red and infra-red windows and (ii) to enhance spatial resolution. Methodologically, the implications of using implicit and explicit solvent, coupled to quantum and thermal samplings of the phase space, will be addressed. Especially, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are explored for a comparison with solely QM methods, in order to address the necessity to consider an accurate description of environmental effects on spectroscopic properties of biological systems.

How and How Much Molecular Conformation Affects Electronic Circular Dichroism: The Case of 1,1-Diarylcarbinols

Chiroptical spectra such as electronic circular dichroism (ECD) are said to be much more sensitive to conformation than their non-chiroptical counterparts, however, it is difficult to demonstrate such a common notion in a clear-cut way. We run DFT and TDDFT calculations on two closely related 1,1-diarylmethanols which show mirror-image ECD spectra for the same absolute configuration. We demonstrate that the main reason for the different chiroptical response of the two compounds lies in different conformational ensembles, caused by a single hydrogen-to-methyl substitution. We conclude that two compounds, having the same configuration but different conformation, may exhibit mirror-image ECD signals, stressing the importance and impact of conformational factors on ECD spectra.

Polarization spectroscopy methods in the determination of interactions of small molecules with nucleic acids - tutorial

The structural characterization of non-covalent complexes between nucleic acids and small molecules (ligands) is of a paramount significance to bioorganic research. Highly informative methods about nucleic acid/ligand complexes such as single crystal X-ray diffraction or NMR spectroscopy cannot be performed under biologically compatible conditions and are extensively time consuming. Therefore, in search for faster methods which can be applied to conditions that are at least similar to the naturally occurring ones, a set of polarization spectroscopy methods has shown highly promising results. Electronic circular dichroism (ECD) is the most commonly used method for the characterization of the helical structure of DNA and RNA and their complexes with ligands. Less common but complementary to ECD, is flow-oriented linear dichroism (LD). Other methods such as vibrational CD (VCD) and emission-based methods (FDCD, CPL), can also be used for suitable samples. Despite the popularity of polarization spectroscopy in biophysics, aside several highly focused reviews on the application of these methods to DNA/RNA research, there is no systematic tutorial covering all mentioned methods as a tool for the characterization of adducts between nucleic acids and small ligands. This tutorial aims to help researchers entering the research field to organize experiments accurately and to interpret the obtained data reliably.

EXAT: EXcitonic analysis tool

We introduce EXcitonic Analysis Tool (EXAT), a program able to compute optical spectra of large excitonic systems directly from the output of quantum mechanical calculations performed with the popular Gaussian 16 package. The software is able to combine in an excitonic scheme the single-chromophore properties and exciton couplings to simulate energies, coefficients, and excitonic spectra (UV-vis, CD, and LD). The effect of the environment can also be included using a Polarizable Continuum Model. EXAT also presents a simple graphical user interface, which shows on-screen both site and exciton properties. To show the potential of the method, we report two applications on a a chiral perturbed BODIPY system and DNA G-quadruplexes, respectively. The program is available online at http://molecolab.dcci.unipi.it/tools/. © 2017 Wiley Periodicals, Inc.

The shape of the electronic circular dichroism spectrum of (2,6-dimethylphenyl)(phenyl)methanol: interplay between conformational equilibria and vibronic effects

Analysis of the interplay between conformational equilibria, solvent effects and vibronic contributions in the ECD spectra.