Red Fluorescence from Organic Microdots: Leveraging Foldamer-Linked Azobenzene for Enhanced Stability and Intensity in Bioimaging Applications.

Azobenzene, while relevant, has faced constraints in biological system applications due to its suboptimal quantum yield and short-wavelength emission. This study presents a pioneering strategy for fabricating organic microdots by coupling foldamer-linked azobenzene, resulting in robust fluorescence intensity and stability, especially in aggregated states, thereby showing promise for bioimaging applications. Comprehensive experimental and computational examinations elucidate the mechanisms underpinning enhanced photostability and fluorescence efficacy. In vitro and in vivo evaluations disclose that the external layer of cis-azo-foldamer microdots performs a self-sacrificial function during photo-bleaching. Consequently, these red-fluorescent microdots demonstrate extraordinary structural and photochemical stabilities over extended periods. The conjugation of a ß-peptide foldamer to the azobenzene chromophore through a glycine linker instigates a blue-shifted and amplified π*-n transition. Molecular dynamics simulations reveal that the aggregated state of cis-azo-foldamers fortifies the stability of cis isomers, thereby augmenting fluorescence efficiency. This investigation furnishes crucial insights into conceptualizing novel, biologically inspired materials, promising stable and enduring imaging applications, and carries implications for diverse arenas such as medical diagnostics, drug delivery, and sensing technologies.

Dynamic embedding of effective harmonic normal mode vibrations in all-atomistic energy gap fluctuations: Case study of light harvesting 2 complex.

Environmental effects in excitation energy transfer have mostly been modeled by baths of harmonic oscillators, but to what extent such modeling provides a reliable description of actual interactions between molecular systems and environments remains an open issue. We address this issue by investigating fluctuations in the excitation energies of the light harvesting 2 complex using a realistic all-atomistic simulation of the potential energy surface. Our analyses reveal that molecular motions exhibit significant anharmonic features, even for underdamped intramolecular vibrations. In particular, we find that the anharmonicity contributes to the broadening of spectral densities and substantial overlaps between neighboring peaks, which complicates the meaning of mode frequencies constituting a bath model. Thus, we develop a strategy to construct a minimally underdamped harmonic bath that has a clear connection to all-atomistic dynamics by utilizing actual normal modes of molecules but optimizing their frequencies such that the resulting bath model can best reproduce the all-atomistic simulation results. By subtracting the underdamped contribution from the entire fluctuations, we also show that identifying a residual spectral density representing all other contributions with overdamped behavior is possible. We find that this can be fitted well with a well-established analytic form of a spectral density function or, alternatively, modeled as explicit time dependent fluctuations with muti-exponential or power law type correlation functions. We provide an assessment and the implications of these possibilities. The approach presented here can also serve as a general strategy to construct a simplified bath model that can effectively represent the underlying all-atomistic bath dynamics.

An Ensemble Docking Approach for Analyzing and Designing Aptamer Heterodimers Targeting VEGF165.

Vascular endothelial growth factor 165 (VEGF165) is a prominent isoform of the VEGF-A protein that plays a crucial role in various angiogenesis-related diseases. It is homodimeric, and each of its monomers is composed of two domains connected by a flexible linker. DNA aptamers, which have emerged as potent therapeutic molecules for many proteins with high specificity and affinity, can also work for VEGF165. A DNA aptamer heterodimer composed of monomers of V7t1 and del5-1 connected by a flexible linker (V7t1:del5-1) exhibits a greater binding affinity with VEGF165 compared to either of the two monomers alone. Although the structure of the complex formed between the aptamer heterodimer and VEGF165 is unknown due to the highly flexible linkers, gaining structural information will still be valuable for future developments. Toward this end of accessing structural information, we adopt an ensemble docking approach here. We first obtain an ensemble of structures for both VEGF165 and the aptamer heterodimer by considering both small- and large-scale motions. We then proceed through an extraction process based on ensemble docking, molecular dynamics simulations, and binding free energy calculations to predict the structures of the VEGF165/V7t1:del5-1 complex. Through the same procedures, we reach a new aptamer heterodimer that bears a locked nucleic acid-modified counterpart of V7t1, namely RNV66:del5-1, which also binds well with VEGF165. We apply the same protocol to the monomeric units V7t1, RNV66, and del5-1 to target VEGF165. We observe that V7t1:del5-1 and RNV66:del5-1 show higher binding affinities with VEGF165 than any of the monomers, consistent with experiments that support the notion that aptamer heterodimers are more effective anti-VEGF165 aptamers than monomeric aptamers. Among the five different aptamers studied here, the newly designed RNV66:del5-1 shows the highest binding affinity with VEGF165. We expect that our ensemble docking approach can help in de novo designs of homo/heterodimeric anti-angiogenic drugs to target the homodimeric VEGF165.

Thermodynamic consequences of stapling side-chains on a peptide ligand using a lactam-bridge: A theoretical study on anti-angiogenic peptides targeting VEGF.

Peptides are promising therapeutic agents for various biological targets due to their high efficacy and low toxicity, and the design of peptide ligands with high binding affinity to the target of interest is of utmost importance in peptide-based drug design. Introducing a conformational constraint to a flexible peptide ligand using a side-chain lactam-bridge is a convenient and efficient method to improve its binding affinity to the target. However, in general, such a small structural modification to a flexible ligand made with the intent of lowering the configurational entropic penalty for binding may have unintended consequences in different components of the binding enthalpy and entropy, including the configurational entropy component, which are still not clearly understood. Toward probing this, we examine different components of the binding enthalpy and entropy as well as the underlying structure and dynamics, for a side-chain lactam-bridged peptide inhibitor and its flexible analog forming complexes with vascular endothelial growth factor (VEGF), using all-atom molecular dynamics simulations. It is found that introducing a side-chain lactam-bridge constraint into the flexible peptide analog led to a gain in configurational entropy change but losses in solvation entropy, solute internal energy, and solvation energy changes upon binding, pinpointing the opportunities and challenges in drug design. The present study features an interplay between configurational and solvation entropy changes, as well as the one between binding enthalpy and entropy, in ligand-target binding upon imposing a conformational constraint into a flexible ligand.

Potential energy interpolation with target-customized weighting coordinates: application to excited-state dynamics of photoactive yellow protein chromophore in water.

Interpolation of potential energy surfaces (PESs) can provide a practical route to performing molecular dynamics simulations with a reliability matching a high-level quantum chemical calculation. An obstacle to its widespread use is perhaps the lack of general and optimal interpolation settings that can be applied in a black-box manner for any given molecular system. How to set up the weights for interpolation is one such task, and we still need to diversify the approaches in order to treat various systems. Here, we develop a new interpolation weighting scheme, which allows us to choose the weighting coordinates in a system-specific manner, by amplifying the contribution from specific internal coordinates. The new weighting scheme with an appropriate selection of coordinates is proved to be effective in reducing the interpolation error along the reaction pathway. As a demonstration, we consider the photoactive yellow protein chromophore system, as it constitutes itself as an interesting target that bears long-standing questions related to excited-state dynamics inside protein environments. We build its two-state diabatic interpolated PES with the new weighting scheme. We indeed see the utility of our scheme by conducting nonadiabatic molecular dynamics simulations with the required semi-global PES based on a limited number of data points.

Variational quantum eigensolver for closed-shell molecules with non-bosonic corrections.

The realization of quantum advantage with noisy-intermediate-scale quantum (NISQ) machines has become one of the major challenges in computational sciences. Maintaining coherence of a physical system with more than ten qubits is a critical challenge that motivates research on compact system representations to reduce algorithm complexity. Toward this end, the variational quantum eigensolver (VQE) used to perform quantum simulations is considered to be one of the most promising algorithms for quantum chemistry in the NISQ era. We investigate reduced mapping of one spatial orbital to a single qubit to analyze the ground state energy in a way that the Pauli operators of qubits are mapped to the creation/annihilation of singlet pairs of electrons. To include the effect of non-bosonic (or non-paired) excitations, we introduce a simple correction scheme in the electron correlation model approximated by the geometrical mean of the bosonic (or paired) terms. Employing it in a VQE algorithm, we assess ground state energies of H2O, N2, and Li2O in good agreement with full configuration interaction (FCI) models respectively, using only 6, 8, and 12 qubits with quantum gate depths proportional to the squares of the qubit counts. With the adopted seniority-zero approximation that uses only one half of the qubit counts of a conventional VQE algorithm, we find that our non-bosonic correction method reaches reliable quantum chemistry simulations at least for the tested systems.

Electron-Binding Dynamics of the Dipole-Bound State: Correlation Effect on the Autodetachment Dynamics.

The nature of the electron-binding forces in the dipole-bound states (DBS) of anions is interrogated through experimental and theoretical means by investigating the autodetachment dynamics from DBS Feshbach resonances of ortho-, meta-, and para-bromophenoxide (BrPhO-). Though the charge-dipole electrostatic potential has been widely regarded to be mainly responsible for the electron binding in DBS, the effect of nonclassical electron correlation has been conceived to be quite significant in terms of its static and/or dynamic contributions toward the binding of the excess electron to the neutral core. State-specific real-time autodetachment dynamics observed by picosecond time-resolved photoelectron velocity-map imaging spectroscopy reveal that the autodetachment processes from the DBS Feshbach resonances of BrPhO- anions cannot indeed be rationalized by the conventional charge-dipole potential. Specifically, the autodetachment lifetime is drastically lengthened depending on differently positioned Br-substitution, and this rate change cannot be explained within the framework of Fermi's golden rule based on the charge-dipole assumption. High-level ab initio quantum chemical calculations with EOM-EA-CCSD, which intrinsically takes into account electron correlations, generate more reasonable predictions on the binding energies than density functional theory (DFT) calculations, and semiclassical quantum dynamics simulations based on the EOM-EA-CCSD data excellently predict the trend in the autodetachment rates. These findings illustrate that static and dynamic properties of the excess electron in the DBS are strongly influenced by correlation interactions among electrons in the nonvalence orbital of the dipole-bound electron and highly polarizable valence orbitals of the bromine atom, which, in turn, dictate the interesting chemical fate of exotic anion species.

Overcoming the Limitation of Spin Statistics in Organic Light Emitting Diodes (OLEDs): Hot Exciton Mechanism and Its Characterization.

Triplet harvesting processes are essential for enhancing efficiencies of fluorescent organic light-emitting diodes. Besides more conventional thermally activated delayed fluorescence and triplet-triplet annihilation, the hot exciton mechanism has been recently noticed because it helps reduce the efficiency roll-off and improve device stability. Hot exciton materials enable the conversion of triplet excitons to singlet ones via reverse inter-system crossing from high-lying triplet states and thereby the depopulation of long-lived triplet excitons that are prone to chemical and/or efficiency degradation. Although their anti-Kasha characteristics have not been clearly explained, numerous molecules with behaviors assigned to the hot exciton mechanism have been reported. Indeed, the related developments appear to have just passed the stage of infancy now, and there will likely be more roles that computational elucidations can play. With this perspective in mind, we review some selected experimental studies on the mechanism and the related designs and then on computational studies. On the computational side, we examine what has been found and what is still missing with regard to properly understanding this interesting mechanism. We further discuss potential future points of computational interests toward aiming for eventually presenting in silico design guides.

Modified Fermi's golden rule rate expressions.

Fermi's golden rule (FGR) serves as the basis for many expressions of spectroscopic observables and quantum transition rates. The utility of FGR has been demonstrated through decades of experimental confirmation. However, there still remain important cases where the evaluation of a FGR rate is ambiguous or ill-defined. Examples are cases where the rate has divergent terms due to the sparsity in the density of final states or time dependent fluctuations of system Hamiltonians. Strictly speaking, assumptions of FGR are no longer valid for such cases. However, it is still possible to define modified FGR rate expressions that are useful as effective rates. The resulting modified FGR rate expressions resolve a long standing ambiguity often encountered in using FGR and offer more reliable ways to model general rate processes. Simple model calculations illustrate the utility and implications of new rate expressions.

Elucidating activation and deactivation dynamics of VEGFR-2 transmembrane domain with coarse-grained molecular dynamics simulations.

The vascular endothelial growth factor receptor 2 (VEGFR-2) is a member of receptor tyrosine kinases (RTKs) and is a dimeric membrane protein that functions as a primary regulator of angiogenesis. As is usual with RTKs, spatial alignment of its transmembrane domain (TMD) is essential toward VEGFR-2 activation. Experimentally, the helix rotations within TMD around their own helical axes are known to participate importantly toward the activation process in VEGFR-2, but the detailed dynamics of the interconversion between the active and inactive TMD forms have not been clearly elucidated at the molecular level. Here, we attempt to elucidate the process by using coarse grained (CG) molecular dynamics (MD) simulations. We observe that inactive dimeric TMD in separation is structurally stable over tens of microseconds, suggesting that TMD itself is passive and does not allow spontaneous signaling of VEGFR-2. By starting from the active conformation, we reveal the mechanism of TMD inactivation through analyzing the CG MD trajectories. We observe that interconversions between a left-handed overlay and a right-handed one are essential for the process of going from an active TMD structure to the inactive form. In addition, our simulations find that the helices can rotate properly when the overlaying structure of the helices interconverts and when the crossing angle of the two helices changes by larger than ~40 degrees. As the activation right after the ligand attachment on VEGFR-2 will take place in the reverse manner of this inactivation process, these structural aspects will also appear importantly for the activation process. The rather large change in helix configuration for activation also explains why VEGFR-2 rarely self-activate and how the activating ligand structurally drive the whole VEGFR-2. This mechanism of TMD activation / inactivation within VEGFR-2 may help in further understanding the overall activation processes of other RTKs.

Molecular mechanism of binding between a therapeutic RNA aptamer and its protein target VEGF: A molecular dynamics study.

Macugen is a therapeutic RNA aptamer against vascular endothelial growth factor (VEGF)-165, the VEGF isoform primarily responsible for angiogenesis. It has been reported that Macugen inhibits angiogenesis by specifically binding to the heparin binding domain (HBD) of VEGF165. The mechanism of the molecular recognition between HBD and Macugen is investigated here using all-atom molecular dynamics simulations. We find that Macugen recognizes HBD by an induced-fit mechanism with major conformational changes in Macugen and almost no changes in the structure of HBD, whereas HBD recognizes Macugen by a conformational selection mechanism.

Mechanism of Ir(ppy)3 Guest Exciton Formation with the Exciplex-Forming TCTA:TPBI Cohost within a Phosphorescent Organic Light-Emitting Diode Environment.

Cohosts based on hole transporting and electron transporting materials often act as exciplexes in the form of intermolecular charge transfer complexes. Indeed, exciplex-forming cohosts have been widely developed as the host materials for efficient phosphorescent organic light-emitting diodes (OLEDs). In host-guest systems of OLEDs, the guest can be excited by two competing mechanisms, namely, excitation energy transfer (EET) and charge transfer (CT). Experimentally, it has been reported that the EET mechanism is dominant and the excitons are primarily formed in the host first and then transferred to the guest in phosphorescent OLEDs based on exciplex-forming cohosts. With this, exciplex-forming cohosts are widely employed for avoiding the formation of trapped charge carriers in the phosphorescent guest. However, theoretical studies are still lacking toward elucidating the relative importance between EET and CT processes in exciting the guest molecules in such systems. Here, we obtain the kinetics of guest excitation processes in a few trimer model systems consisting of an exciplex-forming cohost pair and a phosphorescent guest. We adopt the Förster resonance energy transfer (FRET) rate constants for the electronic transitions between excited states toward solving kinetic master equations. The input parameters for calculating the FRET rate constants are obtained from density functional theory (DFT) and time-dependent DFT. The results show that while the EET mechanism is important, the CT mechanism may still play a significant role in guest excitations. In fact, the relative importance of CT over EET depends strongly on the location of the guest molecule relative to the cohost pair. This is understandable as both the coupling for EET and the interaction energy for CT are strongly influenced by the geometric constraints. Understanding the energy transfer pathways from the exciplex state of cohost to the emissive state of guest may provide insights for improving exciplex-forming materials adopted in OLEDs.

Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes.

Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.

Three States Involving Vibronic Resonance is a Key to Enhancing Reverse Intersystem Crossing Dynamics of an Organoboron-Based Ultrapure Blue Emitter.

The recently developed narrow-band blue-emitting organoboron chromophores based on the multiple-resonance (MR) effect have now become one of the most important components for constructing efficient organic light emitting diodes (OLEDs). While they basically emit through fluorescence, they are also known for showing substantial thermally activated delayed fluorescence (TADF) even with a relatively large singlet-triplet gap (ΔE ST). Indeed, understanding the reverse intersystem crossing (RISC) dynamics behind this peculiar TADF will allow judicious molecular designs toward achieving better performing OLEDs. Explaining the underlying nonadiabatic spin-flip mechanism, however, has often been equivocal, and how the sufficiently fast RISC takes place even with the sizable ΔE ST and vanishingly small spin-orbit coupling is not well understood. Here, we show that a vibronic resonance, namely the frequency matching condition between the vibration and the electronic energy gap, orchestrates three electronic states together and this effect plays a major role in enhancing RISC in a typical organoboron emitter. Interestingly, the mediating upper electronic state is quite high in energy to an extent that its thermal population is vanishingly small. Through semiclassical quantum dynamics simulations, we further show that the geometry dependent non-Condon coupling to the upper triplet state that oscillates with the frequency ΔE ST/ℏ is the main driving force behind the peculiar resonance enhancement. The existence of an array of vibrational modes with strong vibronic rate enhancements provides the ability to sustain efficient RISC over a range of ΔE ST in defiance of the energy gap law, which can render the MR-emitters peculiar in comparison with more conventional donor-acceptor type emitters. Our investigation may provide a new guide for future blue emitting molecule developments.

Analytical gradients for core-excited states in the algebraic diagrammatic construction (ADC) framework.

Expressions for analytical molecular gradients of core-excited states have been derived and implemented for the hierarchy of algebraic diagrammatic construction (ADC) methods up to extended second-order within the core-valence separation (CVS) approximation. We illustrate the use of CVS-ADC gradients by determining relaxed core-excited state potential energy surfaces and optimized geometries for water, formic acid, and benzene. For water, our results show that in the dissociative lowest core-excited state, a linear configuration is preferred. For formic acid, we find that the O K-edge lowest core-excited state is non-planar, a fact that is not captured by the equivalent core approximation where the core-excited atom with its hole is replaced by the "Z + 1" neighboring atom in the periodic table. For benzene, the core-excited state gradients are presented along the Jahn-Teller distorted geometry of the 1s → π* excited state. Our development may pave a new path to studying the dynamics of molecules in their core-excited states.

Uncovering the Conformational Distribution of a Small Protein with Nanoparticle-Aided Cryo-Electron Microscopy Sampling.

Here, we introduce the nanoparticle-aided cryo-electron microscopy sampling (NACS) method to access the conformational distribution of a protein molecule. Two nanogold particles are labeled at two target sites, and the interparticle distance is measured as a structural parameter via cryo-electron microscopy (cryo-EM). The key aspect of NACS is that the projected distance information instead of the global conformational information is extracted from each protein molecule. This is possible because the contrast provided by the nanogold particles is strong enough to provide the projected distance, while the protein itself is invisible due to its low contrast. We successfully demonstrate that various protein conformations, even for small or disordered proteins, which generally cannot be accessed via cryo-EM, can be captured. The demonstrated method with the potential to directly observe the conformational distribution of such systems may open up new possibilities in studying their dynamics at a single-molecule level.

Cooperation between Excitation Energy Transfer and Antisynchronously Coupled Vibrations.

The effects of the environment in energy transfer systems have been continuously studied for decades. Here, we investigate how the energy transfer and the emergence of vibrational correlations cooperate with each other based on simulations with a few numerically approximate mixed quantum classical (MQC) methods. By adopting a two-state system with locally coupled underdamped vibrations that are resonant with the electronic energy gap, we observe prominent energy dissipations from the electronic system to the vibrations, rehighlighting the role of underdamped vibrations as a temporal electronic energy buffer. More importantly, this energy dissipation generates specific phase relations between the two vibrations. Namely, the vibrations become anticorrelated right after the initiation of the energy transfer but then synchronized as the transfer completes. These phase relations are interpreted as a selective activation of an anticorrelated motion of the vibrations and a subsequent deactivation by thermal energy redistribution. Furthermore, we show that a single vibration simultaneously coupled to the two electronic states with opposite phases induces a completely equivalent energy transfer dynamics as the two localized vibrations. Finally, we discuss how the vibrational energy dissipation dynamics is affected by the adopted MQC approaches and warn about the increased subtlety toward properly treating dissipation effects over having reliable population dynamics.

Theoretical Study on the Degree of CO2 Activation in CO2-Coordinated Ni(0) Complexes.

The geometrical characteristic and the degree of CO2 activation of the CO2-coordinated Ni(0) complexes were investigated computationally by quantum chemical means for bidentate and tridentate ligands of PP, PPMeP, and PNP, and sometimes with co-complexing Fe(II) to differently coordinate CO2. We show that the coordination geometry of the central metal is determined by the ligand geometry. The charge and the energy decomposition analyses show that the charge transfer energy through orbital mixing has a strong correlation with CO2 net charge, while the binding energy cannot due to the lack of the coordination number and the deformation energy of the ligand. Among the examined ligands, PNP with negatively charged secondary amine makes Ni(0) an electron-rich atom, which results in an ∼20% higher CO2 activation than those of PP and PPMeP. In particular, Fe(II)-PNP in the CO2-bridged diatomic complex enhances CO2 activation by another ∼20%, partly through the inductive effect of Fe(II), which pulls electron density from Ni-PNP across the CO2-bridge and partly by the backward donation from Fe(II)-PNP. Therefore, the present study encourages us to design a strongly electron-donating ligand and a CO2-bridged diatomic complex to develop more efficient homogeneous catalyst.

Frontiers in Multiscale Modeling of Photoreceptor Proteins.

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.