Computational study of the absorption spectra of green fluorescent protein mutants

In this work, we present a theoretical study of the relationship between molecular structure and the red-shift in absorption spectra of S65G and S65T green fluorescent protein (GFP) mutants. To identify the effects of the protein environment, we combined results from molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics calculations to obtain structural properties, and applied time-dependent density functional theory to calculate the excitation energies. By using results from the MD simulations, we were able to provide a systematic analysis of the structural details that may effect the red-shift in the absorption spectra when taking into account temperature effects. Furthermore, a detailed study of hydrogen bonding during the MD simulations demonstrated differences between S65G and S65T, for example, regarding hydrogen bonding with Glu222. An analysis of the absorption spectra for different forms of the chromophore emphasized the dominance of the anionic forms in solution for the S65G and S65T GFP mutants.     

Backbone dynamics of green fluorescent protein and the effect of histidine 148 substitution

Green fluorescent protein (GFP) and its mutants have become valuable tools in molecular biology. GFP has been regarded as a very stable and rigid protein with the beta-barrel shielding the chromophore from the solvent. Here, we report the 15N nuclear magnetic resonance (NMR) studies on the green fluorescent protein (GFPuv) and its mutant His148Gly. 15N NMR relaxation studies of GFPuv show that most of the beta-barrel of GFP is rigid on the picosecond to nanosecond time scale. For several regions, including the first alpha-helix and beta-sheets 3, 7, 8, and 10, increased hydrogen-deuterium exchange rates suggest a substantial conformational flexibility on the microsecond to millisecond time scales. Mutation of residue 148 located in beta-sheet 7 is known to have a strong impact on the fluorescence properties of GFPs. UV absorption and fluorescence spectra in combination with 1H-15N NMR spectra indicate that the His148Gly mutation not only reduces the absorption of the anionic chromophore state but also affects the conformational stability, leading to the appearance of doubled backbone amide resonances for a number of residues. This suggests the presence of two conformations in slow exchange on the NMR time scale in this mutant.     

Physicochemical properties of blue fluorescent protein determined via molecular dynamics simulation

Blue fluorescent protein (BFP) is a mutant of green fluorescent protein (GFP), where the chromophore has been modified to shift the emitted fluorescence into the blue spectral region. In this study, MD calculations were performed with the GROMACS simulation package and AMBER force field to investigate the dependence of BFPs physicochemical properties on temperature and applied pressure. The MD approach enabled us to calculate the compressibility of protein itself, separately from the nontrivial contribution of the hydration shell, which is difficult to achieve experimentally. The computed compressibility of BFP (3.94 x10(-5) MPa(-1)) is in agreement with experimental values of globular proteins. The center-of-mass diffusion coefficient of BFP and its dependence on temperature and pressure, which plays an important role in its application as a probe for intracellular liquid viscosity measurement, was calculated and found to be in good agreement with photobleaching recovery experimental data. We have shown that decreased temperature as well as applied pressure increases the water viscosity, but the concomitant decrease of the BFP diffusion coefficient behaves differently from Stokes-Einstein formula. It is shown that the number of hydrogen bonds around the protein grows with pressure, which explains the aforementioned deviation. Pressure also reduces root mean square (RMS) fluctuations, especially those of the most flexible residues situated in the loops. The analysis of the RMS fluctuations of the backbone Calpha atoms also reveals that the most rigid part of BFP is the center of the beta-barrel, in accord with temperature B factors obtained from the Protein Data Bank.     

The 2.1A crystal structure of the far-red fluorescent protein HcRed: inherent conformational flexibility of the chromophore

We have determined the crystal structure of HcRed, a far-red fluorescent protein isolated from Heteractis crispa, to 2.1A resolution. HcRed was observed to form a dimer, in contrast to the monomeric form of green fluorescent protein (GFP) or the tetrameric forms of the GFP-like proteins (eqFP611, Rtms5 and DsRed). Unlike the well-defined chromophore conformation observed in GFP and the GFP-like proteins, the HcRed chromophore was observed to be considerably mobile. Within the HcRed structure, the cyclic tripeptide chromophore, Glu(64)-Tyr(65)-Gly(66), was observed to adopt both a cis coplanar and a trans non-coplanar conformation. As a result of these two conformations, the hydroxyphenyl moiety of the chromophore makes distinct interactions within the interior of the beta-can. These data together with a quantum chemical model of the chromophore, suggest the cis coplanar conformation to be consistent with the fluorescent properties of HcRed, and the trans non-coplanar conformation to be consistent with non-fluorescent properties of hcCP, the chromoprotein parent of HcRed. Moreover, within the GFP-like family, it appears that where conformational freedom is permissible then flexibility in the chromophore conformation is possible.     

Structural basis for the fast maturation of Arthropoda green fluorescent protein

Since the cloning of Aequorea victoria green fluorescent protein (GFP) in 1992, a family of known GFP-like proteins has been growing rapidly. Today, it includes more than a hundred proteins with different spectral characteristics cloned from Cnidaria species. For some of these proteins, crystal structures have been solved, showing diversity in chromophore modifications and conformational states. However, we are still far from a complete understanding of the origin, functions and evolution of the GFP family. Novel proteins of the family were recently cloned from evolutionarily distant marine Copepoda species, phylum Arthropoda, demonstrating an extremely rapid generation of fluorescent signal. Here, we have generated a non-aggregating mutant of Copepoda fluorescent protein and solved its high-resolution crystal structure. It was found that the protein beta-barrel contains a pore, leading to the chromophore. Using site-directed mutagenesis, we showed that this feature is critical for the fast maturation of the chromophore.     

Vibrations-determined properties of green fluorescent protein

The physicochemical characteristics of the green fluorescent protein (GFP), including the thermodynamic properties (entropy, enthalpy, Gibbs' free energy, heat capacity), normal mode vibrations, and atomic fluctuations, were investigated. The Gaussian 03 computational chemistry program was employed for normal mode analysis using the AMBER force field. The thermodynamic parameters and atomic fluctuations were then calculated from the vibrational eigenvalues (frequencies) and eigenvectors. The regions of highest rigidity were shown to be the beta-sheet barrel with the central alpha-helix, which bears the chromophore. The most flexible parts of the GFP molecule were the outlying loops that cover the top and bottom of the beta-barrel. This way, the balance between rigidity and flexibility is maintained, which is the optimal relationship for protein stability in terms of Gibbs' free energy. This dual-schemed structure satisfies the requirements for GFP function. In this sense, the structure of GFP resembles a nanoscale drum: a stiff cylinder with flexible vibrating end(s).     

Circularly permuted variants of the green fluorescent protein.

Folding of the green fluorescent protein (GFP) from Aequorea victoria is characterized by autocatalytic formation of its p-hydroxybenzylideneimidazolidone chromophore, which is located in the center of an 11-stranded beta-barrel. We have analyzed the in vivo folding of 20 circularly permuted variants of GFP and find a relatively low tolerance towards disruption of the polypeptide chain by introduction of new termini. All permuted variants with termini in strands of the beta-barrel and about half of the variants with termini in loops lost the ability to form the chromophore. The thermal stability of the permuted GFPs with intact chromophore is very similar to that of the wild-type, indicating that chromophore-side chain interactions strongly contribute to the extraordinary stability of GFP.     

Minimum sequence requirements for the binding of paromomycin to the rRNA decoding site A

We have recently introduced a computational methodology that combines molecular dynamics (MD) simulations, free-energy calculations, and in vitro binding assays to predict the minimum RNA structural requirements for selective, high-affinity RNA binding to small-molecule ligands. Here, we show that this methodology can be applied to the conformationally flexible aminoglycoside antibiotic paromomycin. A RNA consisting of an 11-mer:10-mer duplex that contains one 16S ribosome RNA decoding A-site bound to paromomycin was simulated for 4 ns. The methodology predicts that the 11-mer:10-mer duplex binds to paromomycin with high affinity, whereas smaller RNA duplexes lose complex stability and the ability to bind paromomycin. The predicted high-affinity binding to paromomycin of the 11-mer:10-mer duplex was confirmed experimentally (EC(50) = 0.28 microM), as well as the inability of smaller complexes to bind. Our simulations show good agreement with experiment for dynamic and structural properties of the isolated A-site, including hydrogen-bonding networks and RNA structural rearrangements upon ligand binding. The results suggest that MD simulations can supplement in vitro methods as a tool for predicting minimum RNA-binding motifs for both small, rigid ligands, and large, flexible ligands when structural information is available.     

The 2.0-A crystal structure of eqFP611, a far red fluorescent protein from the sea anemone Entacmaea quadricolor

We have crystallized and subsequently determined to 2.0-A resolution the crystal structure of eqFP611, a far red fluorescent protein from the sea anemone Entacmaea quadricolor. The structure of the protomer, which adopts a beta-can topology, is similar to that of the related monomeric green fluorescent protein (GFP). The quaternary structure of eqFP611, a tetramer exhibiting 222 symmetry, is similar to that observed for the more closely related red fluorescent protein DsRed and the chromoprotein Rtms5. The unique chromophore sequence (Met63-Tyr64-Gly65) of eqFP611, adopts a coplanar and trans conformation within the interior of the beta-can fold. Accordingly, the eqFP611 chromophore adopts a significantly different conformation in comparison to the chromophore conformation observed in GFP, DsRed, and Rtms5. The coplanar chromophore conformation and its immediate environment provide a structural basis for the far red, highly fluorescent nature of eqFP611. The eqFP611 structure extends our knowledge on the range of conformations a chromophore can adopt within closely related members of the green fluorescent protein family.     

Molecular dynamics simulations of enhanced green fluorescent proteins: effects of F64L,S65T and T203Y mutations on the ground-state proton equilibria

Molecular dynamics simulations with the Amber force field are carried out to study two mutants of the green fluorescent protein (GFP), namely EGFP (F64L/S65T) and T203Y-EGFP (E(2)GFP). Those variants display an opposite equilibrium between the structural A and B states, associated with neutral and anionic protonation forms of the chromophore. Configurations of those two states are simulated for each variant and the energetics of their equilibrium in the two mutants is studied by evaluating the change in the relative free energy of A and B states (DeltaG(AB)) upon T203Y mutation. The resulting DeltaDeltaG(AB) agrees with the value inferred from absorption measurements. A comparison of the hydrogen bond network around the chromophore rationalizes the different population of state A and B in EGFP and E(2)GFP. On the basis of structural and energetic considerations, a mechanism for destabilization of the neutral chromophore in S65T mutants is proposed. Simulations of the B state of the S65T variant and of WT GFP are also performed for comparison and to test the force field parameters of the chromophore derived for the present calculations. Possible paths of proton transfer leading to nonfluorescent states of the chromophore are discussed in light of the photodynamical behavior of GFP, as revealed by fluorescence correlation spectroscopy and single-molecule experiments.     

Targeted molecular dynamics simulation studies of binding and conformational changes in E. coli MurD

Enzymes involved in the biosynthesis of bacterial peptidoglycan, an essential cell wall polymer unique to prokaryotic cells, represent a highly interesting target for antibacterial drug design. Structural studies of E. coli MurD, a three-domain ATP hydrolysis driven muramyl ligase revealed two inactive open conformations of the enzyme with a distinct C-terminal domain position. It was hypothesized that the rigid body rotation of this domain brings the enzyme to its closed active conformation, a structure, which was also determined experimentally. Targeted molecular dynamics 1 ns-length simulations were performed in order to examine the substrate binding process and gain insight into structural changes in the enzyme that occur during the conformational transitions into the active conformation. The key interactions essential for the conformational transitions and substrate binding were identified. The results of such studies provide an important step toward more powerful exploitation of experimental protein structures in structure-based inhibitor design.     

MM/PBSA analysis of molecular dynamics simulations of bovine β‐lactoglobulin: Free energy gradients in conformational transitions?

The pH-driven opening and closure of beta-lactoglobulin EF loop, acting as a lid and closing the internal cavity of the protein, has been studied by molecular dynamics (MD) simulations and free energy calculations based on molecular mechanics/Poisson-Boltzmann (PB) solvent-accessible surface area (MM/PBSA) methodology. The forms above and below the transition pH differ presumably only in the protonation state of residue Glu89. MM/PBSA calculations are able to reproduce qualitatively the thermodynamics of the transition. The analysis of MD simulations using a combination of MM/PBSA methodology and the colony energy approach is able to highlight the driving forces implied in the transition. The analysis suggests that global rearrangements take place before the equilibrium local conformation is reached. This conclusion may bear general relevance to conformational transitions in all lipocalins and proteins in general.     

In silico mutations and molecular dynamics studies on a winged bean chymotrypsin inhibitor protein

Winged bean chymotrypsin inhibitor (WCI) has an intruding residue Asn14 that plays a crucial role in stabilizing the reactive site loop conformation. This residue is found to be conserved in the Kunitz (STI) family of serine protease inhibitors. To understand the contribution of this scaffolding residue on the stability of the reactive site loop, it was mutated in silico to Gly, Ala, Ser, Thr, Leu and Val and molecular dynamics (MD) simulations were carried out on the mutants. The results of MD simulations reveal the conformational variability and range of motions possible for the reactive site loop of different mutants. The N-terminus side of the scissile bond, which is close to a beta-barrel, is conformationally less variable, while the C-terminus side, which is relatively far from any such secondary structural element, is more variable and needs stability through hydrogen-bonding interactions. The simulated structures of WCI and the mutants were docked in the peptide-binding groove of the cognate enzyme chymotrypsin and the ability to form standard hydrogen-bonding interactions at P3, P1 and P2' residues were compared. The results of the MD simulations coupled with docking studies indicate that hydrophobic residues like Leu and Val at the 14th position are disruptive for the integrity of the reactive site loop, whereas a residue like Thr, which can stabilize the C-terminus side of the scissile bond, can be predicted at this position. However, the size and charge of the Asn residue made it most suitable for the best maintenance of the integrity of the reactive site loop, explaining its conserved nature in the family.     

In meso crystal structure and docking simulations suggest an alternative proteoglycan binding site in the OpcA outer membrane adhesin

OpcA is an integral outer membrane adhesin protein from Neisseria meningitidis, the causative agent of meningococcal meningitis and septicemia. It binds to sialic acid (SA)-containing polysaccharides on the surface of epithelial cells. The crystal structure of OpcA showed that the protein adopts a 10-stranded beta-barrel structure, with five extensive loop regions on the extracellular side of the membrane. These form a crevice structure, lined with basic residues, which was hypothesized to act as the binding site for polysaccharide ligands. In the current study, a distinctly different OpcA structure has been obtained using crystals grown from a lipidic mesophase. Comparison of the two structures shows that the largest loop (L2), which closes over the end of the beta-barrel in the original crystal form, adopts a much more extended structure by reaching outward and away from the protein. The difference in conformation may be attributable to the absence of zinc ions from the crystallization conditions for the in meso crystal form: in the original structure, two zinc ions were bound to the external loops. Molecular dynamics (MD) simulations performed on the two OpcA models in a lipid bilayer environment demonstrated pronounced loop mobility. These observations support the view that the loop regions of OpcA are capable of a high degree of conformational flexibility. The original binding site for polysaccharide is not present in the in meso crystal form, and is disrupted during MD simulations. Docking analysis suggests a putative alternative location for the SA ligand in the new crystal form of OpcA.     

Interplay between structural rigidity and electrostatic interactions in the ligand binding domain of GluR2

Using molecular dynamics (MD) simulations, computational protein modifications, and a novel theoretical methodology that determines structural rigidity/flexibility (the FIRST algorithm), we investigate how molecular structure and dynamics of the glutamate receptor ligand binding domain (GluR2 S1S2) facilitate its conformational transition. S1S2 is a two-lobe protein, which undergoes a cleft closure conformational transition upon binding an agonist in the cleft between the two lobes; hence it is expected that the mechanism of this conformational transition can be characterized as a hinge-type. However, in the rigidity analysis one lobe of the protein is identified as a single rigid cluster while the other one is structurally flexible, inconsistent with a presumed mechanical hinge mechanism. Instead, we characterize the cleft-closing transition as a load and lock mechanism. We find that when two cross-cleft hydrogen bonds are disrupted the protein undergoes a rapid cleft opening transition. At the same time, the dynamical behavior of the cleft in the presence of the glutamate ligand is only weakly affected by the S652 peptide bond in its flipped conformation observed in the crystal structure. The residue E705 plays significant role in stabilization of the closed conformation via electrostatic interactions. The presence of the E705-K730 salt bridge seems to correlate strongly withthe cleft opening transition in the MD simulations.     

Environment effects on the oscillatory unfolding kinetics of GFP

The chromophore of a green fluorescent protein (GFP) mutant engineered to enhance emission and stability is known to display erratic switchings among a few of its chemical substates and, in particular, between the anionic A and the neutral N substates, whose difference is associated with a proton exchange and a consequent conformation rearrangement. However, when close to unfolding, the A-N switchings suddenly become very regular as shown by fluorescence oscillations that have been recently observed for molecules embedded in wet silica gel. In order to establish whether the matrix hosting the protein is responsible for these oscillations, we investigated the effect of another medium (silanized surfaces), of a different denaturant (urea) and of cosolvents (D(2)O and glycerol). The occurrence of periodic A-N switchings, in the last milliseconds before GFP unfolding, is observed under all investigated conditions, together with three specific frequency values that characterize the pre-unfolding fluorescence. Urea and guanidinium, the denaturants employed in order to unfold GFP, do not lead to appreciable differences in the observed switching parameters, whereas the different media embedding the protein give rise only to frequency shifts that scale with the viscosity of the host. The periodicity of the GFP A-N switchings and their dependence on cosolvents suggest that they could be associated with oscillatory motions between meta-stable conformations of the beta-barrel surrounding the chromophore near protein unfolding.     

The 1.7 A crystal structure of Dronpa: a photoswitchable green fluorescent protein

The green fluorescent protein (GFP), its variants, and the closely related GFP-like proteins possess a wide variety of spectral properties that are of widespread interest as biological tools. One desirable spectral property, termed photoswitching, involves the light-induced alteration of the optical properties of certain GFP members. Although the structural basis of both reversible and irreversible photoswitching events have begun to be unraveled, the mechanisms resulting in reversible photoswitching are less clear. A novel GFP-like protein, Dronpa, was identified to have remarkable light-induced photoswitching properties, maintaining an almost perfect reversible photochromic behavior with a high fluorescence to dark state ratio. We have crystallized and subsequently determined to 1.7 A resolution the crystal structure of the fluorescent state of Dronpa. The chromophore was observed to be in its anionic form, adopting a cis co-planar conformation. Comparative structural analysis of non-photoactivatable and photoactivatable GFPs, together with site-directed mutagenesis of a position (Cys62) within the Dronpa chromophore, has provided a basis for understanding Dronpa photoactivation. Specifically, we propose a model of reversible photoactivation whereby irradiation with light leads to subtle conformational changes within and around the environment of the chromophore that promotes proton transfer along an intricate polar network.     

X-ray structure of Cerulean GFP: a tryptophan-based chromophore useful for fluorescence lifetime imaging

The crystal structure of the cyan-fluorescent Cerulean green fluorescent protein (GFP), a variant of enhanced cyan fluorescent protein (ECFP), has been determined to 2.0 A. Cerulean bears an internal fluorophore composed of an indole moiety derived from Y66W, conjugated to the GFP-like imidazolinone ring via a methylene bridge. Cerulean undergoes highly efficient fluorescence resonance energy transfer (FRET) to yellow acceptor molecules and exhibits significantly reduced excited-state heterogeneity. This feature was rationally engineered in ECFP by substituting His148 with an aspartic acid [Rizzo et al. (2004) Nat. Biotechnol. 22, 445], rendering Cerulean useful for fluorescence lifetime imaging microscopy (FLIM). The X-ray structure is consistent with a single conformation of the chromophore and surrounding residues and may therefore provide a structural rationale for the previously described monoexponential fluorescence decay. Unexpectedly, the carboxyl group of H148D is found in a buried position, directly contacting the indole nitrogen of the chromophore via a bifurcated hydrogen bond. Compared to the similarly constructed ECFP chromophore, the indole group of Cerulean is rotated around the methylene bridge to adopt a cis-coplanar conformation with respect to the imidazolinone ring, resulting in a close edge-to-edge contact of the two ring systems. The double-humped absorbance spectrum persists in single-crystal absorbance measurements, casting doubt on the idea that ground state conformational heterogeneity forms the basis of the two overlapping transitions. At low pH, a blue shift in absorbance of 10-15 nm suggests a pH-induced structural transition that proceeds with a time constant of 47 (+/-2) min and is reversible. Possible interpretations in terms of chromophore isomerization are presented.     

Molecular modeling of formate dehydrogenase: the formation of the Michaelis complex

The formation of the reactive enzyme–substrate complex of formate dehydrogenase has been investigated by molecular dynamics techniques accounting for different conformational states of the enzyme. Simulations revealed that the transport of substrate to the active site through the substrate channel proceeds in the open conformation of enzyme due to the crucial role of the Arg284 residue acting as a vehicle. However, formate binding in the active site of the open conformation leads to the formation of a nonproductive enzyme–substrate complex. The productive Michaelis complex is formed only in the closed enzyme conformation after the substrate and coenzyme have bound, when required rigidity of the binding site and reactive formate orientation due to interactions with Arg284, Asn146, Ile122, and His332 residues is attained. Then, the high occupancy (up to 75%) of the reactive substrate–coenzyme conformation is reached, which was demonstrated by hybrid quantum mechanics/molecular mechanics simulations using various semiempirical Hamiltonians.