Potential Energy Landscape of the Photoinduced Multiple Proton-Transfer Process in the Green Fluorescent Protein: Classical Molecular Dynamics and Multiconfigurational Electronic Structure Calculations

Development of an Iron(II) Complex Exhibiting Thermal- and Photoinduced Double Proton-Transfer-Coupled Spin Transition in a Short Hydrogen Bond

Multiple proton transfer (PT) controllable by external stimuli plays a crucial role in fundamental chemistry, biological activity, and material science. However, in crystalline systems, controlling multiple PT, which results in a distinct protonation state, remains challenging. In this study, we developed a novel tridentate ligand and iron(II) complex with a short hydrogen bond (HB) that exhibits a PT-coupled spin transition (PCST). Single-crystal X-ray and neutron diffraction measurements revealed that the positions of the two protons in the complex can be controlled by temperature and photoirradiation based on the thermal- and photoinduced PCST. The obtained results suggest that designing molecules that form short HBs is a promising approach for developing multiple PT systems in crystals.

Excited State Structural Evolution of a GFP Single-Site Mutant Tracked by Tunable Femtosecond-Stimulated Raman Spectroscopy

Tracking vibrational motions during a photochemical or photophysical process has gained momentum, due to its sensitivity to the progression of reaction and change of environment. In this work, we implemented an advanced ultrafast vibrational technique, femtosecond-stimulated Raman spectroscopy (FSRS), to monitor the excited state structural evolution of an engineered green fluorescent protein (GFP) single-site mutant S205V. This mutation alters the original excited state proton transfer (ESPT) chain. By strategically tuning the Raman pump to different wavelengths (i.e., 801, 539, and 504 nm) to achieve pre-resonance with transient excited state electronic bands, the characteristic Raman modes of the excited protonated (A*) chromophore species and intermediate deprotonated (I*) species can be selectively monitored. The inhomogeneous distribution/population of A* species go through ESPT with a similar ~300 ps time constant, confirming that bridging a water molecule to protein residue T203 in the ESPT chain is the rate-limiting step. Some A* species undergo vibrational cooling through high-frequency motions on the ~190 ps time scale. At early times, a portion of the largely protonated A* species could also undergo vibrational cooling or return to the ground state with a ~80 ps time constant. On the photoproduct side, a ~1330 cm⁻¹ delocalized motion is observed, with dispersive line shapes in both the Stokes and anti-Stokes FSRS with a pre-resonance Raman pump, which indicates strong vibronic coupling, as the mode could facilitate the I* species to reach a relatively stable state (e.g., the main fluorescent state) after conversion from A*. Our findings disentangle the contributions of various vibrational motions active during the ESPT reaction, and offer new structural dynamics insights into the fluorescence mechanisms of engineered GFPs and other analogous autofluorescent proteins.

Mutagenic induction of an ultra-fast water-chain proton wire

Replacement of the hydroxyl group of a hydrophilic sidechain by an H atom in the proton wire of GFP induces formation of a water-chain proton wire. Surprisingly, this "non-native" water chain functions as a proton wire with response times within 10 ps of the wild type protein. This remarkable rate retention is understood as a natural consequence of the well-known Grotthuss mechanism of proton transfer in water.

A new twist in the photophysics of the GFP chromophore: a volume-conserving molecular torsion couple

Dynamics of a nonplanar GFP chromophore are studied experimentally and theoretically. Coupled torsional motion is responsible for the ultrafast decay.

The simple structure of the chromophore of the green fluorescent protein (GFP), a phenol and an imidazolone ring linked by a methyne bridge, supports an exceptionally diverse range of excited state phenomena. Here we describe experimentally and theoretically the photochemistry of a novel sterically crowded nonplanar derivative of the GFP chromophore. It undergoes an excited state isomerization reaction accompanied by an exceptionally fast (sub 100 fs) excited state decay. The decay dynamics are essentially independent of solvent polarity and viscosity. Excited state structural dynamics are probed by high level quantum chemical calculations revealing that the fast decay is due to a conical intersection characterized by a twist of the rings and pyramidalization of the methyne bridge carbon. The intersection can be accessed without a barrier from the pre-twisted Franck–Condon structure, and the lack of viscosity dependence is due to the fact that the rings twist in the same direction, giving rise to a volume-conserving decay coordinate. Moreover, the rotation of the phenyl, methyl and imidazolone groups is coupled in the sterically crowded structure, with the methyl group translating the rotation of one ring to the next. As a consequence, the excited state dynamics can be viewed as a torsional couple, where the absorbed photon energy leads to conversion of the out-of-plane orientation from one ring to the other in a volume conserving fashion. A similar modification of the range of methyne dyes may provide a new family of devices for molecular machines, specifically torsional couples.

The mechanism of oxidation in chromophore maturation of wild-type green fluorescent protein: a theoretical study

DFT calculations suggested that the thermodynamically unfavourable cyclized product was trapped by oxidation.

Photoinduced Chemistry in Fluorescent Proteins: Curse or Blessing?

Photoinduced reactions play an important role in the photocycle of fluorescent proteins from the green fluorescent protein (GFP) family. Among such processes are photoisomerization, photooxidation/photoreduction, breaking and making of covalent bonds, and excited-state proton transfer (ESPT). Many of these transformations are initiated by electron transfer (ET). The quantum yields of these processes vary significantly, from nearly 1 for ESPT to 10^-4-10^-6 for ET. Importantly, even when quantum yields are relatively small, at the conditions of repeated illumination the overall effect is significant. Depending on the task at hand, fluorescent protein photochemistry is regarded either as an asset facilitating new applications or as a nuisance leading to the loss of optical output. The phenomena arising due to phototransformations include (i) large Stokes shifts, (ii) photoconversions, photoactivation, and photoswitching, (iii) phototoxicity, (iv) blinking, (v) permanent bleaching, and (vi) formation of long-lived intermediates. The focus of this review is on the most recent experimental and theoretical work on photoinduced transformations in fluorescent proteins. We also provide an overview of the photophysics of fluorescent proteins, highlighting the interplay between photochemistry and other channels (fluorescence, radiationless relaxation, and intersystem crossing). The similarities and differences with photochemical processes in other biological systems and in dyes are also discussed.

Multiscale Simulations on Spectral Tuning and the Photoisomerization Mechanism in Fluorescent RNA Spinach

Fluorescent RNA aptamer Spinach can bind and activate a green fluorescent protein (GFP)-like chromophore (an anionic DFHBDI chromophore) displaying green fluorescence. Spectroscopic properties, spectral tuning, and the photoisomerization mechanism in the Spinach-DFHBDI complex have been investigated by high-level QM and hybrid ONIOM(QM:AMBER) methods (QM method: (TD)DFT, SF-BHHLYP, SAC-CI, LT-DF-LCC2, CASSCF, or MS-CASPT2), as well as classical molecular dynamics (MD) simulations. First, our benchmark calculations have shown that TD-DFT and spin-flip (SF) TD-DFT (SF-BHHLYP) failed to give a satisfactory description of absorption and emission of the anionic DFHBDI chromophore. Comparatively, SAC-CI, LT-DF-LCC2, and MS-CASPT2 can give more reliable transition energies and are mainly used to further study the spectra of the anionic DFHBDI chromophore in Spinach. The RNA environmental effects on the spectral tuning and the photoisomerization mechanism have been elucidated. Our simulations show that interactions of the anionic cis-DFHBDI chromophore with two G-quadruplexes as well as a UAU base triple suppress photoisomerization of DFHBDI. In addition, strong hydrogen bonds between the anionic cis-DFHBDI chromophore and nearby nucleotides facilitate its binding to Spinach and further inhibit the cis-trans photoisomerization of DFHBDI. Solvent molecules, ions, and loss of key hydrogen bonds with nearby nucleotides could induce dissociation of the anionic trans-DFHBDI chromophore from the binding site. These results provide new insights into fluorescent RNA Spinach and may help rational design of other fluorescent RNAs.

Characterizing the Structures, Spectra, and Energy Landscapes Involved in the Excited-State Proton Transfer Process of Red Fluorescent Protein LSSmKate1

By applying molecular dynamics (MD) simulations and quantum chemical calculations, we have characterized the states and processes involved in the excited-state proton transfer (ESPT) of LSSmKate1. MD simulations identify two stable structures in the electronic ground state of LSSmKate1, one with a protonated chromophore and the other with a deprotonated chromophore, thus leading to two separate low-energy absorption maxima with a large energy spacing, as observed in the calculated and experimentally measured absorption spectra. Proton transfer is induced by electronic excitation. When LSSmKate1 is excited, the electrons in the chromophore are transferred from the phenol ring to the N-acylimine moiety; the acidity of a phenolic hydroxyl group is thus enhanced. The calculated potential energy curves (PECs) exhibit energetic feasibility in the generation of the fluorescent species in LSSmKate1, and the exact agreement between the calculated and experimentally measured values of the large Stokes shift further provides solid theoretical evidence for the ESPT process taking place in photoexcited LSSmKate1. The molecular environments play a significant role in the geometries and absorption/emission energies of the chromophores. Overall, TD-ωB97X-D/molecular mechanics (MM) provides a better description of the optical properties of LSSmKate1 than TD-B3LYP/MM, although it always overestimates the excitation energies.

Transient low-barrier hydrogen bond in the photoactive state of green fluorescent protein

In this paper, we have analyzed the feasibility of spontaneous proton transfer in GFP at the Franck-Condon region directly after photoexcitation. Computation of a sizeable portion of the potential energy surface at the Franck-Condon region of A the structure shows the process of proton transfer to be unfavorable by 3 kcal mol(-1) in S1 if no further structural relaxation is permitted. The ground vibrational state is found to lie above the potential energy barrier of the proton transfer in both S0 and S1. Expectation values of the geometry reveal that the proton shared between the chromophore and W22, and the proton shared between Ser205 and Glu222 are very close to the center of the respective hydrogen bonds, giving support to the claim that the first transient intermediate detected after photoexcitation (I0*) has characteristics similar to those of a low-barrier hydrogen bond [Di Donato et al., Phys. Chem. Chem. Phys., 2012, 13, 16295]. A quantum dynamical calculation of the evolution in the excited state shows an even larger probability of finding those two protons close to the center compared to in the ground state, but no formation of the proton-transferred product is observed. A QM/MM photoactive state geometry optimization, initiated using a configuration obtained by taking the A minimum and moving the protons to the product side, yields a minimum energy structure with the protons transferred and in which the His148 residue is substantially closer to the now anionic chromophore. These results indicate that: (1) proton transfer is not possible if structural relaxation of the surroundings of the chromophore is prevented; (2) protons H1 and H3 especially are found very close to the point halfway between the donor and acceptor after photoexcitation when the zero-point energy is considered; (3) a geometrical parameter exists (the His148-Cro distance) under which the structure with the protons transferred is not a minimum, and that, if included, should lead to the fluorescing I* structure. The existence of an oscillating stationary state between the reactants and products of the triple proton transfer reaction can explain the dual emission reported for the I0* intermediate of wtGFP.

Excited-State Proton-Transfer-Induced Trapping Enhances the Fluorescence Emission of a Locked GFP Chromophore

The chemical locking of the central single bond in core chromophores of green fluorescent proteins (GFPs) influences their excited-state behavior in a distinct manner. Experimentally, it significantly enhances the fluorescence quantum yield of GFP chromophores with an ortho-hydroxyl group, while it has almost no effect on the photophysics of GFP chromophores with a para-hydroxyl group. To unravel the underlying physical reasons for this different behavior, we report static electronic structure calculations and nonadiabatic dynamics simulations on excited-state intramolecular proton transfer, cis–trans isomerization, and excited-state deactivation in a locked ortho-substituted GFP model chromophore (o-LHBI). On the basis of our previous and present results, we find that the S₁ keto species is responsible for the fluorescence emission of the unlocked o-HBI and the locked o-LHBI species. Chemical locking does not change the parts of the S₁ and S₀ potential energy surfaces relevant to enol–keto tautomerization; hence, in both chromophores, there is an ultrafast excited-state intramolecular proton transfer that takes only 35 fs on average. However, the locking effectively hinders the S₁ keto species from approaching the keto S₁/S₀ conical intersections so that most of trajectories are trapped in the S₁ keto region for the entire 2 ps simulation time. Therefore, the fluorescence quantum yield of o-LHBI is enhanced compared with that of unlocked o-HBI, in which the S₁ excited-state decay is efficient and ultrafast. In the case of the para-substituted GFP model chromophores p-HBI and p-LHBI, chemical locking hardly affects their efficient excited-state deactivation via cis–trans isomerization; thus, the fluorescence quantum yields in these chromophores remain very low. The insights gained from the present work may help to guide the design of new GFP chromophores with improved fluorescence emission and brightness.

Emission shaping in fluorescent proteins: role of electrostatics and π-stacking

We obtained the fluorescence spectrum of the GFP with trajectory simulations, and revealed the role of the protein sidechains in emission shifts.

Theoretical Computer-Aided Mutagenic Study on the Triple Green Fluorescent Protein Mutant S65T/H148D/Y145F

Green fluorescent protein (GFP) mutant S65T/H148D has been proposed to host a photocycle that involves an excited-state proton transfer between the chromophore (Cro) and the Asp148 residue and takes place in less than 50 fs without a measurable kinetic isotope effect. It has been suggested that the interaction between the unsuspected Tyr145 residue and the chromophore is needed for the ultrafast sub-50 fs rise in fluorescence. To verify this, we have performed a computer-aided mutagenic study to introduce the additional mutation Y145F, which eliminates this interaction. By means of QM/MM molecular dynamics simulations and time-dependent density functional theory studies, we have assessed the importance of the Cro-Tyr145 interaction and the solvation of Asp148 and shown that in the triple mutant S65T/H148D/Y145F a significant loss in the ultrafast rise of the Stokes-shifted fluorescence should be expected.

Complete Proton Transfer Cycle in GFP and Its T203V and S205V Mutants

Proton transfer is critical in many important biochemical reactions. The unique three-step excited-state proton transfer in avGFP allows observations of protein proton transport in real-time. In this work we exploit femtosecond to microsecond transient IR spectroscopy to record, in D₂O, the complete proton transfer photocycle of avGFP, and two mutants (T203V and S205V) which modify the structure of the proton wire. Striking differences and similarities are observed among the three mutants yielding novel information on proton transfer mechanism, rates, isotope effects, H-bond strength and proton wire stability. These data provide a detailed picture of the dynamics of long-range proton transfer in a protein against which calculations may be compared.

Insight into the structure and the mechanism of the slow proton transfer in the GFP double mutant T203V/S205A

Mutations near the fluorescing chromophore of the green fluorescent protein (GFP) have direct effects on the absorption and emission spectra. Some mutants have significant band shifts and most of the mutants exhibit a loss of fluorescence intensity. In this study we continue our investigation of the factors controlling the excited state proton transfer (PT) process of GFP, in particular to study the effects of modifications to the key side chain Ser205 in wt-GFP, proposed to participate in the proton wire. To this aim we combined mutagenesis, X-ray crystallography, steady-state spectroscopy, time-resolved emission spectroscopy and all-atom explicit molecular dynamics (MD) simulations to study the double mutant T203V/S205A. Our results show that while in the previously described GFP double mutant T203V/S205V the PT process does not occur, in the T203V/S205A mutant the PT process does occur, but with a 350 times slower rate than in wild-type GFP (wt-GFP). Furthermore, the kinetic isotope effect in the GFP double mutant T203V/S205A is twice smaller than in the wt-GFP and in the GFP single mutant S205V, which forms a novel PT pathway. On the other hand, the crystal structure of GFP T203V/S205A does not reveal a viable proton transfer pathway. To explain PT in GFP T203V/S205A, we argue on the basis of the MD simulations for an alternative, novel proton-wire pathway which involves the phenol group of the chromophore and water molecules infrequently entering from the bulk. This alternative pathway may explain the dramatically slow PT in the GFP double mutant T203V/S205A compared to wt-GFP.

Electric field effect on the ground state proton transfer in the H-bonded HBDI complex: an implication of the green fluorescent protein

In this paper, first-principles calculations were performed regarding the electric field effect on the ground state proton transfer (GSPT) in the H-bonded p-hydroxybenzylideneimidazolidinone (HBDI) network that represents the active site of the green fluorescent protein (GFP).

Human cellular retinaldehyde-binding protein has secondary thermal 9-cis-retinal isomerase activity

Cellular retinaldehyde-binding protein (CRALBP) chaperones 11-cis-retinal to convert opsin receptor molecules into photosensitive retinoid pigments of the eye. We report a thermal secondary isomerase activity of CRALBP when bound to 9-cis-retinal. UV/vis and (1)H NMR spectroscopy were used to characterize the product as 9,13-dicis-retinal. The X-ray structure of the CRALBP mutant R234W:9-cis-retinal complex at 1.9 Å resolution revealed a niche in the binding pocket for 9-cis-aldehyde different from that reported for 11-cis-retinal. Combined computational, kinetic, and structural data lead us to propose an isomerization mechanism catalyzed by a network of buried waters. Our findings highlight a specific role of water molecules in both CRALBP-assisted specificity toward 9-cis-retinal and its thermal isomerase activity yielding 9,13-dicis-retinal. Kinetic data from two point mutants of CRALBP support an essential role of Glu202 as the initial proton donor in this isomerization reaction.

Reaction mechanism of photoinduced decarboxylation of the photoactivatable green fluorescent protein: an ONIOM(QM:MM) study

Photoactivatable (PA) fluorescent proteins are a new class of fluorescent proteins in which the intensity of fluorescence is dramatically enhanced through photoinduced decarboxylation process. In the present study, the reaction mechanism of the photoinduced decarboxylation in PA-GFP was investigated by the ONIOM(QM:MM) method. The decarboxylation process starts from the first excited state (IntraCT state) and then proceeds along an InterCT state after the first crossing (or an approximate transition state). Relative to an equilibrium structure in S(0), a barrier of ~94 kcal/mol to reach this approximate transition state is the rate-determining step for the entire decarboxylation process. The InterCT state becomes the open-shell ground state in the product, after the subsequent crossing with a closed-shell state that holds an extra electron on the dissociated CO(2). The present study elucidated for the first time the mechanism of the photoinduced decarboxylation of PA-GFP and supports the widely accepted Kolbe pathway, which could be a common mechanism for the irreversible photoinduced decarboxylation in different fluorescent proteins.

Isomerization mechanism of the HcRed fluorescent protein chromophore

To understand how the protein achieves fluorescence, the isomerization mechanism of the HcRed chromophore is studied both under vacuum and in the solvated red fluorescent protein. Quantum mechanical (QM) and quantum mechanical/molecular mechanical (QM/MM) methods are applied both for the ground and the first excited state. The photoinduced processes in the chromophore mainly involve torsions around the imidazolinone-bridge bond (τ) and the phenoxy-bridge bond (φ). Under vacuum, the isomerization of the cis-trans chromophore essentially proceeds by τ twisting, while the radiationless decay requires φ torsion. By contrast, the isomerization of the cis-trans chromophore in HcRed occurs via simultaneous τ and φ twisting. The protein environment significantly reduces the barrier of this hula twist motion compared with vacuum. The excited-state isomerization barrier via the φ rotation of the cis-coplanar conformer in HcRed is computed to be significantly higher than that of the trans-non-coplanar conformer. This is consistent with the experimental observation that the cis-coplanar-conformation of the chromophore is related to the fluorescent properties of HcRed, while the trans-non-planar conformation is weakly fluorescent or non-fluorescent. Our study shows how the protein modifies the isomerization mechanism, notably by interactions involving the nearby residue Ile197, which keeps the chromophore coplanar and blocks the twisting motion that leads to photoinduced radiationless decay.

Theoretical studies of chromophore maturation in the wild-type green fluorescent protein: ONIOM(DFT:MM) investigation of the mechanism of cyclization

The availability of a gene encoding green fluorescence immediately stimulates interest in the puzzle of autocatalytic formation of the green fluorescent protein (GFP) chromophore. Numerous experimental and theoretical studies have indicated that cyclization is the first and most important step in the maturation process of the GFP. In our previous paper based on cluster models [J. Phys. Chem. B2010, 114, 9698-9705], two possible mechanisms have been investigated with the conclusion that the backbone condensation initiated by deprotonation of the Gly67 amide nitrogen is easier than deprotonation of the Tyr66 α-carbon. However, the impact of the protein environment on the reaction mechanism remains to be explored. In this paper, we investigated the two possible mechanisms with inclusion of protein environmental effects by using molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) calculations. Our calculations reveal no hydrogen bonding network that would facilitate deprotonation of the amide nitrogen of Gly67, although it is the lower energy pathway in the cluster model system. Contrastingly, there is a hydrogen bonding network between Tyr66 α-carbon and Glu222, which is in good agreement with X-ray data. The ONIOM studies show that proton transfer from Tyr66 α-carbon to Glu222 is a long-distance charge transfer process. The charge distribution of the MM region has a significant perturbation to the wave function for the QM region, with the QM energy for the proton transfer product being increased under the influence of the electrostatic protein environment. The barrier for the rate-limiting step in cyclization is quite high, about 40.0 kcal/mol in the case of ONIOM-EE.

The mechanism of cyclization in chromophore maturation of green fluorescent protein: a theoretical study

An intriguing aspect of the green fluorescent protein (GFP) is the autocatalytic post-translational modification that results in the formation of its chromophore. Numerous experimental and theoretical studies indicate that cyclization is the first and the most important step in the maturation process. In this work, two proposed mechanisms for the cyclization were investigated by using the hybrid density functional theory method B3LYP. Cluster models corresponding to the two mechanisms proposed by Wachter et al. [J. Biol. Chem. 2005, 280, 26248-26255] are constructed on the basis of the X-ray crystal structure (PDB entry 2AWJ) and corresponding reaction path potential energy profiles for the two cyclization mechanisms are presented. Our results suggest that the backbone condensation initiated by deprotonation of the Gly67 amide nitrogen is easier than deprotonation of the Tyr66 alpha-carbon. Moreover, Arg96 fulfills the role of stabilizing the enolate moiety, and Glu222 plays the role of a general base. The formation of the cyclized product is found to be 16.0 and 18.6 kcal/mol endothermic with respect to the two models, which is in agreement with experimental observation.

QM/MM studies of structural and energetic properties of the far-red fluorescent protein HcRed

The far-red fluorescent protein HcRed was investigated using molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) calculations. Three models of HcRed (anionic chromophore) were considered, differing in the protonation states of nearby Glu residues (A: Glu214 and Glu146 both protonated; B: Glu214 protonated and Glu146 deprotonated; C: Glu214 and Glu146 both deprotonated). SCC-DFTB/MM MD simulations of model B yield good agreement with the available crystallographic data at ambient pH. Bond lengths in the QM region are well reproduced, with a root mean square (rms) deviation between experimental and average MD data of 0.079 A; the chromophore is almost co-planar, which is consistent with experimental observation; and the five hydrogen bonds involving the chromophore are conserved. QM/MM geometry optimizations were performed on representative snapshot structures from the MD simulations for each model. They confirm the structural features observed in the MD simulations. According to the DFT(B3LYP)/MM results, the cis-conformation of the chromophore is more stable than the trans-form by 9.1-12.9 kcal mol(-1) in model B, and by 12.4-19.9 kcal mol(-1) in model C, consistent with the experimental preference for the cis-isomer. However, in model A when both Glu214 and Glu146 are protonated, the stability is inverted with the trans-form being favored. The different protonation states of the titratable active-site residues Glu214 and Glu146 thus critically influence the manner in which the relative stability and degree of planarity of the cis- and trans-conformers vary with pH. Coupled with the known correlation of chromophore conformation with fluorescence efficiency, this work provides a detailed structural basis for the observed phenomenon that red fluorescent proteins such as HcRed, mKate and Rtms5 show bright fluorescence at high pH.

Multiscale modeling of proteins

The activity within a living cell is based on a complex network of interactions among biomolecules, exchanging information and energy through biochemical processes. These events occur on different scales, from the nano- to the macroscale, spanning about 10 orders of magnitude in the space domain and 15 orders of magnitude in the time domain. Consequently, many different modeling techniques, each proper for a particular time or space scale, are commonly used. In addition, a single process often spans more than a single time or space scale. Thus, the necessity arises for combining the modeling techniques in multiscale approaches. In this Account, I first review the different modeling methods for bio-systems, from quantum mechanics to the coarse-grained and continuum-like descriptions, passing through the atomistic force field simulations. Special attention is devoted to their combination in different possible multiscale approaches and to the questions and problems related to their coherent matching in the space and time domains. These aspects are often considered secondary, but in fact, they have primary relevance when the aim is the coherent and complete description of bioprocesses. Subsequently, applications are illustrated by means of two paradigmatic examples: (i) the green fluorescent protein (GFP) family and (ii) the proteins involved in the human immunodeficiency virus (HIV) replication cycle. The GFPs are currently one of the most frequently used markers for monitoring protein trafficking within living cells; nanobiotechnology and cell biology strongly rely on their use in fluorescence microscopy techniques. A detailed knowledge of the actions of the virus-specific enzymes of HIV (specifically HIV protease and integrase) is necessary to study novel therapeutic strategies against this disease. Thus, the insight accumulated over years of intense study is an excellent framework for this Account. The foremost relevance of these two biomolecular systems was recently confirmed by the assignment of two of the Nobel prizes in 2008: in chemistry for the discovery of GFP and in medicine for the discovery of HIV. Accordingly, these proteins were studied with essentially all of the available modeling techniques, making them ideal examples for studying the details of multiscale approaches in protein modeling.

Quantum-chemical calculations of a long proton wire. Application of a harmonic model to analysis of the structure of an ionic defect in a water chain with an excess proton

Quantum-chemical calculations of molecular complexes (NH(3))(3)Zn(2+)...(H(2)O)(n)...NH(3) (C(n), n = 11, 16, 21, and 30) simulating a proton wire donor-water chain-acceptor were carried out. Earlier found periodicity in the length of the O-H bonds in water chain is explained within the framework of a one-component harmonic model. In complexes C(n), the geometry and electronic structure of ionic defect in water chain with an excess proton were studied. Calculations carried out at ab initio (B3LYP/6-31+G**) and semiempirical (PM3) levels of theory predict different patterns of distribution of the O-H bonds lengths and positive charge on the H-bond hydrogen atoms in the region of ionic defect. The obtained data show how a length of water chain and position of a protonated water link in the chain influence the ionic defect structure. To describe the observed structures of ionic defect, the harmonic model was used and the role of parameters of the H-bonded chain was investigated. The performed analysis explains different mechanisms (concerted and stepwise) of proton transfer along the H-bonded chain derived from ab initio and semiempirical calculation schemes.

Structural and Relaxation Effects in Proton Wire Energetics: Model Studies of the Green Fluorescent Protein Photocycle

We present the results of a systematic series of constrained minimum energy pathway calculations on ground state potential energy surfaces, for a cluster model of the proton chain transfer that mediates the photocycle of the green fluorescent protein, as well as for a model including the solvated protein environment. The calculations vary in terms of the types of modes that are assumed to be capable of relaxing in concert with the movement of the protons and the results demonstrate that the nature and extent of dynamical relaxation has a substantive impact on the activation energy for the proton transfer. We discuss the implications of this in terms of currently available dynamical models and chemical rate theories that might be brought to bear on the kinetics of this important example of proton chain transfer in a biological system.