Search

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.

Photoactivation in the Photoactive Yellow Protein, a bacterial blue light photoreceptor, proceeds via photo-isomerization of the double C=C bond in the covalently attached chromophore. Quantum chemistry calculations, however, have suggested that in addition to double bond photo-isomerization, the isolated chromophore and many of its analogues, can isomerize around a single C-C bond as well. Whereas double bond photo-isomerization has been observed with x-ray crystallography, experimental evidence for single bond photo-isomerization 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 non-adiabatic molecular dynamics simulations to reveal that the locked chromophore isomerizes around the single bond upon photo-activation. Our work thus provides experimental evidence for single bond photo-isomerization in a photoactive yellow protein chromophore analogue and suggests that photo-isomerization is not restricted to the double bonds in conjugated systems. This insight may be useful for designing light-driven molecular switches or motors. peerReviewed

The tetrapyrrole chromophore biliverdin IXα (BV) in the bacteriophytochrome from Deinococcus radiodurans (DrBphP) is usually assumed to be fully protonated, but this assumption has not been systematically validated by experiments or extensive computations. Here, we use force field molecular dynamics simulations and quantum mechanics/molecular mechanics calculations with density functional theory and XMCQDPT2 methods to investigate the effect of the five most probable protonation forms of BV on structural stability, binding pocket interactions, and absorption spectra in the two photochromic states of DrBphP. While agreement with X-ray structural data and measured UV/vis spectra suggest that in both states the protonated form of the chromophore dominates, we also find that a minor population with a deprotonated D-ring could contribute to the red-shifted tail in the absorption spectra.

Proton diffusion (PD) across biological membranes is a fundamental process in many biological systems, and much experimental and theoretical effort has been employed for deciphering it. Here, we report on a spectroscopic probe, which can be tightly tethered to the membrane, for following fast (nanosecond) proton transfer events on the surface of membranes. Our probe is composed of a photoacid that serves as our light-induced proton source for the initiation of the PD process. We use our probe to follow PD, and its pH dependence, on the surface of lipid vesicles composed of a zwitterionic headgroup, a negative headgroup, a headgroup that is composed only from the negative phosphate group, or a positive headgroup without the phosphate group. We reveal that the PD kinetic parameters are highly sensitive to the nature of the lipid headgroup, ranging from a fast lateral diffusion at some membranes to the escape of protons from surface to bulk (and vice versa) at others. By referring to existing theoretical models for membrane PD, we found that while some of our results confirm the quasi-equilibrium model, other results are in line with the nonequilibrium model. peerReviewed

When photoactive molecules interact strongly with confined light modes, new hybrid light–matter states may form: the polaritons. These polaritons are coherent superpositions of excitations of the molecules and of the cavity photon. Recently, polaritons were shown to mediate energy transfer between chromophores at distances beyond the Forster limit. Here we explore the potential of strong coupling for light-harvesting applications by means of atomistic molecular dynamics simulations of mixtures of photoreactive and non-photo-reactive molecules strongly coupled to a single confined light mode. These molecules are spatially separated and present at different concentrations. Our simulations suggest that while the excitation is initially fully delocalized over all molecules and the confined light mode, it very rapidly localizes onto one of the photoreactive molecules, which then undergoes the reaction.

When the enhanced electromagnetic field of a confined light mode interacts with photoactive molecules, the system can be driven into the regime of strong coupling, where new hybrid light-matter states, polaritons, are formed. Polaritons, manifested by the Rabi split in the dispersion, have shown potential for controlling the chemistry of the coupled molecules. Here, we show by angle-resolved steady-state experiments accompanied by multi-scale molecular dynamics simulations that the molecular Stokes shift plays a significant role in the relaxation of polaritons formed by organic molecules embedded in a polymer matrix within metallic Fabry-Pérot cavities. Our results suggest that in the case of Rhodamine 6G, a dye with a significant Stokes shift, excitation of the upper polariton leads to a rapid localization of the energy into the fluorescing state of one of the molecules, from where the energy scatters into the lower polariton (radiative pumping), which then emits. In contrast, for excitonic J-aggregates with a negligible Stokes shift, the fluorescing state does not provide an efficient relaxation gateway. Instead, the relaxation is mediated by exchanging energy quanta matching the energy gap between the dark states and lower polariton into vibrational modes (vibrationally assisted scattering). To understand better how the fluorescing state of a molecule that is not strongly coupled to the cavity can transfer its excitation energy to the lower polariton in the radiative pumping mechanism, we performed multi-scale molecular dynamics simulations. The results of these simulations suggest that non-adiabatic couplings between uncoupled molecules and the polaritons are the driving force for this energy transfer process.

It is included two versions of this item: the Accepted Version which is already Open Access and the Published Version which is under an embargo period till 2022-03-09., 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, This work was supported by the Academy of Finland (Grant No. 323996 to G.G.) as well as the European Research Council (Grant No. ERC-2016-StG-714870 to J.F.) and the Spanish Ministry for Science, Innovation, and Universities—AEI (Grant No. RTI2018-099737-B-I00 to J.F.). We thank the Center for Scientific Computing (CSC-IT Center for Science) for generous computational resources

Ligand-stabilized, atomically precise gold nanoclusters with a metal core of a uniform size of just 1–3 nm constitute an interesting class of nanomaterials with versatile possibilities for applications due to their size-dependent properties and modifiable ligand layers. The key to extending the usability of the clusters in applications is to understand the chemical bonding in the ligand layer as a function of cluster size and ligand structure. Previously, it has been shown that monodispersed gold nanoclusters, stabilized by meta-mercaptobenzoic acid (m-MBA or 3-MBA) ligands and with sizes of 68–144 gold atoms, show ambient stability. Here we show that a combination of nuclear magnetic resonance spectroscopy, UV–vis absorption, infrared spectroscopy, molecular dynamics simulations, and density functional theory calculations reveals a distinct chemistry in the ligand layer, absent in other known thiol-stabilized gold nanoclusters. Our results imply a low-symmetry C1 ligand layer of 3-MBA around the gold cor...

Journal of the American Chemical Society 138(51), 16620 - 16631 (2016). doi:10.1021/jacs.6b05609, There has been considerable debate on the existence of a low-barrier hydrogen bond (LBHB) in the photoactive yellow protein (PYP). The debate was initially triggered by the neutron diffraction study of Yamaguchi et al. ( Proc. Natl. Acad. Sci., U. S. A., 2009, 106, 440−444) who suggested a model in which a neutral Arg52 residue triggers the formation of the LBHB in PYP. Here, we present an alternative model that is consistent within the error margins of the Yamaguchi structure factors. The model explains an increased hydrogen bond length without nuclear quantum effects and for a protonated Arg52. We tested both models by calculations under crystal, solution, and vacuum conditions. Contrary to the common assumption in the field, we found that a single PYP in vacuum does not provide an accurate description of the crystal conditions but instead introduces strong artifacts, which favor a LBHB and a large 1H NMR chemical shift. Our model of the crystal environment was found to stabilize the two Arg52 hydrogen bonds and crystal water positions for the protonated Arg52 residue in free MD simulations and predicted an Arg52 pKa upshift with respect to PYP in solution. The crystal and solution environments resulted in almost identical 1H chemical shifts that agree with NMR solution data. We also calculated the effect of the Arg52 protonation state on the LBHB in 3D nuclear equilibrium density calculations. Only the charged crystal structure in vacuum supports a LBHB if Arg52 is neutral in PYP at the previously reported level of theory ( J. Am. Chem. Soc., 2014, 136, 3542−3552). We attribute the anomalies in the interpretation of the neutron data to a shift of the potential minimum, which does not involve nuclear quantum effects and is transferable beyond the Yamaguchi structure., Published by American Chemical Society, Washington, DC

Nitroxide labels are combined with nucleic acid structures and studied using electron paramagnetic resonance experiments (EPR). As X-ray/NMR structures are unavailable with the nitroxide labels, detailed residue level information, down to atomic resolution, about the effect of these nitroxide labels on local RNA structures is currently lacking. This information is critical to evaluate the choice of spin label. In this study, we compare and contrast the effect of TEMPO-based (NT) and rigid spin (Ç) labels (in both 2’-O methylated and not-methylated forms) on RNA duplexes. We also investigate sequence-dependent effects of NT label on RNA duplex along with the more complex G-quadruplex RNA. Distances measured from molecular dynamics simulations between the two spin labels are in agreement with the EPR experimental data. To understand the effect of labeled oligonucleotides on the structure, we studied the local base pair geometries and global structure in comparison with the unlabeled structures. Based on the structural analysis, we can conclude that TEMPO-based and Ç labels do not significantly perturb the base pair arrangements of the native oligonucleotide. When experimental structures for the spin labelled DNA/RNA molecules are not available, general framework offered by the current study can be used to provide information critical to the choice of spin labels to facilitate future EPR studies.Graphical abstract

In constant pH molecular dynamics simulations, the protonation states of titratable sites can respond to changes of the pH and of their electrostatic environment. Consequently, the number of protons bound to the biomolecule, and therefore the overall charge of the system, fluctuates during the simulation. To avoid artifacts associated with a non-neutral simulation system, we introduce an approach to maintain neutrality of the simulation box in constant pH molecular dynamics simulations, while maintaining an accurate description of all protonation fluctuations. Specifically, we introduce a proton buffer that, like a buffer in experiment, can exchange protons with the biomolecule enabling its charge to fluctuate. To keep the total charge of the system constant, the uptake and release of protons by the buffer are coupled to the titration of the biomolecule with a constraint. We find that, because the fluctuation of the total charge (number of protons) of a typical biomolecule is much smaller than the number of titratable sites of the biomolecule, the number of buffer sites required to maintain overall charge neutrality without compromising the charge fluctuations of the biomolecule, is typically much smaller than the number of titratable sites, implying markedly enhanced simulation and sampling efficiency.

The pKa of the p-mercaptobenzoic acid (pMBA) ligands in the Au102(pMBA)44 nanocluster was measured by using acid–base and IR titration. The observed macroscopic pKa = 6.18 ± 0.05 is significantly more basic than that of free pMBA (pKa = 4.16), and the protonation behavior is anticooperative according to the Hill coefficient n = 0.64 ± 0.04. The cluster is truly water-soluble when more than 22 and insoluble when fewer than 7 ligands are in the deprotonated state. In order to obtain more insight into the anticooperative character, the cluster was modeled at pH ∼6.2 using constant pH molecular dynamics simulations. The pKa values of the individual pMBAs are in the range of 5.18–7.58, depending on the positions of the ligands, generating a nonuniform charge distribution on the cluster surface. The cluster interactions with inorganic and biological environment depend on the surface charge distribution, stressing the importance of understanding the protonation state of the cluster.

Hydration shells around ions are crucial for many fundamental biological and chemical processes. Their local physicochemical properties are quite different from those of bulk water and hard to probe experimentally. We address this problem by combining soft X-ray spectroscopy using a liquid jet and molecular dynamics (MD) simulations together with ab initio electronic structure calculations to elucidate the water–ion interaction in a MgCl2 solution at the molecular level. Our results reveal that salt ions mainly affect the electronic properties of water molecules in close vicinity and that the oxygen K-edge X-ray emission spectrum of water molecules in the first solvation shell differs significantly from that of bulk water. Ion-specific effects are identified by fingerprint features in the water X-ray emission spectra. While Mg2+ ions cause a bathochromic shift of the water lone pair orbital, the 3p orbital of the Cl– ions causes an additional peak in the water emission spectrum at around 528 eV. peerReviewed

We present transferable AMBER-compatible force field parameters for thiolate-protected gold nanoclusters. Five different sized clusters containing both organo-soluble and water-soluble thiolate ligands served as test systems in MD simulations, and parameters were validated against DFT and experimental results. The cluster geometries remain intact during the MD simulations in various solvents, and structural fluctuations and energetics showed agreement with DFT calculations. Experimental diffusion coefficients and crystal structures were also reproduced with sufficient accuracy. The presented parameter set contains the minimum number of cluster-specific parameters enabling the use of these parameters for several different gold nanoclusters. The parameterization of ligands can also be extended to different types of ligands.

The first excited state of the Photoactive Yellow Protein chromophore exhibits a strong charge transfer character and the dipole moments of the excited and ground states differ significantly. Furthermore, the excited state charge distribution changes during the isomerization of this chromophore. These observations suggest that external electric fields can be used to control photo-isomerization, providing a new concept for developing photochromic devices, such as e-paper or optical memory. To test this idea, we performed excited state dynamics simulations and static calculations of a PYP chromophore analogue (pCK − ) in an external electric field. By adjusting direction and strength of the field, we were able to control both selectivity ( i.e., which bond isomerizes) and efficiency (excited state lifetime) of the isomerization process. Simulations like these could potentially be useful in validating whether a given chromophore could be suitable for electrochromic applications.

[Introduction] The concentration of hydronium and hydroxide at the water-air interface has been under debate for a long time. Recent evidence from a range of experiments and theoretical calculations strongly suggests the water surface is somewhat acidic. Using novel polarizable models we have performed potential of mean force calculations of a hydronium ion, a hydroxide ion and a water molecule in a water droplet and a water slab and we were able to rationalize that hydronium, but not hydroxide, is slightly enriched at the surface for two reasons. First, because the hydrogen-bond acceptance capacity of hydronium is weaker than water it is more favorable to have the hydronium oxygen on the surface. Second, hydroxide ions are expelled from the surface of droplets, due to the entropy being lower when a hydroxide ion is hydrated on the surface. As a result, the water dissociation constant pK w increases slightly near the surface. The results are corroborated by calculations of surface tension of NaOH solutions that are in reasonable agreement with experiment. The structural and thermodynamic interpretation of hydronium and hydroxide hydration provided by these calculations opens the route to a better understanding of atmospheric- and surface chemistry. peerReviewed

The complete active space self-consistent field (CASSCF) method is one of the most used approaches to explore theoretically non-adiabatic photochemical processes in medium-size molecular systems. It is also frequently utilized in direct dynamics studies because this method offers a reliable description of the electronic state mixing in the regions of surface crossings (conical intersections). In this context, however, the large computational demand required to perform such on-the-fly molecular dynamics simulations using CASSCF as the electronic structure method means that only relatively small active space can be handled (typically not more than 10 electrons distributed in 10 orbitals). In this article, we present some guidelines to reduce the active space size in CASSCF calculations in order to decrease the computational cost necessary to perform, for example, molecular dynamics simulations. Three important criteria are considered: (i) the degradation of the accuracy in the description of the potential energy surface following the reduction of the active space, (ii) the stability of the reduced active space across the global potential energy surface, and (iii) the cost of the CASSCF calculation using the reduced active space. Two illustrations are given. The first one is the benzene to benzvalene S 1 photochemical pathway, the second one is the photoisomerization of thiomethyl para -coumaric acid, the chromophore in photoactive yellow protein.

Carbonic anhydrases (CAs) are enzymes whose endogenous reaction is the reversible hydration of CO(2) to give HCO(3)(-) and a proton. CA are also known to exhibit weak and promiscuous esterase activity toward activated esters. Here, we report a series of findings obtained with a set of CA inhibitors that showed quite unexpectedly that the compounds were both inhibitors of CO(2) hydration and substrates for the esterase activity of CA. The compounds comprised a monosaccharide core with the C-6 primary hydroxyl group derivatized as a sulfamate (for CA recognition). The remaining four sugar hydroxyl groups were acylated. Using protein X-ray crystallography, the crystal structures of human CA II in complex with four of the sulfamate inhibitors were obtained. As expected, the four structures displayed the canonical CA protein-sulfamate interactions. Unexpectedly, a free hydroxyl group was observed at the anomeric center (C-1) rather than the parent C-1 acyl group. In addition, this hydroxyl group is observed axial to the carbohydrate ring while in the parent structure it is equatorial. A mechanism is proposed that accounts for this inversion of stereochemistry. For three of the inhibitors, the acyl groups at C-2 or at C-2 and C-3 were also absent with hydroxyl groups observed in their place and retention of stereochemistry. With the use of electrospray ionization-Fourier transform ion cyclotron resonance-mass spectrometry (ESI-FTICR-MS), we observed directly the sequential loss of all four acyl groups from one of the carbohydrate-based sulfamates. For this compound, the inhibitor and substrate binding mode were further analyzed using free energy calculations. These calculations suggested that the parent compound binds almost exclusively as a substrate. To conclude, we have demonstrated that acylated carbohydrate-based sulfamates are simultaneously inhibitor and substrate of human CA II. Our results suggest that, initially, the substrate binding mode dominates, but following hydrolysis, the ligand can also bind as a pure inhibitor thereby competing with the substrate binding mode.

International audience; We have performed ab initio CASSCF, CASPT2, and EOM-CCSD calculations on doubly deprotonated p-coumaric acid (pCA2−), the chromophore precursor of the photoactive yellow protein. The results of the calculations demonstrate that pCA2− can undergo only photoisomerization of the double bond. In contrast, the chromophore derivative with the acid replaced by a ketone (p-hydroxybenzylidene acetone, pCK−) undergoes both single- and double-bond photoisomerization, with the single-bond relaxation channel more favorable than the double-bond channel. The substitution alters the nature of the first excited states and the associated potential energy landscape. The calculations show that the electronic nature of the first two (π,π*) excited states are interchanged in vacuo due to the substitution. In pCK−, the first excited state is a charge-transfer (CT π,π*) state, in which the negative charge has migrated from the phenolate ring onto the alkene tail of the chromophore, whereas the locally excited (LE π,π*) state, in which the excitation involves the orbitals on the phenol ring, lies higher in energy and is the fourth excited state. In pCA2−, the CT state is higher in energy due the presence of a negative charge on the tail of the chromophore, and the first excited state is the LE state. In isolated pCA2−, there is a 68 kJ/mol barrier for double-bond photoisomerization on the potential energy surface of this LE state. In water, however, hydrogen bonding with water molecules reduces this barrier to 9 kJ/mol. The barrier separates the local trans minimum near the Franck−Condon region from the global minimum on the excited-state potential energy surface. The lowest energy conical intersection was located near this minimum. In contrast to pCK−, single-bond isomerization is highly unfavorable both in the LE and CT states of pCA2−. These results demonstrate that pCA2− can only decay efficiently in water and exclusively by double-bond photoisomerization. These findings provide a rationale for the experimental observations that pCA2− has both a longer excited-state lifetime and a higher isomerization quantum yield than pCK−.

International audience; We have performed excited-state dynamics simulations of a Photoactive Yellow Protein chromophore analogue in water. The results of the simulations demonstrate that in water the chromophore predominantly undergoes single-bond photoisomerization, rather than double-bond photoisomerization. Despite opposite charge distributions in the chromophore, excited-state decay takes place very efficiently from both single- and double-bond twisted minima in water. Radiationless decay is facilitated by ultrafast solvent reorganization, which stabilizes both minima by specific hydrogen bond interactions. Changing the solvent to the slightly more viscous D(2)O leads to an increase of the excited-state lifetime. Together with previous simulations, the present results provide a complete picture of the effect of the protein on the photoisomerization of the chromophore in PYP: the positive guanidinium group of Arg52 favors double-bond isomerization over single-bond isomerization by lowering the barrier for double-bond isomerization, while the hydrogen bonds with Tyr42 and Glu46 enhance deactivation from the double-bond twisted minimum.

Organisms have evolved a wide variety of mechanisms to utilize and respond to light. In many cases, the biological response is mediated by structural changes that follow photon absorption. These reactions typically occur at femto- to picosecond timescales. As the relevant time and spatial resolutions are notoriously hard to access experimentally, molecular dynamics (MD) simulations are the method of choice to study such ultrafast processes. In the simulations, a multiconfigurational quantum mechanical (QM) description (CASSCF, CASPT2) is required to model the electronic rearrangement of those parts of the system that are involved in the absorption. For the remainder, typically consisting of the apoprotein and the solvent, a simple forcefield model (MM) suffices. QM/MM gradients have to be computed on-the-fly, and surface hopping procedures are needed to model the excited state decay. In this chapter, the computational framework underlying the atomistic simulation of photochemical events is reviewed and a few representative applications are discussed that demonstrate the validity of hybrid QM/MM approaches for photobiological reactions. Keywords: excited state dynamics; force field; molecular mechanics; molecular dynamics; molecular dynamics

We have performed excited-state dynamics simulations of the Arg52Gln (R52Q) mutant of photoactive yellow protein (PYP). The results of these simulations demonstrate that in the mutant the primary events after photoexcitation are different from those in the wild-type. In the mutant, the chromophore predominantly undergoes single bond photoisomerization, whereas in the wild-type, photoisomerization of the double bond occurs. Furthermore, the excited-state lifetime is around three times longer than in wild-type PYP, which agrees well with recent transient absorption measurements. In 20% of the trajectories, we observe the formation of a photoproduct that has the carbonyl oxygen atom of the chromophore flipped by almost 180°, disrupting the hydrogen bond between the chromophore and the backbone amino group of Cys69. This observation, in combination with the fact that the mutant is photoactive, suggests that the break of the hydrogen bond is the key step in the photoactivation process rather than the double bo...

Fluorescent proteins have been widely used as genetically encodable fusion tags for biological imaging. Recently, a new class of fluorescent proteins was discovered that can be reversibly light-switched between a fluorescent and a non-fluorescent state. Such proteins can not only provide nanoscale resolution in far-field fluorescence optical microscopy much below the diffraction limit, but also hold promise for other nanotechnological applications, such as optical data storage. To systematically exploit the potential of such photoswitchable proteins and to enable rational improvements to their properties requires a detailed understanding of the molecular switching mechanism, which is currently unknown. Here, we have studied the photoswitching mechanism of the reversibly switchable fluoroprotein asFP595 at the atomic level by multiconfigurational ab initio (CASSCF) calculations and QM/MM excited state molecular dynamics simulations with explicit surface hopping. Our simulations explain measured quantum yields and excited state lifetimes, and also predict the structures of the hitherto unknown intermediates and of the irreversibly fluorescent state. Further, we find that the proton distribution in the active site of the asFP595 controls the photochemical conversion pathways of the chromophore in the protein matrix. Accordingly, changes in the protonation state of the chromophore and some proximal amino acids lead to different photochemical states, which all turn out to be essential for the photoswitching mechanism. These photochemical states are (i) a neutral chromophore, which can trans-cis photoisomerize, (ii) an anionic chromophore, which rapidly undergoes radiationless decay after excitation, and (iii) a putative fluorescent zwitterionic chromophore. The overall stability of the different protonation states is controlled by the isomeric state of the chromophore. We finally propose that radiation-induced decarboxylation of the glutamic acid Glu215 blocks the proton transfer pathways that enable the deactivation of the zwitterionic chromophore and thus leads to irreversible fluorescence. We have identified the tight coupling of trans-cis isomerization and proton transfers in photoswitchable proteins to be essential for their function and propose a detailed underlying mechanism, which provides a comprehensive picture that explains the available experimental data. The structural similarity between asFP595 and other fluoroproteins of interest for imaging suggests that this coupling is a quite general mechanism for photoswitchable proteins. These insights can guide the rational design and optimization of photoswitchable proteins., Author Summary Proteins whose fluorescence can be reversibly switched on and off hold great promise for applications in high-resolution optical microscopy and nanotechnology. To systematically exploit the potential of such photoswitchable proteins and to enable rational improvements of their properties requires a detailed understanding of the molecular switching mechanism. Here, we have studied the photoswitching mechanism of the reversibly switchable fluoroprotein asFP595 by atomistic molecular dynamics simulations. Our simulations explain measured quantum yields and excited state lifetimes, and also predict the structures of the hitherto unknown intermediates and of the irreversibly fluorescent state. Further, we find that the proton distribution in the active site of the asFP595 controls the photochemical conversion pathways of the chromophore in the protein matrix. Our results show that a tight coupling between trans-cis isomerization of the chromophore and proton transfer is essential for the function of asFP595. The structural similarity between asFP595 and other fluoroproteins suggests that this coupling is a quite general mechanism for photoswitchable proteins. These insights can guide the rational design and optimization of photoswitchable proteins.

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.

Modified and unmodified yeast tRNA(Phe) in solution was simulated to understand the effect of nucleotide modifications on the dynamics of tRNA. High performance computing techniques were employed to obtain a “dynamic picture” at spacial and time resolutions hardly accessible experimentally for these systems. Local flexibility, secondary and tertiary structure, global rearrangements and movements of the whole tRNA were probed by a microsecond of all atom explicit solvent molecular dynamics simulations.The results of our simulations give new insight on experimentally observed biological impact of nucleotide modifications. Amongst other results we find that tRNA modification leads to a decrease in secondary structure and tertiary interactions in the anticodon stem loop. Nucleotides U16 and U17 show different orientations when modified to dihydrouridine. Modification of A58 to 1-methyladenosine causes local rearrangements in the elbow region. Implications of our results on external factor binding are discussed.