Membrane Permeation of Psychedelic Tryptamines by Dynamic Simulations.

Renewed scientific interest in psychedelic compounds represents one of the most promising avenues for addressing the current burden of mental health disorders. Classic psychedelics are a group of compounds that exhibit structural similarities to the naturally occurring neurotransmitter serotonin (5-HT). Acting on the 5-HT type 2A receptors (HT2ARs), psychedelics induce enduring neurophysiological changes that parallel their therapeutic psychological and behavioral effects. Recent preclinical evidence suggests that the ability of psychedelics to exert their action is determined by their ability to permeate the neuronal membrane to target a pool of intracellular 5-HT2ARs. In this computational study, we employ classical molecular dynamics simulations and umbrella sampling techniques to investigate the permeation behavior of 12 selected tryptamines and to characterize the interactions that drive the process. We aim at elucidating the impact of N-alkylation, indole ring substitution and positional modifications, and protonation on their membrane permeability. Dimethylation of the primary amine group and the introduction of a methoxy group at position 5 exhibited an increase in permeability. Moreover, there is a significant influence of positional substitutions on the indole groups, and the protonation of the molecules substantially increases the energy barrier at the center of the bilayer, making the compounds highly impermeable. All the information extracted from the trends predicted by the simulations can be applied in future drug design projects to develop psychedelics with enhanced activity.

Voltage-Gated Ion Channels: Structure, Pharmacology and Photopharmacology.

Voltage-gated ion channels are transmembrane proteins responsible for the generation and propagation of action potentials in excitable cells. Over the last decade, advancements have enabled the elucidation of crystal structures of ion channels. This progress in structural understanding, particularly in identifying the binding sites of local anesthetics, opens avenues for the design of novel compounds capable of modulating ion conduction. However, many traditional drugs lack selectivity and come with adverse side effects. The emergence of photopharmacology has provided an orthogonal way of controlling the activity of compounds, enabling the regulation of ion conduction with light. In this review, we explore the central pore region of voltage-gated sodium and potassium channels, providing insights from both structural and pharmacological perspectives. We discuss the different binding modes of synthetic compounds that can physically occlude the pore and, therefore, block ion conduction. Moreover, we examine recent advances in the photopharmacology of voltage-gated ion channels, introducing molecular approaches aimed at controlling their activity by using photosensitive drugs.

Multiscale Simulations of the Covalent Inhibition of the SARS-CoV-2 Main Protease: Four Compounds and Three Reaction Mechanisms.

We report the results of computational modeling of the reactions of the SARS-CoV-2 main protease (MPro) with four potential covalent inhibitors. Two of them, carmofur and nirmatrelvir, have shown experimentally the ability to inhibit MPro. Two other compounds, X77A and X77C, were designed computationally in this work. They were derived from the structure of X77, a non-covalent inhibitor forming a tight surface complex with MPro. We modified the X77 structure by introducing warheads capable of reacting with the catalytic cysteine residue in the MPro active site. The reaction mechanisms of the four molecules with MPro were investigated by quantum mechanics/molecular mechanics (QM/MM) simulations. The results show that all four compounds form covalent adducts with the catalytic cysteine Cys 145 of MPro. From the chemical perspective, the reactions of these four molecules with MPro follow three distinct mechanisms. The reactions are initiated by a nucleophilic attack of the thiolate group of the deprotonated cysteine residue from the catalytic dyad Cys145-His41 of MPro. In the case of carmofur and X77A, the covalent binding of the thiolate to the ligand is accompanied by the formation of the fluoro-uracil leaving group. The reaction with X77C follows the nucleophilic aromatic substitution SNAr mechanism. The reaction of MPro with nirmatrelvir (which has a reactive nitrile group) leads to the formation of a covalent thioimidate adduct with the thiolate of the Cys145 residue in the enzyme active site. Our results contribute to the ongoing search for efficient inhibitors of the SARS-CoV-2 enzymes.

Effect of stacking interactions on charge transfer states in photoswitches interacting with ion channels.

The activity of ion channels can be reversibly photo-controlled via the binding of molecular photoswitches, often based on an azobenzene scaffold. Those azobenzene derivatives interact with aromatic residues of the protein via stacking interactions. In the present work, the effect of face-to-face and t-shaped stacking interactions on the excited state electronic structure of azobenzene and p-diaminoazobenzene integrated into the NaV1.4 channel is computationally investigated. The formation of a charge transfer state, caused by electron transfer from the protein to the photoswitches, is observed. This state is strongly red shifted when the interaction takes place in a face-to-face orientation and electron donating groups are present on the aromatic ring of the amino acids. The low-energy charge transfer state can interfere with the photoisomerization process after excitation to the bright state by leading to the formation of radical species.

Nature of Hops, Coordinates, and Detailed Balance for Nonadiabatic Simulations in the Condensed Phase.

Photoinduced processes play a crucial role in a multitude of important molecular phenomena. Accurately modeling these processes in an environment other than a vacuum requires a detailed description of the electronic states involved as well as how energy flows are coupled to the surroundings. Nonadiabatic effects must also be included in order to describe the exchange of energy between electronic and nuclear degrees of freedom correctly. In this work, we revisit the ring-opening reaction 1,3-cylohexadiene (CHD) in a solvent environment. Using our newly developed Interface for Non-Adiabatic Quantum mechanics/molecular mechanics in Solvent (INAQS) we trace the evolution of the reaction via hybrid quantum mechanics/molecular mechanics (QM/MM) surface hopping with a focus on the solvent's participation in the nonadiabatic relaxation process and the long-time approach to equilibrium. We explicitly include the MM solvent contribution to the nonadiabatic coupling vector─enabling an accurate approach to equilibrium at long times─and find that in highly multidimensional systems gradients can have little or nothing to do with the nonadiabatic couplings.

Chiral molecule candidates for trapped ion spectroscopy by ab initio calculations: From state preparation to parity violation.

Parity non-conservation (PNC) due to the weak interaction is predicted to give rise to enantiomer dependent vibrational constants in chiral molecules, but the phenomenon has so far eluded experimental observation. The enhanced sensitivity of molecules to physics beyond the Standard Model (BSM) has led to substantial advances in molecular precision spectroscopy, and these may be applied to PNC searches as well. Specifically, trapped molecular ion experiments leverage the universality of trapping charged particles to optimize the molecular ion species studied toward BSM searches, but in searches for PNC, only a few chiral molecular ion candidates have been proposed so far. Importantly, viable candidates need to be internally cold, and their internal state populations should be detectable with high quantum efficiency. To this end, we focus on molecular ions that can be created by near threshold resonant two-photon ionization and detected via state-selective photo-dissociation. Such candidates need to be stable in both charged and neutral chiral versions to be amenable to these methods. Here, we present a collection of suitable chiral molecular ion candidates we have found, including CHDBrI+ and CHCaBrI+, that fulfill these conditions according to our ab initio calculations. We find that organo-metallic species have low ionization energy as neutrals and relatively high dissociation thresholds. Finally, we compute the magnitude of the PNC values for vibrational transitions for some of these candidates. An experimental demonstration of state preparation and readout for these candidates will be an important milestone toward measuring PNC in chiral molecules for the first time.

Providing theoretical insight into the role of symmetry in the photoisomerization mechanism of a non-symmetric dithienylethene photoswitch.

Dithienylethene (DTE) molecular photoswitches have shown to be excellent candidates in the design of efficient optoelectronic devices, due to their high photoisomerization quantum yield (QY), for which symmetry is suggested to play a crucial role. Here, we present a theoretical study on the photochemistry of a non-symmetric dithienylethene photoswitch, with a special emphasis on the effect of asymmetric substitution on the photocyclization and photoreversion mechanisms. We used the Spin-Flip Time Dependent Density Functional Theory (SF-TDDFT) method to locate and characterize the main structures (conical intersections and minima) of the ground state and the first two excited states, S1 and S2, along the ring-opening/closure reaction coordinate of the photocyclization and photoreversion processes, and to identify the important coordinates governing the radiationless decay pathways. Our results suggest that while the main features that characterize the photoisomerization of symmetric DTEs are also present for the photoisomerization of the non-symmetric DTE, the lower energy barrier on S1 along the cycloreversion reaction speaks in favor of a more efficient and therefore a higher cycloreversion QY for the non-symmetric DTEs, making them a better candidate for molecular optoelectronic devices than their symmetric counterparts.

Intramolecular hydrogen transfer in DNA induced by site-selective resonant core excitation.

We present experimental evidence for soft X-ray induced intramolecular hydrogen transfer in the protonated synthetic tri-oligonucleotide d(FUAG) in the gas-phase (FU: fluorouracil). The trinucleotide cations were stored in a cryogenic ion trap and exposed to monochromatic synchrotron radiation. Photoionization and photofragmentation product ion yields were recorded as a function of photon energy. Predominanly glycosidic bond cleavage leading to formation of nucleobase-related fragments is observed. In most cases, glycosidic bond cleavage is accompanied by single or double hydrogen transfer. The combination of absorption-site-sensitive soft X-ray spectroscopy with fragment specific mass spectrometry allows to directly relate X-ray absorption site and fragmentation site. We observe pronounced resonant features in the competition between single and double hydrogen transfer towards nucleobases. A direct comparison of experimental data with time-dependent density functional theory calculations, using short range corrected hybrid functionals, reveal that these hydrogen transfer processes are universal and not limited to population of particular excited states localized at the nucleobases. Instead, hydrogen transfer can occur upon X-ray absorption in any nucleobase and in the DNA backbone. Resonances seem to occur because of site-selective suppression of hydrogen transfer channels. Furthermore, non-covalent interactions of the optimized ground state geometries were investigated to identify intramolecular hydrogen bonds along which hydrogen transfer is most likely.

Influence of the crystal packing in singlet fission: one step beyond the gas phase approximation.

Singlet fission (SF), a multiexciton generation process, has been proposed as an alternative to enhance the performance of solar cells. The gas phase dimer model has shown its utility to study this process, but it does not always cover all the physics and the effect of the surrounding atoms has to be included in such cases. In this contribution, we explore the influence of crystal packing on the electronic couplings, and on the so-called exciton descriptors and electron-hole correlation plots. We have studied three tetracene dimers extracted from the crystal structure, as well as several dimers and trimers of the α and ß polymorphs of 1,3-diphenylisobenzofuran (DPBF). These polymorphs show different SF yields. Our results highlight that the character of the excited states of tetracene depends on both the mutual disposition of molecules and inclusion of the environment. The latter does however not change significantly the interpretation of the SF mechanism in the studied systems. For DPBF, we establish how the excited state analysis is able to pinpoint differences between the polymorphs. We observe strongly bound correlated excitons in the ß polymorph which might hinder the formation of the 1TT state and, consequently, explain its low SF yield.

Binding of azobenzene and p-diaminoazobenzene to the human voltage-gated sodium channel Nav1.4.

The activity of voltage-gated ion channels can be controlled by the binding of photoswitches inside their internal cavity and subsequent light irradiation. We investigated the binding of azobenzene and p-diaminoazobenzene to the human Nav1.4 channel in the inactivated state by means of Gaussian accelerated molecular dynamics simulations and free-energy computations. Three stable binding pockets were identified for each of the two photoswitches. In all the cases, the binding is controlled by the balance between the favorable hydrophobic interactions of the ligands with the nonpolar residues of the protein and the unfavorable polar solvation energy. In addition, electrostatic interactions between the ligand and the polar aminoacids are also relevant for p-diaminoazobenzene due to the presence of the amino groups on the benzene moieties. These groups participate in hydrogen bonding in the most favorable binding pocket and in long-range electrostatic interactions in the other pockets. The thermodinamically preferred binding sites found for both photoswitches are close to the selectivity filter of the channel. Therefore, it is very likely that the binding of these ligands will induce alterations in the ion conduction through the channel.

Site-selective soft X-ray absorption as a tool to study protonation and electronic structure of gas-phase DNA.

The conformation and the electronic structure of gas-phase oligonucleotides depends strongly on the protonation site. 5'-d(FUAG) can either be protonated at the A-N1 or at the G-N7 position. We have stored protonated 5'-d(FUAG) cations in a cryogenic ion trap held at about 20 K. To identify the protonation site and the corresponding electronic structure, we have employed soft X-ray absorption spectroscopy at the nitrogen K-edge. The obtained spectra were interpreted by comparison to time-dependent density functional theory calculations using a short-range exchange correlation functional. Despite the fact that guanine has a significantly higher proton affinity than adenine, the agreement between experiment and theory is better for the A-N1 protonated system. Furthermore, an inverse site sensitivity is observed in which the yield of the nucleobase fragments that contain the absorption site appears substantially reduced, which could be explained by non-statistical fragmentation processes, localized on the photoabsorbing nucleobase.

Theoretical insights into the effect of size and substitution patterns of azobenzene derivatives on the DNA G-quadruplex.

Introducing photoswitches into the DNA G-quadruplex provides excellent opportunities to control folding and unfolding of these assemblies, demonstrating their potential in the development of novel nanodevices with medical and nanotechnology applications. Using a quantum mechanics/molecular mechanics (QM/MM) scheme, we carried out a series of simulations to identify the effect of the size and substitution patterns of three azobenzene derivatives (AZ1, AZ2 and AZ3) on the excitation energies of the two lowest excited states of the smallest photoswitchable G-quadruplex reported to date. We demonstrated that the size and the substitution pattern do not affect the ultrafast cis-trans photoiomerization mechanism of the azobenzene derivatives significantly, in agreement with the experiment. However, molecular dynamics simulations revealed that while AZ2 and AZ3 G-quadruplexes are structurally stable during the simulations, the AZ1 G-quadruplex undergoes larger structural changes and shows two ground state populations that differ in the azobenzene backbone adopting two different conformations. AZ1, with para-para substitution pattern, provides more flexibility to the whole G-quadruplex structure compared to AZ2 and AZ3, and can thus facilitate the photoisomerization reaction between a nonpolymorphic, stacked, tetramolecular G-quadruplex and an unstructured state after trans-cis isomerization occurring in a longer time dynamics, in agreement with the experimental findings. The QM/MM simulations of the absorption spectra indicated that the thermal fluctuation plays a more crucial role in the main absorption band of the azobenzene derivatives than the inclusion of the G-quadruplex, implying that the influence of the G-quadruplex environment is minimal. We propose that the latter is attributed to the position of the azobenzene linkers in the G-quadruplexes, i.e. the edgewise loops containing the azobenzene moieties that are located above the G-quartets, not being fully embedded inside or involved in the stacked structure. Our theoretical findings provide support to a recent study of the photoresponsive formation of photoswitchable G-quadruplex motifs.

Algebraic diagrammatic construction for the polarisation propagator in combination with effective fragment potentials.

The effective fragment potential (EFP) method for the efficient inclusion of solvation effects is combined with the algebraic diagrammatic construction (ADC) scheme for the second- and third-order polarisation propagator. The accuracy of these newly developed EFP-ADC(2) and EFP-ADC(3) methods is tested with respect to supermolecular ADC calculations for a selected set of small solute·solvent complexes. The EFP model for solvation introduces only marginal errors in the excitation energies and oscillator strengths of singlet as well as triplet states, which are strictly localized on the chromophore, significantly below the intrinsic errors of the parent ADC(2) and ADC(3) methods. It is only when delocalization of electron density on the solvent molecules occurs that the error in the excitation energies increases, a well-known behavior of environment models in general. Overall, EFP-ADC schemes prove to be reliable computational approaches to simulate electronic absorption spectra in solution.

Structure and photophysics of indigoids for singlet fission: Cibalackrot.

We report an investigation of structure and photophysics of thin layers of cibalackrot, a sturdy dye derived from indigo by double annulation at the central double bond. Evaporated layers contain up to three phases, two crystalline and one amorphous. Relative amounts of all three have been determined by a combination of X-ray diffraction and FT-IR reflectance spectroscopy. Initially, excited singlet state rapidly produces a high yield of a transient intermediate whose spectral properties are compatible with charge-transfer nature. This intermediate more slowly converts to a significant yield of triplet, which, however, does not exceed 100% and may well be produced by intersystem crossing rather than singlet fission. The yields were determined by transient absorption spectroscopy and corrected for effects of partial sample alignment by a simple generally applicable procedure. Formation of excimers was also observed. In order to obtain guidance for improving molecular packing by a minor structural modification, calculations by a simplified frontier orbital method were used to find all local maxima of singlet fission rate as a function of geometry of a molecular pair. The method was tested at 48 maxima by comparison with the ab initio Frenkel-Davydov exciton model.

Highly Strained Organogold Complexes and Their Gold- or Rhodium-Catalyzed Isomerizations.

The synthesis and characterization of (NHC)(tricyclo[4.1.0.02,7 ]hept-1-yl)gold(I) compounds is reported. Several of these organometallic compounds could be obtained by transmetalation from the organolithium precursor to the corresponding (NHC)gold(I) chlorides. These are by far the most strained stable organogold compounds reported so far. With these new compounds we were able to efficiently and selectively accomplish the first gold-catalyzed rearrangement of organogold compounds, even better yields were obtained by the first rhodium-catalyzed rearrangement of these organogold species. (E/Z)-Mixtures of (NHC)(cyclohex-2-en-1-ylidenemethyl)gold(I) were obtained, none of the other known rearrangement products of tricycloheptanes. Additional insight into the chemical and physical properties of the strained compounds was obtained by spectroscopic and computational studies.

Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework: Theory, implementation, and validation against multi-reference methods.

We report an implementation of non-adiabatic coupling (NAC) forces within the equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) framework via the summed-state approach. Using illustrative examples, we compare NAC forces computed with EOM-CCSD and multi-reference (MR) wave functions (for selected cases, we also consider configuration interaction singles). In addition to the magnitude of the NAC vectors, we analyze their direction, which is important for the calculations of the rate of non-adiabatic transitions. Our benchmark set comprises three doublet radical-cations (hexatriene, cyclohexadiene, and uracil), neutral uracil, and sodium-doped ammonia clusters. When the characters of the states agree among different methods, we observe good agreement between the respective NAC vectors, both in the Franck-Condon region and away. In the cases of large discrepancies between the methods, the disagreement can be attributed to the difference in the states' character, which, in some cases, is very sensitive to electron correlation, both within single-reference and multi-reference frameworks. The numeric results confirm that the accuracy of NAC vectors depends critically on the quality of the underlying wave functions. Within their domain of applicability, EOM-CC methods provide a viable alternative to MR approaches.

Characterization of the Intermediate in and Identification of the Repair Mechanism of (6-4) Photolesions by Photolyases.

Quantum mechanics/molecular mechanics calculations are employed to assign previously recorded experimental spectroscopic signatures of the intermediates occurring during the photo-induced repair of (6-4) photolesions by photolyases to specific molecular structures. Based on this close comparison of experiment and theory it is demonstrated that the acting repair mechanism involves proton transfer from the protonated His365 to the N3' nitrogen of the lesion, which proceeds simultaneously with intramolecular OH transfer along an oxetane-like transition state.

Physicochemical mechanism of light-driven DNA repair by (6-4) photolyases.

DNA photolyases are light-activated enzymes that repair DNA damage induced by ultraviolet (UV) radiation. UV radiation causes two of the most abundant mutagenic and cytotoxic DNA lesions: cyclobutane pyrimidine dimers and 6-4 photolesions. Photolyases selectively bind to DNA and initiate the splitting of mutagenic pyrimidine dimers via photoinduced electron transfer from a flavin adenine dinucleotide anion (FADH(-)) to the lesion triggering its repair. This review discusses the consecutive steps of the repair process, from both experimental and theoretical points of view. It covers the following issues: the process of how photolyases accommodate the lesion into their binding pockets, excitation energy transfer between two involved catalytic cofactors, photoinduced electron transfer to the lesion, the splitting of the pyrimidine dimer radical anion, and the fate of the unstable radical species created after the splitting of the thymine dimer. In particular, mechanisms of the splitting and restoration of the original bases are described in detail, and the most probable repair pathways are outlined.

Emission turn-on and solubility turn-off in conjugated polymers: one- and two-photon-induced removal of fluorescence-quenching solubilizing groups.

The synthesis of highly efficient two-photon uncaging groups and their potential use in functional conjugated polymers for post-polymerization modification are reported. Careful structural design of the employed nitrophenethyl caging groups allows to efficiently induce bond scission by a two-photon process through a combination of exceptionally high two-photon absorption cross-sections and high reaction quantum yields. Furthermore, π-conjugated polyfluorenes are functionalized with these photocleavable side groups and it is possible to alter their emission properties and solubility behavior by simple light irradiation. Cleavage of side groups leads to a turn-on of the fluorescence while solubility of the π-conjugated materials is drastically reduced.

Ultrafast Photoinduced Dynamics in 1,3-Cyclohexadiene: A Comparison of Trajectory Surface Hopping Schemes†.

Photoinduced nonadiabatic processes play a crucial role in a wide range of disciplines, from fundamental steps in biology to modern applications in advanced materials science. A theoretical understanding of these processes is highly desirable, and trajectory surface hopping (TSH) has proven to be a well-suited framework for a wide range of systems. In this work, we present a comprehensive comparison between two TSH algorithms, the conventional Tully's fewest switches surface hopping (FSSH) scheme and the Landau-Zener surface hopping (LZSH), to study the photoinduced ring-opening of 1,3-cyclohexadiene (CHD) to 1,3,5-hexatriene at the spin-flip time-dependent density functional theory (SF-TDDFT) level of theory. Additionally, we compare our results with a literature study at the extended multistate complete active space second-order perturbation theory method (XMS-CASPT2) level of theory. Our results show that the average population and lifetimes estimated with LZSH using SF-TDDFT are closer to the literature (using multireference methods) than those estimated with FSSH using SF-TDDFT. The latter speaks in favor of applying LZSH in combination with the SF-TDDFT method to study larger and more complex systems such as molecular photoswitches where the CHD molecule acts as a backbone. In addition, we present an implementation of Tully's FSSH algorithm as an extension to the PySurf software package.