RESUMEN
We test the performance of the polarizable embedding variational quantum eigensolver self-consistent field (PE-VQE-SCF) model for computing electric field gradients with comparisons to conventional complete active space self-consistent-field (CASSCF) calculations and experimental results. We compute quadrupole coupling constants for ice VIII and ice IX. We find close agreement of the quantum-computing PE-VQE-SCF results with the results from the classical PE-CASSCF calculations and with experiment. Furthermore, we observe that the inclusion of the environment is crucial for obtaining results that match the experimental data. The calculations for ice VIII are within the experimental uncertainty for both CASSCF and VQE-SCF for oxygen and lie close to the experimental value for ice IX as well. With the VQE-SCF, which is based on an adaptive derivative-assembled problem-tailored (ADAPT) ansatz, we find that the inclusion of the environment and the size of the different basis sets do not directly affect the gate counts. However, by including an explicit environment, the wavefunction and therefore the optimization problem become more complicated, which usually results in the need to include more operators from the operator pool, thereby increasing the depth of the circuit.
RESUMEN
Cryptochromes are essential flavoproteins for circadian rhythms and avian magnetoreception. Flavin adenine dinucleotide (FAD), a chromophore within cryptochromes, absorbs blue light, initiating electron transfer processes that lead to a biological signaling cascade. A key step in this cascade is the formation of the FAD semiquinone radical (FADHâ¢), characterized through a specific red-light absorption. The absorption spectra of FADH⢠in cryptochromes are, however, significantly different from those recorded for the cofactor in solution, primarily due to protein-induced shifts in the absorption peaks. This study employs a multiscale approach, combining molecular dynamics (MD) simulations with quantum mechanical/molecular mechanical (QM/MM) methodologies, to investigate the influence of protein dynamics on embedded FADH⢠absorption. We emphasize the role of the protein's polarizable environment in the shaping of the absorption spectrum, crucial for accurate spectral predictions in cryptochromes. Our findings provide valuable insights into the absorption process, advancing our understanding of cryptochrome functioning.
Asunto(s)
Arabidopsis , Criptocromos , Flavina-Adenina Dinucleótido , Simulación de Dinámica Molecular , Teoría Cuántica , Criptocromos/química , Criptocromos/metabolismo , Arabidopsis/metabolismo , Arabidopsis/química , Flavina-Adenina Dinucleótido/química , Flavina-Adenina Dinucleótido/metabolismo , Proteínas de Arabidopsis/química , Proteínas de Arabidopsis/metabolismoRESUMEN
Subcellular membranes have complex lipid and protein compositions, which give rise to organelle-specific membrane packing, fluidity, and permeability. Due to its exquisite solvent sensitivity, the lipophilic fluorescence dye Nile Red has been used extensively to study membrane packing and polarity. Further improvement of Nile Red can be achieved by introducing electron-donating or withdrawing functional groups. Here, we compare the potential of derivatives of Nile Red with such functional substitutions for super-resolution fluorescence microscopy of lipid packing in model membranes and living cells. All studied Nile Red derivatives exhibit cholesterol-dependent fluorescence changes in model membranes, as shown by spectrally resolved stimulated emission depletion (STED) microscopy. STED imaging of Nile Red probes in cells reveals lower membrane packing in fibroblasts from healthy subjects compared to those from patients suffering from Niemann Pick type C1 (NPC1) disease, a lysosomal storage disorder with accumulation of cholesterol and sphingolipids in late endosomes and lysosomes. We also find small but consistent changes in the fluorescence lifetime of the Nile Red derivatives in NPC1 cells, suggesting altered hydrogen-bonding capacity in their membranes. All Nile Red derivatives are essentially non-fluorescent in water but increase their brightness in membranes, allowing for their use in MINFLUX single molecule tracking experiments. Our study uncovers the potential of Nile Red probes with functional substitutions for nanoscopic membrane imaging.
Asunto(s)
Colorantes Fluorescentes , Microscopía Fluorescente , Oxazinas , Oxazinas/química , Humanos , Microscopía Fluorescente/métodos , Colorantes Fluorescentes/química , Colesterol/metabolismo , Fibroblastos/metabolismo , Membrana Celular/metabolismoRESUMEN
It is still a puzzle that has not been entirely solved how migratory birds utilize the Earth's magnetic field for biannual migration. The most consistent explanation thus far is rooted in the modulation of the biological function of the cryptochrome 4 (Cry4) protein by an external magnetic field. This phenomenon is closely linked with the flavin adenine dinucleotide (FAD) cofactor that is noncovalently bound in the protein. Cry4 is activated by blue light, which is absorbed by the FAD cofactor. Subsequent electron and proton transfers trigger radical pair formation in the protein, which is sensitive to the external magnetic field. An important long-lasting redox state of the FAD cofactor is the signaling (FADHâ¢) state, which is present after the transient electron transfer steps have been completed. Recent experimental efforts succeeded in crystallizing the Cry4 protein from Columbia livia (ClCry4) with all of the important residues needed for protein photoreduction. This specific crystallization of Cry4 protein so far is the only avian cryptochrome crystal structure available, which, however, has great similarity to the Cry4 proteins of night migratory birds. The previous experimental studies of the ClCry4 protein included the absorption properties of the protein in its different redox states. The absorption spectrum of the FADH⢠state demonstrated a peculiar red shift compared to the photoabsorption properties of the FAD cofactor in its FADH⢠state in other Cry proteins from other species. The aim of this study is to understand this red shift by employing the tools of computational microscopy and, in particular, a QM/MM approach that relies on the polarizable embedding approximation.
Asunto(s)
Criptocromos , Flavina-Adenina Dinucleótido , Criptocromos/química , Criptocromos/metabolismo , Flavina-Adenina Dinucleótido/química , Flavina-Adenina Dinucleótido/metabolismo , Animales , Proteínas Aviares/química , Proteínas Aviares/metabolismo , Oxidación-ReducciónRESUMEN
We explore Davidson methods for obtaining excitation energies and other linear response properties within the recently developed quantum self-consistent linear response (q-sc-LR) method. Davidson-type methods allow for obtaining only a few selected excitation energies without explicitly constructing the electronic Hessian since they only require the ability to perform Hessian-vector multiplications. We apply the Davidson method to calculate the excitation energies of hydrogen chains (up to H10) and analyze aspects of statistical noise for computing excitation energies on quantum simulators. Additionally, we apply Davidson methods for computing linear response properties such as static polarizabilities for H2, LiH, H2O, OH-, and NH3, and show that unitary coupled cluster outperforms classical projected coupled cluster for molecular systems with strong correlation. Finally, we formulate the Davidson method for damped (complex) linear response, with application to the nitrogen K-edge X-ray absorption of ammonia, and the C6 coefficients of H2, LiH, H2O, OH-, and NH3.
RESUMEN
Determining the properties of molecules and materials is one of the premier applications of quantum computing. A major question in the field is how to use imperfect near-term quantum computers to solve problems of practical value. Inspired by the recently developed variants of the quantum counterpart of the equation-of-motion (qEOM) approach and the orbital-optimized variational quantum eigensolver (oo-VQE), we present a quantum algorithm (oo-VQE-qEOM) for the calculation of molecular properties by computing expectation values on a quantum computer. We perform noise-free quantum simulations of BeH2 in the series of STO-3G/6-31G/6-31G* basis sets and of H4 and H2O in 6-31G using an active space of four electrons and four spatial orbitals (8 qubits) to evaluate excitation energies, electronic absorption, and, for twisted H4, circular dichroism spectra. We demonstrate that the proposed algorithm can reproduce the results of conventional classical CASSCF calculations for these molecular systems.
RESUMEN
The X-ray absorption spectra of aqueous ammonia and ammonium are computed using a combination of coupled cluster singles and doubles (CCSD) with different quantum mechanical and molecular mechanical embedding schemes. Specifically, we compare frozen Hartree-Fock (HF) density embedding, polarizable embedding (PE), and polarizable density embedding (PDE). Integrating CCSD with frozen HF density embedding is possible within the CC-in-HF framework, which circumvents the conventional system-size limitations of standard coupled cluster methods. We reveal similarities between PDE and frozen HF density descriptions, while PE spectra differ significantly. By including approximate triple excitations, we also investigate the effect of improving the electronic structure theory. The spectra computed using this approach show an improved intensity ratio compared to CCSD-in-HF. Charge transfer analysis of the excitations shows the local character of the pre-edge and main-edge, while the post-edge is formed by excitations delocalized over the first solvation shell and beyond.
RESUMEN
Linear response (LR) theory is a powerful tool in classic quantum chemistry crucial to understanding photoinduced processes in chemistry and biology. However, performing simulations for large systems and in the case of strong electron correlation remains challenging. Quantum computers are poised to facilitate the simulation of such systems, and recently, a quantum linear response formulation (qLR) was introduced [Kumar et al., J. Chem. Theory Comput. 2023, 19, 9136-9150]. To apply qLR to near-term quantum computers beyond a minimal basis set, we here introduce a resource-efficient qLR theory, using a truncated active-space version of the multiconfigurational self-consistent field LR ansatz. Therein, we investigate eight different near-term qLR formalisms that utilize novel operator transformations that allow the qLR equations to be performed on near-term hardware. Simulating excited state potential energy curves and absorption spectra for various test cases, we identify two promising candidates, dubbed "proj LRSD" and "all-proj LRSD".
RESUMEN
Core-electron excitations in solvated systems, influenced by solvent geometry and hydrogen bonding, make X-ray absorption spectroscopy (XAS) a valuable tool for assessing solvent-solute interactions. However, calculating XAS spectra with electronic-structure methods has proven challenging due to a delicate interplay between correlation and solvation effects. This study provides a computational procedure for XAS modeling in solvated systems, with water-solvated ammonia and ammonium systems serving as probes. Exploring methodological challenges, we investigate explicit embedding models, specifically the polarizable embedding family, including polarizable density embedding and extended polarizable density embedding. Our linear-response time-dependent density functional theory (LR-TDDFT) XAS calculations reveal the efficiency of this approach, with extended polarizable density embedding emerging as a robust improvement over polarizable density embedding. Contrary to some recent literature, our study challenges the belief that LR-TDDFT cannot accurately describe XAS spectra of ammonia and ammonium solvated in water.
RESUMEN
Bovine rhodopsin is among the most studied proteins in the rhodopsin family. Its primary activation mechanism is the photoisomerization of 11-cis retinal, triggered by the absorption of a UV-visible photon. Different mutants of the same rhodopsin show different absorption wavelengths due to the influence of the specific amino acid residues forming the cavity in which the retinal chromophore is embedded, and rhodopsins activated at different wavelengths are, for example, exploited in the field of optogenetics. In this letter, we present a procedure for systematically investigating color tuning in models of bovine rhodopsin and a set of its mutants embedded in a membrane bilayer. Vertical excitation energy calculations were carried out with the polarizable embedding potential for describing the environment surrounding the chromophore. We show that polarizable embedding outperformed regular electrostatic embedding in determining both the vertical excitation energies and associated oscillator strengths of the systems studied.
Asunto(s)
Retina , Rodopsina , Animales , Bovinos , Rodopsina/química , Retinaldehído , FotonesRESUMEN
We formulate and implement the Variational Quantum Eigensolver Self Consistent Field (VQE-SCF) algorithm in combination with polarizable embedding (PE), thereby extending PE to the regime of quantum computing. We test the resulting algorithm, PE-VQE-SCF, on quantum simulators and demonstrate that the computational stress on the quantum device is only slightly increased in terms of gate counts compared to regular VQE-SCF. On the other hand, no increase in shot noise was observed. We illustrate how PE-VQE-SCF may lead to the modeling of real chemical systems using a simulation of the reaction barrier of the Diels-Alder reaction between furan and ethene as an example.
RESUMEN
The recently developed extended polarizable density embedding (PDE-X) model is evaluated for the spectroscopic properties of organic chromophores solvated in water, including both one- and two-photon absorption properties. The PDE-X embedding model systematically improves vertical excitation energies over the preceding polarizable density embedding model (PDE). PDE-X shows more modest improvements over existing embedding models for oscillator strengths and two-photon absorption cross-sections, which are more sensitive properties. We argue that the origin of these discrepancies is related to the description of polarization effects, suggesting directions for future development of the embedding model.
RESUMEN
Natamycin is an antifungal polyene macrolide that is used as a food preservative but also to treat fungal keratitis and other yeast infections. In contrast to other polyene antimycotics, natamycin does not form ion pores in the plasma membrane, but its mode of action is poorly understood. Using nuclear magnetic resonance (NMR) spectroscopy of deuterated sterols, we find that natamycin slows the mobility of ergosterol and cholesterol in liquid-ordered (Lo) membranes to a similar extent. This is supported by molecular dynamics (MD) simulations, which additionally reveal a strong impact of natamycin dimers on sterol dynamics and water permeability. Interference with sterol-dependent lipid packing is also reflected in a natamycin-mediated increase in membrane accessibility for dithionite, particularly in bilayers containing ergosterol. NMR experiments with deuterated sphingomyelin (SM) in sterol-containing membranes reveal that natamycin reduces phase separation and increases lipid exchange in bilayers with ergosterol. In ternary lipid mixtures containing monounsaturated phosphatidylcholine, saturated SM, and either ergosterol or cholesterol, natamycin interferes with phase separation into Lo and liquid-disordered (Ld) domains, as shown by NMR spectroscopy. Employing the intrinsic fluorescence of natamycin in ultraviolet-sensitive microscopy, we can visualize the binding of natamycin to giant unilamellar vesicles (GUVs) and find that it has the highest affinity for the Lo phase in GUVs containing ergosterol. Our results suggest that natamycin specifically interacts with the sterol-induced ordered phase, in which it disrupts lipid packing and increases solvent accessibility. This property is particularly pronounced in ergosterol containing membranes, which could underlie the selective antifungal activity of natamycin.
RESUMEN
Quantum-mechanical (QM) and classical embedding models approximate a supermolecular quantum-chemical calculation. This is particularly useful when the supermolecular calculation has a size that is out of reach for present QM models. Although QM and classical embedding methods share the same goal, they approach this goal from different starting points. In this study, we compare the polarizable embedding (PE) and frozen-density embedding (FDE) models. The former is a classical embedding model, whereas the latter is a density-based QM embedding model. Our comparison focuses on solvent effects on optical spectra of solutes. This is a typical scenario where super-system calculations including the solvent environment become prohibitively large. We formulate a common theoretical framework for PE and FDE models and systematically investigate how PE and FDE approximate solvent effects. Generally, differences are found to be small, except in cases where electron spill-out becomes problematic in the classical frameworks. In these cases, however, atomic pseudopotentials can reduce the electron-spill-out issue.
RESUMEN
The polarizable density embedding (PDE) model is a focused QM/QM fragment-based embedding model designed to model solvation effects on molecular properties. We extend the PDE model to include exchange and nonadditive exchange-correlation (for DFT) in the embedding potential in addition to the existing electrostatic, polarization, and nonelectrostatic effects already present. The resulting model, termed PDE-X, yields localized electronic excitation energies that accurately capture the range dependence of the solvent interaction and gives close agreement with full quantum mechanical (QM) results, even when using minimal QM regions. We show that the PDE-X embedding description consistently improves the accuracy of excitation energies for a diverse set of organic chromophores. The improved embedding description leads to systematic solvent effects that do not average out when applying configurational sampling.
RESUMEN
We extend the polarizable density embedding (PDE) model to support the calculation of nuclear magnetic resonance (NMR) shielding constants using gauge-including atomic orbitals (GIAOs) within a density functional theory (DFT) framework. The PDE model divides the total system into fragments, describing some by quantum mechanics (QM) and the others through an embedding model. The PDE model uses anisotropic polarizabilities, inter-fragment two-electron Coulomb integrals, and a non-local repulsion operator to emulate the QM effects. The terms involving Coulomb integrals are straightforwardly extended with GIAOs. In contrast, we consider two approaches to handle the gauge dependency of the non-local operator, employing either simple symmetrization or a gauge transformation. We find the latter approach to be most stable with respect to increasing the basis set size of the QM region. We examine the accuracy of the PDE model for calculating NMR shielding constants on several solutes in a water solution. The performance is compared with the classical polarizable embedding (PE) model in addition to supermolecular reference calculations. Based on these systems, we address the basis set convergence characteristics and the QM region size requirements. Furthermore, we investigate the performance of the PDE model for a system with significant electron spill-out. In many cases, we find that the PDE model outperforms the PE model, especially regarding the accuracy of nuclear shielding constants when using small QM region sizes and in systems with significant electron spill-out.
RESUMEN
The polarizable density embedding model is combined with the multiconfigurational self-consistent field and MC-srDFT electronic structure methods to calculate solvatochromic shifts of the n-π* absorption of acrolein and the π-π* absorption of the para-nitrophenolate anion in aqueous solution. Differences between linear-response (LR) and state-specific (SS) solvent shifts are analyzed by assessing the contributions of different terms in the solvent potential. This comparison shows that the differences are not only due to the intrinsically different response of LR and SS excitation energies to the polarizability of the environment but also due to a different response to the static part of the environment potential. These observations show that even in nonpolarizable environments, LR and SS calculations based on SCF (orbital optimization) methods do not necessarily agree on the spectral shift. The difference can be as large as, or even dominate, the difference due to dynamical polarization.
Asunto(s)
Teoría Cuántica , Agua , Acroleína , Solventes/química , Agua/químicaRESUMEN
Natamycin is a polyene macrolide, widely employed to treat fungal keratitis and other yeast infections as well as to protect food products against fungal molds. In contrast to other polyene macrolides, such as nystatin or amphotericin B, natamycin does not form pores in yeast membranes, and its mode of action is not well understood. Here, we have employed a variety of spectroscopic methods, computational modeling, and membrane reconstitution to study the molecular interactions of natamycin underlying its antifungal activity. We find that natamycin forms aggregates in an aqueous solution with strongly altered optical properties compared to monomeric natamycin. Interaction of natamycin with model membranes results in a concentration-dependent fluorescence increase which is more pronounced for ergosterol- compared to cholesterol-containing membranes up to 20 mol% sterol. Evidence for formation of specific ergosterol-natamycin complexes in the bilayer is provided. Using nuclear magnetic resonance (NMR) and electron spin resonance (ESR) spectroscopy, we find that natamycin sequesters sterols, thereby interfering with their well-known ability to order acyl chains in lipid bilayers. This effect is more pronounced for membranes containing the sterol of fungi, ergosterol, compared to those containing mammalian cholesterol. Natamycin interferes with ergosterol-dependent transport of lysine by the yeast transporter Lyp1, which we propose to be due to the sequestering of ergosterol, a mechanism that also affects other plasma membrane proteins. Our results provide a mechanistic explanation for the selective antifungal activity of natamycin, which can set the stage for rational design of novel polyenes in the future.
Asunto(s)
Natamicina , Proteínas de Saccharomyces cerevisiae , Sistemas de Transporte de Aminoácidos Básicos/metabolismo , Antibacterianos/metabolismo , Antifúngicos/química , Colesterol/química , Ergosterol/química , Lisina/metabolismo , Natamicina/metabolismo , Natamicina/farmacología , Polienos , Saccharomyces cerevisiae/metabolismo , Proteínas de Saccharomyces cerevisiae/metabolismo , Esteroles/metabolismoRESUMEN
Much effort has been invested in the search for modulators of membrane transport proteins such as P-glycoprotein (P-gp) to improve drug bioavailability and reverse multidrug resistance in cancer. Nonionic surfactants, a class of pharmaceutical excipients, are known to inhibit such proteins, but knowledge about the exact mechanism of this inhibition is scarce. Here, we perform multiscale molecular dynamics simulations of one of these surfactants, polysorbate 20 (PS20), to reveal the behavior of such compounds on the molecular level and thereby discover the molecular mechanism of the P-gp inhibition. We show that the amphiphilic headgroup of PS20 is too hydrophobic to partition in the water phase, which drives the binding of PS20 to the amphiphilic drug-binding domain of P-gp and thereby causes the inhibition of the protein. Based on our findings, we conclude that PS20 primarily inhibits P-gp through direct binding to the drug-binding domain (DBD) from the extracellular leaflet.
Asunto(s)
Miembro 1 de la Subfamilia B de Casetes de Unión a ATP , Polisorbatos , Subfamilia B de Transportador de Casetes de Unión a ATP , Excipientes/química , Polisorbatos/química , Tensoactivos/químicaRESUMEN
The coupling between the quantum (QM) and classical (MM) regions is often one of the computational bottlenecks when applying polarizable QM/MM to computational spectroscopy. In this Article, we explore three strategies to approximate the QM/MM coupling based on multipole expansion techniques. The implementations of these approximations are benchmarked in terms of both accuracy and computational efficiency and are furthermore applied to the calculation of spectroscopic properties including both one- and two-photon absorption strengths. We show that the proposed strategies provide significant computational savings without compromising the accuracy of the calculated spectroscopic properties.