Excited State Electronic Interconversion and Structural Transformation

Feb 21, 2019 - George Augustine , sriram Srinivasa raghavan , Kamini NumbiRamudu , Shanmugam Easwaramoorthi , Ganesh Shanmugam , Jaimohan ...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF TEXAS DALLAS

B: Biophysics; Physical Chemistry of Biological Systems and Biomolecules

Excited State Electronic Interconversion and Structural Transformation of Engineered Red Emitting Green Fluorescent Protein Mutant George Augustine, sriram Srinivasa raghavan, Kamini NumbiRamudu, Shanmugam Easwaramoorthi, Ganesh Shanmugam, Jaimohan Seetharani Murugaiyan, Krishnasamy Gunasekaran, Chinju Govind, Venugopal Karunakaran, and Niraikulam Ayyadurai J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.8b10516 • Publication Date (Web): 21 Feb 2019 Downloaded from http://pubs.acs.org on February 21, 2019

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 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 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.

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 28 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

1

Excited State Electronic Interconversion and Structural Transformation of

2

Engineered Red Emitting Green Fluorescent Protein Mutant

3

George Augustine ‡* Sriram Raghavan, †* Kamini NumbiRamudu, ‡ Shanmugam Easwaramoorthi, § Ganesh

4

Shanmugam,⊥ Jaimohan Seetharani Murugaiyan║ Krishnasamy Gunasekaran,† Chinju Govind∇ Venugopal

5

Karunakaran, ∇ and Niraikulam Ayyadurai, ‡

6

AUTHOR ADDRESS

7

‡Department

8

Leather Research Institute (CSIR-CLRI), Chennai, India, 600 020.

9

† Department

of Biochemistry and Biotechnology, Council of Scientific and Industrial Research-Central

of Crystallography and Biophysics, University of Madras, Chennai, India, 600 025.

10

§ Inorganic

11

⊥ Organic

12

║ Division

13



14

National Institute for Interdisciplinary Science and Technology, Thiruvananthapuram 695 019,

15

Kerala, India

16

AUTHOR INFORMATION

17

Corresponding Author

18

‡ Dr.

19

Chennai, India. [email protected]

& Physical Chemistry Laboratory, CSIR-CLRI, India, 600 020.

& Bio-organic Chemistry, CSIR-CLRI, India, 600 020. of Advanced Materials, CSIR-CLRI, Chennai, India, 600 020.

Photosciences and Photonics Section, Chemical Sciences and Technology Division, CSIR-

N. Ayyadurai, Division of Biochemistry and Biotechnology, CSIR-CLRI, Adyar,

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry 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 2 of 28

1

ABSTRACT

2

Red-emitting fluorescent protein (RFP) with large Stokes shift offer limited auto-fluorescence background

3

used in deep tissue imaging. Here, by introducing the free amino group in Aequorea victoria, the

4

electrostatic charges of p-hydroxybenzylidene imidazolinone (p-HBI) chromophore of green fluorescent

5

protein (GFP) has been altered resulting in an unusual, 85 nm red shifted fluorescence. The structural and

6

biophysical analysis suggested that the red shift is due to positional shift occupancy of Glu222 and Arg96

7

resulted along with extended conjugation and a relaxed chromophore. Femtosecond transient absorption

8

spectra exhibited that the excited state relaxation dynamics of rGFP (τ4 = 234 ps) are faster compared to the

9

avGFP (τ4 = 3.0 ns). The nanosecond time-resolved emission spectra of rGFP reveal the continuous spectral

10

shift during emission by a solvent reorientation in the chromophore. Finally, the molecular dynamics

11

simulations revealed the rearrangement of the hydrogen bond interactions in the chromophore vicinity,

12

reshaped the symmetric distribution of Van der Waal space to fine tune the GFP structure resulting with

13

highly red shifted rGFP.

14 15 16 17 18 19 20 21

ACS Paragon Plus Environment

2

Page 3 of 28 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

1

1. Introduction

2

Aequorea victoria green fluorescent proteins (avGFP) have remarkably revolutionized the area of cell

3

biology research attributed to the inherent fluorescence quantum yields.1-4 Later, the red fluorescent proteins

4

(RFP) have encountered considerable attention than GFPs, as they have enhanced fluorescence emission

5

above 600 nm, which is more advantageous in bio-imaging owing to their reduced autofluoresence, lower

6

light scattering, and higher tissue penetration.5 In general, the RFP variants share similar chromophore core

7

motif of GFPs namely p-hydroxybenzylidene-imidazolinone (p-HBI) except for the additional acylimine

8

(C=N-C=O) group at the C1 atom of the Cα-Cβ bond of tyrosine (Tyr66), which is formed through

9

oxidation and have an extended π-conjugated network.6 This extended conjugation network induces

10

additional resonance states thereby decreasing the electronic energy gap leading to red shifted fluorescence.

11

Several research groups have aimed to mimic the structure in avGFP with red-shifted excitation/emission

12

spectra either through functional group modification of the chromophore or by changing the protein

13

environment perturbing specific interactions. For instance, Yellow fluorescent protein (YFP) (Thr203)

14

evolved from avGFP, forms π-stacking interactions with an p-HBI chromophore ensuing in a notable red-

15

shifted spectrum (20 nm).7 Moreover, Mishin research group attempted in creating red chromophore in

16

avGFP through the non-oxidative mechanism of Ser65 dehydration resulting in a double bond formation

17

between Cα and Cβ atoms of Asn68 and Glu222, respectively.8 Likewise, Budisa and colleagues introduced

18

the tryptophan analogue in the cyan fluorescent protein leading to a golden variant with a large Stokes

19

shift.9 Primarily, RFP was originated from Anthozoa species10 after which several red emitting RFPs (620

20

nm) were constructed by the combination of rational design and stepwise random mutagenesis, directed

21

evolution and somatic hyper mutation.11-12 Though GFP and RFP share 99% sequence similarity, the

22

folding ability and spectroscopical properties of both proteins are different. However, it is still challenging

23

to develop a RFP with large Stokes shift values by mimicking avGFP variant. Here, we hypothesized that

24

the reorganization of cis-trans isomerization of chromophore Tyr66 of the non-intrinsic photoswitch

25

forming GFP-HS variant which may result with the large red shift. For that, we genetically introduced the

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry 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 4 of 28

1

amine chemistry in the ortho-position of Tyr66 to mimic the electrostatic charge similar to RFP. As

2

expected, the site-specific and global incorporation of 3-amino-L-tyrosine (3-NH2-L-Tyr) in the avGFP

3

resulted in a mimic of RFP (rGFP) chromophore containing protein. Herein, we report that the positional

4

reorientation of the chromophore by introducing additional amine interacting moiety (Tyr66) in GFP which

5

yielded 85 nm red-shifted green fluorescent protein (rGFP) compared with the wild counterpart. The amino

6

substitution also has effects on physiological properties such as oligomeric state, excited state energy, pH

7

response and etc. Comparative investigations of FP’s were also evaluated through small angle X-ray

8

scattering (SAXs), steady-state and time-resolved fluorescence studies. Femtosecond time-resolved

9

transient absorption spectra and molecular dynamics simulation were carried out to identify the altered

10

excited state relaxation path and isomeric effects (cis-trans) on rGFP. The effect of water hydration channel

11

and solvent occupiable spheres were studied to understand the altered space in avGFP and rGFP,

12

respectively. Geometrical contact alteration upon addition of NH2 group in the chromophore was compared

13

with GFP, which was studied using structural, biophysical and molecular dynamics simulations. Finally,

14

nanosecond time-resolved emission spectral (TRES) studies of rGFP and its native GFP reveals the lacking

15

of critical inter-molecular interaction with the chromophore.

16

2. Methods

17

2.1 Site Specific incorporation of Non-canonical amino acids

18

Site specific incorporation of 3-NH2-L-Tyr was performed in E. coli JW2581 containing pAD-GFPY66am

19

and pAC-DOPA- 6TRN (non-canonical tRNA/synthetase) as mentioned in Ayyadurai et al.13, 14 (2011) The

20

cells having both the plasmids for site specific incorporation were grown in LB broth supplemented with

21

ampicillin (50 mg/mL) and tetracycline (7.5 mg/mL) then was incubated overnight at 37 °C. The cells in

22

the overnight cultures were collected, washed thoroughly in phosphate buffer saline (PBS), and inoculated

23

into glycerol minimal media (MM) supplemented with 0.03 mM tyrosine (Tyr), ampicillin (50 mg/mL) and

24

tetracycline (7.5 mg/mL). The cultures were allowed to grow at 37°C until attaining the mid log phase (A600

ACS Paragon Plus Environment

4

Page 5 of 28 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

1

0.6 to 0.8) and 1 mM of 3-NH2-L-Tyr was added before induction. The optimum expression of site

2

specifically incorporated 3-NH2-L-Tyr was achieved upon 0.2% arabinose induction for 7 h at 37 °C.

3

Interestingly, rGFP was obtained from both site-specific and global incorporation methods. Then, the cells

4

were harvested and frozen at – 80 °C for purification.

5

2.2 Protein expression and Purification

6 7

The pET28-avGFP plasmid containing the avGFP gene was successfully cloned into pQE80L using

8

BamH1 and HindIII restriction enzyme, as described by George et al. (2017).15 The avGFP

9

expression was induced in Escherichia coli tyrosine auxotroph JW2581 by the addition of 1 mM

10

IPTG in Luria Bertani (LB) broth and incubated at 37 oC for 6 h. At the same time, 3-NH2-L-Tyr

11

labeled protein expression was also carried out in the same E. coli auxotroph cells harboring

12

pQE80L-avGFP. The cells were first grown in minimal medium consisting of 0.3 mM tyrosine as

13

a natural analogue, until the depletion of tyrosine in cells during the mid-log phase (A600 0.6 to 0.8

14

OD). Consequently, the tyrosine substrate analogue 3-NH2-L-Tyr was added and induced for the

15

target gene expression by adding 1 mM IPTG. The optimum level of rGFP (3-NH2-L-Tyr

16

incorporated) expression was obtained under these conditions. After the expression, the cells were

17

ruptured in a high-pressure homogenizer (Stansted, Essex, UK) and the chromophore-containing

18

avGFP (green) and rGFP (red) proteins were purified with the Ni-NTA affinity column. Further,

19

the His-tag purified protein was subjected to the Gel permeation chromatography (GPC) by AKTA

20

Explorer FPLC system containing Sephadex G25 HR and Superdex 200 HR column at 37 oC.

21

2.3 Crystallization, Data collection, structure solution statistics.

22

The avGFP was crystallized using vapor diffusion, hanging and sitting drop method at 293 K with

23

a protein concentration of about 10 mg/ml in a well buffer containing 10 % PEG 8000, 100 mM

24

HEPES pH 8.0 and 100 mM MgCl2. Spiky needles like crystal grew within a day. Micro seeding

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry 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 6 of 28

1

was carried out to increase the size of the crystal and hexagonal prism type crystal morphology was

2

observed. The data were collected at 100 K with double crystal single wavelength monochromatic

3

beam source with an experimental wavelength of 0.98 Å. The data was collected using pixel

4

detector Dectris Pilatus 2 M at Elettra’s beam line 5.2 R light source (Italy), 3.2 Å was integrated

5

using XDS and consequently scaled and merged using SCALA.16-18 The initial model was accessed

6

by PHENIX

7

presence of non-crystallographic or crystal twinning present in the data. Initially, all anisotropic

8

factors B were evaluated for occurrence of anomalies due to poor data strength and low-resolution

9

completion, based on bijvoet difference. No pseudo-translation was found in the model; hence, the

10

data was set to isotopic for further refinement. PHENIX phaser using maximum probability was

11

used to solve the structure, with a model template of PDB ID: 3V3D used for Molecular

12

replacement at 1.95 Å.19 The 3V3D model was fixed by removing all the water molecules. Initially,

13

rigid body refinement procedure was utilized employing real-space refinement. The Chromophore

14

was fixed and restrained with phoenix- elbow refinement. The model was validated using

15

MolProbity analysis for any statistical deviance.20

19

The xtriage routine was conducted to evaluate the obtained model, to know the

16 17

3. Results and discussion

18

3.1 Development of congener rGFP and its spectral characterization

19

To mimic RFP, we introduced an electron donating amino group in the ortho-position of the Tyr residue of

20

anionic avGFP, by employing misaminoacylation21 and orthogonal translational machinery.22 avGFP and

21

congener rGFP were purified with Ni-NTA affinity column, to obtain a final concentration around ~28 and

22

20 mgL-1, respectively. Furthermore, the quality of NH2-Tyr incorporation was ensured by SDS-PAGE

23

(Supporting Information, Figure S1), and Matrix Assisted Laser Desorption-Time of Flight (MALDI-TOF)

24

analysis. A good correlation between calculated mass (28266.7 Da) and determined mass (28266.1 Da) was

ACS Paragon Plus Environment

6

Page 7 of 28 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

1

found, which are shown in Supporting Information, Figure S2. The MALDI-TOF analysis confirmed that

2

the purified GFP variants with Tyr analogues possessed an incorporation efficiency of >90%, which

3

resulted in “Red” GFP (rGFP). The “rGFP” obtained from the site and residue-specific incorporation

4

method exhibited absorption and emission maximum at

5

respectively (Supporting Information, Figure S9A&B). The emission maxima of rGFP were observed at

6

610 nm which is red-shifted to the tune of 85 nm (0.34 eV) when compared to that of avGFP (525 nm, 2.40

7

eV) (Supporting Information, Figure S9A&B). Previously, several attempts have been made to create the

8

most-red shifted proteins from anionic GFP chromophore. To our knowledge, the developed rGFP is the

9

most red shifted variant of avGFP characterized to date apart from golden FP 9 (λabs = 466 nm, λem = 574

10

nm) and EYFP23 (λabs = 514 nm, λem = 527 nm). The emission maxima of rGFP was about 27 nm (0.082

11

eV) red shifted from commercially available red fluorescent protein (DsRed, λem = 583 nm). 24 It is similar

12

to DsRed derived mKeima (λem = 620 nm), mKate variant (LSSmKate-2 λem = 605 nm), mcherry (λem = 611

13

nm) and mRuby (λem = 605 nm) variants.25-27 The extended red shifts of DsRed derivatives are not only

14

attributed to the mutations introduced in the chromophore, but also to the rearrangements in the interactions

15

(hydrogen bonding, electrostatic and hydrophobic) of the chromophore with surrounding key amino acid

16

residues. For example, the large Stokes shift has been explained by various mechanisms. mPlum had been

17

created with extended Stokes shift due to solvent reorganization28 and excited state relaxation of the

18

chromophore, enabled by 2.56 Å flexible hydrogen bond between the Ile 65 (Chromophore residue) and

19

Glu16.29 On the other hand, Tag RFP675’s large Stokes shift was due to the re-establishment of dipole-

20

dipole interaction between the chromophore and protein matrix.30 In case of the mKeima,25 LSS mKate,27

21

mBeRFP

22

(ESPT). Based on the above genetic incorporation, spectroscopic and theoretical studies, we assumed that

23

the enhanced red shift caused by mutations (NH2-Tyr) could result in DsRed mimic of avGFP variant.

24 25

3.2 X-ray crystallographic investigations and rGFP modelling.

31

and mCherry,

32

525 nm (2.37eV) and at

610 nm (2.03eV),

large Stokes shift were provided by engineered excited state proton transfer

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry 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 8 of 28

1

The structural basis for the extended Stokes shift, the 3D structure of avGFP were obtained using X-ray

2

crystallography (Supporting Information, Table S1), which was solved at physiological pH 7.4 with a

3

resolution of 3.2 Å (Protein Data Bank, PDB: 5WWK). From crystallographic structural investigation, we

4

determined the hydrogen bond distance between chromophore and surrounding key residues, which were

5

as follows: OG Ser205–OE1 Glu222 was 4.79 Å, OG Ser205–HOH was 2.31 Å, O Asn146–HOH was 2.77

6

Å, OG Thr 203-HOH was 3 Å and O Asn 146–ND1 His148 was 4.17 Å. The hydrogen bonds generally

7

have energies between 0.2–40 kcal mol-1 with a bond length ranging between 2.2–2.5 Å, which are

8

considered as strong, 2.5-3.2 Å as moderate and 3.2-4.0 Å as weak bond.33 Dynamic bond positioning

9

average may vary, but there was a proton contact between His148-Asn146 and Ser205-Glu222, which was

10

contributing to the weak electrostatic interaction. Electron density map with 1.80-sigma cut-off of 2Fo-Fc

11

form represented continuous electron density between O2 CR2 – NH1, NH2-Arg 96, which had strong

12

hydrogen bond of 2.46 Å and 2.80 Å (Supporting Information Figure S6), respectively. The solvent voids

13

of overall avGFP and its active site region are shown in Figure 1A (1-2). The rGFP was also diffracted to

14

a resolution of 5 Å and was resilient to higher quality crystallographic trials. Hence, SAXs profiling and

15

protein modeling were carried out for rGFP in comparison with avGFP. Through molecular dynamics

16

approach, we computationally constructed large red shifted rGFP by adding NH2 group in the ortho position

17

of Tyr residue in the GFP chromophore. For full computational details, see the supporting information. As

18

a result, the 3D structure of the rGFP was modeled with chromophore consisting of o-amino p-

19

hydroxybenzylidene-imidazolinone (oApHBI). To discern the structural basis for red shifted absorption

20

and emission in rGFP, avGFP was superimposed with rGFP. Even though the backbone key residues Val60,

21

Thr61, Gln91, and Gln219 of avGFP are superimposable with rGFP backbone, residues Arg96, ASN146

22

and Glu222 were significantly deviated from their position, to accommodate the oApHBI chromophore

23

within the compact β barrel (Figure 1B3). Also, the electrostatic landscape of rGFP was slightly altered in

24

comparison with GFP, which is shown in Figure 1C. These residues deviation caused an alteration in the H

25

bond rearrangement (Supporting Information, Figure. S20) and electrostatic interaction between the

26

chromophore and key residues in rGFP (Figure 1C).

ACS Paragon Plus Environment

8

Page 9 of 28 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

1

3.3 SAX investigation

2

The rGFP was resilient to crystallization and is oligomeric in nature, resembling most tetrameric red shifted

3

FPs (tetrameric E2 orange and E2-Crimson).34, 35 Before the SAXs modeling, the proteins were subjected

4

to dynamic light scattering (DLS) to resolve the protein state (Supporting Information, Figure. S7). The

5

rGFP showed high polydispersity with a broad peak (750 nm) indicating the predominance of oligomeric

6

species (Supporting Information, Figure. S7B). In contrast, smaller sizes for avGFP were obtained

7

reflecting the formation of a monomeric form of the protein in the solution (Supporting Information, Figure

8

S7A). Hence, SAXs profiles were acquired to determine the dynamic behavior of the rGFP and avGFP, to

9

understand their corresponding oligomeric state, monodispersity, dynamic motions (determined by the

10

average radius of gyration, Rg) 36- 37 (Supporting Information, Figure S4 & Table S2). The intensity profile

11

of avGFP and rGFP were plotted as double logarithmic scales (I (Q) vs. Log Q), to compare the monomeric

12

to oligomeric states of the protein (Supporting Information, Figure S4.3). Guinier’s linear fit was plotted to

13

understand the monodispersity of the proteins (Supporting Information, Figure. S4.1). Radius of gyration

14

(Rg) and Radius of cross section (Rc) values of rGFP were determined as 34.5 nm and 4.9 nm, respectively

15

(Supporting Information, Table S2). Persistent length (L) or compactness of the scattering biomolecule was

16

found to be 30 nm. Kratky plot (Fig. 4-2) [I (Q)*Q2 vs. Q] of SAXs data showed bell-shaped curves that

17

exhibited prolate nature of the protein with defined particle size. Similar to Kratky plot, pair wise distance

18

distribution P(r) analysis of the protein was also used to calculate Dmax value of the full-length protein using

19

indirect Fourier transformation analysis (Supporting Information, Figure S4.4). To obtain a better insight

20

of the structural nature of the protein, the monomeric GFP having p2 symmetry was fitted to dummy atom

21

model of rGFP having p32 symmetry, which demonstrated oligomeric propensity along with χ2 scalar fitting

22

values with agreeable deviations (Supporting Information, Table S2). The superimposed dummy scattered

23

model of rGFP and avGFP are shown in Supporting information, Figure (4.5).

24

3.4 Molecular dynamic simulation and QM/MM

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry 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 10 of 28

1

To comprehend, the excited state absorption at molecular level, the chromophore of GFP and rGFP in the

2

gas phase has extensively been studied using molecular dynamics simulations (QM/MM). During the

3

current study, it was found that the excited state spectra displayed a major peak at around 418 nm (2.96 eV)

4

for avGFP and red shifted band at 498 nm (2.48 eV) for rGFP. Despite variation in the computed values,

5

the absorption energy mentioned above for rGFP is red shifted by 0.48 eV (80 nm) relative to the avGFP.

6

We noticed acceptable correlations between computed excitation energies and corresponding experimental

7

absorption band maxima of avGFP (2.37 eV) and rGFP (2.46 eV), within the established standard

8

deviations of QM/MM protocol. Indeed, this pronounced 0.48 eV energy difference prompted us to

9

investigate the nature of the rGFP’s chromophore conformation (cis-trans) at the excited state through

10

molecular simulation. It has been reported that the spectral perturbation arising from torsional and dihedral

11

angle variations around the chromophore is proposed to play a vital role in spectral deviations.38,

12

Substitution of Asp148 in GFP variant (S65T/H148D) had led to an increased trans-isomer conformational

13

probability known as asynchronous hula twist, attributed to the altered excited state.40 In our case, we found

14

that cis and trans conformations of the chromophore are prevalent in wild GFP and rGFP, respectively.

15

Conformational densities of dihedral angle (ø, τ) variations obtained by MD for cis-trans isomerism

16

concerning spatial occupancy are shown in Figure 2A and Figure S11-13. From these results, we observed

17

that the ethylene linker between the imidazolone and benzyl ring of Tyr in the core chromophore did not

18

show any major distortions and are within limited standard deviations. Main spectral perturbation (extended

19

red shift) can only be attributed to the incorporation of NH2-Tyr along with acquired trans-isomers in rGFP.

20

As per the previous report, excited state dynamics of Wt GFP chromophores (p-HBI) are governed by

21

Excited State Proton Transfer (ESPT), in which proton shuttles via connected networks of His148, Asp146,

22

Thr203, Ser205, and Glu222. But for the targeted proton transfer (PT) to occur, the H-bond distance

23

between the p-hydroxyphenyl group of the chromophore (donor) and the carbonyl group of Glu (acceptor)

24

should be very short .41 In particular, the central core (p-HBI) consisting of ionizable phenoxy group is

25

directly attributed to the formation of a hydrogen bond with His148 (gatekeeper), which is responsible for

39

ACS Paragon Plus Environment

10

Page 11 of 28 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

1

proton transfer only in Wt GFP but not in GFP variants (sf GFP and EGFP).42,43 X-ray crystallography

2

studies also confirmed this hypothesis. Based on these observations, we also anticipated the related

3

mechanism of proton transfer exhibiting a red shift in rGFP in comparison with avGFP. In contrast, through

4

crystallographic study, we found that the hydrogen bond distances between the chromophore and

5

surrounding residues are long, which can be classified in the range of moderate to weak categories

6

concerning bond strengths, as described above. Interestingly, according to our GFP crystal structure, we

7

deduced that His148 in avGFP and rGFP do not display direct Hydrogen bond interaction with the phenoxy

8

group of p-HBI. In detail, upon excitation, the hydroxyl group in GFP’s aromatic chromophore gets ionized

9

and tends to lose the proton leading to the formation of the anionic state. At this stage, the anionic phenolate

10

group is stabilized (balanced) by His148 residue through the water-mediated hydrogen bond network.

11

When hydroxyl group acquires partial negativity, nitrogen in imidazolone attains partial positive charge

12

and is stabilized by Glu222 network through ESPT. The system works like an electrostatic assembly similar

13

to the one given in Figure 3D, where the entire electrostatic charge network is counter balanced.

14

Introducing additional amine group in anionic avGFP tends to form self zwitterionic assembly in its

15

chromophore shown in Figure 3C. In detail, partial charge species association on NH2 group is more

16

possible than the imine in imidazolone, which is due to the proximity of amine group to phenolate moiety

17

of chromophore. Here, the distance between the oxygen of hydroxyl and nitrogen in the amine is around

18

2.6 Å, which lies in the strong hydrogen bond region. The radial hydrogen bond distance is of 1.4 Å.

19

Incorporation of amine in ortho position not only led to the spectral change but the chromophore

20

benzylidene moiety considerably shifted its locus towards polar Thr203 and Ser205. Additional pH

21

dependent spectral studies were also carried out to confirm the charge species of rGFP. The UV-visible

22

absorption spectra of rGFP between the pH 5 and pH 7 distinguish the spectral features corresponding to

23

the deprotonated and zwitterionic form (Supporting information, Figure S9E). At pH 5, the absorption

24

observed with the maximum at ~400 nm is can be attributed to the zwitterionic form and the emission

25

spectra was measured by exciting the sample at 360 nm ensures the significant protonated state population

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry 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 12 of 28

1

(Supporting information, Figure S9F). The observation of emission spectra of rGFP with similar spectral

2

features at different pH and irrespective of the excitation wavelength suggests that no excited state proton

3

transfer occurring in rGFP (Supporting information, Figure S9F-G).

4

The hydrogen bond networks stabilized the symmetric charge distribution around the chromophore residues

5

of avGFP (Arg96 & Glu222). While the incorporation of amine decreased the hydrogen bond interactions

6

in rGFP as compared to wild avGFP (Supporting Information, Figure S20A). The amine incorporation can

7

be interpreted based on it being self-zwitterionic and its stabilizing nature and thus does not require

8

additional hydrogen bonds to stabilize charge equilibration. Also, we were able to identify symmetrical Van

9

der Waals distribution around the chromophore for avGFP, and its positional deviation for the amino

10

derived variants for rGFP (Figure 3A). The Van der Waals sphere in rGFP shifts from a mean average

11

position involving a process of charge equilibration i.e., protonated amine being neutralized by a carboxylic

12

group of Glu222. Though Glu222 is known to be a substantial residue, Arg96 can also be characterized as

13

an important residue in maintaining the charge equilibration. Along with residual networks, water hydration

14

pathway also plays a vital role in relaxation dynamics. Solvent accessible sphere near chromophore is given

15

in Figure 3B. The solvent channel is highly adapted such that water channelizes through a particular path

16

to solvent accessible voids near the chromophore. It is known that Glu222 and Arg96 not only play a vital

17

role in chromophore formation but also in maintaining charge equilibrium. Based on the residual occupancy

18

distance plot it was identified that arginine had conserved contact but Glu222 changed with respect to time

19

(Figure 2B & Supporting Information, Figure S14-19).

20

3.5 rGFP maturation analysis

21

Initial spectroscopic characterization and detailed understanding of protein structure by X-ray

22

crystallographic studies and QM/MM approach emphasized investigation of red chromophore

23

transformation and maturation of rGFP, which was compared with rapidly maturing mono red fluorescent

24

protein (mRFP) variant by fluorescent measurements (Supporting Information, Figure S8). Both mRFP and

ACS Paragon Plus Environment

12

Page 13 of 28 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

1

rGFP yielded a total maturation time and a half-life (t1/2) of about ~6 and ~3 h, respectively. These results

2

demonstrated that rGFP matures faster than the Wt DsRed (t1/2~12 h) and have a maturation time very close

3

to mRuby27 (2.8 h) and mScarlet44 (2.9 h) variants of DsRed, which render the protein suitable for bio-

4

imaging applications. Hence, it can be deduced that the substitutions of NH2-Tyr in avGFP do not distort

5

the protein secondary structure (Supporting Information, Figure S3) and instead, it accelerates the protein

6

folding and maturation by intermolecular rearrangements of key residues concerning the interaction with

7

the chromophore.

8

3.6 Effect of Cryogenic temperature on rGFP by Steady state fluorescent analysis

9

The steady-state fluorescence spectra of rGFP and avGFP were measured at 100 to 400 K (Figure 4A).

10

While decreasing the temperature (310-245 K) the fluorescent intensity of avGFP was gradually diminished

11

(Figure 4A). In contrast, the fluorescence intensity of rGFP were considerably increased at their respective

12

emissions (λem = 610, 525 nm). Particularly, rGFP displayed three distinct minor peaks with the maximum

13

at 491, 521, and 570 nm upon exciting the sample at 440 nm, 100 K without suppressing the major peak at

14

610 nm. Similar phenomenon had been observed for BFP in which the dissociation of the Hydrogen bond

15

between the chromophore and His148 and strong water-mediated linkage have been attributed to the

16

enhanced fluorescent intensity at a lower temperature.45 Supporting this, rGFP also displayed strong

17

dissociation of a hydrogen bond with chromophore as compared to GFP and had strong water-mediated

18

linkage near chromophore (Figure 3B). Consequently, it enhanced the rigidity of chromophore leading to

19

increased fluorescent intensity at freezing temperature.

20

3.7 ESPT and Fluorescence Life time analysis

21

These charge equilibration and proton transfer mechanism of rGFP chromophore was demonstrated using

22

molecular simulation studies. To support this, we measured absorption and steady-state fluorescence

23

spectra of both the proteins in D2O and H2O (Figure 4B-C & Supporting Information, Figure S9). As a

24

result, we observed that the deuteration (D2O-τ1-0.97 ns, τ2-2.89 ns) of rGFP did not extend the lifetime,

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry 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 28

1

which was more or less similar to rGFP in H2O (H2O-τ1-0.955 ns, τ2-2.87 ns) (Supporting Information,

2

Table S3). To our knowledge, among the variants of DsRed, rGFP protein significantly extended the

3

average life time (3.86 ns) of red shifted protein at 610 nm as compared to DsRed variants, like mKate2

4

(2.5 ns), mRuby2 (2.5 ns), TagRFP-T (2.3 ns) and mCherry (1.5 ns).44 Similarly, mRFP consistently showed

5

biexponential decay behavior with the average fluorescence lifetime of 4.22 ns. This indicated that these

6

proteins could have acquired nearly identical but two different molecular conformations of the chromophore

7

present in the excited state, as observed in other biexponential decay fitted RFPs (eqFP611 and HcRed). To

8

explore the spectral states and electron transfer of the proteins in the presence of an electron acceptor, we

9

subjected the proteins to the fluorescent spectrum analysis in Benzoquinone (BQ) (Figure 4E-F). Increasing

10

the concentration of BQ (1.25 mM) quenched the emission of both the proteins instead of promoting the

11

electron transfer.

12

3.8 Femtosecond time-resolved transient absorption spectra

13

The femtosecond time-resolved transient absorption spectra of avGFP and rGFP in 0.1 M phosphate buffer

14

at pH 7.4 were measured by exciting the sample at 500 nm. In Figure 5A1, the transient absorption spectra

15

of avGFP consists of a positive excited-state absorption (ESA) band below 450 nm and negative band from

16

450 to 600 nm due to the overlapping ground-state bleaching (GSB) and stimulated emission (SE)

17

contributions. The transient absorption spectra decayed with the increase of spectral delay time and did not

18

attain the equilibrium even after 1.5 ns. Similarly, transient absorption spectra of rGFP were measured,

19

which have been shown in Figure 5A2. The ESA maximum was centered at around 414 nm, and a negative

20

band of GSB was observed from 450 to 600 nm. Even though, the spectral features were the same as that

21

of avGFP, in rGFP the excited state relaxation dynamics were faster and came back to equilibrium at around

22

400 ps. The kinetic decays of avGFP and rGFP at probing wavelength of 420 and 520 nm were obtained

23

by exciting at 500 nm, which are shown in Figure 5B1. The analysis of the femtosecond transient absorption

24

spectra consisting of a three-dimensional data set (wavelength, time and change in absorbance) was

25

performed with the global analysis program GLOTARAN46. Four exponential components were optimally

ACS Paragon Plus Environment

14

Page 15 of 28 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

1

fitted to the transient absorption spectra to describe the excited state dynamics completely. The transient

2

absorption spectra of avGFP were fitted with four exponentials globally and the resulting time constants

3

were τ1 = 384 fs, τ2 = 1.66 ps, τ3 = 13.26 ps, and τ4 = 3.0 ns (fixed). The time constant of 3.0 ns was fixed

4

explicitly to represent one of the components of fluorescence lifetime of avGFP obtained from TCSPC.

5

Similarly, the time constants for rGFP were τ1 = 296 fs, τ2 = 4.25 ps, τ3 = 16.06 ps, and τ4 = 234 ps. It was

6

observed that the excited state relaxation dynamics of rGFP was faster when compared to avGFP. The

7

faster dynamics in rGFP could be attributed to the changes in the conformational structure and

8

alteration in the hydrogen bonding network of the chromophore resulting from the incorporation

9

of the NH2 group compared to the avGFP and needs further investigation.

10 11

3.9 Time-resolved emission spectra

12

Time-resolved emission spectra with nanosecond time resolution were measured to investigate the excited-

13

state chemistry of avGFP and rGFP. For rGFP, the early time spectra are dominated by a rapidly decaying

14

red-shifted emission. The emission was broad, which has also been observed in the corresponding studies

15

of Discosoma sp. variant mPlum or mRFP emission (Figure 6A). For the wild-type (avGFP) there is no

16

evidence of red shift and a single band is detected with a spectrum very similar to that of the Aequorea

17

victoria emission which decays as a function of time after excitation (Supporting Information, Figure S10

18

A). Detailed analysis showed that the spectrum narrowed as a function of time after excitation in rGFP. The

19

narrowing is most apparent as subsidence of the fluorescence on the red edge at the time scale of 1.65

20

nanoseconds. However, in rGFP, no isoemissive point was observed. The TRES and time-resolved area

21

normalized emission spectra (TRANES) of rGFP, and avGFP spectrum are shown in (Figure 6B &

22

Supporting Information, Figure S10). The spectra observed for rGFP (425 and 610 nm) at 2.20 ns closely

23

matches the steady-state emission spectrum of which peaks at 610 nm. A significant observation is that the

24

TRANES intersect was not observed in GFP and rGFP variants. As reported earlier, an isoemissive point

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry 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 16 of 28

1

can only be detected in a fluorescent variant with discrete emitting states, and it disappears in a system that

2

undergoes a continuous spectral shift during emission.47,48 For rGFP, the observation of non-isoemissive

3

point indicated that the system would undergo a continuous spectral shift from native GFP to rGFP (610

4

nm) emission. Genetic incorporation of ortho NH2 induced a red shift in the fluorescence due to

5

chromophore structural modification. This is because of the interaction between the dipole of the

6

chromophore and that of the surrounding solvent molecules. The data suggest that NH2-Tyr could provoke

7

the positional shift occupancy of Glu222 and Arg96 resulting in extended conjugation and relaxed

8

chromophore, where higher solvent polarity in GFP resulted in the greater red shift.

9

4. Conclusion

10

In conclusion, this study imparts new insights into the dynamic behavior of rGFP with extended Stokes

11

shift. Earlier reports proved that the unusual red shift of RFP (TagRFP 675) attributed to the multiple

12

hydrogen bonds with the chromophore. Interestingly, in contrast to the multiple hydrogen bonding, we

13

observed a reduced number of hydrogen bonds between the chromophore and surrounding key residues.

14

Upon incorporation of the amino group, the chromophore (o-NH2p-HBI) has transformed from cis to trans

15

isomer. We observed the enhanced fluorescence lifetime and rapid, excited state relaxation in rGFP by the

16

incorporation of NH2-Tyr, which might be effective modalities for extended Stokes shift. Besides, the

17

rearrangement of the hydrogen bonding network displayed an alteration in the symmetrical Van der Waals

18

distribution around the oApHBI of rGFP as compared to the avGFP. Above all, the enhanced fluorescence

19

intensity at cryogenic temperature unraveled new avenues for the application of rGFP in newly emerging

20

technique, specifically cryogenic single molecule super-resolution imaging. These results elucidated the

21

extended Stokes shift in detail, along with designed computational and experimental methods for

22

developing FPs with new properties accounting for desired applications.

23

ASSOCIATED CONTENT

24

Supporting Information.

ACS Paragon Plus Environment

16

Page 17 of 28 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

1

General details on the material and methods, full characterization of the proteins, optimization of the

2

conditions of avGFP and rGFP biosynthesis, Crystal structure details, Biophysical characterization,

3

Bioinformatics analysis.

4

Notes

5

There are no conflicts to declare

6

ACKNOWLEDGMENT

7

Department of Science and Technology initially supported this research under Fast Track Grant

8

(SB/YS/LS-217/2013). The authors gratefully acknowledge “Science and Technology Revolution in

9

Leather with a Green Touch” (STRAIT-1.2.3) and Zero Emission Research Initiative for Solid Waste

10

(ZERIS), 12th five-year plan project under Council of Scientific and Industrial Research (CSIR). The first

11

author George A thankfully acknowledges CSIR for the award of Senior Research Fellowship. We thank

12

the Director, CSIR-CLRI for his support during the project and allowing the work to be conducted at CSIR-

13

CLRI.

14

REFERENCES

15 16

(1) Misteli, T.; Spector, D. L. Applications of the green fluorescent protein in cell biology and biotechnology. Nat. Biotechnol. 1997, 15, 961–963.

17 18

(2) Zimmer, M. Green Fluorescent Protein (GFP):  Applications, Structure, and Related Photophysical Behavior. Chem. Rev. 2002, 102, 759–781.

19 20

(3) Lei, C.; Huang, Y.; Nie, Z.; Hu, J.; Li, L.; Lu, G.; Han, Y.; Yao, S. A Supercharged Fluorescent Protein as a Versatile Probe for Homogeneous DNA Detection and Methylation Analysis.

21 22 23 24

Angew. Chem. Int. Ed. 2014, 53, 8358–8362. (4) Acharya, A.; Bogdanov, A. M.; Grigorenko, B. L.; Bravaya, K. B.; Nemukhin, A. V.; Lukyanov, K. A.; Krylov, A. I. Photoinduced Chemistry in Fluorescent Proteins: Curse or Blessing? Chem. Rev. 2017, 117, 758-795.

25 26

(5) Shcherbakova, D. M.; Subach, O. M.; Verkhusha, V. V. Red fluorescent proteins: advanced imaging applications and future design. Angew. Chem. Int. Ed. 2012, 51, 10724–10738.

27 28

(6) Miyawaki, A.; Shcherbakova, D. M.; Verkhusha, V. V. Red fluorescent proteins: chromophore formation and cellular applications. Curr. Opin. Struct. Biol. 2012, 22, 679-688.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry 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

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

Page 18 of 28

(7) Tsien, R. Y. The green fluorescent protein. Annu. Rev. Biochem. 1998, 67, 509 – 544. (8) Mishin, A. S.; Subach, F.V.; Yampolsky, I.V.; King, W.; Lukyanov, K. A.; Verkhusha, V. V. The first mutant of the Aequorea victoria green Fluorescent protein that forms a red chromophore. Biochemistry. 2008, 47, 4666-4673. (9) Bae, J. H.; Rubini, M.; Jung, G.; Wiegand, G.; Seifert, M. H.; Azim, M. K.; Kim, J. S.; Zumbusch, A.; Holak, T. A.; Moroder, L. et al., Expansion of the genetic code enables design of a novel "gold" class of green fluorescent proteins. J. Mol. Biol. 2003, 328, 1071-1081. (10) Matz, M. V.; Fradkov, A. F.; Labas, Y. A.; Savitsky, A. P.; Zaraisky, A. G.; Markelov, M. L.; Lukyanov, S. A. Fluorescent proteins from nonbioluminescent Anthozoa species. Nat Biotechnol. 1999, 17, 969-973. (11) Campbell, R. E.; Tour, O.; Palmer, A. E.; Steinbach, P. A.; Baird, G. S.; Zacharias, D. A.; Tsien, R. Y. A monomeric red fluorescent protein, Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 7877– 7882. (12) Wang, L.; Jackson, W. C.; Steinbach, P. A.; Tsien, R. Y. Evolution of new nonantibody proteins via iterative somatic hypermutation. Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 16745– 16749.

17 18

(13) Seyedsayamdost, M. R.; Chan, C. T.; Mugnaini, V.; Stubbe, J.; Bennati, M. PELDOR spectroscopy with DOPA-β2 and NH2Y-α2s: distance measurements between residues involved

19

in the radical propagation pathway of E. coli ribonucleotide reductase. J. Am. Chem. Soc. 2007,

20

129, 15748−15749.

21 22 23 24 25 26 27 28 29

(14) Ayyadurai, N.; Deepankumar, K.; Saravanan Prabhu, N.; Lee, S.; Yun, H. A facile and efficient method for the incorporation of multiple unnatural amino acids into a single protein. Chem. Comm, 2011, 47, 3430-3432. (15) George, A.; Krishna Priya, G.; Ilamaran, M.; Kamini, N.R.; Ganesh, S.; Easwaramoorthi, S.; Ayyadurai, N. Accelerated Strain‐Promoted and Oxidation‐Controlled Cyclooctyne‐Quinone Cycloaddition for Cell Labeling. Chemistry Select, 2017, 24, 7117-7122. (16) Rebuffi, L.; Sánchez del Río, M.; Busetto, E.; Scardi, P. Understanding the instrumental profile of synchrotron radiation X-ray powder diffraction beamlines. J. Synchrotron Radiat 2017, 24, 622-635.

30 31

(17) Computational Collaborative project Number 4, The CC4 Suite: The program for Protein crystallography. Acta Crystallogr. Sect. D 1994, 50, 760-763.

32 33

(18) Adams, P. D.; Afonine, P. V.; Bunkoczi, G.; Chen, V. B.; Davis, I.W.; Echols, N.; Headd, J.J.; Hung, L.W.; Kapral, G. J.; Grosse-Kunstleve, R. W.; McCoy, A. J.; Moriarty, N. W.;

34

Oeffner, R.; Read, R. J.; Richardson, D. C.; Richardson, J. S.; Terwilliger, T. C.; Zwart, P.H.

35

PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta

36

Crystallogr. Sect. D Biol. Crystallogr. 2010, 66, 213-221.

ACS Paragon Plus Environment

18

Page 19 of 28 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

1 2

(19) De Meulenaere, E.; Nguyen Bich, N.; De Wergifosse, M.; Van Hecke, K.; Van Meervelt, L.; Vanderleyden, J.; Champagne, B.; Clays, K. Improving the Second-Order Nonlinear Optical

3

Response of Fluorescent Proteins: The Symmetry Argument. J. Am. Chem. Soc. 2013, 135,

4

4061-4069.

5 6

(20) Chen, V. B.; Arendall, W. B.; Headd, J. J.; Keedy, D. A.; Immormino, M.; Kapral, G. J.; Murray, L. W.; Richardson, J. S.; Richardson, D. C. MolProbity: all-atom structure validation

7

for macromolecular crystallography. Acta Crystallogr. Sect. D Biol. Crystallogr. 2010, 66, 12-

8

21.

9 10 11

(21) Moor N.; Klipcan, L.; Safro, M. G. Bacterial and eukaryotic phenylalanyl-tRNA synthetases catalyze misaminoacylation of tRNA (Phe) with 3, 4-dihydroxy-L-phenylalanine.Chem. Biol. 2011, 18, 1221-1229.

12 13

(22) Kowal, A. K.; Kohrer, C.; RajBhandary, U. L. Twenty-first aminoacyl-tRNA synthetasesuppressor tRNA pairs for possible use in site-specific incorporation of amino acid analogues

14

into proteins in eukaryotes and in eubacteria. Proc. Natl. Acad. Sci. USA. 2001, 98, 2268-2273.

15 16

(23) Nagai, R.; Ibata, K.; Park, E. S.; Kubota, M.; Mikoshiba, K.; Miyawaki, A. A variant of yellow fluorescent protein with fast and efficient maturation for cell-biological applications. Nature.

17 18 19 20 21 22 23 24 25 26 27 28 29

2002, 20, 1585-1588. (24) Verkhusha, V. V.; Chudakov, D. M.; Gurskaya, N. G.; Lukyanov, S.; Lukyanov, K. A. Common pathway for the red chromophore formation in fluorescent proteins and chromoproteins. Chem. Biol. 2004, 11, 845-854. (25) Henderson, J. N.; Osborn, M. F.; Koon, N.; Gepshtein, R.; Huppert, D.; Remington, S. J. Excited state proton transfer in the red fluorescent protein mKeima. J. Am. Chem. Soc. 2009, 131, 13212-13213. (26) 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. (27) Kredel, S.; Oswald, F. Nienhaus, K.; Deuschle, K.; Röcker, C.; Wolff, M.; Heilker, R; Nienhaus, G. U.; Wiedenmann, mRuby, a Bright Monomeric Red Fluorescent Protein for Labeling of Subcellular Structures J. PLoS One. 2009, 4, 1–7.

30 31

(28) Abbyad, P.; Childs, W.; Shi, X.; Boxer, S. G. Dynamic Stokes shift in green fluorescent protein variants. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 20189-20194.

32 33

(29) Faraji, S.; Krylov, A. I. On the Nature of an Extended Stokes Shift in the mPlum Fluorescent Protein. J. Phys. Chem. B. 2015, 119, 13052-13062

34 35

(30) Piatkevich, K. D.; Malashkevich, V. N.; Almo, S. C.; Verkhusha, V. V. Engineering ESPT Pathways Based on Structural Analysis of LSSmKate Red Fluorescent Proteins with Large

36

Stokes Shift. J. Am. Chem. Soc. 2010, 132, 10762-10770.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry 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 20 of 28

1 2

(31) Yang, J.; Wang, L.; Yang, F.; Luo, H.; Xu, L.; Lu, J.; Zeng, S.; Zhang, Z. mBeRFP, an improved large stokes shift red fluorescent protein. PLoS One. 2013, 8, 6-11.

3 4

(32) Shen, Y.; Chen, Y.; Wu, J.; Shaner, N. C.; Campbell, R. E. Engineering of mCherry variants with long Stokes shift, red-shifted fluorescence, and low cytotoxicity. PLoS One. 2017, 12, 1-

5

14.

6 7

(33) Pauling, L.; Corey, R. B.; Branson, H. R. Two Hydrogen-Bonded Helical Configurations of the Polypeptide Chain. Proc. Natl. Acad. Sci. U.S.A. 1951, 37, 205-211.

8 9

(34) Strack, R. L.; Hein, B.; Bhattacharyya, D.; Hell, S. W.; Keenan, R. J.; Glick, B. S. Biochemistry. A rapidly maturing far-red derivative of DsRed-Express2 for whole-cell labeling.

10

2009, 48, 8279-8281.

11 12

(35) Strack, R. L.; Bhattacharyya, D.; Glick B. S.; Keenan, R. J. Noncytotoxic orange and red/green derivatives of DsRed-Express2 for whole-cell labelling. BMC Biotechnol. 2009, 32, 1-10.

13 14

(36) Petoukhov, M. V.; Franke, D.; Shkumatov, A. V.; Tria, G.; Kikhney, A. G.; Gajda, M.; Gorba, C.; Mertens, H. D. T.; Konarev, P. V., Svergun, D. I. New developments in the ATSAS program

15 16 17 18

package for small-angle scattering data analysis. J. Appl. Crystallogr. 2012, 45, 342-350. (37) Singh, J.; Khan, M. I.; Singh Yadav, S. P.; Srivastava, A.; Sinha, K. K.; Ashish.; Das, P.; Kundu, B. L-Asparaginase of Leishmania donovani: Metabolic target and its role in Amphotericin B resistance. Int. J. Parasitol. Drugs Drug Resist. 2017, 7, 337–349.

19 20

(38) Grigorenko, B. L.; Nemukhin, A. V.; Polyakov, I. V.; Morozov, D. I.; Krylov A. I. FirstPrinciples Characterization of the Energy Landscape and Optical Spectra of Green Fluorescent

21

Protein along the A→I→B Proton Transfer Route. J. Am. Chem. Soc. 2013, 135, 11541-11549;

22 23

(39) Kramer, Z. C.; Takahashi, K.; Skodje, R. T. Water Catalysis and Anticatalysis in Photochemical Reactions: Observation of a Delayed Threshold Effect in the Reaction Quantum

24 25 26 27 28 29 30 31 32 33 34 35

Yield. J. Am. Chem. Soc. 2010, 132, 15154-15157. (40) Zhang, Q.; Chen, X.; Cui, G.; Fang, W.; Thiel, W. Concerted Asynchronous Hula‐Twist Photoisomerization in the S65T/H148D Mutant of Green Fluorescent Protein. Angew. Chem. Int. Ed. 2014, 53, 8649-8653. (41) Kennis, J. T. M.; Larsen, D. S.; van Stokkum I. H. M.; Vengris, M.; van Thor J. Uncovering the hidden ground state of green fluorescent protein J. Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 17988–17993. (42) Seifert, M. H. J.; Georgescu, J.; Ksiazek, D.; Smialowski, P.; Rehm, T.; Steipe, B.; Holak, T. A. Backbone dynamics of green fluorescent protein and the effect of histidine 148 substitution. Biochemistry. 2003, 42, 2500-2512. (43) Oltrogge, L. M.; Wang, Q.; Boxer, S. G. Ground-State Proton Transfer Kinetics in Green Fluorescent Protein Biochemistry. 2014, 53, 5947-5957

ACS Paragon Plus Environment

20

Page 21 of 28 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

1 2 3 4 5 6

The Journal of Physical Chemistry

(44) Bindels, D. S.; Haarbosch, L.; Van Weeren, L.; Postma, M.; Wiese, K. E.; Mastop, M.; Aumonier, S.; Gotthard, G.; Royant, A.; Hink, M. A.; et al., mScarlet: a bright monomeric red fluorescent protein for cellular imaging. Nat. Methods. 2016, 14, 53-56. (45) Mauring, K.; Deich, J.; Rosell, F. I.; Mcananey, T. B.; Moerner, W. E.; Boxer, S. G. Enhancement of the Fluorescence of the Blue Fluorescent Proteins by High Pressure or Low Temperature. J. Phys. Chem. B. 2005, 109, 12976-12981.

7 8

(46) Snellenburg, J.; Laptenok, S.; Seger, R.; Mullen, K. M.; Van Stokkum, I. H. M. Glotaran: A Java-based graphical user interface for the R package TIMP J. Stat. Software, 2012, 49, 1-22.

9 10

(47) Koti, A. S. R.; Periasamy, N. Time-Resolved Area-Normalized Emission Spectroscopy (TRANES):  A Novel Method for Confirming Emission from Two Excited States. J.Phys. Chem.

11 12 13 14 15 16

A, 2001, 115, 7094-7099. (48) Yoon, E.; Konold, P. E.; Lee, J.; Joo, T.; Jimenez, R. Far-Red Emission of mPlum Fluorescent Protein Results from Excited-State Interconversion between Chromophore Hydrogen-Bonding States. J.Phys. Chem. Lett. 2016, 7, 2170−2174.

17 18 19 20 21 22 23 24 25 26 27 28 29 30

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry 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 22 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14

Figure 1. Structural representation of avGFP and rGFP proteins. (A) Surface model (1) of the

15

avGFP where sphere represents the solvent accessible spheres (SAS). (2) SAS near the

16

chromophore surface view. (B) Molecular structure of rGFP (1), avGFP (2) chromophore and

17

selected residues in their vicinity are shown; Superimposed model (3) of rGFP, an avGFP

18

chromophore and selected residues in the vicinity of the chromophores are shown. Hydrogen bonds

19

are represented as a dashed black line. (C) Electrostatic surface map of avGFP and rGFP.

20 21 22 23 24

ACS Paragon Plus Environment

22

Page 23 of 28 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

1 2 3 4 5 6 7

Figure 2. Dihedral angle and a radial plot representing residual variation. Torsional density

8

variation of ethaline linker moiety in avGFP and rGFP (A1-4) with dihedral angle (ø,τ) density

9

representing cis-trans isomerism probability in avGFP and rGFP, respectively. (B1-4) distance pair

10

principal component analysis (PCA) of the residual interatomic vector between Chromophore

11

Arg96 and Glu222, respectively.

12 13 14 15 16

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry 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 24 of 28

1 2 3

4 5 6 7 8 9 10

Figure 3. Graphical representation of avGFP and rGFP to evaluate the ESPT function. (A) Van der Waal

11

space of avGFP (1) and rGFP (2, 3) site-specific and global substitution and (B) Solvent sphere near the

12

chromophore. (C) Predicted protonation pathway of rGFP chromophore in trans-form (1) and avGFP in cis

13

form (2) significant interacting residues around the chromophore. (D) Mechanism of proposed

14

chromophore electrostatic stability.

ACS Paragon Plus Environment

24

Page 25 of 28 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

1 2 3 4 5 6 7 8 9 10 11 12 13 14

Figure 4. (A) Temperature dependent steady-state fluorescence of GFP and rGFP at pH 7.4 from 225 to

15

325 K. (B-C) Normalized time-resolved fluorescence decay for avGFP (B) and rGFP (C) equilibrated with

16

H2O (black) and D2O (red) at 515 nm and 610 nm , respectively as obtained by fluorescence upconversion

17

spectroscopy to understand the ESPT. Excitation occurs at 1 ps on the time axis. (D) Time resolved

18

fluorescence spectroscopic analysis of avGFP, rGFP and mRFP dissolved in 50 mM Phosphate buffer pH

19

7.4. (E-F) The avGFP and rGFP were prepared at a similar concentration, equilibrated with Benzoquinone

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry 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 26 of 28

1

(BQ), excited at the 400 nm and 460 nm and monitored the fluorescence emission spectra of avGFP (E)

2

and rGFP (F) at 525 nm and 610 nm to evaluate the ESPT function.

3

4

5

6

7

8

9

10

11

12

13

14

Figure 5. (A) Femtosecond transient absorption spectra of avGFP (1) and rGFP (2) in 0.1M phosphate

15

buffer at pH 7.4 obtained by exciting at 500 nm at RT. The grey dashed line represents the steady-state

ACS Paragon Plus Environment

26

Page 27 of 28 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

1

absorption spectra. (B) Kinetic profile at 420 and 520 nm obtained from femtosecond transient absorption

2

spectra of avGFP (red) and rGFP (black) at pH 7.4. Inset of panel B2 shows the expanded kinetic trace at

3

520 nm for both avGFP and rGFP.

4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

ACS Paragon Plus Environment

27

The Journal of Physical Chemistry 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 28 of 28

1

Figure 6. Time-resolved intensity (A) and (B) area normalized emission spectra (TRANES) of

2

rGFP excited at 375 nm. Time increases in the direction of the arrow.

3 4 5 6 7 8 9 10 11 12 13

TOC GRAPHICS

14 15 16 17

ACS Paragon Plus Environment

28