Molecular Dynamic Indicators of the ... - ACS Publications

Subscriber access provided by University of Texas Libraries

Article

Molecular Dynamic Indicators of the Photoswitching Properties of GFPs Daryna Smyrnova, Benjamien Moeyaert, Servaas Michielssens, Johan Hofkens, Peter Dedecker, and Arnout Ceulemans J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.5b04826 • Publication Date (Web): 25 Aug 2015 Downloaded from http://pubs.acs.org on August 31, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Molecular Dynamic Indicators of the Photoswitching Properties of GFPs Daryna Smyrnova,† Benjamien Moeyaert,‡ Servaas Michielssens,† Johan Hofkens,‡ Peter Dedecker,‡ and Arnout Ceulemans∗,† Quantum Chemistry and Physical Chemistry Division, and Molecular Imaging and Photonics Division, Department of Chemistry, KU Leuven, Celestijnenenlaan 200F, 3001 Heverlee, Belgium E-mail: [email protected]

∗

To whom correspondence should be addressed Quantum Chemistry and Physical Chemistry Division, Department of Chemistry , KU Leuven, Celestijnenenlaan 200F, 3001 Heverlee, Belgium ‡ Molecular Imaging and Photonics Division, Department of Chemistry, KU Leuven, Celestijnenenlaan 200F, 3001 Heverlee, Belgium †

1

ACS Paragon Plus Environment


1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Reversibly photoswitchable fluorescent proteins (RSFPs) are highly useful probes for a range of applications including diffraction-unlimited fluorescence microscopy. It was previously shown that reversible photoswitching involves not only cis-trans isomerization and protonation-deprotonation of the chromophore, but also results in a marked difference in β-barrel flexibility. In this work, we performed flexibility profiling and functional mode analysis (FMA) using molecular dynamics calculations to study how the flexibility of the RSFP β-barrel influences the photoswitching properties of several fluorescent proteins. We also used Partial Least-Squared (PLS) FMA to detect promising mutation sites for the modulation of photoswitching properties of RSFPs. Our results show that the flexibility of RSFP does depend on its state with a systematically higher flexibility in the darkstate compared to the bright-state. In particular our method highlights the importance of Val157 in Dronpa, which upon mutation yields striking difference in the collective motions of the two mutants. Overall, we show that PLS-FMA yields information, complementary to static structures, that can guide rational design of fluorescent proteins.

Keywords RSFP, PLS, FMA

Introduction Fluorescent proteins (FPs) are 30kDa proteins with a β-barrel structure. After translation and folding, they autocatalytically form a chromophore, requiring only the presence of molecular oxygen. 1,2 This unique feature has made them particularly useful for in vivo fluorescence imaging. Their biophysical and spectroscopic properties have been modulated by mutating the amino acid sequence, resulting in brighter and more photostable fluorescent proteins, ranging over almost the entire visible spectrum. 3–9 2


Page 2 of 32

Page 3 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


An interesting class of FPs is the photoactivatable FPs (PAFPs). These labels can change their emissive properties on demand. For instance, reversibly photoswitchable FPs (RSFPs) such as Dronpa 10–12 are originally in a bright state, and can be repeatedly turned off with blue and back on with purple light. Other PAFPs such as the irreversibly photoconvertible FPs (PCFPs) originally exist in a bright green-emissive state, and can by application of a high flux of purple light be irreversibly converted to a red-emissive state. Examples of PCFPs are EosFP 13 and Dendra2. 14 Some recently developed proteins combine both reversible photoswitching and irreversible photoconversion, such as IrisFP, 15 NijiFP 16 and pcDronpa. 17 The interest in PAFPs is mainly driven by the development of diffraction-unlimited fluorescence microscopy techniques like PALM, RESOLFT and pcSOFI 18–24 which all heavily rely on photoswitching or photoconversion of FPs. With these superresolution techniques, an arbitrarily high resolution can in principle be achieved, although the results depend on both instrumental parameters and the properties of the labels that are used. 25 Ideally, RSFPs have two thermostable states that are spectrally well separated, have a high resistance to photobleaching and reliable switching rates. 26,27 The molecular mechanism underlying photoswitching has been the subject of intense study. In most cases the absorption band of the bright fluorescent form of RSFPs is consistent with the deprotonated, anionic state of the chromophore, while the dark, non-fluorescent state corresponds to the neutral form of the chromophore. From crystal structures of RSFPs in both states, it is also seen that the bright state chromophore is in a cis-conformation, while the dark state assumes a trans conformation due to a so-called hula-twist around the chromophore’s methylene bridge 12,28–34 or a rotation around a double bond 35–38 (see Fig. 1). Additionally, NMR studies of the bright and dark states of Dronpa suggest a different flexibility of these two states. 39 In this work we have investigated whether this difference in flexibility observed in Dronpa is a recurring theme in other switchable proteins, by investigating this effect in pairs of

3



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

closely-related switchable proteins that differ by only a few mutations yet exhibit different switching properties. We also have analyzed whether this could serve as a simple indicator of the photoswitching ability of the protein.

Figure 1: Structural and photophysical characteristics of Dronpa. (A) Photoswitching scheme of the Dronpa protein. Upon irradiation with 405 and 488 nm light, Dronpa switches between dark and bright states. (B) Representative views of the chromophore and its nearby residues. Structures of the dark and bright state are shown in gray [Protein Data Bank (PDB) ID code 2POX] and green (PDB ID code 2Z1O), respectively. 12,40 The CYG chromophore in Dronpa adapts distinct cis and trans conformations in the bright and dark state, respectively. (C) Numbering of the β-sheets in barrel-fold of the FP. In particular, we studied two sets of closely related proteins with different switching properties. First, we studied Dronpa and its mutants rsFastLime, 12 rsKame 41 and pcDronpa. 17 Although differing in only 1 to 5 residues, these proteins have markedly different photoswitching properties (see Table 1). In particular rsKame is slow-switching (has lower photoisomerization quantum yield in comparison to Dronpa) and rsFastLime is fast-switching (has 4


Page 4 of 32

Page 5 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


higher photo-isomerization quantum yield in comparison to Dronpa). In the case of pcDronpa only the photoswitching green species was considered for the study. Photo-isomerization quantum yields for all of the proteins presented in this study are summarized in the Table S1. We additionally studied EosFP and IrisFP (like in the case of pcDronpa we considered only the green photoswitchable species before the irreversible green-to-red photoconversion). These two proteins have a 76% sequence similarity compared to Dronpa and differ by only two residues (IrisFP = EosFP_F173S_F191L). The 191 residue is pointing out of the βbarrel and the mutation at this position is said to have no influence on the spectroscopy. 15 However, while EosFP is considered to be non-photoswitchable, IrisFP does photoswitch; its photoswitching properties are attributed to a higher spatial freedom of the chromophore in its pocket and the disruption of a Van der Waals connection between Met159 and hydroxybenzilidene moiety. 15,42 Moreover, the HYG chromophore of EosFP/IrisFP, pivotal for the green-to-red photoconversion property, can also be found in pcDronpa. Here, we examine the β-barrel flexibility of the dark and bright state of Dronpa and some of its slower and faster switching variants. By means of Functional Mode Analysis (FMA) we elucidate the positions of the aminoacids contributing the most to the difference in flexibility between the on- and off-state and then see how this contribution changes in the closest mutants. Table 1: Effect of several mutations on the photoswitching properties of Dronpa. Key mutation(s) V157G V157I V60A+C62H+N94S+N102I+E218G

Mutant rsFastLime 12 rsKame 41 pcDronpa 17

5

Properties Fast photoswitching Slow photoswitching Slow photoswitching,green-to-red photoconvertible



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Methods Chromophore parametrization For the parametrization of the chromophores in the cis-state, the input structures were taken from the PDB (Dronpa bright state: 2Z1O and IrisFP bright state: 2VVH). The trans-state structures were created manually using Avogadro. 43 This was done by rotating the hydroxyphenyl moiety of the chromophore around the chromophore’s methylene bridge. Further each chromophore was optimized in vacuum by high level QM method. Also each chromophore was capped by ACE (acetylated) and NME (methylamidated) residues following the standard charge derivation procedure for the amberFF9X. 44 For structure preparation, the parameters were determined from quantum chemical calculations at HF/6-31G* level using Gaussian09. 45 The ESP charges were determined with the R.E.D. tool, 46 which aided in automatizing the optimization of the structures and later on to map the atomic charges. The dihedral angle parameters for φ and θ were taken from literature 29 where they were on the B3LYP/6-31G* level. Parameter files obtained with AmberTools 47 were converted to Gromacs format using the Acepype script. 48 The structures of the chromophore and sidechains with the atomtypes and the atomic point charges can be found in Fig. S1.

Molecular dynamics simulations For the first part of the study, proteins in both the bright and dark state were taken from the PDB (bright structures as mentioned in the previous paragraph, Dronpa dark state: 2POX, IrisFP dark state: 2VVJ, pcDronpa green-on state: 4HQ8, pcDronpa green-off state: 4HQ9, EosFP green state: 1ZUX). Taking into account the importance of the crystallographic waters we have carefully examined the PDB files used in the simulation. In particular for Dronpa several of them are available. The X-ray structure with the highest resolution is 2Z1O, also we compared it to 2Z6X, 39 2IOV, 12 2IE2 49 structures and found that it includes all of the important crystallographic waters present in the other structures. It looks like the 6


Page 6 of 32

Page 7 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


2IOV chain A structure contains an extra water inside of the chromophore’s cavity, however this water is absent in chains B,C and D, hence we have mapped it as of lesser importance and avoided including it to our simulations. Comparison of these four structures (only chain A) can be found in Fig. S2. In order to examine transient states, we replaced, in the bright state structure, the bright state chromophore by the dark state chromophore and changed the protonation of His193 in the bright state’s crystal structure. All of the mutations were made using the Rosetta package. 50 MD simulations were performed with Gromacs 4.6.5 51 and the Amber99sb forcefield, 47 which was shown to perform well in the prediction of NMR results. 52 Particle mesh Ewald was used to include long-range interactions with a non-bonded cut-off set to 13.5 Å. Hydrogen bonds were constrained using the SHAKE algorithm. The time step was 2 fs and trajectory files were written out every 20 ps. The temperature was controlled using the velocity-rescale algorithm. In the bright state, His193 was protonated and the chromophore was in the anionic form, while in the dark state the chromophore was turned into the neutral form and His193 was neutral as well (with the proton on Nδ). 28 For the simulation, all of the crystallographic waters were maintained, since they play an important role in the chromophore surrounding. The charged protein was neutralized by Na+ and Cl- ions at a concentration of 0.15M and solved in a cubic TIP3P water box, which was 10 Å bigger (in each direction) than the dimensions of the protein. Each chromophore transformation and the mutations were made on the original crystal structure. Next, the system was minimized (steepest descent method followed by L-FBGS) and equilibrated with all the heavy atoms being restrained, at T=300K and 1bar (100ps). It was then passed to the production run with no constraints for 100 ns. For each protein, 5 independent equilibration and production runs were performed, resulting in 5×100 ns of statistically independent simulations.

7



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Flexibility estimation RMSF calculation For the estimation of protein flexibility, the RMSF of the backbone was measured. The first 20 ns were cut out of each trajectory, and the trajectories were aligned on the first frame. The average RMSF was calculated over 400 ns. PLS-based FMA calculations To probe which residues have an influence on the photoswitching properties, we performed a Partial Least-Squares based Functional Mode Analysis (PLS-FMA). 53 These calculations detect the principal collective mode which correlates with a given functional property of the protein, in our case photoswitching. Since the RSFPs assume either a bright or a dark state, we label these states as +1 and -1, respectively. PLS-FMA works by taking the trajectory of the simulation and using half of it as a training set and half as a model set. Hence we mixed the trajectory for the bright and dark protein making sure the training set contains both trajectories (i.e. first 200 ns from the bright protein trajectory, then 200 ns from the dark protein, then again bright and dark). An example of the correlation between the training and cross-validation model set can be seen in Fig. S3. This calculation resulted in ensembleweighted maximum correlated motions (ewMCM) trajectories for different pairs of proteins, which are discussed in the Results section. To avoid the problem of overfitting, we trained the model to up to 50 components (see Fig S2). Using the Pearson correlation coefficient we concluded that the best fit occurs for 5 components with a correlation coefficient of almost 94%, after which the correlation on the test set decreases due to overfitting.

8


Page 8 of 32

Page 9 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Results and Discussion RMSF as a criterion for the photoswitching property It has been reported that the flexibility of the FP β-barrel is dramatically affected by the photoswitching. 39 This finding, based upon NMR studies, states that the dark-state Dronpa has higher flexibility around the hydroxyphenyl ring of the chromophore and the neighboring β-strand. In particular, the aminoacids in front of the chromophore’s hydroxyphenyl ring, namely in β-sheets 7, 8, 9, 10 and 11 (see Fig. 1C), show broadening of the NMR signal due to a higher flexibility of this region, occurring on a timescale that is short compared to the NMR relaxation time. In the bright state, the chromophore is tethered to the β-barrel via a hydrogen bond with Ser142 and a stacking interaction with protonated His193. In the dark state, however, this connection is lost, allowing those and neighbouring aminoacids to gain more flexibility. Also, the hydrogen bonding between β strands 7 and 10 is weak, further adding to the dark state flexibility. Broadening of the spectral signatures of aminoacids 65-69, in the α-helix holding the chromophore and adjacent to the loop 3, was also observed. The results of the MD simulations and the RMSF flexibility analysis are represented in Fig. 2 for Dronpa and Fig. S4 and S5 for IrisFP and pcDronpa, respectively. From these data, it is clear that for all three proteins the RMSF of the dark state on average is indeed higher than that of the bright state. Taking a closer look, we can see that the increased flexibility of the dark state is situated in the β-sheets close to the chromophore’s hydroxyphenyl moiety, especially around the β-strands 7 and 10 and the α-helix (see Fig. 1C for the arrangement of the β-sheets). Furthermore, the Dronpa trajectories show that the higher flexibility is mostly due to the rearrangement of the His193, Arg66 and also sterical hindrance of the chromophore by the protein environment around residues Val157, Phe173 and Glu140. A noteworthy feature is the higher flexibility that is observed at position Tyr116 and Asn105. If we look into the structural rearrangement between the bright and dark state we arrive 9



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


to the conclusions observed earlier: 40,54 particularly strong H-bond between Ser142 and the chromophore in Dronpa bright and absence of such a bond in Dronpa dark. Also, Dronpa bright has a stable H-bond network formed around the chromophore through residues Glu212 – His193 – Glu144 – Arg66 (Fig. S6). In Dronpa dark this network is lost and the chromophore possesses higher flexibility, due to the isomerization of the chromophore and deprotonation of His193 28 (Fig. S7). Also due to the transient repulsion during the isomerization of the chromophore, Val 157 flips by 100 thus giving way to accommodate deprotonated ◦

His193. The latter serves as a H-bond acceptor for Ser142, although this bond is not very stable since Ser142 forms a second H-bond with the backbone atom of Glu140 as well, thus giving rise to flexibility in the β-sheets 7 and 10. Next, we compared EosFP, considered non-switchable, to IrisFP. In Fig S3 it can be seen that the average flexibility of the dark-state of IrisFP is only slightly higher compared to the bright state. In β-sheets 4, 8, and 9, adjacent to the Ser173 which introduced the switching compared to EosFP, the flexibility of the dark-state is comparable to the flexibility of the bright-state. This is in agreement with the explanation for the effect of the Phe173Ser mutation on the photoswitching properties, 15 namely the replacement of Phe173 by Ser, resulting in a conformational change of the Met159 residue and the creation of a larger cavity. Also less tight packing of the cavity around the chromophore allows higher flexibility for Ser173. That leads to higher flexibility of the chromophore which is in line with the observed lowering of the QY 55 (also see the comparison of the chromophore’s flexibilities in Fig. S8A). Correlation between the Phe173 conformation and the chromophore state in Dronpa was also observed by Moors et al. 29 This observation about the flexibility in the β-sheet 7 and 10 for Dronpa also holds for IrisFP, which has structural rearrangement very similar to Dronpa. We also studied pcDronpa and compared its dark state flexibility to Dronpa and IrisFP. The dark-state of pcDronpa has an even higher flexibility than it does in Dronpa and IrisFP (see Fig. S4). This can be explained by the combination of features of both the IrisFP and

11



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Dronpa environment. pcDronpa has a chromophore which is identical to IrisFP’s and the sequence and structure are very similar to Dronpa’s. From this, we can conclude that in Dronpa-like photoswitching proteins (where excitation light also brings the FP to the dark state) 56 the dark-state indeed possesses higher flexibility in β7-10 sheets. This most probably results from a lower stabilization of the trans-state by the protein environment (which also possibly contributes to a higher energy in comparison to the cis-state). This observation is in line with the conclusions made about IrisFP and EosFP. 15,55 In this study authors assign higher enthalpy to the dark state with the difference of 68 kJ/mol for IrisFP and 114 kJ/mol for EosFP. Another study on asFP595 (a protein which contrary to Dronpa is formed in the dark state and excitation light brings it to the bright state) by Mironov et al. presents the free energy difference upon thermal isomerization of the chromophore. The stable dark-state has a free energy difference of ≈ 25 kJ/mol in comparison to the bright-state. 38 IrisFP has a higher flexibility than Dronpa, but also a lower molecular brightness (27.9 mM-1 cm-1 and 72.2 mM-1 cm-1 respectively). 15,40 The same applies to pcDronpa which has a flexibility in the bright state similar to that of Dronpa as well as a similarly high molecular brightness (97.8 mM-1 cm-1 ). 17 We notice increased RMSF in rsKame in the bright-state and accordingly it also has lower brightness than Dronpa. 41 Hence there indeed seems to be a correlation between brightness and flexibility in the bright state (see Fig. S8B). Ultimately high flexibility could be one of the reasons of fluorescence absence in the dark-state.

Mapping collective protein motions associated with the photoswitching Although our RMSF calculations clearly show that the flexibility of the RSFP’s β-barrel depends on their state, these techniques lack the sensitivity to distinguish which aminoacids are crucial for triggering these differences. One way to identify the aminoacids contributing the most to the inflated flexibility in 12


Page 12 of 32

Page 13 of 32




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 32

most promising candidates to either improve the photoswitching property or to inhibit it. Our PLS analysis clearly distinguishes two separate regions of interest: the first one consists of β-sheets 7,8, 9 and 10 which are situated adjacent to the chromophore’s hydroxyphenyl ring, the second region being the α-helix holding the chromophore. Particularly in IrisFP, Ser173, responsible for the photoswitching properties, has a high ewMCM, just like Ile157 and to lesser extent His193. The high ewMCM around Gly155 is most likely the result of glycine’s lack of side chain. In Dronpa, the higher ewMCM values are found at residues Val157, Glu144 and Tyr116 while in pcDronpa Ser142, Met159 and His193 show high ewMCM values. These residues are known to be involved into the process of photoswitching and as can be seen from Fig. 3, are mostly located in β8. 56 Also, it is known that mutations in the α-helix influence the photophysical properties of RSFPs. In particular the A69T mutation in pcDronpa significantly lowers its brightness and alters its pKa, 17 while the V60A mutation in Dronpa produces a "fast-forming" variant with an increased brightness. 17 In our case aminoacids 65-68 also contribute to the difference in promoting motions of the protein in bright and dark state. The results obtained by the PLS-based FMA confirm that there is a significant difference between pcDronpa, Dronpa and IrisFP in loop 3 (see the loop numbering in the Fig.

1)

which seems to contribute the most to the photoswitching property change in Dronpa and pcDronpa, while in IrisFP this movement seems to be distributed more over the two other loops. This can be explained by the fact that in loop 2 in IrisFP, there is one extra Gly residue (Gly180) that is not present in in Dronpa and pcDronpa. This makes loop 2 longer and more flexible. Also, at position 73 (which is in loop 3, following the Dronpa numbering) IrisFP has Gln, while in Dronpa and pcDronpa it is Val, which does not form a stable Hbond with Asp184 (corresponding to Lys183 in case of IrisFP). It is interesting to note how conformational changes in the loops are connected with the chromophore’s conformational change. Loop 3 in particular is connected with the α-helix and residues 65-69 adjacent to the chromophore, so its higher flexibility in the trans-state may contribute to the higher

14


Page 15 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


displacement of the chromophore. The same aminoacids in the α-helix show broadening of the NMR signal in the dark state Dronpa. 39 Another aminoacid of interest in Dronpa is Tyr116, which is Asn in IrisFP and pcDronpa. Remarkably, this aminoacid was mutated from Tyr in pcDronpa, yielding pcDronpa2 which has an improved green-to-red photoconversion rate. The interaction of the aminoacid at position 116 with the chromophore is indeed different between Dronpa and pcDronpa/IrisFP. In case of Dronpa the CYG chromophore forms H-bonds with Tyr116 through a water molecule, while in pcDronpa and IrisFP the HYG chromophore interacts with Asn116 sidechain through two coordinated water molecules. The PLS-FMA clearly identifies several of the known aminoacids participating in the modulation of the photoswitching properties. Also it has added a new dynamic level to the study of the structural properties of the RSFPs. We elaborated on this concept by testing whether our technique can probe the effect of point mutations on the flexibility of RSFPs. The results are presented in the following section.

Single-point mutation effect on the flexibility of RSFP Although the previous section has proved that PLS-FMA can contribute to identifying the aminoacids which are determining the photoswitching properties we would like to fine-tune the interpretation of the results which we obtained from it. The easiest way to do so is to test the RMSF/FMA protocol on a single-point mutant set. Also by doing so we have systematized the difference in flexibility which was noted in the previous section and connected it to the photo-physical properties. From the previous experiments, we found that in the Dronpa-like switching RSFPs the dark state has a higher flexibility while the FMA analysis pointed out the key aminoacids which promote the ewMCM in dark or bright state, i.e. promote the difference in the motions in the two states. We now investigate how the flexibility of Dronpa changes in its fastswitching mutant rsFastLime and a slow-switching mutant rsKame. The basic explanation 15



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of the difference in the QY of the photo-isomerization seems straightforward: Gly157 is a small aminoacid hence it makes a bigger cavity around the chromophore and in such a way decreases the steric hindrance for a cis-trans isomerization of the chromophore. 12 In the case of rsKame, Leu157 is a slightly more bulky aminoacid in comparison to Val, hence it increases the sterical hindrance when the chromophore undergoes the isomerization. 41 However, comparison of the effect of the single V157G/L mutation on the flexibility in Dronpa, rsFastLime and rsKame in their bright and dark state does not indicate which of those 3 proteins would be a faster photo-switcher. Both "fast" and "slow" switchers have increased flexibility in the dark state in comparison to Dronpa (see results in Fig. S7), and rsKame shows higher flexibility in the bright-state as well. Apparently the switching mechanism itself may play a role and we therefore investigate if the speed of dark-to-bright switching is hindered by sterical arrangement of the aminoacids around the chromophore. We have checked whether the rate of isomerization may be dependent on the intermediate state when the protein environment is still not fully adapted to the chromophore. Recent experimental studies on the photoswitching mechanism of Dronpa by van Thor et al. as well as Yadav et al. 30,58 indicate that the first step for the photoswitching process is the isomerization of the chromophore which is followed by ground-state proton transfer. Here, following the notion of Andresen et al. about the sterical hindrance of the chromophore by Val157, 40 we have tried to mimic this situation which we will denote as a transient state. In the transient state, the chromophore is already in its neutral trans-state, while the protein is not yet fully adjusted and still corresponds to the bright-state. Also the His193 side chain is in its neutral form (proton on Nδ). Then we run the MD simulation where during the first 10-20 ns (see Fig. S9) the most important aminoacids like Arg66 and His193 change their conformation to the one which is closer to that of a dark-state protein. So although all of the events related to the photoswitching have already happened the protein is still not fully relaxed (here we are talking about the relaxation times taking more than microseconds) in relation to the crystal dark state.

16


Page 16 of 32

Page 17 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Mimicking a transient dark state will allow us to compare these results to the nonswitching proteins where the dark-state simply doesn’t exist so that we cannot make use of the available crystal structure (i.e. in EosFP). First, we checked whether the simulated transient structure adjusts to the dark crystal structure and we found that RMSD from the Dronpa dark crystal structure is ≈0.5 nm2 , while for original dark state the RMSD is ≈0.2 nm2 . We noticed that the concerted motion of His193 and Arg66 takes place during the first 10-15 ns of the simulation for Dronpa and rsFastLime, while for rsKame it takes more than 20 ns (see Fig S8). The flexibility of the backbone in the simulated transient states is represented in Fig. 4A. The results show us that in the transient dark-state, the mutated proteins’ (rsKame and rsFastLime) flexibility does differ in comparison to Dronpa. Even more, we observe a higher flexibility for rsKame (slow-switcher) in the transient dark state, while for rsFastLime the difference is almost not present, since the protein quickly converges to the natural dark-state structure. The resulting difference between the flexibility in the bright and transient dark state can be found in the Fig. 5. It is plotted against the on-off photo-isomerization quantum yield (QY), although this observation does not necessarily imply a direct causal relationship. Indeed rsFastLime has the highest QY but the smallest difference in the flexibility, while for rsKame the opposite holds place. Also if we look at the off-state relaxation half-times then for Dronpa it is 840 min and for rsFastLime – 8 min. 12 Although not being reported explicitly, Rosenbloom et al show that rsKame has a longer dwell time in the dark state. 41 All of this points to a higher energy difference between the dark and bright state in the slow-switching RSFP. The results of the FMA analysis on rsKame and rsFastLime can be found in Fig. 6. In the case of rsFastLime, the collective motions of the bright and dark states have very few differences, while in case of rsKame the amplitude of ewMCM is much higher, meaning that the bright and dark states have very different functional modes. The differences between collective motions arise from the difference in the energy landscape of the two states. The

17




Page 18 of 32

Page 19 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Figure 5: The quantum yield of the on-off isomerization plotted vs. the difference in the flexibility between the transient dark and the bright state. The quantum yields are taken from literature and the flexibility difference is calculated for the chromophore and the aminoacids situated in the β-sheets in front of it, namely 135-144, 150-159, 166-174, 186-194 and 206-212 following Dronpa numbering. transition state and fall into the second state. Hence we can say that if two states are close, the photoswitching will occur faster. This leads us to conclusion that there is a link between the difference in the collective motions and the speed of the photoswitching. However their connection would be rather phenomenological than direct. Combining the FMA and RMSF results in both the natural and transient dark states we come to two conclusions. First of all, in rsFastLime both the bright and dark state are very close energetically and are separated only by a low energy barrier. Secondly, for rsKame the bright state is energetically less favourable than in Dronpa, while its dark state is separated by a high energy barrier. Most probably in case of rsKame there is more than one transition state each corresponding to a change in a principal mode differentiating the bright and dark states. Another pair of close mutants we have looked at is IrisFP and EosFP. Although IrisFP is a mutant of EosFP, we have also tried to compare EosFP to the results of back-mutation S173F introduced into IrisFP. The combinations used for the analysis can be found in Fig. 19




Page 20 of 32

Page 21 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


similar to that of EosFP. Hence we calculated the flexibility of the transient state of IrisFP S173F. Its flexibility is higher than in the case of dark state IrisFP S173F and shows a behavior more similar to EosFP. If we look at the flexibility of EosFP in the bright and dark state then we can say that the latter one is substantially higher in the β-sheet 7, 8, 9 and 10, located around Phe173, Ile157 and Glu140, His193. All these residues are able to sterically hinder the chromophore’s isomerization. Also, the average flexibility of EosFP in the dark state is higher than that for any other RSFP. The difference between the bright and transient dark state for IrisFP and EosFP is shown on the Fig. 5. Like in the case with Dronpa it can be noticed that the difference is much higher for non-switching EosFP than for IrisFP. Combining the results of RMSF and FMA for single-point mutants of EosFP and Dronpa it becomes apparent that the effect of a single mutation can barely be seen when comparing it only in the existing crystal structures. PDB structures provide static picture where the protein environment can easily incorporate the mutation, but photo-isomerization is a highly dynamic process involving many intermediate states on the way from the bright to dark state. That is why transient states can play a role when trying to predict the effect of the mutation. Although they correspond only to the ground state, there might be a phenomenological link between transient-bright state difference and the QY of the photo-isomerization. In fact thermal relaxation is a ground-state process and the transient state can be contributing to the reduced photo-isomerization QY by constantly repopulating the bright state.

Conclusions In this paper we have investigated flexibility and functional mode sensitivity as dynamic indicators of photoswitching properties in Dronpa, IrisFP and several of their mutants. We found that for Dronpa and its photoconvertible mutant pcDronpa, the flexibility, as quantified by RMSF analysis, is considerably higher in the dark-state, while for IrisFP this effect

21



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

is limited to only two β-sheets close to the chromophore. Furthermore, we can see that the flexibility of EosFP in the dark-state is high as well, in fact much higher than in Dronpa, IrisFP or pcDronpa. This can be explained by the fact that normally the dark-state in Dronpa-like RSFPs has a higher energy in comparison with the bright-state. As a result, in non-photoswitchable FPs the dark state is simply not or barely populated due to a too high energy difference separating two states. Increased flexibility can relate to the lower stabilization of the H-bond network between the aminoacids, so although the dark-state is virtually possible in such protein as EosFP, it is not long-lived and the chromophore quickly relaxes to the bright-state structure. This hypothesis is supported by the fact that the final dark states of the RSFPs have lower flexibility than the transient states connecting the bright and dark states. We observed a correlation between the brightness of fluorescence and the flexibility of the bright-state, not only for Dronpa, pcDronpa and IrisFP, but also for the Dronpa mutants rsFastLime and rsKame. The FMA proves to be a sensitive tool to probe mutation sites, which would affect the most the photoswitching property of the protein. Specifically for Dronpa the sites of interest are the aminoacids situated in the β-sheets 7-10 and also the ones in the loop connected to the α-helix.

PLS-based FMA analysis of the effect of a single mutation

(Dronpa/rsFastLime/rsKame) shows a striking difference in the results. In particular in rsFastLime the difference in the collective motions between bright and dark state is barely noticeable and involves mostly the α-helix. In rsKame it is much higher and is brought about by the β-sheets 7-10. Dronpa is situated somewhere in between these two and shows high difference in the motions of Val157. Thus the lower difference in the internal collective motions of the bright and darks state of RSFP could contribute to the faster-switching property. In conclusion, flexibility analysis of RSFPs correlates with photoswitching properties, but as an indicator provides only a low resolution measure for protein design. On the other

22


Page 22 of 32

Page 23 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


hand FMA yields complementary information to predict mutation sites that have a high probability to modulate photoswitching properties. Further studies involving the study of the free energy difference and the detailed study of the mechanism involved into the process of photoswitching both thermally and photochemically are needed to understand the whole picture.

Acknowledgement Financial support from the Flemish Government through the KU Leuven concerted action scheme is gratefully acknowledged. The computational resources and services used in this work were provided by the VSC (Flemish Supercomputer Center), funded by the Hercules Foundation and the Flemish Government. B.M. is funded by a Ph.D. grant from the Agency for Innovation by Science and Technology (IWT) Flanders. P.D. thanks the Research Foundation-Flanders (FWO Vlaanderen) for a postdoctoral fellowship and the KULeuven for a research professorship.

Supporting Information Available Chromophore’s parameters and flexibility anlysis for all the proteins which were discussed above. This material is available free of charge via the Internet at http://pubs.acs.org/.

References (1) Shimomura, O.; Johnson, F. H.; Saiga, Y. Extraction, Purification and Properties of Aequorin, a Bioluminescent Protein from the Luminous Hydromedusan,Aequorea. J. Cell. Comp. Physiol. 1962, 59, 223–239. (2) Prasher, D. C.; Eckenrode, V. K.; Ward, W. W.; Prendergast, F. G.; Cormier, M. J. 23



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Primary Structure of the Aequorea Victoria Green-Gluorescent Protein. Gene 1992, 111, 229–233. (3) Shaner, N. C.; Steinbach, P. A.; Tsien, R. Y. A Guide to Choosing Fluorescent Proteins. Nat. Methods 2005, 2, 905–909. (4) Heim, R.; Cubitt, A. B.; Tsien, R. Y. Improved Green Fluorescence. Nature 1995, 373, 663–664. (5) Cormack, B. P.; Valdivia, R. H.; Falkow, S. FACS-Optimized Mutants of the Green Fluorescent Protein (GFP). Gene 1996, 173, 33–38. (6) Yang, T.-T. Improved Fluorescence and Dual Color Detection with Enhanced Blue and Green Variants of the Green Fluorescent Protein. J. Biol. Chem. 1998, 273, 8212–8216. (7) Ormö, M.; Cubitt, A. B.; Kallio, K.; Gross, L. A.; Tsien, R. Y.; Remington, S. J. Crystal Structure of the Aequorea Victoria Green Fluorescent Protein. Science 1996, 273, 1392–1395. (8) Shaner, N. C.; Campbell, R. E.; Steinbach, P. A.; Giepmans, B. N. G.; Palmer, A. E.; Tsien, R. Y. Improved Monomeric Red, Orange and Yellow Fluorescent Proteins Derived from Discosoma sp. Red Fluorescent Protein. Nat. Biotechnol. 2004, 22, 1567– 1572. (9) Lelimousin, M.; Noirclerc-Savoye, M.; Lazareno-Saez, C.; Paetzold, B.; Le Vot, S.; Chazal, R.; Macheboeuf, P.; Field, M. J.; Bourgeois, D.; Royant, A. Intrinsic Dynamics in ECFP and Cerulean Control Fluorescence Quantum Yield. Biochemistry (Mosc.) 2009, 48, 10038–10046. (10) Ando, R. Regulated Fast Nucleocytoplasmic Shuttling Observed by Reversible Protein Highlighting. Science 2004, 306, 1370–1373.

24


Page 24 of 32

Page 25 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


(11) Habuchi, S.; Ando, R.; Dedecker, P.; Verheijen, W.; Mizuno, H.; Miyawaki, A.; Hofkens, J. From The Cover: Reversible Single-Molecule Photoswitching in the GFPLike Fluorescent Protein Dronpa. Proc. Natl. Acad. Sci. 2005, 102, 9511–9516. (12) Stiel, A. C.; Trowitzsch, S.; Weber, G.; Andresen, M.; Eggeling, C.; Hell, S. W.; Jakobs, S.; Wahl, M. C. 1.8 Å Bright-State Structure of the Reversibly Switchable Fluorescent Protein Dronpa Guides the Generation of Fast Switching variants. Biochem. J. 2007, 402, 35–42. (13) Wiedenmann, J.; Ivanchenko, S.; Oswald, F.; Schmitt, F.; R´’ocker, C.; Salih, A.; Spindler, K.-D.; Nienhaus, G. U. EosFP, a Fluorescent Marker Protein with UVInducible Green-to-Red Fluorescence Conversion. Proc. Natl. Acad. Sci. 2004, 101, 15905–15910. (14) Gurskaya, N. G.; Verkhusha, V. V.; Shcheglov, A. S.; Staroverov, D. B.; Chepurnykh, T. V.; Fradkov, A. F.; Lukyanov, S.; Lukyanov, K. A. Engineering of a Monomeric Green-to-Red Photoactivatable Fluorescent Protein Induced by Blue Light. Nat. Biotechnol. 2006, 24, 461–465. (15) Adam, V.; Lelimousin, M.; Boehme, S.; Desfonds, G.; Nienhaus, K.; Field, M. J.; Wiedenmann, J.; McSweeney, S.; Nienhaus, G. U.; Bourgeois, D. Structural Characterization of IrisFP, an Optical Highlighter Undergoing Multiple Photo-Induced Transformations. Proc. Natl. Acad. Sci. 2008, 105, 18343–18348. (16) Adam, V.; Moeyaert, B.; David, C. C.; Mizuno, H.; Lelimousin, M.; Dedecker, P.; Ando, R.; Miyawaki, A.; Michiels, J.; Engelborghs, Y. Rational Design of Photoconvertible and Biphotochromic Fluorescent Proteins for Advanced Microscopy Applications. Chem. Biol. 2011, 18, 1241–1251. (17) Moeyaert, B.; Nguyen Bich, N.; De Zitter, E.; Rocha, S.; Clays, K.; Mizuno, H.; van Meervelt, L.; Hofkens, J.; Dedecker, P. Green-to-Red Photoconvertible Dronpa Mutant 25



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

for Multimodal Super-resolution Fluorescence Microscopy. ACS Nano 2014, 8, 1664– 1673. (18) Betzig, E.; Patterson, G. H.; Sougrat, R.; Lindwasser, O. W.; Olenych, S.; Bonifacino, J. S.; Davidson, M. W.; Lippincott-Schwartz, J.; Hess, H. F. Imaging Intracellular Fluorescent Proteins at Nanometer Resolution. Science 2006, 313, 1642–1645. (19) Rust, M. J.; Bates, M.; Zhuang, X. Sub-Diffraction-Limit Imaging by Stochastic Optical Reconstruction Microscopy (STORM). Nat. Methods 2006, 3, 793–796. (20) Hofmann, M.; Eggeling, C.; Jakobs, S.; Hell, S. W. Breaking the Diffraction Barrier in Fluorescence Microscopy at low Light Intensities by Using Reversibly Photoswitchable Proteins. Proc. Natl. Acad. Sci. 2005, 102, 17565–17569. (21) Dedecker, P.; Hotta, J.-i.; Flors, C.; Sliwa, M.; Uji-i, H.; Roeffaers, M. B. J.; Ando, R.; Mizuno, H.; Miyawaki, A.; Hofkens, J. Subdiffraction Imaging Through the Selective Donut-Mode Depletion of Thermally Stable Photoswitchable Fluorophores: Numerical Analysis and Application to the Fluorescent Protein Dronpa. J. Am. Chem. Soc. 2007, 129, 16132–16141. (22) Grotjohann, T.; Testa, I.; Reuss, M.; Brakemann, T.; Eggeling, C.; Hell, S. W.; Jakobs, S. rsEGFP2 Enables Fast RESOLFT Nanoscopy of Living Cells. eLife 2012, 1, e00248. (23) Moeyaert, B.; Dedecker, P. In Photoswitch. Proteins; Cambridge, S., Ed.; Springer New York, 2014; Vol. 1148; pp 261–276. (24) Duwé, S.; Moeyaert, B.; Dedecker, P. In Curr. Protoc. Chem. Biol.; Mahal, L., Romesberg, F., Shah, K., Shamu, C., Strano, M. S., Thomas, C., Eds.; John Wiley & Sons, Inc., 2015; pp 27–41.

26


Page 26 of 32

Page 27 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


(25) Vandenberg, W.; Leutenegger, M.; Lasser, T.; Hofkens, J.; Dedecker, P. DiffractionUnlimited Imaging: From Pretty Pictures to Hard Numbers. Cell Tissue Res. 2015, 360, 151–178. (26) Heilemann, M.; Dedecker, P.; Hofkens, J.; Sauer, M. Photoswitches: Key Molecules for Subdiffraction-Resolution Fluorescence Imaging and Molecular Quantification. Laser Photonics Rev. 2009, 3, 180–202. (27) Dedecker, P.; De Schryver, F. C.; Hofkens, J. Fluorescent Proteins: Shine on, You Crazy Diamond. J. Am. Chem. Soc. 2013, 135, 2387–2402. (28) Li, X.; Chung, L. W.; Mizuno, H.; Miyawaki, A.; Morokuma, K. A Theoretical Study on the Nature of On- and Off-States of Reversibly Photoswitching Fluorescent Protein Dronpa: Absorption, Emission, Protonation, and Raman. J. Phys. Chem. B 2010, 114, 1114–1126. (29) Moors, S. L. C.; Michielssens, S.; Flors, C.; Dedecker, P.; Hofkens, J.; Ceulemans, A. How Is cis-trans Isomerization Controlled in Dronpa Mutants? A Replica Exchange Molecular Dynamics Study. J. Chem. Theory Comput. 2008, 4, 1012–1020. (30) Warren, M. M.; Kaucikas, M.; Fitzpatrick, A.; Champion, P.; Timothy Sage, J.; van Thor, J. J. Ground-State Proton Transfer in the Photoswitching Reactions of the Fluorescent Protein Dronpa. Nat. Commun. 2013, 4, 1461. (31) Gayda, S.; Nienhaus, K.; Nienhaus, G. U. Mechanistic Insights into Reversible Photoactivation in Proteins of the GFP Family. Biophys. J. 2012, 103, 2521–2531. (32) Habuchi, S.; Dedecker, P.; Hotta, J.-i.; Flors, C.; Ando, R.; Mizuno, H.; Miyawaki, A.; Hofkens, J. Photo-Induced Protonation/Deprotonation in the GFP-like Fluorescent Protein Dronpa: Mechanism Responsible for the Reversible Photoswitching. Photochem. Photobiol. Sci. 2006, 5, 567–576.

27



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(33) Schäfer, L. V.; Groenhof, G.; Boggio-Pasqua, M.; Robb, M. A.; Grubmüller, H. Chromophore Protonation State Controls Photoswitching of the Fluoroprotein asFP595. PLoS Comput. Biol. 2008, 4, e1000034. (34) Kao, Y.-T.; Zhu, X.; Min, W. Protein-flexibility Mediated Coupling Between Photoswitching Kinetics and Surrounding Viscosity of a Photochromic Fluorescent Protein. Proc. Natl. Acad. Sci. 2012, 109, 3220–3225. (35) Jonasson, G.; Teuler, J.-M.; Vallverdu, G.; Mérola, F.; Ridard, J.; Lévy, B.; Demachy, I. Excited State Dynamics of the Green Fluorescent Protein on the Nanosecond Time Scale. J. Chem. Theory. Comput. 2011, 7, 1990–1997. (36) Li, X.; Chung, L. W.; Mizuno, H.; Miyawaki, A.; Morokuma, K. Primary Events of Photodynamics in Reversible Photoswitching Fluorescent Protein Dronpa. J. Phys. Chem. Lett. 2010, 1, 3328–3333. (37) Schäfer, L.; Groenhof, G.; Klingen, A.; Ullmann, G.; Boggio-Pasqua, M.; Robb, M.; Grubmüller, H. Photoswitching of the Fluorescent Protein asFP595: Mechanism, Proton Pathways, and Absorption Spectra. Angew. Chem. Int. Ed. 2007, 46, 530–536. (38) Mironov, V. A.; Khrenova, M. G.; Grigorenko, B. L.; Savitsky, A. P.; Nemukhin, A. V. Thermal Isomerization of the Chromoprotein asFP595 and Its Kindling Mutant A143G: QM/MM Molecular Dynamics Simulations. J. Phys. Chem. B 2013, 117, 13507–13514. (39) Mizuno, H.; Mal, T. K.; Walchli, M.; Kikuchi, A.; Fukano, T.; Ando, R.; Jeyakanthan, J.; Taka, J.; Shiro, Y.; Ikura, M. Light-dependent Regulation of Structural Flexibility in a Photochromic Fluorescent Protein. Proc. Natl. Acad. Sci. 2008, 105, 9227– 9232. (40) Andresen, M.; Stiel, A. C.; Trowitzsch, S.; Weber, G.; Eggeling, C.; Wahl, M. C.; Hell, S. W.; Jakobs, S. Structural Basis for Reversible Photoswitching in Dronpa. Proc. Natl. Acad. Sci. 2007, 104, 13005–13009. 28


Page 28 of 32

Page 29 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


(41) Rosenbloom, A. B.; Lee, S.-H.; To, M.; Lee, A.; Shin, J. Y.; Bustamante, C. Optimized Two-Color Super Resolution Imaging of Drp1 During Mitochondrial Fission with a Slow-Switching Dronpa Variant. Proc. Natl. Acad. Sci. 2014, 111, 13093–13098. (42) Duan, C.; Adam, V.; Byrdin, M.; Ridard, J.; Kieffer-Jaquinod, S.; Morlot, C.; Arcizet, D.; Demachy, I.; Bourgeois, D. Structural Evidence for a Two-Regime Photobleaching Mechanism in a Reversibly Switchable Fluorescent Protein. J. Am. Chem. Soc. 2013, 135, 15841–15850. (43) Hanwell, M. D.; Curtis, D. E.; Lonie, D. C.; Vandermeersch, T.; Zurek, E.; Hutchison, G. R. Avogadro: an Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J. Cheminformatics 2012, 4, 17. (44) Cornell, W. D.; Cieplak, P.; Bayly, C. I.; Gould, I. R.; Merz, K. M.; Ferguson, D. M.; Spellmeyer, D. C.; Fox, T.; Caldwell, J. W.; Kollman, P. A. A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules. J. Am. Chem. Soc. 1995, 117, 5179–5197. (45) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A. et al. Gaussian 09. 2009. (46) Dupradeau, F.-Y.; Pigache, A.; Zaffran, T.; Savineau, C.; Lelong, R.; Grivel, N.; Lelong, D.; Rosanski, W.; Cieplak, P. The R.E.D. tools: advances in RESP and ESP Charge Derivation and Force Field Library Building. Phys. Chem. Chem. Phys. 2010, 12, 7821–7839. (47) Case, D. A.; Babin, V.; Berryman, J. T.; Betz, R. M.; Cai, Q.; Cerutti, D. S.; Cheatham, T. E.; III, T. A. D.; Duke, R. E.; Gohlke, H. et al. AMBER 14. 2014. (48) Sousa da Silva, A. W.; Vranken, W. F. ACPYPE - AnteChamber PYthon Parser interfacE. BMC Res. Notes 2012, 5, 367. 29



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(49) Wilmann, P. G.; Turcic, K.; Battad, J. M.; Wilce, M. C.; Devenish, R. J.; Prescott, M.; Rossjohn, J. The 1.7 Å Crystal Structure of Dronpa: A Photoswitchable Green Fluorescent Protein. J. Mol. Biol. 2006, 364, 213 – 224. (50) Leaver-Fay, A.; Tyka, M.; Lewis, S. M.; Lange, O. F.; Thompson, J.; Jacak, R.; Kaufman, K. W.; Renfrew, P. D.; Smith, C. A.; Sheffler, W. et al. Methods Enzymol.; Elsevier, 2011; Vol. 487; pp 545–574. (51) Pronk, S.; Pall, S.; Schulz, R.; Larsson, P.; Bjelkmar, P.; Apostolov, R.; Shirts, M. R.; Smith, J. C.; Kasson, P. M.; van der Spoel, D. GROMACS 4.5: Hhigh-Throughput and Highly Parallel Open Source Molecular Simulation Toolkit. Bioinformatics 2013, 29, 845–854. (52) Lange, O. F.; van der Spoel, D.; de Groot, B. L. Scrutinizing Molecular Mechanics Force Fields on the Submicrosecond Timescale with NMR Data. Biophys. J. 2010, 99, 647–655. (53) Krivobokova, T.; Briones, R.; Hub, J. S.; Munk, A.; de Groot, B. L. Partial LeastSquares Functional Mode Analysis: Application to the Membrane Proteins AQP1, Aqy1, and CLC-ec1. Biophys. J. 2012, 103, 786–796. (54) Nam, K.-H.; Kwon, O. Y.; Sugiyama, K.; Lee, W.-H.; Kim, Y. K.; Song, H. K.; Kim, E. E.; Park, S.-Y.; Jeon, H.; Hwang, K. Y. Structural Characterization of the Photoswitchable Fluorescent Protein Dronpa-C62S. Biochem. Biophys. Res. Commun. 2007, 354, 962–967. (55) Adam, V.; Berardozzi, R.; Byrdin, M.; Bourgeois, D. Phototransformable Fluorescent Proteins: Future Challenges. Curr. Opin. Chem. Biol. 2014, 20, 92–102. (56) Bourgeois, D.; Adam, V. Reversible Photoswitching in Fluorescent Proteins: A mechanistic View. IUBMB Life 2012, 64, 482–491.

30


Page 30 of 32

Page 31 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


(57) Hub, J. S.; de Groot, B. L. Detection of Functional Modes in Protein Dynamics. PLoS Comput. Biol. 2009, 5, e1000480. (58) Yadav, D.; Lacombat, F.; Dozova, N.; Rappaport, F.; Plaza, P.; Espagne, A. Real-Time Monitoring of Chromophore Isomerization and Deprotonation during the Photoactivation of the Fluorescent Protein Dronpa. J. Phys. Chem. B 2014, 1601–1602.

31




Page 32 of 32

Molecular Dynamic Indicators of the ... - ACS Publications

Recommend Documents