Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Stalling chromophore synthesis of the fluorescent protein Venus reveals the molecular basis of the final oxidation step

Husam Sabah Auhim ab, Bella L. Grigorenko cd, Tessa K. Harris a, Ozan E. Aksakal a, Igor V. Polyakov cd, Colin Berry a, Gabriel dos Passos Gomes ef, Igor V. Alabugin *g, Pierre J. Rizkallah h, Alexander V. Nemukhin *cd and D. Dafydd Jones *a
aSchool of Biosciences, Molecular Biosciences Division, Cardiff University, Sir Martin Evans Building, Cardiff, CF10 3AX, UK. E-mail: jonesdd@cardiff.ac.uk; Tel: +44 (0)29 2087 4290
bDepartment of Biology, College of Science, University of Baghdad, Baghdad, Iraq
cChemistry Department, Lomonosov Moscow State University, Leninskie Gory, 1-3, Moscow, Russian Federation. E-mail: anem@lcc.chem.msu.ru; Tel: +7 495 939 1096
dEmanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russian Federation
eDepartment of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON M5S 3H6, Canada
fDepartment of Computer Science, University of Toronto, 214 College St., Toronto, Ontario M5T 3A1, Canada
gDepartment of Chemistry and Biochemistry, Florida State University, Tallahassee, Fl 32306, USA. E-mail: alabugin@chem.fsu.edu; Tel: +1 850 644 5795
hSchool of Medicine, Cardiff University, CF14 4XN, UK

Received 7th December 2020 , Accepted 26th March 2021

First published on 31st March 2021


Abstract

Fluorescent proteins (FPs) have revolutionised the life sciences, but the mechanism of chromophore maturation is still not fully understood. Here we show that incorporation of a photo-responsive non-canonical amino acid within the chromophore stalls maturation of Venus, a yellow FP, at an intermediate stage; a crystal structure indicates the presence of O2 located above a dehydrated enolate form of the imidazolone ring, close to the strictly conserved Gly67 that occupies a twisted conformation. His148 adopts an “open” conformation so forming a channel that allows O2 access to the immature chromophore. Absorbance spectroscopy supported by QM/MM simulations suggests that the first oxidation step involves formation of a hydroperoxyl intermediate in conjunction with dehydrogenation of the methylene bridge. A fully conjugated mature chromophore is formed through release of H2O2, both in vitro and in vivo. The possibility of interrupting and photochemically restarting chromophore maturation and the mechanistic insights open up new approaches for engineering optically controlled fluorescent proteins.


Introduction

Fluorescent proteins (FPs) represent an important family of proteins that emit light in the visible region of the spectrum without the requirement of any additional cofactor.1–4 Their unique fluorescence properties have revolutionised the life sciences through their use as genetically encoded imaging probes and sensors.5,6In situ synthesis of the chromophore is critical and occurs through the stepwise covalent arrangement of three contiguous residues Xaa-Tyr-Gly (where Xaa can be various amino acid residues), with only molecular O2 needed (Fig. 1a). In the case of Venus, these three residues are 65-Gly-Tyr-Gly-67 (Fig. 1a). Venus is a yellow fluorescent derivative of the classical Aequorea victoria GFP protein.7 Amongst the various mutations introduced to generate this yellow variant, two key changes are the introduction of glycine at residue 65 (from serine), which is thought to influence chromophore flexibility,8 and tyrosine at residue 203 (from threonine); Tyr203 π-stacks with the chromophore causing a red shift in excitation and emission into the yellow region.7,9 The mature chromophore has an extended conjugated bond network comprised of the phenolic P ring and the imidazolone I ring linked via a methylene bridge (see Fig. 1a for ring and atom notation), which is located within the core of the β-barrel structure.1,4,10,11 The first step in chromophore formation is cyclisation through linkage of the backbone amine of the strictly conserved Gly67 residue to the carbonyl carbon of the Xaa 65 residue. Cyclisation is followed by the dehydration/O2-dependent oxidation steps to generate a fully conjugated system. There has been much discussion about the exact mechanism of chromophore formation, especially the order of the dehydration and dehydrogenation/oxidation steps,4,8,10–12 and the role O2 plays in the rate-limiting oxidation step.13 The availably of structures representing trapped intermediates has helped shed some light on potential mechanisms14–19 but the exact role O2 plays in the maturation process still needs to be fully addressed. This is in part hampered by neither the crucial O2-bound form nor the oxygenated intermediates being experimentally observed. The former is especially critical as it can act as a realistic starting point for theoretical simulations. Gly67 is another key residue that is critical to chromophore synthesis and is strictly conserved amongst fluorescent proteins,17,20 but its role in maturation is unclear.
image file: d0sc06693a-f1.tif
Fig. 1 Effect of azF incorporation at residue 66 of Venus. (a) Scheme outlining the basic maturation of Venus. The bond and ring nomenclature are described in the lower structure. (b) Solution colours of VenusWT and Venus66azF. (c) Solution colours of Venus66azF before (dark) and after (UV) illumination with UV light. (d) Absorbance (solid line) and fluorescence (dashed line) of Venus66azF. Red, light blue, grey, orange, green, purple, dark blue represent spectra after 0, 1, 5, 10, 15, 30, 45 and 60 min. Full spectral properties are shown in ESI Table S1. Full fluorescence emission time course and in vivo imaging are shown in ESI Fig. S1.

Manipulating the chemical properties of the chromophore by protein engineering, either directly through changes to two of the three chromophore residues, or indirectly through changing the chromophore environment, have generated a range of new fluorescent proteins, including Venus itself,7 with properties suited to their particular application.6 One of the most important FP class for super-resolution imaging is the photo-controllable FPs, whereby fluorescence is either switched on/off, or spectral properties significantly shifted in response to light.21–24 Mechanisms of action involves chemical modifications such as decarboxylation of Glu222 (e.g. PA-GFP25), backbone cleavage (e.g. Kaede26) and chromophore hydration (e.g. Dreiklang27), or conformational changes such as chromophore cis/trans isomerisation (e.g. Dronpa,28 rsEGFP29). The use of photochemically active non-canonical amino acids (ncAA) has further expanded optical control approaches.30 Phenyl azide photochemistry is particularly useful as we have used it previously to turn on, off, or switch the fluorescence properties of green31,32 and red33 FP types. Replacement of the chromophore forming tyrosine residue in super-folding GFP (sfGFP)32 and mCherry33 is known to impede fluorescence until irradiated; the azide group is thought to act as an excited-state quencher until converted to the amine.

Here we use the photochemical properties of genetically encoded phenyl azide to stall Venus chromophore maturation at an immature non-fluorescent intermediate (termed im-Venus66azF) state before UV irradiation instigates maturation to a final fluorescent form. The structure of the intermediate reveals the protein has undergone the dehydration but not the oxidation step. Additional new structural features add further new insights, including an essential role for the strictly conserved Gly67 and, for the first time, experimental observation of a putative molecular O2 in direct proximity to an enolate form of the dehydrated immature chromophore. The combination of experimental spectroscopy with quantum mechanics/molecular mechanics (QM/MM) simulations allowed us to propose a mechanism for the O2 dependent oxidation step whereby a hydroperoxyl intermediate is formed as part of the oxidation mechanism.

Results

Incorporation of phenyl azide chemistry within the Venus chromophore

To incorporate a photochemical switch into Venus, we replaced the chromophore residue Tyr66 with the aromatic non-canonical amino acid p-azido-L-phenylalanine (azF);34,35 effectively the hydroxyl group of tyrosine is replaced with an azide group. The variant termed Venus66azF is produced as a colourless protein (Fig. 1b) with no inherent fluorescence (Fig. 1b–d and S1), and acts an excellent starting point for an optically controlled FP. Mass spectrometry confirmed that azF was incorporated and that the chromophore had at least gone through the cyclisation–dehydration step (Fig. S2). The absorbance spectrum of the colourless non-fluorescent Venus66azF reveals the protein absorbs in the UV region. Absorbance peaks at ∼350 nm (∼3.5 eV) with a significant shoulder at ∼360–370 nm (∼3.44–3.35 eV); no fluorescence is observed from this species (Fig. 1d and S1a, b). Irradiation of Venus66azF with low intensity (6 W from a handheld lamp) near-UV light converts the protein to a coloured and fluorescent form both in vitro (Fig. 1c, d and S1a, b) and in vivo (Fig. S1c and d), with an absorbance λmax shifting to a broad 440 nm peak (2.82 eV) (ESI Table S1), similar to that observed previously for the amine chromophore product in sfGFP.32 Fluorescence emission increased over 100-fold on activation, allowing imaging of bacterial cells expressing the activated Venus66azF (Fig. S1d). Spectroscopic analysis of the time course of activation reveals a potential two-step process. Within 1 min of UV exposure, the double hump absorbance spectrum feature disappears, with the remainder of the time points forming a clear isosbestic point at ∼390 nm (3.18 eV) suggesting conversion directly from one form to another. As the initial dark sample diverges away from the isosbestic point, it suggests a primary conversion to an intermediate species, which is also non-fluorescent, followed by a slower conversion to the activated form. The final product of UV irradiation was the phenyl amine version of the chromophore, as observed previously for sfGFP66azF,32 but it appeared sensitive to fragmentation on denaturation due to cleavage close to the chromophore (see ESI Fig. S2 and associated discussion); fragmentation within the chromophore on denaturation is commonly observed for GFP-like proteins.24,36 SDS-PAGE analysis of the equivalent UV timepoints suggests that UV irradiation does not generate a completely fragmented protein (ESI Fig. 2); after 5 min fragmentation remains relatively constant (20–25%), in line with previous denaturation-dependent fragmentation observations.24,36 The mass spectrometry analysis did not indicate a secondary decarboxylation (e.g. Glu222) event that can sometimes be associated with UV irradiation of fluorescent proteins.37 UV irradiation of VenusWT did not result in any spectral shifts (ESI Fig. S1e and f).

Structure of im-Venus66azF

The structure of the pre-photoactivated Venus66azF (from herein referred to as immature or im-Venus66azF) was determined to 1.9 Å resolution (see ESI Table S2 for structural statistics) to gain an insight into the molecular nature of the non-fluorescent, immature chromophore. The general structure of im-Venus66azF is similar to mature wild-type Venus (henceforth termed VenusWT), with a root-mean-squared deviation (RMSD) over the backbone of 0.39 Å (ESI Fig. S3a). The most significant differences occur in and around the chromophore (termed CRO). Electron density fits well to a cyclised chromophore containing the intact azide group present on the P ring (Fig. 2 and S3b). The lack of the electron density protruding from the I ring at the C65 position suggests the dehydration step in chromophore maturation has occurred by this point (ESI Fig. S3b and c). Compared to VenusWT, the I ring component shifts position in im-Venus66azF due to the more acute angle in the methylene bridge (Fig. S4b). A comparison of the residues surrounding the CRO shows that several residues exist in different conformations compared to Venus (Fig. S4). His148, a critical H-bonding residue to the mature CRO, exists fully in the open-gate conformation. Glu222, another residue necessary for function, and Tyr203 (π-stacks with the CRO in VenusWT) also shift position with respect to the CRO (Fig. 2a). The newly introduced phenyl azide group occupies a similar position as the tyrosyl residue in Venus (Fig. 2a and S4). The azide moiety itself is located between Gln204-Ser205 and Tyr145-Ser147, displacing a water molecule normally present in Venus, with Ser205 hydroxyl group making potential polar interactions with the azide.
image file: d0sc06693a-f2.tif
Fig. 2 Structure of im-Venus66azF proximal to the chromophore. (a) Comparison of im-Venus66azF (grey; PDB 6sm0) with Venus (yellow; PDB 1mwy9). CRO is the chromophore (Gly65-Tyr/AzF66-Gly67). (b) Position of the O2 molecule in im-Venus66azF. The top panel shows the electron density (2Fo-Fc, 1.0σ) for the CRO, Y203 and O2 together with an additional water molecule. The lower panel removes the electron density for clarity. Relevant distances are shown in Å.

Further analysis of the chromophore, including comparative estimates of bond distances between the refined structure of im-Venus66azF with the refined Venus structure, reveals several novel features that provide us with insights into CRO maturation and the role of molecular O2 in the process. In im-Venus66azF, the backbone carbonyl of the strictly conserved Gly67 is twisted ∼180° out of position compared to that observed in other FPs (Fig. 2a, b and ESI Fig. 3d, e). To our knowledge, the only other time this conformation has been observed is in the unpublished structure of an immature chromophore of a GFP maturation disabling mutant determined by the Getzoff group (PDB 2qt2; ESI Fig. S5). There is additional electron density sandwiched between the Gly67 carbonyl group and the hydroxyl group of Tyr203 (Fig. 2b and ESI Fig. 3d and e), which we have assigned to molecular O2 after attempting to refine the structure with either one or two H2O. We found the presence of Tyr203 effectively blocks the ability of two water molecules to occupy this position whereas one water molecule left an elongated tail of density. We found that molecular O2 fitted best; the difference map is featureless around the elongated density fitted to O2 with no positive (no atoms unaccounted for) nor negative difference (atoms that should not be where they are modelled), as can be seen in Fig. 2b. The O2 molecule lies between the twisted carbonyl oxygen of the Gly67 and the hydroxyl group of Tyr203 above the plane of I ring element (Fig. 2b). O2 has been postulated to be positioned either above the plane of the chromophore facing Glu222/Tyr203 or below the chromophore plane facing Arg96;18,38 here O2 is above the plane of the chromophore on the Glu222/Tyr203 face (Fig. 2b).

The structure also provides an insight into the nature of the trapped chromophore intermediate (Fig. 3). During maturation, O2 is thought to be involved in generating the final C[double bond, length as m-dash]C that links the P-ring and I ring (C66β–C66α to C66β[double bond, length as m-dash]C66α). Comparative analysis of the C66β–C66α bond lengths in the refined structure suggests it could be a single bond in im-Venus66azF and a double bond in Venus (Fig. 3). The bond angle between C66γ–C66β–C66α is also more acute for im-Venus66azF (118° versus 131° for Venus; Fig. 3). In the I-ring, the C66–O66 bond is 1.48 Å for im-Venus66azF, which is longer than would be expected for a keto-carbonyl C[double bond, length as m-dash]O bond (1.20 Å) as observed in Venus (Fig. 3b). This longer C–O bond also makes polar contacts with the critical maturation residue Arg96, which occupies a near-identical position in VenusWT (Fig. 2a). Thus, we predict that the enolate is the most likely form of C66–O66 (Fig. 3a) with Arg96 stabilising the negative charge. The negative charge on the enolate can be offset by a positive charge delocalised around the I ring, as proposed in Fig. 3a.


image file: d0sc06693a-f3.tif
Fig. 3 Chromophore structures of im-Venus66azF and mature VenusWT fluorescent proteins. (a) im-Venus66azF with proposed chemical structure; (b) VenusWT (PDB code 1myw9) with the chemical structure. A comparative analysis of the selected chemical bond lengths is shown.

Simulation of the oxidative step of chromophore maturation

To correlate the structural and functional observations, modelling of reaction intermediates was performed using the atomic coordinates of im-Venus66azF obtained in this work as a starting point. Structures of potential intermediates were assumed on the basis of previous experience in modelling chromophore maturation in GFP38 and optimized in QM/MM calculations. Analysis of the obtained results for Venus66azF shows that mechanistic steps during the oxidation step of chromophore maturation are essentially the same as in the wild-type GFP.38

The overall scheme based on the simulations is shown in Fig. 4, and the corresponding structures were associated with the observed absorbance data. The simulations show that the first step involves Gly67 switching to its energetically more favourable (∼7 kcal mol−1) canonical configuration. Thus, it appears we were fortunate that the crystalline form of im-Venus66azF was trapped in the observed conformation shown above. Triplet state oxygen can now access the I ring with concomitant protonation of Glu222, which acts as a general acid/base in the maturation scheme. The oxidation steps then proceed starting from partial negative charge transfer to O2, which switches from the triplet to singlet state. Glu222 is protonated with N65 donating the proton (Fig. S6). The O2 then attacks C65 (and not C66α) generating a peroxy intermediate (Fig. 4 and S5). The peroxy anion then abstracts a proton from C66β to form the stable hydroperoxyl intermediate with Glu222 protonating N65 (Fig. 4 and S6); the hydroperoxyl species has a theoretical absorbance of 360 nm (3.44 eV). As well as Gly67 converting through to its canonical conformation, the formation of the methylene bridge between the I and P rings results in a shift to a configuration similar to that observed for the mature chromophore observed in VenusWT; the methylene bridge bond angle is now 134° (Fig. 5a), similar to that observed for VenusWT (Fig. 3b).


image file: d0sc06693a-f4.tif
Fig. 4 Reaction scheme for oxygenation of the chromophore based on modelling results. The panels show relevant segments of the assumed protein structures optimized in QM/MM calculations. Triplet and singlet systems are differentiated by black and green lines, as indicated in the diagram. The ΔE between each state is shown in red with units of kcal mol−1. Significant changes between each state are outlined in the diagram.

image file: d0sc06693a-f5.tif
Fig. 5 Models of the hydroperoxyl intermediate. (a) Model structure as determined by QM/MM simulations with the hydroperoxyl group bound to C65 as indicated. On the right-hand side is a comparison of the modelled hydroperoxyl intermediate (cyan) with the crystal structure of im-Venus66azF. (b and c) Stereoelectronic interactions in hydroperoxyl intermediate (hydrogen atoms omitted for clarity with the numbers related to energies in kcal mol−1). Both the hydroperoxyl moiety acting as a σ-acceptor (b) and p-donor (c) are shown.

We also undertook natural bond orbital (NBO) analysis to evaluate the hydroperoxyl intermediate. NBO analysis transforms electron density from DFT calculations into localized orbitals that are closely tied to the chemical bonding concepts. In particular, NBO analysis is commonly used to evaluate hyperconjugative stabilising interactions.39,40 The lone pair of the hydroperoxyl group (–OOH) in the model has near-ideal alignment with the two C–N bonds (Fig. 5b and c). This favourable stereoelectronic arrangement activates stabilizing hyperconjugative image file: d0sc06693a-t1.tif interactions, which can partially compensate for the loss of aromatic stabilisation in the I ring. This is complemented by two strong image file: d0sc06693a-t2.tif interactions. The importance of the latter effect is expected to grow further in the transition state for the final C–O bond scission where it provides an important transition state stabilization effect that can significantly assist the final aromatizing step of the cascade.41–43

The final step occurring over the irradiation period is the conversion to a fully mature fluorescent chromophore. Two events need to be considered: reduction of the azide and full conjugation of the chromophore through a loss of the hydroperoxyl moiety. The spectral properties suggest that the final end product is likely to be the phenyl-amine form of the mature chromophore as has been observed before in superfolder version of GFP (sfGFP; Table S1 and ref. 32); this was confirmed by mass spectrometry (ESI Fig. S2). Simulations concur with this with the final product having a theoretical absorbance max at 444 nm (2.79 eV), close to the 440 nm (2.82 eV) observed in Fig. 1d. Full chromophore conjugation with the azide left intact will generate a species less stable than the preceding step and has a theoretical absorbance maximum at 451 nm (2.75 eV) (Fig. 6a). The alternative route appears more likely: reduction to the phenylamine followed by loss of the hydroperoxyl group (generating H2O2). The phenylamine version of hydroperoxyl intermediate has a theoretical absorbance maximum of 367 nm (3.38 eV). Based on the experimental time course observed in the Fig. 1d, the initial species mix dominated by the hydroperoxyl intermediate (computed λmax = 360 nm, 3.44 eV) within 1 min of irradiation converts to the dominant hydroperoxyl phenylamine (computed λmax = 367 nm, 3.38 eV) that then forms the mature phenylamine chromophore.


image file: d0sc06693a-f6.tif
Fig. 6 Conversion of phenyl azide intermediate to mature chromophore. (a) Potential routes to mature phenyl amine chromophore from the hydroperoxyl intermediate. The two routes are azide reduction followed by full chromophore conjugation (black route) and vice versa (grey route). The theoretical λmax values are shown. It should be pointed out that the presumed species without the hydroperoxyl and phenyl azide present (λmax = 451 nm, 2.75 eV) is considerably higher in energy than the hydroperoxyl intermediate. Therefore, the upper grey route is unlikely, as also clarified in (b). (b) Mechanistic analysis highlights the contrasting effects of azide and amine substituents on the elimination of hydroperoxide in the final aromatisation step. The decelerating effect of the acceptor azide on the OOH elimination is deactivated once the azide is converted into a strongly donating amine moiety.

Oxygen access to the chromophore and the twisted Gly67

The structure of im-Venus66azF revealed the presence of O2 within the core of the protein close to the I ring of the chromophore. There are two potential ways in which O2 becomes located in such a position: during the folding process or diffusion into the core after folding to the β-can structure and initial chromophore cyclisation. Simulations, whereby O2 is replaced by a water molecule, reveal that the twisted conformation of the immature chromophore is more stable by 18.1 kcal mol−1 than the canonical form (Fig. 7a). The presence of water is unsurprising given that O2 is not required for the preceding cyclisation and dehydration events. Thus, water helps stabilise the twisted form of Gly67 present in im-Venus66azF. If this is the case, O2 is likely to displace the water molecule as part of the maturation mechanism, meaning it has to gain access to the core of the protein. Analysis of the im-Venus66azF structure reveals that the “open” His148 conformation may play a role (Fig. 7b and c). A channel leads to the chromophore in im-Venus66azF that is only available if His148 occupies the “open” state; when His148 occupies the “closed” state normally observed in mature fluorescent proteins, the tunnel becomes blocked.
image file: d0sc06693a-f7.tif
Fig. 7 Role of water and internal tunnels in chromophore maturation. (a) Simulation of the twisted conformation stabilised by a water molecule in place of O2 (and for reference, the canonical form usually observed in mature VenusWT conformation). The twisted form is more stable by 18.1 kcal mol−1 compared to the canonical form. The water molecule of interest is highlighted. Chromophore accessibility, as shown by (b) spheres and (c) CAVER tunnel analysis.44 The alternative conformation for His148 (yellow sticks) in VenusWT is shown.

Discussion

Photo-controllable FPs have become an important tool in modern super-resolution cell imaging.45 The use of phenyl azide chemistry to control fluorescence, here and in green31,32 and red33 fluorescent proteins, provides a simple and general common mechanism to implement photocontrol across a broad range of the FP colour palette. Indeed, the use of a ncAA in conjunction with a reprogrammed amber stop codon allows secondary control: production of the fusion protein constructs with and without the FP adjunct.46 While the demonstration of the conversion of an essentially colourless immature fluorescent protein to an active form through phenylazide photo-decaging confirms the approach as a means to photo-control FPs, insights into the chromophore maturation process are arguably the most important aspect of this work.

Venus is an engineered version of the original Aequorea victoria GFP, that includes mutations that comprise the chromophore (S65G) and directly interact with the chromophore (e.g. T203Y) so contributing towards its red-shifted fluorescence properties. There are currently two models that describe the overarching process by which chromophore maturation in A. victoria derived fluorescent proteins occur:4 “cyclisation–oxidation–dehydration”12,15 and “cyclisation–dehydration–oxidation”13,18,47 both of which are supported by experimental data albeit through using different approaches. Our data supports the “cyclisation–dehydration–oxidation model, at least in the context of Venus66azF and under the conditions we used. As has been pointed out previously, both maturation processes may indeed occur in parallel with factors such as oxygen concentration and local mutational events influencing the order.4 In our study, the experimentally obtained trapped intermediate structure of im-Venus66azF provided evidence that dehydration had already taken place (Fig. S3) but oxidation had yet to occur (Fig. 3). Coupled with the simulation data, we propose a chromophore synthesis route for Venus66azF outlined in Scheme 1.


image file: d0sc06693a-s1.tif
Scheme 1 Proposed mechanism for the chromophore maturation, including details of the final oxidation step.

The crystal structure of the im-Venus66azF provides evidence for the nature of an intermediate prior to the final oxidation step: the I ring in the enolate form with the C66α and C66β forming a single bond (Fig. 3a). The long C66–O66 bond length is indicative of the enolate whose negative charge is stabilised by critical chromophore maturation residue, Arg96 (Fig. 2a).16,20,48 The crystal structure also provides evidence of the location of O2. Two relative positions with respect to the chromophore plane have been proposed: on the Arg96 (ref. 18) or the Glu222 face.38 Here we show that in our system molecular oxygen is placed on the Glu222 face, directly above the I ring (Fig. 2), in agreement with recent simulation data for GFP.38 Arg96 has been suggested as the oxygen activator through the positively charged side-chain18 but this is unlikely to be the case here and may instead play a role in stabilising the enolate form of the I ring.

The twisted Gly67 configuration is clearly observed in the crystal structure and differs from the canonical position normally present in FPs (Fig. 2). Glycine has a less restricted ψ angle range. In Venus, the Gly67 ψ dihedral angle is −23° compared to 167° in im-Venus66azF, with twisted conformation energetically less favourable (by 7 kcal mol−1) when O2 is present and also hinders access of oxygen to the I ring (Fig. 2). The twisted conformation could be a legacy of the cyclisation reaction whereby the nucleophilic attack of N67 on C65 will require rotation of the Gly67 ψ angle leading to its observed placement in im-Venus66azF. While we were fortunate to have trapped this immediate form in the crystalline state, the configuration could be stabilised by the O67 H-bonding with the C65 hydroxyl group (observed by Getzoff and colleagues17) before condensation. The simulation data also revealed that a single spatially defined water molecule can stabilise the twisted conformation over the canonical form more so than the O2 bound form (Fig. 5a) that our data suggests is likely be present in the crystal structure. This could suggest that a hydrated form precedes the oxygen bound step. It is interesting to speculate that the origin of the water molecule may be the product of the cyclisation/dehydration reaction that precedes oxidation. Indeed, while our crystal structure of im-Venus66azF strongly indicates O2 is present, we cannot rule out that a population of the structures has a water molecule present in the same position: fitting the observed electron density to a dynamic water over two sites suggests this could be feasible. However, we feel the evidence for O2 being the main species is: (1) if the O in water takes the position of one or other O atoms in the O2 molecule, this will generate steric issues due to the presence of the H atoms; (2) the position of O2 in the structure correlates with the QM/MM work here and elsewhere38 and available tunnels through to the chromophore with H148 in the open configuration. The most likely scenario is that there is mixture of O2 and water occupying the site, with O2 occupancy likely being the dominant population observed in our crystal structure. It is interesting to speculate that the alternative conformation of G67 may induce a dipole in O2 due to the close proximity of the δ on O67, resulting in an interaction network spanning to Tyr203.

If a water molecule originally stabilised the twisted Gly67 conformation, this suggests that O2 needs to access the protein core for the final oxidation step to take place. We propose that His148 plays a key role in this process. His148 is dynamic49 and has been observed in both the “open” and “closed” conformations with the former not normally reported in the crystal structure as it is a minor component when observed (for example Arpino et al.50 Reddington et al.32 and Brejc et al.51); the closed conformation is the major form observed in mature FPs as His148 in this configuration H-bonds to the chromophore and plays a critical role in function.1,10,52,53 Im-Venus66azF almost exclusively exists in the open conformation (Fig. 2a and S4) that generates a channel through to the chromophore (Fig. 7b and c). Such a tunnel at a similar position has been observed previously for GFP-like proteins.54,55 In Dreikling, a reversible photo-switchable close relative of Venus, H148 exists predominantly in its “open” conformation that can potentially allow water access to the chromophore as part of the hydration event that underlies photochemical control.27,56 Thus, His148 may acts as a “gatekeeper” residue, so determining access to the chromophore as well as its functional role (H-bond to chromophore phenol group).49 Given that oxidation is the rate-limiting step in maturation,4 it is interesting to speculate that the exchange rate between the two His148 conformations may play a role in defining this rate.

Simulations suggest that the next dominant form is hydroperoxyl intermediate attached to C65 and not C66α, as suggested by others.18 The presence of a hydroperoxyl intermediate attached to C65 has been proposed previously for GFP.38 Computational modelling of GFP suggests that O2 may form a bridge between the C65–C66α with a concerted proton transfer mechanism initially from C66β resulting in the formation of the C65 attached hydroperoxyl intermediate and thus the β-methylene bridge.38 The argument against attachment to C66α as the intermediate comes from the observed spectra data (Fig. 2b), whereby the dominant 350–360 nm peaks for the intermediate suggests some extension of the conjugated double bond system (here proposed to be from the phenyl azide to C66α–C66β; Fig. 5). The formation of a hydroperoxyl intermediate at C66α would prevent the formation of the double bond with C66β. The electron lone pair on the C65 hydroperoxyl moiety also aligns perfectly with the two C–N bonds that, in turn, helps stabilise the five-membered I ring. We propose that the nearby Glu222 plays a vital role acting as a general acid/base during the formation of the hydroperoxyl intermediate through the first abstraction and then the donation of a proton to N66 (Fig. S6). The importance of Glu222 to maturation has been observed previously, with the E222Q mutation in EGFP considerably slowing maturation.20,48 During this process we propose that a peroxyanion is formed, which abstracts the proton from the activated C66β to form hydroperoxyl intermediate (Fig. 4 and S6), as suggested for GFP.38 Thus, the formation of the C66α[double bond, length as m-dash]C66β double bond occurs before the generation of H2O2 and not concurrently (Scheme 1).

The proposed final and rate-limiting step in the process is the formation of the fully conjugated fluorescent chromophore. In our system, we believe this is a UV induced phase that happens in two steps due to the presence of the azido group: (1) conversion of the phenyl azide to the phenylamine; (2) loss of the hydroperoxyl group so generating a fully conjugated system. In our proposed model, we suggest that the reduction of the azide to an amine occurs first (Scheme 1). This is based on the computed theoretical absorbance of each species (Fig. 6a) and on the expectation that azide conversion to strongly donating amine group would significantly help with loss of the hydroperoxyl moiety (see Fig. 6b for mechanistic details). Furthermore, the departure of the OOH group needed for the conversion of hydroperoxyl intermediate to the final maturated chromophore is expected to be greatly facilitated when azide, a mild acceptor (Hammett parameter σp = +0.08), is changed to NH2, a strong donor (σp = −0.66). In the azide-substituted peroxide, the lone pair of N66 is not able to fully assist in the departure of the OOH group as its electron density is partially delocalized in the other direction, towards the azide. This stereo-electronic tug-of-war is removed once the amine is formed. As the NH2 group is a powerful donor, the electron density is no longer shifted from N66 to the aryl ring; the lone pair on N66 is now free to stabilize the transition state for heterolytic C⋯OOH bond scission. A fully conjugated chromophore with the azide group attached has a computed λmax of 451 nm, whereas the amine version of hydroperoxyl intermediate is 367 nm. Given the observed absorbance time course in Fig. 2 goes from a mixed species with two peaks between 340–360 nm to a single species at 360 nm that directly converts to the single species at 440 nm, the logical progression of the computed spectra are 360 nm (phenyl azide/hydroperoxyl form) to 367 nm (phenylamine/hydroperoxyl forms) to 444 nm (mature amine chromophore). UV light may also play a role beyond the initial conversion from the azide to the amine as prolonged exposure is needed to fully develop fluorescence (Fig. S1 and ref. 30). From comparison of the absorbance and fluorescence time course, the initial intermediate observed after 1 min is not itself significantly fluorescent nor is the initial species absorbing at 351 nm suggesting additional UV-induced events are needed for full fluorescence. It is not clear why the further input of UV light is needed but the final step from the hydroperoxyl intermediate requires energy to overcome the last reaction barrier involving proton transfer from the I ring to E222 as part of the process of H2O2 generation.38 Given that Venus66azF remains largely intact on UV exposure (Fig. S2b) and does not have a HYG chromophore motif, backbone fragmentation is unlikely to be the main end-point of activation as observed in green-to-red photoactive FPs such as Kaede26 and EosFP.57 While simulations provide strong evidence (Fig. 6) that conversion of the azide to the amine is the initial step in the process we cannot completely rule out that this step may happen later in the scheme.

Conclusion

For the first time, our X-ray structure of im-Venus66azF provides an experimental evidence concerning insight into the role played by oxygen in fluorescent protein chromophore maturation. The fortuitous trapping of oxygen between Tyr203 and an alternative conformation of the strictly conserved Gly67 before photo-decaging of azF66 was critical. From our structure and subsequent modelling, the cyclisation–dehydration–oxidation model mechanism is the likely route to chromophore maturation, at least for Venus66azF, but with new insights concerning the nature of the enolate intermediate together with mechanistic details of the O2-dependent steps. These new insights not only improve our fundamental understanding of chromophore maturation but may help aid our ability to generate improved fluorescent protein for future applications.

Methods and materials

Engineering, production and structure of Venus66azF

The generation of the mutant and subsequent recombinant production of Venus66azF is outlined in the ESI and based on a previously published procedure.55 Protein production and purification was carried out in the dark to prevent photolysis of azF. Crystallisation, diffraction data collection and structure determination of im-Venus66azF was performed as outlined in the ESI. A long, spiny crystal of im-Venus66azF was observed after a few days in the A12 condition of the PACT premier TM HT-96 screen (0.01 M zinc chloride, 0.1 M sodium acetate, pH 5 and 20% w/v PEG 6000). Data were collected at the Diamond Light Source (Harwell, UK) at beamline I04-1. Data reduction was completed with XIA2 (ref. 58) using XDS. Data scaling, merging and analysis were completed with AIMLESS59 and TRUNCATE60 in the CCP4 Package. Structure solution was made via molecular replacement with PHASER61 using PDB entry 4J88, the structure of sfGFP66azF. The structure was refined with REFMAC5,62 using the standard restraints supplied with the program. The geometry restraints for the chromophore were derived from ideal geometry with JLIGAND,63 allowing REFMAC to determine weights automatically. COOT64 was used for graphical sessions to adjust the model to match the electron density maps. The crystallographic statistics are provided in ESI Table 2.

Absorbance and fluorescence spectroscopy

Absorbance spectra were recorded using a Cary 60 spectrophotometer (Agilent) in a 1 cm pathlength quartz cuvette. The molar absorbance coefficients were calculated by recording absorbance spectra of Venus samples with a known concentration (5–10 μM) and then extrapolated to 1 M using the Beer–Lambert law. Emission spectra were recorded using a Cary Eclipse fluorimeter (Varian) using a 5 mm × 5 mm QS quartz cuvette. Protein concentration of 0.5 μM for VenusWT or 10 μM for Venus66azF was used with a fixed scan rate of 120 nm min−1 with a 5 nm slit width. The excitation wavelength for VenusWT and Venus66azF were determined from their absorbance spectrum. Bacterial live cell imaging was performed as outlined in the ESI.

Photolysis of Venus66azF

Photolysis experiments were carried out using a UVM-57 Handheld UV lamp (6 W; 302 nm UV, UVP Cambridge, UK) and 1 cm pathlength quartz cuvette (Hellma), essentially as described previously.31,54 A maximum of 500 μl of protein sample (10 μM) was pipetted into a cuvette and exposed to the UV (302 nm) for the indicated periods of time at a distance of 1 cm. The absorbance spectra and emission spectra were recorded immediately afterwards, as described above.

QM/MM simulations

Coordinates of heavy atoms of the im-Venus66azF crystal structure were used to construct a full atom three-dimensional model system. Previous examples in quantum mechanics/molecular mechanics (QM/MM) simulations of the mechanism of chromophore maturation included the wild-type GFP,38 its Gly65–Gly66–Gly67 mutant,65 and the recent modeling of Dreiklang,56,66 a photo-switchable protein with the chromophore formed from the same amino acid residues Gly65, Tyr66, Gly67, as in Venus. Hydrogen atoms were added manually using molecular mechanics tools; the side chains of Arg and Lys were assumed as positively charged, the side chains of Glu and Asp as negatively charged. The model protein molecule was fully surrounded by explicit water molecules.

Structures of possible intermediates in the maturation reaction were optimized in QM/MM calculations. A large fraction of the chromophore-containing pocket was assigned to the QM-part. The Gly65–azF66–Gly67 fragment of the immature chromophore (CRO), the side chains from Arg96, Tyr203, Ser205, Glu222 and 4 water molecules were included. This initial composition was considered to model structures without the oxygen molecule. In majority of calculations, the O2 species was inserted to the cavity near CRO. Calculations of energies and energy gradients in QM were carried out using Kohn–Sham DFT with the PBE0 functional67 and the cc-pVDZ basis set. The AMBER force field was used in MM. The NWChem software package68 was applied to scan fragments of potential energy surface. These scans along with the previous experience in modelling chromophore maturation in GFP35 allowed us to construct protein structures of potential intermediates, which were optimized in QM/MM calculations. To model the system in the triplet electronic state, the unrestricted DFT approach was used.

Vertical excitation energies at selected points on the ground state potential energy surface were computed using the extended multiconfigurational quasi-degenerate perturbation theory in the second order (XMCQDPT2)69 the protocol that we verified earlier and used extensively in studies of the photoreceptor proteins.70 Here, the perturbation theory calculations were based on the complete active space self-consistent field (CASSCF) wavefunctions obtained by distributing 16 electrons over 12 orbitals and using density averaging over 15 states. To perform these calculations using the Firefly quantum chemistry package,71 large molecular clusters including the QM parts of the system were selected. Natural Bond Orbital (NBO) analysis was used to evaluate stereoelectronic interactions.39,40 Geometry optimizations for NBO evaluations were performed with SMD72 for solvation corrections and the unrestricted wB97X DFT functional73 (with an integration grid of pruned 175[thin space (1/6-em)]974 for first-row atoms and 250[thin space (1/6-em)]974 for atoms in the second and later rows) with the 6-311++G(2d,p) basis set for all atoms. Grimme's D2 version for empirical dispersion74,75 was also included. Natural Bond Orbital (NBO) analyses were performed with NBO6 linked to Gaussian 16. They were used to gauge the magnitude of the hyperconjugative interactions in the presented systems.

Author contributions

All authors contributed to the writing of the paper and analysing data. HSA undertook structural biology and functional characterisation of the Venus proteins; TH helped prepare mutant proteins. OEA prepared mutant protein and generated mass spectra. PJR collected structural data and helped with structure determination and refinement, and with preparing the manuscript. BLG, IVP and AVN developed computational models and carried out computer simulations. IA and GPG provided analysis of stereoelectronic factors involved in the formation and reactivity of hydroperoxide intermediates. CB analysed data. DDJ conceived and directed the project, and contributed to data analysis.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

We would like to thank the staff at the Diamond Light Source (Harwell, UK) for the supply of facilities and beam time, especially beamline I04-1 staff, under beamtime code mx18812. We thank BBSRC (BB/H003746/1 and BB/M000249/1), EPSRC (EP/J015318/1). H. S. A. was supported by the Higher Committee for Education Development in Iraq. OEA was supported by an EPSRC DTP IDTH Physics for Life studentship. We would like to thank the Protein Technology Hub, School of Biosciences, Cardiff University for use of facilities and Cardiff University Department of Chemistry Analytical Unit for performing mass spectroscopy. A. V. N., B. L. G. and I. V. P. thank the Russian Science Foundation (project 19-73-20032) for support of the modelling part of this work. The calculations were carried out using the equipment of the shared research facilities of HPC computing resources at Lomonosov Moscow State University. The use of supercomputer resources of the Joint Supercomputer Center of the Russian Academy of Sciences is also acknowledged.

References

  1. R. Y. Tsien, Annu. Rev. Biochem., 1998, 67, 509–544 CrossRef CAS PubMed.
  2. M. Zimmer, Chem. Soc. Rev., 2009, 38, 2823–2832 RSC.
  3. A. A. Pakhomov and V. I. Martynov, Chem. Biol., 2008, 15, 755–764 CrossRef CAS PubMed.
  4. T. D. Craggs, Chem. Soc. Rev., 2009, 38, 2865–2875 RSC.
  5. N. C. Shaner, P. A. Steinbach and R. Y. Tsien, Nat. Methods, 2005, 2, 905–909 CrossRef CAS PubMed.
  6. E. A. Rodriguez, R. E. Campbell, J. Y. Lin, M. Z. Lin, A. Miyawaki, A. E. Palmer, X. Shu, J. Zhang and R. Y. Tsien, Trends Biochem. Sci., 2017, 42, 111–129 CrossRef CAS PubMed.
  7. T. Nagai, K. Ibata, E. S. Park, M. Kubota, K. Mikoshiba and A. Miyawaki, Nat. Biotechnol., 2002, 20, 87–90 CrossRef CAS PubMed.
  8. T. Sen, A. V. Mamontova, A. V. Titelmayer, A. M. Shakhov, A. A. Astafiev, A. Acharya, K. A. Lukyanov, A. I. Krylov and A. M. Bogdanov, Int. J. Mol. Sci., 2019, 20, 5229 CrossRef CAS PubMed.
  9. A. Rekas, J. R. Alattia, T. Nagai, A. Miyawaki and M. Ikura, J. Biol. Chem., 2002, 277, 50573–50578 CrossRef CAS PubMed.
  10. M. Zimmer, Chem. Rev., 2002, 102, 759–781 CrossRef CAS PubMed.
  11. R. M. Wachter, Acc. Chem. Res., 2007, 40, 120–127 CrossRef CAS PubMed.
  12. L. Zhang, H. N. Patel, J. W. Lappe and R. M. Wachter, J. Am. Chem. Soc., 2006, 128, 4766–4772 CrossRef CAS PubMed.
  13. B. G. Reid and G. C. Flynn, Biochemistry, 1997, 36, 6786–6791 CrossRef CAS PubMed.
  14. N. V. Pletneva, V. Z. Pletnev, K. A. Lukyanov, N. G. Gurskaya, E. A. Goryacheva, V. I. Martynov, A. Wlodawer, Z. Dauter and S. Pletnev, J. Biol. Chem., 2010, 285, 15978–15984 CrossRef CAS PubMed.
  15. M. A. Rosenow, H. A. Huffman, M. E. Phail and R. M. Wachter, Biochemistry, 2004, 43, 4464–4472 CrossRef CAS PubMed.
  16. T. I. Wood, D. P. Barondeau, C. Hitomi, C. J. Kassmann, J. A. Tainer and E. D. Getzoff, Biochemistry, 2005, 44, 16211–16220 CrossRef CAS PubMed.
  17. D. P. Barondeau, C. J. Kassmann, J. A. Tainer and E. D. Getzoff, Biochemistry, 2005, 44, 1960–1970 CrossRef CAS PubMed.
  18. D. P. Barondeau, J. A. Tainer and E. D. Getzoff, J. Am. Chem. Soc., 2006, 128, 3166–3168 CrossRef CAS PubMed.
  19. M. A. Rosenow, H. N. Patel and R. M. Wachter, Biochemistry, 2005, 44, 8303–8311 CrossRef CAS PubMed.
  20. J. A. Sniegowski, M. E. Phail and R. M. Wachter, Biochem. Biophys. Res. Commun., 2005, 332, 657–663 CrossRef CAS PubMed.
  21. D. M. Shcherbakova and V. V. Verkhusha, Curr. Opin. Chem. Biol., 2014, 20, 60–68 CrossRef CAS PubMed.
  22. V. Adam, R. Berardozzi, M. Byrdin and D. Bourgeois, Curr. Opin. Chem. Biol., 2014, 20, 92–102 CrossRef CAS PubMed.
  23. J. Lippincott-Schwartz and G. H. Patterson, Trends Cell Biol., 2009, 19, 555–565 CrossRef CAS PubMed.
  24. A. Acharya, A. M. Bogdanov, B. L. Grigorenko, K. B. Bravaya, A. V. Nemukhin, K. A. Lukyanov and A. I. Krylov, Chem. Rev., 2017, 117, 758–795 CrossRef CAS PubMed.
  25. J. N. Henderson, R. Gepshtein, J. R. Heenan, K. Kallio, D. Huppert and S. J. Remington, J. Am. Chem. Soc., 2009, 131, 4176–4177 CrossRef CAS PubMed.
  26. R. Ando, H. Hama, M. Yamamoto-Hino, H. Mizuno and A. Miyawaki, Proc. Natl. Acad. Sci. U. S. A., 2002, 99, 12651–12656 CrossRef CAS PubMed.
  27. T. Brakemann, A. C. Stiel, G. Weber, M. Andresen, I. Testa, T. Grotjohann, M. Leutenegger, U. Plessmann, H. Urlaub, C. Eggeling, M. C. Wahl, S. W. Hell and S. Jakobs, Nat. Biotechnol., 2011, 29, 942–947 CrossRef CAS PubMed.
  28. R. Ando, H. Mizuno and A. Miyawaki, Science, 2004, 306, 1370–1373 CrossRef CAS PubMed.
  29. T. Grotjohann, I. Testa, M. Leutenegger, H. Bock, N. T. Urban, F. Lavoie-Cardinal, K. I. Willig, C. Eggeling, S. Jakobs and S. W. Hell, Nature, 2011, 478, 204–208 CrossRef CAS PubMed.
  30. W. Niu and J. Guo, Mol. BioSyst., 2013, 9, 2961–2970 RSC.
  31. S. C. Reddington, A. J. Baldwin, R. Thompson, A. Brancale, E. M. Tippmann and D. D. Jones, Chem. Sci., 2015, 6, 1159–1166 RSC.
  32. S. C. Reddington, P. J. Rizkallah, P. D. Watson, R. Pearson, E. M. Tippmann and D. D. Jones, Angew. Chem., Int. Ed., 2013, 52, 5974–5977 CrossRef CAS PubMed.
  33. S. C. Reddington, S. Driezis, A. M. Hartley, P. D. Watson, P. J. Rizkallah and D. D. Jones, RSC Adv., 2015, 5, 77734–77738 RSC.
  34. J. Chin, S. Santoro, A. Martin, D. King, L. Wang and P. Schultz, J. Am. Chem. Soc., 2002, 124, 9026–9027 CrossRef CAS PubMed.
  35. S. Reddington, P. Watson, P. Rizkallah, E. Tippmann and D. D. Jones, Biochem. Soc. Trans., 2013, 41, 1177–1182 CrossRef CAS PubMed.
  36. J. Wei, J. S. Gibbs, H. D. Hickman, S. S. Cush, J. R. Bennink and J. W. Yewdell, J. Biol. Chem., 2015, 290, 16431–16439 CrossRef CAS PubMed.
  37. J. J. van Thor, T. Gensch, K. J. Hellingwerf and L. N. Johnson, Nat. Struct. Biol., 2002, 9, 37–41 CrossRef CAS PubMed.
  38. B. L. Grigorenko, A. I. Krylov and A. V. Nemukhin, J. Am. Chem. Soc., 2017, 139, 10239–10249 CrossRef CAS PubMed.
  39. G. D. P. Gomes and I. V. Alabugin, in Applied Theoretical Organic Chemistry, ed. D. J. Tantillo, 2018, pp. 451–502,  DOI:10.1142/q0119.
  40. F. Weinhold, C. R. Landis and E. D. Glendening, Int. Rev. Phys. Chem., 2016, 35, 399–440 Search PubMed.
  41. G. D. P. Gomes, V. Vil, A. Terent'ev and I. V. Alabugin, Chem. Sci., 2015, 6, 6783–6791 RSC.
  42. E. Juaristi, G. Dos Passos Gomes, A. O. Terent'ev, R. Notario and I. V. Alabugin, J. Am. Chem. Soc., 2017, 139, 10799–10813 CrossRef CAS PubMed.
  43. V. A. Vil, Y. A. Barsegyan, L. Kuhn, M. V. Ekimova, E. A. Semenov, A. A. Korlyukov, A. O. Terent'ev and I. V. Alabugin, Chem. Sci., 2020, 11, 5313–5322 RSC.
  44. E. Chovancova, A. Pavelka, P. Benes, O. Strnad, J. Brezovsky, B. Kozlikova, A. Gora, V. Sustr, M. Klvana, P. Medek, L. Biedermannova, J. Sochor and J. Damborsky, PLoS Comput. Biol., 2012, 8, e1002708 CrossRef CAS PubMed.
  45. D. M. Shcherbakova, P. Sengupta, J. Lippincott-Schwartz and V. V. Verkhusha, Annu. Rev. Biophys., 2014, 43, 303–329 CrossRef CAS PubMed.
  46. J. A. Arpino, A. J. Baldwin, A. R. McGarrity, E. M. Tippmann and D. D. Jones, PLoS One, 2015, 10, e0127504 CrossRef PubMed.
  47. R. Heim, D. C. Prasher and R. Y. Tsien, Proc. Natl. Acad. Sci. U. S. A., 1994, 91, 12501–12504 CrossRef CAS PubMed.
  48. J. A. Sniegowski, J. W. Lappe, H. N. Patel, H. A. Huffman and R. M. Wachter, J. Biol. Chem., 2005, 280, 26248–26255 CrossRef CAS PubMed.
  49. M. H. Seifert, D. Ksiazek, M. K. Azim, P. Smialowski, N. Budisa and T. A. Holak, J. Am. Chem. Soc., 2002, 124, 7932–7942 CrossRef CAS PubMed.
  50. J. A. Arpino, P. J. Rizkallah and D. D. Jones, Acta Crystallogr., Sect. D: Biol. Crystallogr., 2014, 70, 2152–2162 CrossRef CAS PubMed.
  51. K. Brejc, T. K. Sixma, P. A. Kitts, S. R. Kain, R. Y. Tsien, M. Ormo and S. J. Remington, Proc. Natl. Acad. Sci. U. S. A., 1997, 94, 2306–2311 CrossRef CAS PubMed.
  52. S. T. Hsu, G. Blaser and S. E. Jackson, Chem. Soc. Rev., 2009, 38, 2951–2965 RSC.
  53. A. M. Hartley, A. J. Zaki, A. R. McGarrity, C. Robert-Ansart, A. V. Moskalenko, G. F. Jones, M. F. Craciun, S. Russo, M. Elliott, J. E. Macdonald and D. D. Jones, Chem. Sci., 2015, 6, 3712–3717 RSC.
  54. A. M. Hartley, H. L. Worthy, S. C. Reddington, P. J. Rizkallah and D. D. Jones, Chem. Sci., 2016, 7, 6484–6491 RSC.
  55. H. L. Worthy, H. S. Auhim, W. D. Jamieson, J. R. Pope, A. Wall, R. Batchelor, R. L. Johnson, D. W. Watkins, P. Rizkallah, O. K. Castell and D. D. Jones, Commun. Chem., 2019, 2, 83 CrossRef.
  56. T. Sen, Y. Ma, I. V. Polyakov, B. L. Grigorenko, A. V. Nemukhin and A. I. Krylov, J. Phys. Chem. B, 2021, 125, 757–770 CrossRef CAS PubMed.
  57. J. Wiedenmann, S. Ivanchenko, F. Oswald, F. Schmitt, C. Rocker, A. Salih, K. D. Spindler and G. U. Nienhaus, Proc. Natl. Acad. Sci. U. S. A., 2004, 101, 15905–15910 CrossRef CAS PubMed.
  58. G. Winter, J. Appl. Crystallogr., 2009, 43, 186–190 CrossRef.
  59. P. Evans, Acta Crystallogr., Sect. D: Biol. Crystallogr., 2006, 62, 72–82 CrossRef PubMed.
  60. S. Bailey, Acta Crystallogr., Sect. D: Biol. Crystallogr., 1994, 50, 760–763 CrossRef PubMed.
  61. A. J. McCoy, R. W. Grosse-Kunstleve, P. D. Adams, M. D. Winn, L. C. Storoni and R. J. Read, J. Appl. Crystallogr., 2007, 40, 658–674 CrossRef CAS PubMed.
  62. G. N. Murshudov, P. Skubak, A. A. Lebedev, N. S. Pannu, R. A. Steiner, R. A. Nicholls, M. D. Winn, F. Long and A. A. Vagin, Acta Crystallogr., Sect. D: Biol. Crystallogr., 2011, 67, 355–367 CrossRef CAS PubMed.
  63. A. A. Lebedev, P. Young, M. N. Isupov, O. V. Moroz, A. A. Vagin and G. N. Murshudov, Acta Crystallogr., Sect. D: Biol. Crystallogr., 2012, 68, 431–440 CrossRef CAS PubMed.
  64. P. Emsley and K. Cowtan, Acta Crystallogr., Sect. D: Biol. Crystallogr., 2004, 60, 2126–2132 CrossRef PubMed.
  65. B. L. Grigorenko, E. D. Kots, A. I. Krylov and A. V. Nemukhin, Mendeleev Commun., 2019, 29, 187–189 CrossRef CAS.
  66. B. L. Grigorenko, I. V. Polyakov, A. I. Krylov and A. V. Nemukhin, J. Phys. Chem. B, 2019, 123, 8901–8909 CrossRef CAS PubMed.
  67. C. Adamo and V. Barone, J. Chem. Phys., 1999, 110, 6158–6170 CrossRef CAS.
  68. E. Apra, E. J. Bylaska, W. A. de Jong, N. Govind, K. Kowalski, T. P. Straatsma, M. Valiev, H. J. J. van Dam, Y. Alexeev, J. Anchell, V. Anisimov, F. W. Aquino, R. Atta-Fynn, J. Autschbach, N. P. Bauman, J. C. Becca, D. E. Bernholdt, K. Bhaskaran-Nair, S. Bogatko, P. Borowski, J. Boschen, J. Brabec, A. Bruner, E. Cauet, Y. Chen, G. N. Chuev, C. J. Cramer, J. Daily, M. J. O. Deegan, T. H. Dunning Jr, M. Dupuis, K. G. Dyall, G. I. Fann, S. A. Fischer, A. Fonari, H. Fruchtl, L. Gagliardi, J. Garza, N. Gawande, S. Ghosh, K. Glaesemann, A. W. Gotz, J. Hammond, V. Helms, E. D. Hermes, K. Hirao, S. Hirata, M. Jacquelin, L. Jensen, B. G. Johnson, H. Jonsson, R. A. Kendall, M. Klemm, R. Kobayashi, V. Konkov, S. Krishnamoorthy, M. Krishnan, Z. Lin, R. D. Lins, R. J. Littlefield, A. J. Logsdail, K. Lopata, W. Ma, A. V. Marenich, J. Martin Del Campo, D. Mejia-Rodriguez, J. E. Moore, J. M. Mullin, T. Nakajima, D. R. Nascimento, J. A. Nichols, P. J. Nichols, J. Nieplocha, A. Otero-de-la-Roza, B. Palmer, A. Panyala, T. Pirojsirikul, B. Peng, R. Peverati, J. Pittner, L. Pollack, R. M. Richard, P. Sadayappan, G. C. Schatz, W. A. Shelton, D. W. Silverstein, D. M. A. Smith, T. A. Soares, D. Song, M. Swart, H. L. Taylor, G. S. Thomas, V. Tipparaju, D. G. Truhlar, K. Tsemekhman, T. Van Voorhis, A. Vazquez-Mayagoitia, P. Verma, O. Villa, A. Vishnu, K. D. Vogiatzis, D. Wang, J. H. Weare, M. J. Williamson, T. L. Windus, K. Wolinski, A. T. Wong, Q. Wu, C. Yang, Q. Yu, M. Zacharias, Z. Zhang, Y. Zhao and R. J. Harrison, J. Chem. Phys., 2020, 152, 184102 CrossRef CAS PubMed.
  69. A. A. Granovsky, J. Chem. Phys., 2011, 134, 214113 CrossRef PubMed.
  70. A. V. Nemukhin, B. L. Grigorenko, M. G. Khrenova and A. I. Krylov, J. Phys. Chem. B, 2019, 123, 6133–6149 CrossRef CAS PubMed.
  71. A. A. Granovsky, Firefly version 8, http://classic.chem.msu.su/gran/firefly/index.html Search PubMed.
  72. A. V. Marenich, C. J. Cramer and D. G. Truhlar, J. Phys. Chem. B, 2009, 113, 6378–6396 CrossRef CAS PubMed.
  73. J. D. Chai and M. Head-Gordon, Phys. Chem. Chem. Phys., 2008, 10, 6615–6620 RSC.
  74. S. Grimme, J. Antony, S. Ehrlich and H. Krieg, J. Chem. Phys., 2010, 132, 154104 CrossRef PubMed.
  75. S. Grimme, S. Ehrlich and L. Goerigk, J. Comput. Chem., 2011, 32, 1456–1465 CrossRef CAS PubMed.

Footnotes

Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sc06693a
These authors contributed equally to the work.

This journal is © The Royal Society of Chemistry 2021
Click here to see how this site uses Cookies. View our privacy policy here.