The role of microenvironments on computed vibrationally-resolved emission spectra: The case of oxazines.

Oxazine dyes act as reporters of their near environment by the response of their fluorescence spectra. At the same time, their fluorescence spectra exhibit a pronounced vibrational progression. In this work, we computationally investigate the impact of near-environment models consisting of aggregated water as well as betaine molecules on the vibrational profile of fluorescence spectra of different oxazine derivatives. For aggregated betaine and a water molecule located above the plane of the dyes, we observe a distinct modification of the vibrational profile, which is more pronounced than the effect of a continuum description of a solvent environment. Our analysis shows that this effect cannot be explained by a pure change in the electronic structure, but that also vibrational degrees of freedom of the environment can be decisive for the vibrational profile and should, hence, not generally be neglected.

Vibrational embedding theory.

We suggest a consistent framework for the embedding of reduced-space correlated vibrational wave functions in a potential of the remaining modes and generalize this concept to arbitrary many subspaces. We present an implementation of this framework for vibrational coupled-cluster theory and response treatments. For C=O stretches of small molecules, we show that the embedded treatment accelerates convergence for enlarging subsets. For the water dimer and trimer as well as a water wire in bacteriorhodopsin, we investigate different partitioning schemes for the embedding approach: In the local partitioning of the vibrations, the modes dominated by motions in the same spatial region are correlated, whereas in the energy-based partitioning, modes of similar fundamental frequencies are correlated. In most cases, we obtain better agreement with superset reference results for the local partitioning than for energy-based partitioning. This work represents an important step toward multi-level methodologies in vibrational-structure theory required for its application to sizable (bio-)molecular systems.

Theoretical and Numerical Comparison of Quantum- and Classical Embedding Models for Optical Spectra.

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.

Fragment-based approach for the efficient calculation of the refractive index of metal-organic frameworks.

Increasing demands on materials in the field of optical applications require novel materials. Metal-organic frameworks (MOFs) are a prominent class of hybrid inorganic-organic materials with a modular layout. This allows the fine-tuning of their optical properties and the tailored design of optical systems. In the present theoretical study, an efficient method to calculate the refractive index (RI) of MOFs is introduced. For this purpose, the MOF is split into disjoint fragments, the linkers and the inorganic building units. The latter are disassembled until metal ions are obtained. The static polarizabilities are calculated individually using molecular density functional theory (DFT). From these, the MOF's RI is calculated. To obtain suitable polarizabilities, an exchange-correlation functional benchmark was performed first. Subsequently, this fragment-based approach was applied to a set of 24 MOFs including Zr-based MOFs and ZIFs. The calculated RI values were compared to the experimental values and validated using HSE06 hybrid functional DFT calculations with periodic boundary conditions. The examination of the MOF set revealed a speed up of the RI calculations by the fragment-based approach of up to 600 times with an estimated maximal deviation from the periodic DFT results below 4%.

Binding of a Pyrene-Based Fluorescent Amyloid Ligand to Transthyretin: A Combined Crystallographic and Molecular Dynamics Study.

Misfolding and aggregation of transthyretin (TTR) cause several amyloid diseases. Besides being an amyloidogenic protein, TTR has an affinity for bicyclic small-molecule ligands in its thyroxine (T4) binding site. One class of TTR ligands are trans-stilbenes. The trans-stilbene scaffold is also widely applied for amyloid fibril-specific ligands used as fluorescence probes and as positron emission tomography tracers for amyloid detection and diagnosis of amyloidosis. We have shown that native tetrameric TTR binds to amyloid ligands based on the trans-stilbene scaffold providing a platform for the determination of high-resolution structures of these important molecules bound to protein. In this study, we provide spectroscopic evidence of binding and X-ray crystallographic structure data on tetrameric TTR complex with the fluorescent salicylic acid-based pyrene amyloid ligand (Py1SA), an analogue of the Congo red analogue X-34. The ambiguous electron density from the X-ray diffraction, however, did not permit Py1SA placement with enough confidence likely due to partial ligand occupancy. Instead, the preferred orientation of the Py1SA ligand in the binding pocket was determined by molecular dynamics and umbrella sampling approaches. We find a distinct preference for the binding modes with the salicylic acid group pointing into the pocket and the pyrene moiety outward to the opening of the T4 binding site. Our work provides insight into TTR binding mode preference for trans-stilbene salicylic acid derivatives as well as a framework for determining structures of TTR-ligand complexes.

Fragmentation-Based Decomposition of a Metalloenzyme-Substrate Interaction: A Case Study for a Lytic Polysaccharide Monooxygenase.

We present a novel decomposition scheme for electronic interaction energies based on the flexible formulation of fragmentation schemes through fragment combination ranges (FCRs; J. Chem. Phys., 2021, 155, 164105). We devise a clear additive decomposition with contribution of nondisjoint fragments and correction terms for overlapping fragments and apply this scheme to the metalloenzyme-substrate complex of a lytic polysaccharide monooxygenase (LPMO) with an oligosaccharide. By this, we further illustrate the straightforward adaptability of the FCR-based schemes to novel systems. Our calculations suggest that the description of the electronic structure is a larger error source than the fragmentation scheme. In particular, we find a large impact of the basis set size on the interaction energies. Still, the introduction of three-body interaction terms in the fragmentation setup improves the agreement to the supermolecular reference. Yet, the qualitative results for the decomposition scheme with two-body terms only largely agree within the investigated electronic-structure approaches and basis sets, which are B97-3c, DFT (TPSS and B3LYP), and MP2 methods. The overlap contributions are found to be small, allowing analysis of the interaction energy into individual amino acid residues: We find a particularly strong interaction between the substrate and the LPMO copper active site.

Tailored anharmonic-harmonic vibrational profiles for fluorescent biomarkers.

We propose a hybrid anharmonic-harmonic scheme for vibrational broadenings, which embeds a reduced-space vibrational configuration interaction (VCI) anharmonic wave function treatment in the independent-mode displaced harmonic oscillator (IMDHO) model. The resulting systematically-improvable VCI-in-IMDHO model allows including the vibronic effects of all vibrational degrees of freedom, while focusing the effort on the important degrees of freedom with minimal extra computational effort compared to a reduced-space VCI treatment. We show for oligothiophene examples that the VCI-in-IMDHO approach can yield accurate vibrational profiles employing smaller vibrational spaces in the VCI part than the reduced-space VCI approach. By this, the VCI-in-IMDHO model enables accurate calculation of vibrational profiles of common fluorescent dyes with more than 100 vibrational degrees of freedom. We illustrate this for three examples of fluorescent biomarkers of current interest. These are the oligothiophene-based fluorescent dye called HS84, 1,4-diphenylbutadiene, and an anthracene diimide. For all examples, we assess the impact of the anharmonic treatment on the vibrational broadening, which we find to be more pronounced for the intensities than for the peak positions.

Epigenetic Priming of Bladder Cancer Cells With Decitabine Increases Cytotoxicity of Human EGFR and CD44v6 CAR Engineered T-Cells.

Background: Treatment of B-cell malignancies with CD19-directed chimeric antigen receptor (CAR) T-cells marked a new era in immunotherapy, which yet has to be successfully adopted to solid cancers. Epigenetic inhibitors of DNA methyltransferases (DNMTi) and histone deacetylases (HDACi) can induce broad changes in gene expression of malignant cells, thus making these inhibitors interesting combination partners for immunotherapeutic approaches. Methods: Urothelial carcinoma cell lines (UCC) and benign uroepithelial HBLAK cells pretreated with the DNMTi decitabine or the HDACi romidepsin were co-incubated with CAR T-cells directed against EGFR or CD44v6, and subsequent cytotoxicity assays were performed. Effects on T-cell cytotoxicity and surface antigen expression on UCC were determined by flow cytometry. We also performed next-generation mRNA sequencing of inhibitor-treated UCC and siRNA-mediated knockdown of potential regulators of CAR T-cell killing. Results: Exposure to decitabine but not romidepsin enhanced CAR T-cell cytotoxicity towards all UCC lines, but not towards the benign HBLAK cells. Increased killing could neither be attributed to enhanced target antigen expression (EGFR and CD44v6) nor fully explained by changes in the T-cell ligands PD-L1, PD-L2, ICAM-1, or CD95. Instead, gene expression analysis suggested that regulators of cell survival and apoptosis were differentially induced by the treatment. Decitabine altered the balance between survival and apoptosis factors towards an apoptosis-sensitive state associated with increased CAR T-cell killing, while romidepsin, at least partially, tilted this balance in the opposite direction. Knockdown experiments with siRNA in UCC confirmed BID and BCL2L1/BCLX as two key factors for the altered susceptibility of the UCC. Conclusion: Our data suggest that the combination of decitabine with CAR T-cell therapy is an attractive novel therapeutic approach to enhance tumor-specific killing of bladder cancer. Since BID and BCL2L1 are essential determinants for the susceptibility of a wide variety of malignant cells, their targeting might be additionally suitable for combination with immunotherapies, e.g., CAR T-cells or checkpoint inhibitors in other malignancies.

Tau Protein Binding Modes in Alzheimer's Disease for Cationic Luminescent Ligands.

The bi-thiophene-vinylene-benzothiazole (bTVBT4) ligand developed for Alzheimer's disease (AD)-specific detection of amyloid tau has been studied by a combination of several theoretical methods and experimental spectroscopies. With reference to the cryo-EM tau structure of the tau protofilament ( Nature 2017, 547, 185), a periodic model system of the fibril was created, and the interactions between this fibril and bTVBT4 were studied with nonbiased molecular dynamics simulations. Several binding sites and binding modes were identified and analyzed, and the results for the most prevailing fibril site and ligand modes are presented. A key validation of the simulation work is provided by the favorable comparison of the theoretical and experimental absorption spectra of bTVBT4 in solution and bound to the protein. It is conclusively shown that the ligand-protein binding occurs at the hydrophobic pocket defined by the residues Ile360, Thr361, and His362. This binding site is not accessible in the Pick's disease (PiD) fold, and fluorescence imaging of bTVBT4-stained brain tissue samples from patients diagnosed with AD and PiD provides strong support for the proposed tau binding site.

A unified and flexible formulation of molecular fragmentation schemes.

We present a flexible formulation for energy-based molecular fragmentation schemes. This framework does not only incorporate the majority of existing fragmentation expansions but also allows for flexible formulation of novel schemes. We further illustrate its application in multi-level approaches and for electronic interaction energies. For the examples of small water clusters, a small protein, and protein-protein interaction energies, we show how this flexible setup can be exploited to generate a well-suited multi-level fragmentation expansion for the given case. With such a setup, we reproduce the electronic protein-protein interaction energy of ten different structures of a neurotensin and an extracellular loop of its receptor with a mean absolute deviation to the respective super-system calculations below 1 kJ/mol.

Toward Accurate Theoretical Vibrational Spectra: A Case Study for Maleimide.

We employ and combine a number of recent developments in vibrational structure methods to push their current size limitations toward molecules with tens of modes and showcase their availability for the maleimide molecule. In particular, we assess the use of different rectilinear vibrational coordinates, namely, normal coordinates, hybrid optimized and localized coordinates, and flexible adaptation of local coordinates of nuclei coordinates. These different coordinate parameterizations are employed in conjunction with the adaptive density-guided approach to generate potential energy surfaces (PESs). A screening procedure is furthermore introduced, which provides estimates of the importance of individual terms in the PES, resulting in significant reductions in the computational cost of the PES construction. We find that all three sets of coordinates provide approximately the same level of accuracy in vibrational structure calculations and report fundamental excitation energies with a mean absolute deviation of less than 12 cm-1 when compared to experimental data. We expect that similar accuracy in vibrational structure calculations can be achieved for molecules of larger size using the proposed procedures.

Adaptive density-guided approach to double incremental potential energy surface construction.

We present a combination of the recently developed double incremental expansion of potential energy surfaces with the well-established adaptive density-guided approach to grid construction. This unique methodology is based on the use of an incremental expansion for potential energy surfaces, known as n-mode expansion; an incremental many-body representation of the electronic energy; and an efficient vibrational density-guided approach to automated determination of grid dimensions and granularity. The reliability of the method is validated calculating potential energy surfaces and obtaining fundamental excitation energies for three moderate-size chain-like molecular systems. The use of our methodology leads to considerable computational savings for potential energy surface construction compared to standard approaches while maintaining a high level of accuracy in the resulting potential energy surfaces. Additional investigations indicate that our method can be applied to covalently bound and strongly interacting molecular systems, even though these cases are known to be very unfavorable for fragmentation schemes. We therefore conclude that the presented methodology is a robust and flexible approach to potential energy surface construction, which introduces considerable computational savings without compromising the accuracy of vibrational spectra calculations.

Vibrationally resolved emission spectra of luminescent conjugated oligothiophenes from anharmonic calculations.

We report on accurate and efficient calculations of vibrationally resolved emission spectra for oligothiophenes from anharmonic vibrational configuration interaction wave-function calculations in reduced vibrational spaces. These reduced spaces are chosen based on the independent mode displaced harmonic oscillator model. Good agreement with experiment is obtained for all-trans oligothiophenes with two to five rings also when employing only a few active modes. Vibrational modes incorporating inter-ring carbon-carbon stretches and a ring breathing mode are found to be the main players in the vibrational progression for the emission from the first excited electronic state for all investigated oligothiophene derivatives. The presented framework is here illustrated for oligothiophenes, but we have made no underlying system-dependent assumptions and believe it to become a valuable tool for the rational design of fluorescence biomarkers.

Binding sites for luminescent amyloid biomarkers from non-biased molecular dynamics simulations.

A very stable binding site for the interaction between a pentameric oligothiophene and an amyloid-ß(1-42) fibril has been identified by means of non-biased molecular dynamics simulations. In this site, the probe is locked in an all-trans conformation with a Coulombic binding energy of 1200 kJ mol-1 due to the interactions between the anionic carboxyl groups of the probe and the cationic ε-amino groups in the lysine side chain. Upon binding, the conformationally restricted probes show a pronounced increase in molecular planarity. This is in line with the observed changes in luminescence properties that serve as the foundation for their use as biomarkers.

Anharmonic vibrational spectra from double incremental potential energy and dipole surfaces.

We extend the fragmentation-based double incremental expansion in FALCON coordinates (DIF) and its linear-scaling analogue [C. König and O. Christiansen, J. Chem. Phys., 2016, 145, 064105] to dipole surfaces. Thereby, we enable the calculation of intensities in vibrational absorption spectra from these cost-efficient property surfaces. We validate the obtained potential energy and dipole surfaces by vibrational spectra calculations employing damped response theory for correlated vibrational coupled cluster wave functions. Our largest calculation on a hexa-phenyl includes all 180 vibrational degrees of freedom of the system, which illustrates the potential of both the DIF schemes for property surface generation and the use of damped response theory from high-dimensional correlated vibrational wave functions. Generally, we obtain good agreement between the spectra calculated from the DIF property surfaces and the non-fragmented analogues. Moreover, when adopting suitable electronic structure methods, good agreement with respect to the experiment can be obtained, as shown for the example of 5-methylfurfural and RI-MP2. In conclusion, our results illustrate that the presented scheme with linearly scaling surfaces enables high quality spectra, as long as reasonably sized fragments can be defined. With this work, we push the realistic limits of vibrational spectra calculations from vibrational wave function methods and accurate electronic structure calculations to significantly larger systems than currently accessible.

FALCON: A method for flexible adaptation of local coordinates of nuclei.

We present a flexible scheme for calculating vibrational rectilinear coordinates with well-defined strict locality on a certain set of atoms. Introducing a method for Flexible Adaption of Local COordinates of Nuclei (FALCON) we show how vibrational subspaces can be "grown" in an adaptive manner. Subspace Hessian matrices are set up and used to calculate and analyze vibrational modes and frequencies. FALCON coordinates can more generally be used to construct vibrational coordinates for describing local and (semi-local) interacting modes with desired features. For instance, spatially local vibrations can be approximately described as internal motion within only a group of atoms and delocalized modes can be approximately expressed as relative motions of rigid groups of atoms. The FALCON method can support efficiency in the calculation and analysis of vibrational coordinates and energies in the context of harmonic and anharmonic calculations. The features of this method are demonstrated on a few small molecules, i.e., formylglycine, coumarin, and dimethylether as well as for the amide-I band and low-frequency modes of alanine oligomers and alpha conotoxin.

Biomechanical characteristics of bioabsorbable magnesium-based (MgYREZr-alloy) interference screws with different threads.

PURPOSE: Degradable magnesium implants have received increasing interest in recent years. In anterior cruciate ligament reconstruction surgery, the well-known osteoconductive effects of biodegradable magnesium alloys may be useful. The aim of this study was to examine whether interference screws made of MgYREZr have comparable biomechanical properties to commonly used biodegradable screws and whether a different thread on the magnesium screw has an influence on the fixation strength. METHODS: Five magnesium (MgYREZr-alloy) screws were tested per group. Three different groups with variable thread designs (Designs 1, 2, and 3) were produced and compared with the commercially available bioabsorbable Bioacryl rapid polylactic-co-glycolic acid screw Milagro®. In vitro testing was performed in synthetic bone using artificial ligament fixed by an interference screw. The constructs were pretensioned with a constant load of 60 N for 30 s followed by 500 cycles between 60 N and 250 N at 1 Hz. Construct displacements between the 1st and 20th and the 21st and 500th cycles were recorded. After a 30 s break, a maximum load to failure test was performed at 1 mm/s measuring the maximum pull-out force. RESULTS: The maximum loads to failure of all three types of magnesium interference screws (Design 1: 1,092 ± 133.7 N; Design 2: 1,014 ± 103.3 N; Design 3: 1,001 ± 124 N) were significantly larger than that of the bioabsorbable Milagro® interference screw (786.8 ± 62.5 N) (p < 0.05). However, the greatest maximum load was found with magnesium screw Design 1. Except for a significant difference between Designs 1 and 2, there were no further significant differences among the four groups in displacement after the 20th cycle. CONCLUSIONS: Biomechanical testing showed higher pull-out forces for magnesium compared with a commercial polymer screw. Hence, they suggest better stability and are a potential alternative. The thread geometry does not significantly influence the stability provided by the magnesium implants. This study shows the first promising results of a degradable material, which may be a clinical alternative in the future.

Hybrid Optimized and Localized Vibrational Coordinates.

We present a new type of vibrational coordinates denoted hybrid optimized and localized coordinates (HOLCs) aiming at a good set of rectilinear vibrational coordinates supporting fast convergence in vibrational stucture calculations. The HOLCs are obtained as a compromise between the recently promoted optimized coordinates (OCs) and localized coordinates (LCs). The three sets of coordinates are generally different from each other and differ from standard normal coordinates (NCs) as well. In determining the HOLCs, we optimize the vibrational self-consistent field (VSCF) energy with respect to orthogonal transformation of the coordinates, which is similar to determining OCs but for HOLCs we additionally introduce a penalty for delocalization, by using a measure of localization similar to that employed in determining LCs. The same theory and implementation covers OCs, LCs, and HOLCs. It is shown that varying one penalty parameter allows for connecting OCs and LCs. The HOLCs are compared to NCs, OCs, and LCs in their nature and performance as basis for vibrational coupled cluster (VCC) response calculations of vibrational anharmonic energies for a small set of simple systems comprising water, formaldehyde, and ethylene. It is found that surprisingly good results can be obtained with HOLCs by using potential energy surfaces as simple as quadratic Taylor expansions. Quite similar coordinates are found for the already established OCs but obtaining these OCs requires much more elaborate and expensive potential energy surfaces and localization is generally not guaranteed. The ability to compute HOLCs for somewhat larger systems is demonstrated for coumarin and the alanine quadramer. The good agreement between HOLCs and OCs, together with the much easier applicability of HOLCs for larger systems, suggests that HOLCs may be a pragmatically very interesting option for anharmonic calculations on medium to large molecular systems.

Automatic determination of important mode-mode correlations in many-mode vibrational wave functions.

We introduce new automatic procedures for parameterizing vibrational coupled cluster (VCC) and vibrational configuration interaction wave functions. Importance measures for individual mode combinations in the wave function are derived based on upper bounds to Hamiltonian matrix elements and/or the size of perturbative corrections derived in the framework of VCC. With a threshold, this enables an automatic, system-adapted way of choosing which mode-mode correlations are explicitly parameterized in the many-mode wave function. The effect of different importance measures and thresholds is investigated for zero-point energies and infrared spectra for formaldehyde and furan. Furthermore, the direct link between important mode-mode correlations and coordinates is illustrated employing water clusters as examples: Using optimized coordinates, a larger number of mode combinations can be neglected in the correlated many-mode vibrational wave function than with normal coordinates for the same accuracy. Moreover, the fraction of important mode-mode correlations compared to the total number of correlations decreases with system size. This underlines the potential gain in efficiency when using optimized coordinates in combination with a flexible scheme for choosing the mode-mode correlations included in the parameterization of the correlated many-mode vibrational wave function. All in all, it is found that the introduced schemes for parameterizing correlated many-mode vibrational wave functions lead to at least as systematic and accurate calculations as those using more standard and straightforward excitation level definitions. This new way of defining approximate calculations offers potential for future calculations on larger systems.

Vibronic-structure tracking: a shortcut for vibrationally resolved UV/Vis-spectra calculations.

The vibrational coarse structure and the band shapes of electronic absorption spectra are often dominated by just a few molecular vibrations. By contrast, the simulation of the vibronic structure even in the simplest theoretical models usually requires the calculation of the entire set of normal modes of vibration. Here, we exploit the idea of the mode-tracking protocol [M. Reiher and J. Neugebauer, J. Chem. Phys. 118, 1634 (2003)] in order to directly target and selectively calculate those normal modes which have the largest effect on the vibronic band shape for a certain electronic excitation. This is achieved by defining a criterion for the importance of a normal mode to the vibrational progressions in the absorption band within the so-called "independent mode, displaced harmonic oscillator" (IMDHO) model. We use this approach for a vibronic-structure investigation for several small test molecules as well as for a comparison of the vibronic absorption spectra of a truncated chlorophyll a model and the full chlorophyll a molecule. We show that the method allows to go beyond the often-used strategy to simulate absorption spectra based on broadened vertical excitation peaks with just a minimum of computational effort, which in case of chlorophyll a corresponds to about 10% of the cost for a full simulation within the IMDHO approach.