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