Alternant Hydrocarbon Diradicals as Optically Addressable Molecular Qubits.

High-spin molecules allow for bottom-up qubit design and are promising platforms for magnetic sensing and quantum information science. Optical addressability of molecular electron spins has also been proposed in first-row transition-metal complexes via optically detected magnetic resonance (ODMR) mechanisms analogous to the diamond-nitrogen-vacancy color center. However, significantly less progress has been made on the front of metal-free molecules, which can deliver lower costs and milder environmental impacts. At present, most luminescent open-shell organic molecules are π-diradicals, but such systems often suffer from poor ground-state open-shell characters necessary to realize a stable ground-state molecular qubit. In this work, we use alternancy symmetry to selectively minimize radical-radical interactions in the ground state, generating π-systems with high diradical characters. We call them m-dimers, referencing the need to covalently link two benzylic radicals at their meta carbon atoms for the desired symmetry. Through a detailed electronic structure analysis, we find that the excited states of alternant hydrocarbon m-diradicals contain important symmetries that can be used to construct ODMR mechanisms leading to ground-state spin polarization. The molecular parameters are set in the context of a tris(2,4,6-trichlorophenyl)methyl (TTM) radical dimer covalently tethered at the meta position, demonstrating the feasibility of alternant m-diradicals as molecular color centers.

Non-Hermitian molecular dynamics simulations of exciton-polaritons in lossy cavities.

The observation that materials can change their properties when placed inside or near an optical resonator has sparked a fervid interest in understanding the effects of strong light-matter coupling on molecular dynamics, and several approaches have been proposed to extend the methods of computational chemistry into this regime. Whereas the majority of these approaches have focused on modeling a single molecule coupled to a single cavity mode, changes to chemistry have so far only been observed experimentally when very many molecules are coupled collectively to multiple modes with short lifetimes. While atomistic simulations of many molecules coupled to multiple cavity modes have been performed with semi-classical molecular dynamics, an explicit description of cavity losses has so far been restricted to simulations in which only a very few molecular degrees of freedom were considered. Here, we have implemented an effective non-Hermitian Hamiltonian to explicitly treat cavity losses in large-scale semi-classical molecular dynamics simulations of organic polaritons and used it to perform both mean-field and surface hopping simulations of polariton relaxation, propagation, and energy transfer.

phbuilder: A Tool for Efficiently Setting up Constant pH Molecular Dynamics Simulations in GROMACS.

Constant pH molecular dynamics (MD) is a powerful technique that allows the protonation state of residues to change dynamically, thereby enabling the study of pH dependence in a manner that has not been possible before. Recently, a constant pH implementation was incorporated into the GROMACS MD package. Although this implementation provides good accuracy and performance, manual modification and the preparation of simulation input files are required, which can be complicated, tedious, and prone to errors. To simplify and automate the setup process, we present phbuilder, a tool that automatically prepares constant pH MD simulations for GROMACS by modifying the input structure and topology as well as generating the necessary parameter files. phbuilder can prepare constant pH simulations from both initial structures and existing simulation systems, and it also provides functionality for performing titrations and single-site parametrizations of new titratable group types. The tool is freely available at www.gitlab.com/gromacs-constantph. We anticipate that phbuilder will make constant pH simulations easier to set up, thereby making them more accessible to the GROMACS user community.

Tuning the Coherent Propagation of Organic Exciton-Polaritons through the Cavity Q-factor.

Transport of excitons in organic materials can be enhanced through polariton formation when the interaction strength between these excitons and the confined light modes of an optical resonator exceeds their decay rates. While the polariton lifetime is determined by the Q(uality)-factor of the optical resonator, the polariton group velocity is not. Instead, the latter is solely determined by the polariton dispersion. Yet, experiments suggest that the Q-factor also controls the polariton propagation velocity. To understand this observation, the authors perform molecular dynamics simulations of Rhodamine chromophores strongly coupled to Fabry-Pérot cavities with various Q-factors. The results suggest that propagation in the aforementioned experiments is initially dominated by ballistic motion of upper polariton states at their group velocities, which leads to a rapid expansion of the wavepacket. Cavity decay in combination with non-adiabatic population transfer into dark states, rapidly depletes these bright states, causing the wavepacket to contract. However, because population transfer is reversible, propagation continues, but as a diffusion process, at lower velocity. By controlling the lifetime of bright states, the Q-factor determines the duration of the ballistic phase and the diffusion coefficient in the diffusive regime. Thus, polariton propagation in organic microcavities can be effectively tuned through the Q-factor.

Multi-scale molecular dynamics simulations of enhanced energy transfer in organic molecules under strong coupling.

Exciton transport can be enhanced in the strong coupling regime where excitons hybridize with confined light modes to form polaritons. Because polaritons have group velocity, their propagation should be ballistic and long-ranged. However, experiments indicate that organic polaritons propagate in a diffusive manner and more slowly than their group velocity. Here, we resolve this controversy by means of molecular dynamics simulations of Rhodamine molecules in a Fabry-Pérot cavity. Our results suggest that polariton propagation is limited by the cavity lifetime and appears diffusive due to reversible population transfers between polaritonic states that propagate ballistically at their group velocity, and dark states that are stationary. Furthermore, because long-lived dark states transiently trap the excitation, propagation is observed on timescales beyond the intrinsic polariton lifetime. These insights not only help to better understand and interpret experimental observations, but also pave the way towards rational design of molecule-cavity systems for coherent exciton transport.

Optical control of ultrafast structural dynamics in a fluorescent protein.

The photoisomerization reaction of a fluorescent protein chromophore occurs on the ultrafast timescale. The structural dynamics that result from femtosecond optical excitation have contributions from vibrational and electronic processes and from reaction dynamics that involve the crossing through a conical intersection. The creation and progression of the ultrafast structural dynamics strongly depends on optical and molecular parameters. When using X-ray crystallography as a probe of ultrafast dynamics, the origin of the observed nuclear motions is not known. Now, high-resolution pump-probe X-ray crystallography reveals complex sub-ångström, ultrafast motions and hydrogen-bonding rearrangements in the active site of a fluorescent protein. However, we demonstrate that the measured motions are not part of the photoisomerization reaction but instead arise from impulsively driven coherent vibrational processes in the electronic ground state. A coherent-control experiment using a two-colour and two-pulse optical excitation strongly amplifies the X-ray crystallographic difference density, while it fully depletes the photoisomerization process. A coherent control mechanism was tested and confirmed the wave packets assignment.

gmXtal: Cooking Crystals with GROMACS.

Molecular dynamics (MD) simulations are routinely performed of biomolecules in solution, because this is their native environment. However, the structures used in such simulations are often obtained with X-ray crystallography, which provides the atomic coordinates of the biomolecule in a crystal environment. With the advent of free electron lasers and time-resolved techniques, X-ray crystallography can now also access metastable states that are intermediates in a biochemical process. Such experiments provide additional data, which can be used, for example, to optimize MD force fields. Doing so requires that the simulation of the biomolecule is also performed in the crystal environment. However, in contrast to simulations of biomolecules in solution, setting up a crystal is challenging. In particular, because not all solvent molecules are resolved in X-ray crystallography, adding a suitable number of solvent molecules, such that the properties of the crystallographic unit cell are preserved in the simulation, can be difficult and typically is a trial-and-error based procedure requiring manual interventions. Such interventions preclude high throughput applications. To overcome this bottleneck, we introduce gmXtal, a tool for setting up crystal simulations for MD simulations with GROMACS. With the information from the protein data bank (rcsb.org) gmXtal automatically (i) builds the crystallographic unit cell; (ii) sets the protonation of titratable residues; (iii) builds missing residues that were not resolved experimentally; and (iv) adds an appropriate number of solvent molecules to the system. gmXtal is available as a standalone tool https://gitlab.com/pbuslaev/gmxtal .

Serial Femtosecond Crystallography Reveals that Photoactivation in a Fluorescent Protein Proceeds via the Hula Twist Mechanism.

Chromophore cis/trans photoisomerization is a fundamental process in chemistry and in the activation of many photosensitive proteins. A major task is understanding the effect of the protein environment on the efficiency and direction of this reaction compared to what is observed in the gas and solution phases. In this study, we set out to visualize the hula twist (HT) mechanism in a fluorescent protein, which is hypothesized to be the preferred mechanism in a spatially constrained binding pocket. We use a chlorine substituent to break the twofold symmetry of the embedded phenolic group of the chromophore and unambiguously identify the HT primary photoproduct. Through serial femtosecond crystallography, we then track the photoreaction from femtoseconds to the microsecond regime. We observe signals for the photoisomerization of the chromophore as early as 300 fs, obtaining the first experimental structural evidence of the HT mechanism in a protein on its femtosecond-to-picosecond timescale. We are then able to follow how chromophore isomerization and twisting lead to secondary structure rearrangements of the protein ß-barrel across the time window of our measurements.

Best Practices in Constant pH MD Simulations: Accuracy and Sampling.

Various approaches have been proposed to include the effect of pH in molecular dynamics (MD) simulations. Among these, the λ-dynamics approach proposed by Brooks and co-workers [Kong, X.; Brooks III, C. L. J. Chem. Phys. 1996, 105, 2414-2423] can be performed with little computational overhead and hfor each typeence be used to routinely perform MD simulations at microsecond time scales, as shown in the accompanying paper [Aho, N. et al. J. Chem. Theory Comput. 2022, DOI: 10.1021/acs.jctc.2c00516]. At such time scales, however, the accuracy of the molecular mechanics force field and the parametrization becomes critical. Here, we address these issues and provide the community with guidelines on how to set up and perform long time scale constant pH MD simulations. We found that barriers associated with the torsions of side chains in the CHARMM36m force field are too high for reaching convergence in constant pH MD simulations on microsecond time scales. To avoid the high computational cost of extending the sampling, we propose small modifications to the force field to selectively reduce the torsional barriers. We demonstrate that with such modifications we obtain converged distributions of both protonation and torsional degrees of freedom and hence consistent pKa estimates, while the sampling of the overall configurational space accessible to proteins is unaffected as compared to normal MD simulations. We also show that the results of constant pH MD depend on the accuracy of the correction potentials. While these potentials are typically obtained by fitting a low-order polynomial to calculated free energy profiles, we find that higher order fits are essential to provide accurate and consistent results. By resolving problems in accuracy and sampling, the work described in this and the accompanying paper paves the way to the widespread application of constant pH MD beyond pKa prediction.

Scalable Constant pH Molecular Dynamics in GROMACS.

Molecular dynamics (MD) computer simulations are used routinely to compute atomistic trajectories of complex systems. Systems are simulated in various ensembles, depending on the experimental conditions one aims to mimic. While constant energy, temperature, volume, and pressure are rather straightforward to model, pH, which is an equally important parameter in experiments, is more difficult to account for in simulations. Although a constant pH algorithm based on the λ-dynamics approach by Brooks and co-workers [Kong, X.; Brooks III, C. L. J. Chem. Phys.1996, 105, 2414-2423] was implemented in a fork of the GROMACS molecular dynamics program, uptake has been rather limited, presumably due to the poor scaling of that code with respect to the number of titratable sites. To overcome this limitation, we implemented an alternative scheme for interpolating the Hamiltonians of the protonation states that makes the constant pH molecular dynamics simulations almost as fast as a normal MD simulation with GROMACS. In addition, we implemented a simpler scheme, called multisite representation, for modeling side chains with multiple titratable sites, such as imidazole rings. This scheme, which is based on constraining the sum of the λ-coordinates, not only reduces the complexity associated with parametrizing the intramolecular interactions between the sites but also is easily extendable to other molecules with multiple titratable sites. With the combination of a more efficient interpolation scheme and multisite representation of titratable groups, we anticipate a rapid uptake of constant pH molecular dynamics simulations within the GROMACS user community.

Controlling Exciton Propagation in Organic Crystals through Strong Coupling to Plasmonic Nanoparticle Arrays.

Exciton transport in most organic materials is based on an incoherent hopping process between neighboring molecules. This process is very slow, setting a limit to the performance of organic optoelectronic devices. In this Article, we overcome the incoherent exciton transport by strongly coupling localized singlet excitations in a tetracene crystal to confined light modes in an array of plasmonic nanoparticles. We image the transport of the resulting exciton-polaritons in Fourier space at various distances from the excitation to directly probe their propagation length as a function of the exciton to photon fraction. Exciton-polaritons with an exciton fraction of 50% show a propagation length of 4.4 µm, which is an increase by 2 orders of magnitude compared to the singlet exciton diffusion length. This remarkable increase has been qualitatively confirmed with both finite-difference time-domain simulations and atomistic multiscale molecular dynamics simulations. Furthermore, we observe that the propagation length is modified when the dipole moment of the exciton transition is either parallel or perpendicular to the cavity field, which opens a new avenue for controlling the anisotropy of the exciton flow in organic crystals. The enhanced exciton-polariton transport reported here may contribute to the development of organic devices with lower recombination losses and improved performance.

Identifying Vibrations that Control Non-adiabatic Relaxation of Polaritons in Strongly Coupled Molecule-Cavity Systems.

The strong light-matter coupling regime, in which excitations of materials hybridize with excitations of confined light modes into polaritons, holds great promise in various areas of science and technology. A key aspect for all applications of polaritonic chemistry is the relaxation into the lower polaritonic states. Polariton relaxation is speculated to involve two separate processes: vibrationally assisted scattering (VAS) and radiative pumping (RP), but the driving forces underlying these two mechanisms are not fully understood. To provide mechanistic insights, we performed multiscale molecular dynamics simulations of tetracene molecules strongly coupled to the confined light modes of an optical cavity. The results suggest that both mechanisms are driven by the same molecular vibrations that induce relaxation through nonadiabatic coupling between dark states and polaritonic states. Identifying these vibrational modes provides a rationale for enhanced relaxation into the lower polariton when the cavity detuning is resonant with specific vibrational transitions.

The Photocycle of Bacteriophytochrome Is Initiated by Counterclockwise Chromophore Isomerization.

Photoactivation of bacteriophytochrome involves a cis-trans photoisomerization of a biliverdin chromophore, but neither the precise sequence of events nor the direction of the isomerization is known. Here, we used nonadiabatic molecular dynamics simulations on the photosensory protein dimer to resolve the isomerization mechanism in atomic detail. In our simulations the photoisomerization of the D ring occurs in the counterclockwise direction. On a subpicosecond time scale, the photoexcited chromophore adopts a short-lived intermediate with a highly twisted configuration stabilized by an extended hydrogen-bonding network. Within tens of picoseconds, these hydrogen bonds break, allowing the chromophore to adopt a more planar configuration, which we assign to the early Lumi-R state. The isomerization process is completed via helix inversion of the biliverdin chromophore to form the late Lumi-R state. The mechanistic insights into the photoisomerization process are essential to understand how bacteriophytochrome has evolved to mediate photoactivation and to engineer this protein for new applications.

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers.

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. This study aims to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

Multi-scale dynamics simulations of molecular polaritons: The effect of multiple cavity modes on polariton relaxation.

Coupling molecules to the confined light modes of an optical cavity is showing great promise for manipulating chemical reactions. However, to fully exploit this principle and use cavities as a new tool for controlling chemistry, a complete understanding of the effects of strong light-matter coupling on molecular dynamics and reactivity is required. While quantum chemistry can provide atomistic insight into the reactivity of uncoupled molecules, the possibilities to also explore strongly coupled systems are still rather limited due to the challenges associated with an accurate description of the cavity in such calculations. Despite recent progress in introducing strong coupling effects into quantum chemistry calculations, applications are mostly restricted to single or simplified molecules in ideal lossless cavities that support a single light mode only. However, even if commonly used planar mirror micro-cavities are characterized by a fundamental mode with a frequency determined by the distance between the mirrors, the cavity energy also depends on the wave vector of the incident light rays. To account for this dependency, called cavity dispersion, in atomistic simulations of molecules in optical cavities, we have extended our multi-scale molecular dynamics model for strongly coupled molecular ensembles to include multiple confined light modes. To validate the new model, we have performed simulations of up to 512 Rhodamine molecules in red-detuned Fabry-Pérot cavities. The results of our simulations suggest that after resonant excitation into the upper polariton at a fixed wave vector, or incidence angle, the coupled cavity-molecule system rapidly decays into dark states that lack dispersion. Slower relaxation from the dark state manifold into both the upper and lower bright polaritons causes observable photo-luminescence from the molecule-cavity system along the two polariton dispersion branches that ultimately evolves toward the bottom of the lower polariton branch, in line with experimental observations. We anticipate that the more realistic cavity description in our approach will help to better understand and predict how cavities can modify molecular properties.

Covalent and non-covalent coupling of a Au102 nanocluster with a fluorophore: energy transfer, quenching and intracellular pH sensing.

Interactions between an atomically precise gold nanocluster Au102(p-MBA)44 (p-MBA = para mercaptobenzoic acid) and a fluorescent organic dye molecule (KU, azadioxatriangulenium) are studied. In solution, the constituents form spontaneously a weakly bound complex leading to quenching of fluorescence of the KU dye via energy transfer. The KU can be separated from the complex by lowering pH, leading to recovery of fluorescence, which forms a basis for an optical reversible pH sensor. However, the sensor is not a stable entity, which could be delivered inside cells. For this purpose, a covalently bound hybrid is synthesized by linking the KU dye to the ligand layer of the cluster via an ester bond. Covalent linking facilitates entry of the cluster-dye hybrids into cells via endocytosis. Inside cells, the hybrids accumulate in endosomes where Au102 releases its cargo via hydrolysis of the ester bond. Changes of the local pH inside endosomes regulate reversible fluorescence due to variations in the interactions between the Au102 cluster and the dye. This work presents a concept for delivering reporter molecules into cells by using atomically precise gold nanoclusters as carriers and paves the way for future developments of cluster-reporter sensors for in vivo measurements of e.g. absolute pH values or ion concentrations.

Ultrafast structural changes within a photosynthetic reaction centre.

Photosynthetic reaction centres harvest the energy content of sunlight by transporting electrons across an energy-transducing biological membrane. Here we use time-resolved serial femtosecond crystallography1 using an X-ray free-electron laser2 to observe light-induced structural changes in the photosynthetic reaction centre of Blastochloris viridis on a timescale of picoseconds. Structural perturbations first occur at the special pair of chlorophyll molecules of the photosynthetic reaction centre that are photo-oxidized by light. Electron transfer to the menaquinone acceptor on the opposite side of the membrane induces a movement of this cofactor together with lower amplitude protein rearrangements. These observations reveal how proteins use conformational dynamics to stabilize the charge-separation steps of electron-transfer reactions.

The primary structural photoresponse of phytochrome proteins captured by a femtosecond X-ray laser.

Phytochrome proteins control the growth, reproduction, and photosynthesis of plants, fungi, and bacteria. Light is detected by a bilin cofactor, but it remains elusive how this leads to activation of the protein through structural changes. We present serial femtosecond X-ray crystallographic data of the chromophore-binding domains of a bacterial phytochrome at delay times of 1 ps and 10 ps after photoexcitation. The data reveal a twist of the D-ring, which leads to partial detachment of the chromophore from the protein. Unexpectedly, the conserved so-called pyrrole water is photodissociated from the chromophore, concomitant with movement of the A-ring and a key signaling aspartate. The changes are wired together by ultrafast backbone and water movements around the chromophore, channeling them into signal transduction towards the output domains. We suggest that the observed collective changes are important for the phytochrome photoresponse, explaining the earliest steps of how plants, fungi and bacteria sense red light.

Plants adapt to the availability of light throughout their lives because it regulates so many aspects of their growth and reproduction. To detect the level of light, plant cells use proteins called phytochromes, which are also found in some bacteria and fungi. Phytochrome proteins change shape when they are exposed to red light, and this change alters the behaviour of the cell. The red light is absorbed by a molecule known as chromophore, which is connected to a region of the phytochrome called the PHY-tongue. This region undergoes one of the key structural changes that occur when the phytochrome protein absorbs light, turning from a flat sheet into a helix. Claesson, Wahlgren, Takala et al. studied the structure of a bacterial phytochrome protein almost immediately after shining a very brief flash of red light using a laser. The experiments revealed that the structure of the protein begins to change within a trillionth of a second: specifically, the chromophore twists, which disrupts its attachment to the protein, freeing the protein to change shape. Claesson, Wahlgren, Takala et al. note that this structure is likely a very short-lived intermediate state, which however triggers more changes in the overall shape change of the protein. One feature of the rearrangement is the disappearance of a particular water molecule. This molecule can be found at the core of many different phytochrome structures and interacts with several parts of the chromophore and the phytochrome protein. It is unclear why the water molecule is lost, but given how quickly this happens after the red light is applied it is likely that this disappearance is an integral part of the reshaping process. Together these events disrupt the interactions between the chromophore and the PHY-tongue, enabling the PHY-tongue to change shape and alter the structure of the phytochrome protein. Understanding and controlling this process could allow scientists to alter growth patterns in plants, such as crops or weeds.

Photoactive Yellow Protein Chromophore Photoisomerizes around a Single Bond if the Double Bond Is Locked.

Photoactivation in the Photoactive Yellow Protein, a bacterial blue-light photoreceptor, proceeds via photoisomerization of the double C═C bond in the covalently attached chromophore. Quantum chemistry calculations, however, have suggested that in addition to double-bond photoisomerization, the isolated chromophore and many of its analogues can isomerize around a single C-C bond as well. Whereas double-bond photoisomerization has been observed with X-ray crystallography, experimental evidence of single-bond photoisomerization is currently lacking. Therefore, we have synthesized a chromophore analogue, in which the formal double bond is covalently locked in a cyclopentenone ring, and carried out transient absorption spectroscopy experiments in combination with nonadiabatic molecular dynamics simulations to reveal that the locked chromophore isomerizes around the single bond upon photoactivation. Our work thus provides experimental evidence of single-bond photoisomerization in a photoactive yellow protein chromophore analogue and suggests that photoisomerization is not restricted to the double bonds in conjugated systems. This insight may be useful for designing light-driven molecular switches or motors.