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.

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.

Modeling One-Electron Oxidation Potentials and Hole Delocalization in Double-Stranded DNA by Multilayer and Dynamic Approaches.

The number of innovative applications for DNA nowadays is growing quickly. Its use as a nanowire or electrochemical biosensor leads to the need for a deep understanding of the charge-transfer process along the strand, as well as its redox properties. These features are computationally simulated and analyzed in detail throughout this work by combining molecular dynamics, multilayer schemes, and the Marcus theory. One-electron oxidation potential and hole delocalization have been analyzed for six DNA double strands that cover all possible binary combinations of nucleotides. The results have revealed that the one-electron oxidation potential decreases with respect to the single-stranded DNA, giving evidence that the greater rigidity of a double helix induces an increase in the capacity of storing the positive charge generated upon oxidation. In addition, the hole is mainly stored in nucleobases with large reducer character, i.e., purines, especially when those are arranged in a stacked configuration in the same strand. From the computational point of view, the sampling needed to describe biological systems implies a significant computational cost. Here, we show that a small number of representative conformations generated by clustering analysis provides accurate results when compared with those obtained from sampling, reducing considerably the computational cost.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Computational Photophysics in the Presence of an Environment.

Most processes triggered by ultraviolet (UV) or visible (vis) light in nature take place in complex biological environments. The first step in these photophysical events is the excitation of the absorbing system or chromophore to an electronically excited state. Such an excitation can be monitored by the UV-vis absorption spectrum. A precise calculation of the UV-vis spectrum of a chromophore embedded in an environment is a challenging task that requires the consideration of several ingredients, besides an accurate electronic-structure method for the excited states. Two of the most important are an appropriate description of the interactions between the chromophore and the environment and accounting for the vibrational motion of the whole system. In this contribution, we review the most common theoretical methodologies to describe the environment (including quantum mechanics/continuum and quantum mechanics/molecular mechanics models) and to account for vibrational sampling (including Wigner sampling and molecular dynamics). Further, we illustrate in a series of examples how the lack of these ingredients can lead to a wrong interpretation of the electronic features behind the UV-vis absorption spectrum.

Solvent reorganization triggers photo-induced solvated electron generation in phenol.

The analysis of the absorption spectrum and density of states of a cluster of phenol solvated with 15 water molecules indicates that the reorganization of the water molecules, facilitating the formation of solvated electrons, is a plausible mechanism in the photodissociation of phenol. Using quantitative wavefunction analysis, we demonstrate that while charge-transfer states involving electron transfer from phenol to water are mainly dark, a considerable number of them exists below the maximum of the ππ* absorption band and could be populated by internal conversion. These low-lying charge-transfer states do not show extended O-H distances, but are found for large electron-hole separations at which several water molecules can solvate and stabilize the transferred electron. Thus, charge-transfer states in solvated phenol can be stabilized by two factors: (i) elongation of the O-H bond, as was extensively discussed in the past, and (ii) reorganization of solvent molecules, as it is shown here.

Finite-temperature Wigner phase-space sampling and temperature effects on the excited-state dynamics of 2-nitronaphthalene.

The concept of finite temperature Wigner phase-space sampling allowing the population of vibrationally excited states is introduced and employed to study temperature effects on the absorption spectrum of 2-nitronaphtalene (2NN) and its relaxation dynamics. It is found that, despite the fact that the general deactivation mechanism of 2NN after light irradiation does not change with increasing temperature, i.e., after excitation to the singlet manifold, 2NN deactivates via internal conversion in less than 100 fs and then undergoes intersystem crossing to the triplet manifold before decaying to the lowest triplet state in less than 150 fs, the intersystem crossing rate increases at higher temperatures. This increase is attributed to the appearance of a new deactivation pathway, which is not operative at lower temperatures. The present example illustrates that appropriate initial conditions beyond the idealized temperature of 0 K are indispensable to obtain reliable excited state mechanisms, which can be then related to experimental conditions.

Mechanism of Ultrafast Intersystem Crossing in 2-Nitronaphthalene.

Nitronaphthalene derivatives efficiently populate their electronically excited triplet states upon photoexcitation through ultrafast intersystem crossing (ISC). Despite having been studied extensively by time-resolved spectroscopy, the reasons behind their ultrafast ISC remain unknown. Herein, we present the first ab initio nonadiabatic molecular dynamics study of a nitronaphthalene derivative, 2-nitronaphthalene, including singlet and triplet states. We find that there are two distinct ISC reaction pathways involving different electronic states at distinct nuclear configurations. The high ISC efficiency is explained by the very small electronic and nuclear alterations that the chromophore needs to undergo during the singlet-triplet transition in the dominating ISC pathway after initial dynamics in the singlet manifold. The insights gained in this work are expected to shed new light on the photochemistry of other nitro polycyclic aromatic hydrocarbons that exhibit ultrafast intersystem crossing.