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.

One-Electron Oxidation Potentials and Hole Delocalization in Heterogeneous Single-Stranded DNA.

The study of DNA processes is essential to understand not only its intrinsic biological functions but also its role in many innovative applications. The use of DNA as a nanowire or electrochemical biosensor leads to the need for a deep investigation of the charge transfer process along the strand as well as of the redox properties. In this contribution, the one-electron oxidation potential and the charge delocalization of the hole formed after oxidation are computationally investigated for different heterogeneous single-stranded DNA strands. We have established a two-step protocol: (i) molecular dynamics simulations in the frame of quantum mechanics/molecular mechanics (QM/MM) were performed to sample the conformational space; (ii) energetic properties were then obtained within a QM1/QM2/continuum approach in combination with the Marcus theory over an ensemble of selected geometries. The results reveal that the one-electron oxidation potential in the heterogeneous strands can be seen as a linear combination of that property within the homogeneous strands. In addition, the hole delocalization between different nucleobases is, in general, small, supporting the conclusion of a hopping mechanism for charge transport along the strands. However, charge delocalization becomes more important, and so does the tunneling mechanism contribution, when the reducing power of the nucleobases forming the strand is similar. Moreover, charge delocalization is slightly enhanced when there is a correlation between pairs of some of the interbase coordinates of the strand: twist/shift, twist/slide, shift/slide, and rise/tilt. However, the internal structure of the strand is not the predominant factor for hole delocalization but the specific sequence of nucleotides that compose the strand.

The protein environment restricts the intramolecular charge transfer character of the luciferine/luciferase complex.

The electronic characterization of the luciferine/luciferase complex is fundamental to tune its photophysical properties and develop more efficient devices based on this luminiscent system. Here, we apply molecular dynamics simulations, hybrid quantum mechanics/molecular mechanics (QM/MM) calculations and transition density analysis to compute the absorption and emission spectra of luciferine/luciferase and analyze the nature of the relevant electronic state and its behaviour with the intramolecular and intermolecular degrees of freedom. It is found that the torsional motion of the chromophore is hampered by the presence of the enzyme, reducing the intramolecular charge transfer nature of the absorbing and emitting state. In addition, such a reduced charge transfer character does not correlate in a strong way neither with the intramolecular motion of the chromophore nor with the chromophore/amino-acid distances. However, the presence of a polar environment around the oxygen atom of the thiazole ring of the oxyluciferin, coming from both the protein and the solvent, enhances the charge transfer character of the emitting state.

Intramolecular and intermolecular hole delocalization rules the reducer character of isolated nucleobases and homogeneous single-stranded DNA.

The use of DNA strands as nanowires or electrochemical biosensors requires a deep understanding of charge transfer processes along the strand, as well as of the redox properties. These properties are computationally assessed in detail throughout this study. By applying molecular dynamics and hybrid QM/continuum and QM/QM/continuum schemes, the vertical ionization energies, adiabatic ionization energies, vertical attachment energies, one-electron oxidation potentials, and delocalization of the hole generated upon oxidation have been determined for nucleobases in their free form and as part of a pure single-stranded DNA. We show that the reducer ability of the isolated nucleobases is explained by the intramolecular delocalization of the positively charged hole, while the enhancement of the reducer character when going from aqueous solution to the strand correlates very well with the intermolecular hole delocalization. Our simulations suggest that the redox properties of DNA strands can be tuned by playing with the balance between intramolecular and intermolecular charge delocalization.

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.

An Efficient Multilayer Approach to Model DNA-Based Nanobiosensors.

In this work, we present a full computational protocol to successfully obtain the one-electron reduction potential of nanobiosensors based on a self-assembled monolayer of DNA nucleobases linked to a gold substrate. The model is able to account for conformational sampling and environmental effects at a quantum mechanical (QM) level efficiently, by combining molecular mechanics (MM) molecular dynamics and multilayer QM/MM/continuum calculations within the framework of Marcus theory. The theoretical model shows that a guanine-based biosensor is more prone to be oxidized than the isolated nucleobase in water due to the electrostatic interactions between the assembled guanine molecules. In addition, the redox properties of the biosensor can be tuned by modifying the nature of the linker that anchor the nucleobases to the metal support.

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis.

Herein, an Energy Decomposition Analysis (EDA) scheme extended to the framework of QM/MM calculations in the context of electrostatic embeddings (QM/MM-EDA) including atomic charges and dipoles is applied to assess the effect of the QM region size on the convergence of the different interaction energy components, namely, electrostatic, Pauli, and polarization, for cationic, anionic, and neutral systems interacting with a strong polar environment (water). Significant improvements are found when the bulk solvent environment is described by a MM potential in the EDA scheme as compared to pure QM calculations that neglect bulk solvation. The predominant electrostatic interaction requires sizable QM regions. The results reported here show that it is necessary to include a surprisingly large number of water molecules in the QM region to obtain converged values for this energy term, contrary to most cluster models often employed in the literature. Both the improvement of the QM wave function by means of a larger basis set and the introduction of polarization into the MM region through a polarizable force field do not translate to a faster convergence with the QM region size, but they lead to better results for the different interaction energy components. The results obtained in this work provide insight into the effect of each energy component on the convergence of the solute-solvent interaction energy with the QM region size. This information can be used to improve the MM FFs and embedding schemes employed in QM/MM calculations of solvated systems.

MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations.

We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.

Computation of Oxidation Potentials of Solvated Nucleobases by Static and Dynamic Multilayer Approaches.

The determination of the redox properties of nucleobases is of paramount importance to get insight into the charge-transfer processes in which they are involved, such as those occurring in DNA-inspired biosensors. Although many theoretical and experimental studies have been conducted, the value of the one-electron oxidation potentials of nucleobases is not well-defined. Moreover, the most appropriate theoretical protocol to model the redox properties has not been established yet. In this work, we have implemented and evaluated different static and dynamic approaches to compute the one-electron oxidation potentials of solvated nucleobases. In the static framework, two thermodynamic cycles have been tested to assess their accuracy against the direct determination of oxidation potentials from the adiabatic ionization energies. Then, the introduction of vibrational sampling, the effect of implicit and explicit solvation models, and the application of the Marcus theory have been analyzed through dynamic methods. The results revealed that the static direct determination provides more accurate results than thermodynamic cycles. Moreover, the effect of sampling has not shown to be relevant, and the results are improved within the dynamic framework when the Marcus theory is applied, especially in explicit solvent, with respect to the direct approach. Finally, the presence of different tautomers in water does not affect significantly the one-electron oxidation potentials.

New challenges for photopharmacology and computational modeling.

Mous et al. recently reported the molecular mechanism of chloride transport through a light-activated pumping rhodopsin, a key process involved in a range of cellular functions. Their results open exciting new challenges for photopharmacology and computational modeling that should be addressed in the coming years.

On the Permeation of Polychlorinated Dibenzodioxins and Dibenzofurans through Lipid Membranes: Classical MD and Hybrid QM/MM-EDA Analysis.

The permeation of dioxin-like pollutants, namely, chlorinated dibenzodioxins and dibenzofurans, through lipid membranes has been simulated using classic molecular dynamics (CMD) combined with the umbrella sampling approach. The most toxic forms of chlorinated dibenzodioxin and dibenzofuran, 2,3,7,8-tetrachloro-p-dibenzodioxin (TCDD) and 2,3,7,8-tetrachlorodibenzofuran (TCDF), and a dioleyl-phosphatidylcholine (DOPC) lipid membrane of 50 Å wide have been chosen for our study. The free energy profile shows the penetration process is largely favoured thermodynamically (ΔG ≈ -12 kcal/mol), with a progressively decrease of the free energy until reaching the energy minima at distances of 8 Å and 9.5 Å from the centre of the membrane for, respectively, TCDD and TCDF. At the centre of the membrane, both molecules display subtle local maxima with free energy differences of 0.5 and 1 kcal/mol with respect to the energy minima for TCDD and TCDF, respectively. Furthermore, the intermolecular interactions between the molecules and the lipid membrane have been characterized at the minima and the local maxima using hybrid quantum mechanics/molecular mechanics energy decomposition analysis (QM/MM-EDA). Total interaction energies of -17.5 and -16.5 kcal/mol have been found at the energy minima for TCDD and TCDF, respectively. In both cases, the dispersion forces govern the molecule-membrane interactions, no significant changes have been found at the local maxima, in agreement with the classical free energy profile. The small differences found in the results obtained for TCDD and TCDF point out that the adsorption and diffusion processes through the cell membrane are not related to the different toxicity shown by these pollutants.

Cosolvent and Dynamic Effects in Binding Pocket Search by Docking Simulations.

The lack of conformational sampling in virtual screening projects can lead to inefficient results because many of the potential drugs may not be able to bind to the target protein during the static docking simulations. Here, we performed ensemble docking for around 2000 United States Food and Drug Administration (FDA)-approved drugs with the RNA-dependent RNA polymerase (RdRp) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as a target. The representative protein structures were generated by clustering classical molecular dynamics trajectories, which were evolved using three solvent scenarios, namely, pure water, benzene/water and phenol/water mixtures. The introduction of dynamic effects in the theoretical model showed improvement in docking results in terms of the number of strong binders and binding sites in the protein. Some of the discovered pockets were found only for the cosolvent simulations, where the nonpolar probes induced local conformational changes in the protein that lead to the opening of transient pockets. In addition, the selection of the ligands based on a combination of the binding free energy and binding free energy gap between the best two poses for each ligand provided more suitable binders than the selection of ligands based solely on one of the criteria. The application of cosolvent molecular dynamics to enhance the sampling of the configurational space is expected to improve the efficacy of virtual screening campaigns of future drug discovery projects.

Activation by oxidation and ligand exchange in a molecular manganese vanadium oxide water oxidation catalyst.

Despite their technological importance for water splitting, the reaction mechanisms of most water oxidation catalysts (WOCs) are poorly understood. This paper combines theoretical and experimental methods to reveal mechanistic insights into the reactivity of the highly active molecular manganese vanadium oxide WOC [Mn4V4O17(OAc)3]3- in aqueous acetonitrile solutions. Using density functional theory together with electrochemistry and IR-spectroscopy, we propose a sequential three-step activation mechanism including a one-electron oxidation of the catalyst from [Mn2 3+Mn2 4+] to [Mn3+Mn3 4+], acetate-to-water ligand exchange, and a second one-electron oxidation from [Mn3+Mn3 4+] to [Mn4 4+]. Analysis of several plausible ligand exchange pathways shows that nucleophilic attack of water molecules along the Jahn-Teller axis of the Mn3+ centers leads to significantly lower activation barriers compared with attack at Mn4+ centers. Deprotonation of one water ligand by the leaving acetate group leads to the formation of the activated species [Mn4V4O17(OAc)2(H2O)(OH)]- featuring one H2O and one OH ligand. Redox potentials based on the computed intermediates are in excellent agreement with electrochemical measurements at various solvent compositions. This intricate interplay between redox chemistry and ligand exchange controls the formation of the catalytically active species. These results provide key reactivity information essential to further study bio-inspired molecular WOCs and solid-state manganese oxide catalysts.

Characterization of cisplatin/membrane interactions by QM/MM energy decomposition analysis.

We extend for the first time a quantum mechanical energy decomposition analysis scheme based on deformation electron densities to a hybrid electrostatic embedding quantum mechanics/molecular mechanics framework. The implemented approach is applied to characterize the interactions between cisplatin and a dioleyl-phosphatidylcholine membrane, which play a key role in the permeation mechanism of the drug inside the cells. The interaction energy decomposition into electrostatic, induction, dispersion and Pauli repulsion contributions is performed for ensembles of geometries to account for conformational sampling. It is evidenced that the electrostatic and repulsive components are predominant in both polar and non-polar regions of the bilayer.

The Permeation Mechanism of Cisplatin Through a Dioleoylphosphocholine Bilayer*.

The investigation of the intermolecular interactions between platinum-based anticancer drugs and lipid bilayers is of special relevance to unveil the mechanisms involved in different steps of the anticancer mode of action of these drugs. We have simulated the permeation of cisplatin through a model membrane composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine lipids by means of umbrella sampling classical molecular dynamics simulations. The initial physisorption of cisplatin into the polar region of the lipid membrane is controlled by long-range electrostatic interactions with the choline groups in a first step and, in a second step, by long-range electrostatic and hydrogen bond interactions with the phosphate groups. The second half of the permeation pathway, in which cisplatin diffuses through the nonpolar region of the bilayer, is characterized by the drop of the interactions with the polar heads and the rise of attractive interactions with the non-polar tails, which are dominated by van der Waals contributions. The permeation free-energy profile is explained by a complex balance between the drug/lipid interactions and the energy and entropy contributions associated with the dehydration of the drug along the permeation pathway and with the decrease and increase of the membrane ordering along the first and second half of the mechanism, respectively.

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.

Stacking Effects on Anthraquinone/DNA Charge-Transfer Electronically Excited States.

The design of more efficient photosensitizers is a matter of great importance in the field of cancer treatment by means of photodynamic therapy. One of the main processes involved in the activation of apoptosis in cancer cells is the oxidative stress on DNA once a photosensitizer is excited by light. As a consequence, it is very relevant to investigate in detail the binding modes of the chromophore with DNA, and the nature of the electronically excited states that participate in the induction of DNA damage, for example, charge-transfer states. In this work, we investigate the electronic structure of the anthraquinone photosensitizer intercalated into a double-stranded poly(dG-dC) decamer model of DNA. First, the different geometric configurations are analyzed by means of classical molecular dynamics simulations. Then, the excited states for the most relevant poses of anthraquinone inside the binding pocket are computed by an electrostatic-embedding quantum mechanics/molecular mechanics approach, where anthraquinone and one of the nearby guanine residues are described quantum mechanically to take into account intermolecular charge-transfer states. The excited states are characterized as monomer, exciton, excimer, and charge-transfer states based on the analysis of the transition density matrix, and each of these contributions to the total density of states and absorption spectrum is discussed in terms of the stacking interactions. These results are relevant as they represent the footing for future studies on the reactivity of anthraquinone derivatives with DNA and give insights on possible geometrical configurations that potentially favor the oxidative stress of DNA.

Orbital-free photophysical descriptors to predict directional excitations in metal-based photosensitizers.

The development of dye-sensitized solar cells, metalloenzyme photocatalysis or biological labeling heavily relies on the design of metal-based photosensitizes with directional excitations. Directionality is most often predicted by characterizing the excitations manually via canonical frontier orbitals. Although widespread, this traditional approach is, at the very least, cumbersome and subject to personal bias, as well as limited in many cases. Here, we demonstrate how two orbital-free photophysical descriptors allow an easy and straightforward quantification of the degree of directionality in electron excitations using chemical fragments. As proof of concept we scrutinize the effect of 22 chemical modifications on the archetype [Ru(bpy)3]2+ with a new descriptor coined "substituent-induced exciton localization" (SIEL), together with the concept of "excited-electron delocalization length" (EEDL n ). Applied to quantum ensembles of initially excited singlet and the relaxed triplet metal-to-ligand charge-transfer states, the SIEL descriptor allows quantifying how much and whereto the exciton is promoted, as well as anticipating the effect of single modifications, e.g. on C-4 atoms of bpy units of [Ru(bpy)3]2+. The general applicability of SIEL and EEDL n is further established by rationalizing experimental trends through quantification of the directionality of the photoexcitation. We thus demonstrate that SIEL and EEDL descriptors can be synergistically employed to design improved photosensitizers with highly directional and localized electron-transfer transitions.

Excimer Intermediates en Route to Long-Lived Charge-Transfer States in Single-Stranded Adenine DNA as Revealed by Nonadiabatic Dynamics.

The ultrafast time evolution of a single-stranded adenine DNA is studied using a hybrid multiscale quantum mechanics/molecular mechanics (QM/MM) scheme coupled to nonadiabatic surface hopping dynamics. As a model, we use (dA)20 where a stacked adenine tetramer is treated quantum chemically. The dynamical simulations combined with on-the-fly quantitative wave function analysis evidence the nature of the long-lived electronically excited states formed upon absorption of UV light. After a rapid decrease of the initially excited excitons, relaxation to monomer-like states and excimers occurs within 100 fs. The former monomeric states then relax into additional excimer states en route to forming stabilized charge-transfer states on a longer timescale of hundreds of femtoseconds. The different electronic-state characters is reflected on the spatial separation between the adenines: excimers and charge-transfer states show a much smaller spatial separation than the monomer-like states and the initially formed excitons.

Biological evaluation of novel thiomaltol-based organometallic complexes as topoisomerase IIα inhibitors.

Topoisomerase IIα (topo2α) is an essential nuclear enzyme involved in DNA replication, transcription, recombination, chromosome condensation, and highly expressed in many tumors. Thus, topo2α-targeting has become a very efficient and well-established anticancer strategy. Herein, we investigate the cytotoxic and DNA-damaging activity of thiomaltol-containing ruthenium-, osmium-, rhodium- and iridium-based organometallic complexes in human mammary carcinoma cell lines by means of several biological assays, including knockdown of topo2α expression levels by RNA interference. Results suggest that inhibition of topo2α is a key process in the cytotoxic mechanism for some of the compounds, whereas direct induction of DNA double-strand breaks or other DNA damage is mostly rather minor. In addition, molecular modeling studies performed for two of the compounds (with Ru(II) as the metal center) evinces that these complexes are able to access the DNA-binding pocket of the enzyme, where the hydrophilic environment favors the interaction with highly polar complexes. These findings substantiate the potential of these compounds for application as antitumor metallopharmaceuticals.