Search

Photoswitches are becoming increasingly popular in pharmacology due to the possibility of modifying their activity with light. Hence, it is crucial to understand the photophysics of these compounds to identify promising light-activated drugs. We focused our study on DAD, an azobenzene derivative that, according to a previous experimental investigation, can restore visual function in blind mice due to trans - cis photoisomerization upon light absorption. With the present computational study, we aim to characterize the absorption spectrum of DAD, and to understand its photoisomerization mechanism by means of conformational search analysis, quantum mechanical (QM) and hybrid QM/continuum calculations, and classical molecular dynamics simulations. Moreover, we explored the effect of the derivation (DAD vs azobenzene), the protonation (DAD vs DADH 2 2+ , the two possible protonation states) and the solvation (vacuum vs water) on the photoisomerization. Similarly to azobenzene, we showed that the photoisomerization of both protonation states of DAD begin with the population of the bright S 2 state. Then, it crosses to the S 1 surface and relaxes along the rotation of the azo dihedral to a S 1 / S 0 crossing point. The latter is close to a transition state that connects the trans and cis geometries on the ground state. Finally, our results suggested that amino derivation, nonprotonation and water solvation could improve the quantum yield of the photoisomerization.

The modeling of the bioluminescent system of fireflies is key to understand the binding mode of the oxyluciferin/luciferase complex and its photophysical properties with the aim of developing high-efficiency devices and techniques. In this work, we present the OLUFF force field, which is able to describe the interactions to sample the conformational space of the four possible oxyluciferin emitters in both ground and excited state. This force field has been parameterized to reproduce quantum mechanical (QM) energies and geometrical parameters. Moreover, it has been validated by comparing probability distribution functions, minimized structures, infrared spectra and normal mode analysis obtained from OLUFF-based molecular dynamic (MD) simulations with their QM counterparts. Additionally, ground state simulations have also been performed using the general amber force field (GAFF) and compared with the OLUFF. It has been demonstrated that the OLUFF not only reproduces well the QM properties, but also improves the results from the GAFF., Competing Interests: Declarations Conflict of interest The authors declare no conflict of interest., (© 2024. The Author(s), under exclusive licence to the European Photochemistry Association, European Society for Photobiology.)

The interactions between π-systems in dimers of aromatic molecules lead to particularly stable conformations within the relative orientations of the monomers. Extensive research has been conducted on the properties of these complexes in the neutral state. However, in recent decades, there has been a significant surge in applications harnessing these structures for electrical purposes. Therefore, this study places particular emphasis on a deeper understanding of the redox properties of these compounds and how to modify them. To achieve this, we have focused on modeling the effect of a wide range of functional groups on the redox properties of benzene derivatives, observing a correlation between these properties and the change in the molecular dipole moment. Then, we investigated the effect of π-stacking interactions on these properties in dimers formed by either identical or different monomers. In both cases, there is an enhancement of the reducing character of the systems due to these interactions. Upon oxidation, the charge is distributed proportionally to the redox potential of each monomer. Therefore, if there is heterogeneity in these potentials, the properties of the complete cationic system will be influenced by the monomer with a greater tendency to undergo oxidation. The considered models serve as an excellent example for studying the behavior of nucleobases in DNA or aromatic amino acids, among others., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We report the measurement of impulsive stimulated x-ray Raman scattering in neutral liquid water. An attosecond pulse drives the excitations of an electronic wavepacket in water molecules. The process comprises two steps: a transition to core-excited states near the oxygen atoms accompanied by transition to valence-excited states. Thus, the wavepacket is impulsively created at a specific atomic site within a few hundred attoseconds through a nonlinear interaction between the water and the x-ray pulse. We observe this nonlinear signature in an intensity-dependent Stokes Raman sideband at 526 eV. Our measurements are supported by our state-of-the-art calculations based on the polarization response of water dimers in bulk solvation and propagation of attosecond x-ray pulses at liquid density.

G-quadruplexes (G4s) are important therapeutic and photopharmacological targets in cancer research. Small-molecule ligands targeting G4s offer a promising strategy to block DNA transactions and induce genetic instability in cancer cells. While numerous G4-ligands have been reported, relatively few examples exist of compounds whose G4-interactive binding properties can be modulated using light. Herein, we report the photophysical characterization of a novel ortho -fluoroazobenzene derivative, Py-Azo4F-3N , that undergoes reversible two-way isomerization upon visible light exposure. Using a combination of biophysical techniques, including affinity and selectivity assays, structural and computational analysis, and cytotoxicity experiments in cancer cell lines, we carefully characterized the G4-interactive binding properties of both isomers. We identify the trans isomer as the most promising form of interacting and stabilizing G4s, enhancing their ablation capability in cancer cells. Our research highlights the importance of light-responsive molecules in achieving precise control over G4 structures, demonstrating their potential in innovative anticancer strategies.

G-quadruplex (G4) DNA structures are increasingly acknowledged as promising targets in cancer research, and the development of G4-specific stabilizing compounds may lay a fundamental foundation in precision medicine for cancer treatment. Here, we propose a light-responsive G4-binder for precise modulation of drug activation, providing dynamic and spatiotemporal control over G4-associated biological processes contributing to cancer cell death. We developed a specialized fluorinated azobenzene (AB) switch equipped with a quinoline unit and a positively charged carboxamide side chain, Q-Azo4F-C, designed for targeted binding to G4 structures within cells. Biophysical studies, combined with molecular dynamics simulations, provide insights into the unique coordination modes of the photoswitchable ligand in its trans and cis configurations when interacting with G4s. The observed variations in complexation processes between the two isomeric states in different cancer cell lines manifest in more than 25-fold reversible cytotoxic activity. Immunostaining conducted with the structure-specific G4 antibody (BG4), establishes a direct correlation between cytotoxicity and the varying extent of G4 induction regulated by the two isoforms. Finally, we demonstrate the photo-driven reversible regulation of G4 structures in lung cancer cells by Q-Azo4F-C. Our findings highlight the potential of light-responsive G4-binders in advancing precision cancer therapy through dynamic control of G4-mediated pathways., (© 2024 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

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., (© 2024 The Authors. ChemPhysChem published by Wiley-VCH GmbH.)

In this work, we have theoretically determined the one-electron oxidation potentials and charge transfer mechanisms in complex systems based on a self-assembled monolayer of guanine molecules adsorbed on a gold surface through different organic linkers. Classical molecular dynamics simulations were carried out to sample the conformational space of both the neutral and the cationic species. Thus, the redox potentials were determined for the ensembles of geometries through multiscale quantum-mechanics/molecular-mechanics/continuum solvation model calculations in the framework of the Marcus theory and in combination with an additive scheme previously developed. In this context, conformational sampling, description of the environment, and effects caused by the linker have been considered. Applying this methodology, we unravel the phenomena of electric current transport by evaluating the different stages in which charge transfer could occur. The results revealed how the positive charge migrates from the organic layer to the gold surface. Specifically, the transport mechanism seems to take place mainly along a single ligand and driven with the help of the electrostatic interactions of the surrounding molecules. Aside, several self-assembled monolayers with different linkers have been analyzed to understand how the nature of that moiety can tune the redox properties and the efficiency of the transport. We have found that the conjugation between the guanine and the linker, at the same time conjugated to the gold surface, gives rise to a more efficient transport. In conclusion, the established computational protocol sheds light on the mechanism behind charge transport in electrochemical DNA-based biosensor nanodevices.

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.

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 (HT 2A Rs), 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-HT 2A Rs. 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.

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 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.

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.

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 Na V 1.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.

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.

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.

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., (© 2022 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

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.

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., Competing Interests: Declaration of interests None declared by author., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

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 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.

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.

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 [Mn 4 V 4 O 17 (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 [Mn 2 3+ Mn 2 4+ ] to [Mn 3+ Mn 3 4+ ], acetate-to-water ligand exchange, and a second one-electron oxidation from [Mn 3+ Mn 3 4+ ] to [Mn 4 4+ ]. Analysis of several plausible ligand exchange pathways shows that nucleophilic attack of water molecules along the Jahn-Teller axis of the Mn 3+ centers leads to significantly lower activation barriers compared with attack at Mn 4+ centers. Deprotonation of one water ligand by the leaving acetate group leads to the formation of the activated species [Mn 4 V 4 O 17 (OAc) 2 (H 2 O)(OH)] - featuring one H 2 O 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., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

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., (© 2021 Wiley-VCH GmbH.)

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 Na v 1.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.

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.

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.

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., (This journal is © The Royal Society of Chemistry 2020.)

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.

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.

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.

The absorption properties of Temoporfin, a second-generation photosensitizer employed in photodynamic therapy, are calculated with an electrostatic-embedding quantum mechanics/molecular mechanics (QM/MM) scheme in methanol. The suitability of several ensembles of geometries generated by different sampling techniques, namely classical-molecular-dynamics (MD) and QM/MM-MD thermal sampling, Wigner quantum sampling and a hybrid protocol, which combines the thermal and quantum approaches, is assessed. It is found that a QM description of the chromophore during the sampling is needed in order to achieve a good agreement with respect to the experimental spectrum. Such a good agreement is obtained with both QM/MM-MD and Wigner samplings, demonstrating that a proper description of the anharmonic motions of the chromophore is not relevant in the computation of the absorption properties. In addition, it is also found that solvent organization is a rather fast process and a long sampling is not required. Finally, it is also demonstrated that the same exchange-correlation functional should be employed in the sampling and in the computation of the excited states properties to avoid unphysical triplet states with relative energies close or below 0 eV.

Interactions between chromophores and biological environments may alter the electronic properties of the chromophores. A three-layered QM/QM/MM ONIOM scheme with electrostatic embedding is implemented to investigate the influence of an additional QM layer on excited-state calculations with respect to a standard QM/MM description. The implemented ONIOM scheme is employed to compute the electronic excitations of the photosensitizer methylene blue interacting with a solvated DNA double strand. It is shown that the additional quantum mechanical description of several nucleobases in the vertical energy calculations induces energy shifts in the excited states of methylene blue, compared to the energies of a traditional QM/MM scheme, where the solvated double strand is described fully classically. The energy shifts present an electrostatic component, caused by a charge redistribution of the environment, and an electronic-coupling component, originated by the mixing between the electronic bright state of methylene blue and a charge-transfer state between methylene blue and guanine. In addition, hydrogen bonding and stacking interactions are stronger when the environment is described quantum mechanically during the geometry optimizations than when it is fully described by molecular mechanics. These larger intermolecular interactions produce further energy shifts in the excitation energies of the photosensitizer.

Liposomal formulations facilitate the administration of hydrophobic drugs, avoiding precipitation and aggregation phenomena when injected in polar solvents. The integration of the photosensitizer into the liposome may alter the fluidity of the system and, thus, modify the delivery process of the drug. Such a change has been observed for the liposomal formulation of Temoporfin, which is one of the most potent chlorin photosensitizers employed in photodynamic therapy. Here, all-atom molecular dynamics simulations have been performed to identify the nature of the intermolecular interactions that might be responsible of the different lipids freedom of motion when the drug is introduced in the bilayer. It is found that Temoporfin participates as a hydrogen donor in strong hydrogen-bonding interactions with the polar groups of the phospholipids. The theoretical analysis suggests that the rigidity of drug/liposome complexes can be modulated by considering the different hydrogen-bond ability of the photosensitizer and the carrier material.

The influence of vibrational motion on electronic excited state properties is investigated for the organic chromophore 2-nitronaphtalene in methanol. Specifically, the performance of two vibrational sampling techniques - Wigner sampling and sampling from an ab initio molecular dynamics trajectory- is assessed, in combination with implicit and explicit solvent models. The effects of the different sampling/solvent combinations on the energy and electronic character of the absorption bands are analyzed in terms of charge transfer and exciton size, computed from the electronic transition density. The absorption spectra obtained using sampling techniques and its underlying properties are compared to those of the electronic excited states calculated at the Franck-Condon equilibrium geometry. It is found that the absorption bands of the vibrational ensembles are red-shifted compared to the Franck-Condon bright states, and this red-shift scales with the displacement from the equilibrium geometry. Such displacements are found larger and better described when using ensembles from the harmonic Wigner distribution than snapshots from the molecular dynamics trajectory. Particularly relevant is the torsional motion of the nitro group that quenches the charge transfer character of some of the absorption bands. This motion, however, is better described in the molecular dynamics trajectory. Thus, none of the vibrational sampling approaches can satisfactorily capture all important aspects of the nuclear motion. The inclusion of solvent also red-shifts the absorption bands with respect to the gas phase. This red-shift scales with the charge-transfer character of the bands and is found larger for the implicit than for the explicit solvent model. The advantages and drawbacks of the different sampling and solvent models are discussed to guide future research on the calculation of UV-vis spectra of nitroaromatic compounds.

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.

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., (© 2018 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

The inclusion of solvent effects in the calculation of excited states is vital to obtain reliable absorption spectra and density of states of solvated chromophores. Here we analyze the performance of three classical approaches to describe aqueous solvent in the calculation of the absorption spectra and density of states of pyridine, tropone, and tropothione. Specifically, we compare the results obtained from quantum mechanics/polarizable continuum model (QM/PCM) versus quantum mechanics/molecular mechanics (QM/MM) in its electrostatic-embedding (QM/MMee) and polarizable-embedding (QM/MMpol) fashions, against full-QM computations, in which the solvent is described at the same level of theory as the chromophore. We show that QM/PCM provides very accurate results describing the excitation energies of ππ* and nπ* transitions, the last ones dominated by strong hydrogen-bonding effects, for the three chromophores. The QM/MMee approach also performs very well for both types of electronic transitions, although the description of the ππ* ones is slightly worse than that obtained from QM/PCM. The QM/MMpol approach performs as well as QM/PCM for describing the energy of ππ* states, but it is not able to provide a satisfactory description of hydrogen-bonding effects on the nπ* states of pyridine and tropone. The relative intensity of the absorption bands is better accounted for by the explicit-solvent models than by the continuum-solvent approach.

The characterization of electronic properties of metal complexes embedded in membrane environments is of paramount importance to develop efficient photosensitizers in optogenetic applications. Molecular dynamics and QM/MM simulations together with quantitative wave function analysis reveal a directional electronic redistribution of the exciton formed upon excitation of [Ru(bpy) 2 (bpy-C17)] 2+ when going from water to a lipid bilayer, despite the fact that the media influence neither the metal-to-ligand charge-transfer character nor the excitation energy of the absorption spectra. When the photosensitizer is embedded into the DOPC lipid membrane, exciton population is mainly located in the bypyridyl sites proximal to the positively charged surface of the bilayer due to electrostatic interactions. This behavior shows that the electronic structure of metal complexes can be controlled through the binding to external species, underscoring the crucial role of the environment in directing the electronic flow upon excitation and thus helping rational tuning of optogenetic agents.

Abstract: Cyclobutane thymine dimerization is the most prominent DNA photoinduced damage. While the ultrafast mechanism that proceeds in the singlet manifold is nowadays well established, the triplet-state pathway is not completely understood. Here we report the underlying mechanism of the photosensitized dimerization process in the triplet state. Quantum chemical calculations, combined with wavefunction analysis, and nonadiabatic molecular dynamics simulations demonstrate that this is a stepwise reaction, traversing a long-lived triplet biradical intermediate, which is characterized as a Frenkel exciton with very small charge-transfer character. The low yield of the reaction is regulated by two factors: (i) a relatively large energy barrier that needs to be overcome to form the exciton intermediate, and (ii) a bifurcation of the ground-state potential-energy surface that mostly leads back to the Franck-Condon region because dimerization requires a very restricted combination of coordinates and velocities at the event of non-radiative decay to the ground state.

The characterization of the electronically excited states of DNA strands populated upon solar UV light absorption is essential to unveil light-induced DNA damage and repair processes. We report a comprehensive analysis of the electronic properties of the UV spectrum of single-stranded polyadenine based on theoretical calculations that include excitations over eight nucleobases of the DNA strand and environmental effects by a multiscale quantum mechanics/molecular mechanics scheme, conformational sampling by molecular dynamics, and a meaningful interpretation of the electronic structure by quantitative wavefunction analysis. We show that electronic excitations are extended mainly over two nucleobases with additional important contributions of monomer-like excitations and excitons delocalized over three monomers. Half of the spectral intensity derives from locally excited and Frenkel exciton states, while states with partial charge-transfer character account for the other half and pure charge-transfer states represent only a minor contribution. The hypochromism observed when going from the isolated monomer to the strand occurs independently from delocalization and charge transfer and is instead explained by long-range environmental perturbations of the monomer states.

Skin photoprotection is commonly believed to rely on the photochemistry of 5,6-dihydroxyindole (DHI)- and 5,6-dihydroxyindole-2-carboxylic acid (DHICA)-based eumelanin building blocks. Attempts to elucidate the underlying excited-state relaxation mechanisms have been partly unsuccessful due to the marked instability to oxidation. We report a study of the excited-state deactivation of DHI using steady-state and time-resolved fluorescence accompanied by high-level quantum-chemistry calculations including solvent effects. Spectroscopic data show that deactivation of the lowest excited state of DHI in aqueous buffer proceeds on the 100 ps time scale and is 20 times faster than in methanol. Quantum-chemical calculations reveal that the excited-state decay mechanism is a sequential proton-coupled electron transfer, which involves the initial formation of a solvated electron from DHI, followed by the transfer of a proton to the solvent. This unexpected finding would prompt a revision of current notions about eumelanin photophysics and photobiology.

The lack of molecular mechanistic understanding of the interaction between metal complexes and biomolecules hampers their potential medical use. Herein we present a robust procedure combining resonant X-ray emission spectroscopy and multiscale molecular dynamics simulations, which allows for straightforward elucidation of the precise interaction mechanism at the atomic level. The report unveils an unforeseen hydrolysis process and DNA binding of [Pt{N(p-HC 6 F 4 )CH 2 } 2 py 2 ] (Pt103), which showed potential cytotoxic activity in the past. Pt103 preferentially coordinates to adjacent adenine sites, instead of guanine sites as in cisplatin, because of its hydrogen bond ability. Comparison with previous research on cisplatin suggests that selective binding to guanine or adenine may be achieved by controlling the acidity of the compound.

Multi-configurational second order perturbation theory (CASPT2) has become a very popular method for describing excited-state properties since its development in 1990. To account for systematic errors found in the calculation of dissociation energies, an empirical correction applied to the zeroth-order Hamiltonian, called the IPEA shift, was introduced in 2004. The errors were attributed to an unbalanced description of open-shell versus closed-shell electronic states and is believed to also lead to an underestimation of excitation energies. Here we show that the use of the IPEA shift is not justified and the IPEA should not be used to calculate excited states, at least for organic chromophores. This conclusion is the result of three extensive analyses. Firstly, we survey the literature for excitation energies of organic molecules that have been calculated with the unmodified CASPT2 method. We find that the excitation energies of 356 reference values are negligibly underestimated by 0.02 eV. This value is an order of magnitude smaller than the expected error based on the calculation of dissociation energies. Secondly, we perform benchmark full configuration interaction calculations on 137 states of 13 di- and triatomic molecules and compare the results with CASPT2. Also in this case, the excited states are underestimated by only 0.05 eV. Finally, we perform CASPT2 calculations with different IPEA shift values on 309 excited states of 28 organic small and medium-sized organic chromophores. We demonstrate that the size of the IPEA correction scales with the amount of dynamical correlation energy (and thus with the size of the system), and gets immoderate already for the molecules considered here, leading to an overestimation of the excitation energies. It is also found that the IPEA correction strongly depends on the size of the basis set. The dependency on both the size of the system and of the basis set, contradicts the idea of a universal IPEA shift which is able to compensate for systematic CASPT2 errors in the calculation of excited states.

The formation of cyclobutane thymine dimers is one of the most important DNA carcinogenic photolesions induced by ultraviolet irradiation. The long debated question whether thymine dimerization after direct light excitation involves singlet or triplet states is investigated here for the first time using nonadiabatic molecular dynamics simulations. We find that the precursor of this [2 + 2] cycloaddition reaction is the singlet doubly π 2 π* 2 excited state, which is spectroscopically rather dark. Excitation to the bright 1 ππ* or dark 1 nπ* excited states does not lead to thymine dimer formation. In all cases, intersystem crossing to the triplet states is not observed during the simulated time, indicating that ultrafast dimerization occurs in the singlet manifold. The dynamics simulations also show that dimerization takes place only when conformational control happens in the doubly excited state.

Although nitrobenzene is the smallest nitro-aromatic molecule, the nature of its electronic structure is still unclear. Most notably, the lowest-energy absorption band was assessed in numerous studies providing conflicting results regarding its charge-transfer character. In this study, we employ a combination of molecular dynamics and quantum chemical methods to disentangle the nature of the lowest-energy absorption band of nitrobenzene. Surprisingly, the charge-transfer transition from the benzene moiety to the nitro group is found to be quenched by a flow of charge into the opposite direction induced by vibrational motion. Beyond clarifying the charge-transfer character of nitrobenzene, we show that the widely used approach of analyzing the ground-state minimum-energy geometry provides completely wrong conclusions.

The fact that drugs currently used in the treatment of Leishmania are highly toxic and associated with acquired resistance has promoted the search for new therapies for treating American tegumentary leishmaniasis (ATL). In this study, BALB/c mice were injected in the hind paw with Leishmania (Leishmania) amazonensis and subsequently treated with a combination of nitric oxide (NO) donor (cis-[Ru(bpy) 2imN(NO)](PF6)3) (Ru-NO), given by intraperitoneal injection, and oral Brazilian propolis for 30 days. Ru-NO reached the center of the lesion and increased the NO level in the injured hind paw without lesion exacerbation. Histological and immunological parameters of chronic inflammation showed that this combined treatment increased the efficacy of macrophages, determined by the decrease in the number of parasitized cells, leading to reduced expression of proinflammatory and tissue damage markers. In addition, these drugs in combination fostered wound healing, enhanced the number of fibroblasts, pro-healing cytokines and induced collagen synthesis at the lesion site. Overall, our findings suggest that the combination of the NO donor Ru-NO and Brazilian propolis alleviates experimental ATL lesions, highlighting a new therapeutic option that can be considered for further in vivo investigations as a candidate for the treatment of cutaneous leishmaniasis.

Phenothiazinium dyes are used as photosensitizers in photodynamic therapy. Their mode of action is related to the generation of triplet excited states by intersystem crossing. Therefore, rationalizing the factors that influence intersystem crossing is crucial to improve the efficacy of photodynamic therapy. Here we employ quantum mechanics/molecular mechanics calculations to investigate the effect of aqueous and nucleic acid environments on the intersystem crossing mechanism in methylene blue. We find that the mechanism by which the triplet states are generated depends strongly on the environment. While intersystem crossing in water is mediated exclusively by vibronic spin-orbit coupling, it is enhanced in DNA due to a second pathway driven by electronic spin-orbit coupling. Competing charge-transfer processes, which are also possible in the presence of DNA, can therefore be suppressed by a suitable structural functionalization, thereby increasing the efficacy of photodynamic therapy., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The understanding of interactions between small molecules and DNA is crucial to design new anticancer drugs targeted to DNA. Methylene blue (MB) is a phenothiazinium dye that has shown promising results in photodynamic therapy treatment. The noncovalent binding of methylene blue to DNA was experimentally and theoretically analyzed in the past, but certain features of the binding mode are still not clear. In this work, force field molecular dynamics simulations were performed to simulate the binding of methylene blue to alternating GC and AT sequences at two different ionic strengths. External, intercalative, minor groove, and major groove binding modes are discussed based on energetic and structural analyses. External and major groove complexes were found to be unstable structures, although for poly(dA-dT) the major groove binding mode cannot be discarded, especially at high ionic strengths. Minor groove and intercalative binding leads to stable adducts. The most energetically favorable orientation of the dye inside the minor groove is different for the two DNA sequences because of the different balances between the DNA deformation energy and the dye/DNA interaction energy. The intercalative binding is the most important interaction mode. The dye undergoes rotational transitions inside the intercalative pocket for both DNA sequences giving rise to three dye/DNA adducts that have different energetic and structural features. This rotational motion explains the different behavior found in experiments for the GC and AT nucleic acids at different ionic strengths.

Paracoccidioidomycosis (PCM) is a granulomatous disease caused by a dimorphic fungus, Paracoccidioides brasiliensis (Pb). To determine the influence of nitric oxide (NO) on this disease, we tested cis-[Ru(bpy)2(NO)SO3](PF6), ruthenium nitrosyl, which releases NO when activated by biological reducing agents, in BALB/c mice infected intravenously with Pb 18 isolate. In a previous study by our group, the fungicidal activity of ruthenium nitrosyl was evaluated in a mouse model of acute PCM, by measuring the immune cellular response (DTH), histopathological characteristics of the granulomatous lesions (and numbers), cytokines, and NO production. We found that cis-[Ru(bpy)2(NO)SO3](PF6)-treated mice were more resistant to infection, since they exhibited higher survival when compared with the control group. Furthermore, we observed a decreased influx of inflammatory cells in the lung and liver tissue of treated mice, possibly because of a minor reduction in fungal cell numbers. Moreover, an increased production of IL-10 and a decrease in TNF-α levels were detected in lung tissues of infected mice treated with cis-[Ru(bpy)2(NO)SO3](PF6). Immunohistochemistry showed that there was no difference in the number of VEGF- expressing cells. The animals treated with cis-[Ru(bpy)2(NO)SO3](PF6) showed high NO levels at 40 days after infection. These results show that NO is effectively involved in the mechanism that regulates the immune response in lung of Pb-infected mice. These data suggest that NO is a resistance factor during paracoccidioidomycosis by controlling fungal proliferation, influencing cytokine production, and consequently moderating the development of a strong inflammatory response.