Inspecting fluctuation and coordination around chromophore inside green fluorescent protein from water to nonpolar solvent

Green fluorescent protein (GFP) is a widely used biomarker that demands systematical rational approaches to its structure function redesign. In this work, we mainly utilized atomistic molecular dynamics simulations to inspect and visualize internal fluctuation and coordination around chromophore inside GFP, from water to nonpolar octane solvent. We found that GFP not only maintains its β-barrel structure well into the octane, but also sustains internal residue and water coordination to position the chromophore stably while suppress dihedral fluctuations of the chromophore, so that functional robustness of GFP is achieved. Our accompanied fluorescence microscope measurements accordingly confirmed the GFP functioning into the octane. Furthermore, we identified that crucial water sites inside GFP along with permeable pores on the β-barrel of the protein are largely preserved from the water to the octane solvent, which allows sufficiently fast exchanges of internal water with the bulk or with the water layer kept on the surface of the protein. By additionally pulling GFP from bulk water to octane, we suggest that the GFP function can be well maintained into the nonpolar solvent as long as, first, the protein does not denature in the nonpolar solvent nor across the polar-nonpolar solvent interface; second, a minimal set of water molecules are in accompany with the protein; third, the nonpolar solvent molecules may need to be large enough to be nonpermeable via the water pores on the β-barrel.     

Ligand binding to heme proteins: connection between dynamics and function

Ligand binding to heme proteins is studied by using flash photolysis over wide ranges in time (100 ns-1 ks) and temperature (10-320 K). Below about 200 K in 75% glycerol/water solvent, ligand rebinding occurs from the heme pocket and is nonexponential in time. The kinetics is explained by a distribution, g(H), of the enthalpic barrier of height H between the pocket and the bound state. Above 170 K rebinding slows markedly. Previously we interpreted the slowing as a "matrix process" resulting from the ligand entering the protein matrix before rebinding. Experiments on band III, an inhomogeneously broadened charge-transfer band near 760 nm (approximately 13,000 cm-1) in the photolyzed state (Mb*) of (carbonmonoxy)myoglobin (MbCO), force us to reinterpret the data. Kinetic hole-burning measurements on band III in Mb* establish a relation between the position of a homogeneous component of band III and the barrier H. Since band III is red-shifted by 116 cm-1 in Mb* compared with Mb, the relation implies that the barrier in relaxed Mb is 12 kJ/mol higher than in Mb*. The slowing of the rebinding kinetics above 170 K hence is caused by the relaxation Mb*----Mb, as suggested by Agmon and Hopfield [(1983) J. Chem. Phys. 79, 2042-2053]. This conclusion is supported by a fit to the rebinding data between 160 and 290 K which indicates that the entire distribution g(H) shifts. Above about 200 K, equilibrium fluctuations among conformational substates open pathways for the ligands through the protein matrix and also narrow the rate distribution. The protein relaxations and fluctuations are nonexponential in time and non-Arrhenius in temperature, suggesting a collective nature for these protein motions. The relaxation Mb*----Mb is essentially independent of the solvent viscosity, implying that this motion involves internal parts of the protein. The protein fluctuations responsible for the opening of the pathways, however, depend strongly on the solvent viscosity, suggesting that a large part of the protein participates. While the detailed studies concern MbCO, similar data have been obtained for MbO2 and CO binding to the beta chains of human hemoglobin and hemoglobin Zürich. The results show that protein dynamics is essential for protein function and that the association coefficient for binding from the solvent at physiological temperatures in all these heme proteins is governed by the barrier at the heme.     

Dynamics of ligand rebinding to unfolded MbCO by guanidine HCl

Femtosecond vibrational spectroscopy was used to probe a functionally important dynamics and residual structure of myoglobin unfolded by 4 M guanidine HCl. The spectra of the dissociated CO indicated that the residual structure of unfolded myoglobin (Mb) forms a few hydrophobic cavities that could accommodate the dissociated ligand. Geminate rebinding (GR) of CO to the unfolded Mb is three-orders-of-magnitude faster and more efficient than the native Mb but similar to a model heme in a viscous solvent, suggesting that the GR of CO to heme is accelerated by the longer retention of the dissociated ligand near the Fe atom by the poorly-structured protein matrix of the unfolded Mb or viscous solvent. The inefficient GR of CO in native Mb, while dissociated CO is trapped in the primary heme pocket located near the active binding site, indicates that the tertiary structure of the pocket in native Mb plays a functionally significant role.     

转基因斑马鱼分析胰岛β-细胞发育情况

斑马鱼的个体小、高产和体外受精等特点使其已经迅速成为研究脊椎动物器官发育和人类疾病的模式生物之一。我们建立了一个转基因斑马鱼动物模型来研究胰岛β-细胞的发育。首先,构建了斑马鱼胰岛素（Insulin ,INS） 启动子与绿色荧光蛋白(GFP) 组成的表达载体, 命名为INS:GFP。其次,将质粒在斑马鱼1-细胞期注射到细胞质内。最后我们成功获得了生殖系稳定遗传胰岛素转基因斑马鱼,在成鱼和幼鱼期均可以通过GFP标记β-细胞。通过方便的荧光筛选,我们观察到胰岛在受精后18h开始形成,1－5d后由初始的脊索中线两侧向右迁移。从我们构建的胰岛素转基因斑马鱼,可以直观判断胰岛的发育情况,为研究胰岛的发育、损伤和再生提供了一个简便和直观的新型工具。     

Structural dynamics of distal histidine replaced mutants of myoglobin accompanied with the photodissociation reaction of the ligand

Protein dynamics observed by the transient grating (TG) method are studied for some site-directed mutants at the distal histidine of myoglobin (H64L, H64Q, H64V). The time profiles of the TG signals are very sensitive to the amino acid residue of the 64 position. It was found that the sensitivity is mostly caused by the different rates of the ligand escape from the protein to solvent and the magnitude of the molecular volume change. Several molecular origins of the volume difference between MbCO and Mb, such as the electrostatic interaction in the distal pocket, movement of helices, and distal water, are proposed. Interestingly, the volume difference between the CO-trapped Mb inside the protein interior and Mb is similar to that of the partial molar volume of CO in organic solvent. The effect of mutation on the nature of the CO trapped site is discussed.     

Myosin-induced volume increase of the hyper-mobile water surrounding actin filaments

Microwave dielectric spectroscopy can measure the rotational mobility of water molecules that hydrate proteins and the hydration-shell volume. Using this technique, we have recently shown that apart from typical hydrating water molecules with lowered mobility there are other water molecules around the actin filaments (F-actin) which have a much higher mobility than that of bulk water [Biophys. J. 85 (2003) 3154]. We report here that the volume of this water component (hyper-mobile water) markedly increases without significant change of the volume of the ordinary hydration shell when the myosin motor-domain (S1, myosin subfragment-1) binds to F-actin. No hyper-mobile component was found in the hydration shell of S1 itself. The present results strongly suggest that the solvent space around S1 bound to F-actin is diffusionally asymmetric, which supports our model of force generation by actomyosin proposed previously [op. cit.].     

Neutron Spin-Echo Studies of Hemoglobin and Myoglobin: Multiscale Internal Dynamics

Neutron spin-echo spectroscopy was used to study structural fluctuations that occur in hemoglobin (Hb) and myoglobin (Mb) in solution. Using neutron spin-echo data up to a very high momentum transfer q (∼ 0.62 Å^− 1), we characterized the internal dynamics of these proteins at the levels of dynamic pair correlation function and self-correlation function in the time range of several picoseconds to a few nanoseconds. In the same protein solution, data transition from pair correlation motion to self-correlation motion as the momentum transfer q increases. At low q, coherent scattering dominates; at high q, observations are largely due to incoherent scattering. The low q data were interpreted in terms of an effective diffusion coefficient; on the other hand, the high q data were interpreted in terms of mean square displacements. Comparison of data from the two homologous proteins collected at different temperatures and protein concentrations was used to assess the contributions made by translational and rotational diffusion and internal modes of motion to the data. The temperature dependence of decay times can be attributed to changes in the viscosity and temperature of the solvent, as predicted by the Stokes-Einstein relationship. This is true for contributions from both diffusive and internal modes of motion, indicating an intimate relationship between the internal dynamics of the proteins and the viscosity of the solvent. Viscosity change associated with protein concentration can account for changes in diffusion observed at different concentrations, but is apparently not the only factor involved in the changes in internal dynamics observed with change in protein concentration. Data collected at high q indicate that internal modes in Mb are generally faster than those in Hb, perhaps due to the greater surface-to-volume ratio of Mb and the fact that surface groups tend to exhibit faster motion than buried groups. Comparison of data from Hb and data from Mb at low q indicates an unexpectedly rapid motion of Hb αβ dimers relative to one another. Dynamic motion of subunits is increasingly perceived as important to the allosteric behavior of Hb. Our data demonstrate that this motion is highly sensitive to protein concentration, temperature, and solvent viscosity, indicating that great care needs to be exercised in interpreting its effect on protein function.     

Time resolved thermodynamics associated with ligand photorelease in heme peroxidases and globins: Open access channels versus gated ligand release

Heme proteins represent a diverse class of biomolecules responsible for an extremely diverse array of physiological functions including electron transport, monooxygenation, ligand transport and storage, cellular signaling, respiration, etc. An intriguing aspect of these proteins is that such functional diversity is accomplished using a single type of heme macrocycle based upon iron protoporphyrin IX. The functional diversity originates from a delicate balance of inter-molecular interactions within the protein matrix together with well choreographed dynamics that modulate the heme electronic structure as well as ligand entry/exit pathways from the bulk solvent to the active site. Of particular interest are the dynamics and energetics associated with the entry/exit of ligands as this process plays a significant role in regulating the rates of heme protein activity. Time-resolved photoacoustic calorimetry (PAC) has emerged as a powerful tool through which to probe the underlying energetics associated with small molecule dissociation and release to the bulk solvent in heme proteins on time scales from tens of nanoseconds to several microseconds. In this review, the results of PAC studies on various classes of heme proteins are summarized highlighting how different protein structures affect the thermodynamics of ligand migration. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

The influence of solvent composition on global dynamics of human butyrylcholinesterase powders: a neutron-scattering study

A major result of incoherent elastic neutron-scattering experiments on protein powders is the strong dependence of the intramolecular dynamics on the sample environment. We performed a series of incoherent elastic neutron-scattering experiments on lyophilized human butyrylcholinesterase (HuBChE) powders under different conditions (solvent composition and hydration degree) in the temperature range from 20 to 285 K to elucidate the effect of the environment on the enzyme atomic mean-square displacements. Comparing D(2)O- with H(2)O-hydrated samples, we were able to investigate protein as well as hydration water molecular dynamics. HuBChE lyophilized from three distinct buffers showed completely different atomic mean-square displacements at temperatures above approximately 200 K: a salt-free sample and a sample containing Tris-HCl showed identical small-amplitude motions. A third sample, containing sodium phosphate, displayed highly reduced mean-square displacements at ambient temperature with respect to the other two samples. Below 200 K, all samples displayed similar mean-square displacements. We draw the conclusion that the reduction of intramolecular protein mean-square displacements on an Angstrom-nanosecond scale by the solvent depends not only on the presence of salt ions but also on their type.     

Viscosity dependence of the solute quenching of the tryptophanyl fluorescence of proteins

We have studied the viscosity dependence of the acrylamide quenching of the fluorescence on the internal tryptophan residues in cod parvalbumin and ribonuclease T1, as well as the model systems. N-acetyl-L-tryptophanamide and glucagon. For the latter systems, the apparent rate constant, kq(app), for acrylamide quenching shows a typical diffusion-limited behavior. For parvalbumin and ribonuclease T1, however, the viscosity dependence of kq(app) is quite different. There is little change in the kq(app) values on increasing the bulk viscosity from 1 to 10 cP (by addition of glycerol), but a further increase from 10 to 100 cP results in a significant reduction in the kq(app). Both an unfolding mechanism and a quencher penetration mechanism are considered to explain the results. Only the penetration mechanism is found to be consistent, and our data are interpreted as indicating that the rate-limiting step for quenching goes from being that of diffusion through the protein matrix, at low viscosity, to diffusion through the bulk solvent, at high viscosity. By also considering the Kramers' relationship in fitting our data, we are able to obtain insight regarding the coupling between internal fluctuations in the structure of the protein and motion of the bulk solvent. For parvalbumin and ribonuclease T1, the internal dynamics are found to be very weakly coupled to the bulk.     

Hydration energy landscape of the active site cavity in cytochrome P450cam

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein–ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol⁻¹) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable. Proteins 32:381–396, 1998. © 1998 Wiley-Liss, Inc.     

Geminate rebinding in trehalose-glass embedded myoglobins reveals residue-specific control of intramolecular trajectories.

It is becoming increasingly apparent that hydrophobic cavities (also referred to as xenon cavities) within proteins have significant functional implications. The potential functional role of these cavities in modulating the internal dynamics of carbon monoxide in myoglobin (Mb) is explored in the present study by using glassy matrices derived from trehalose to limit protein dynamics and to eliminate ligand exchange between the solvent and the protein. By varying the temperature (-15 to 65 degrees C) and humidity for samples of carbonmonoxy myoglobin embedded in trehalose-glass, it is possible to observe a hierarchy of distinct geminate recombination phases that extend from nanosecond to almost seconds that can be directly associated with rebinding from specific hydrophobic cavities. The use of mutant forms of Mb reveals the role of key residues in modulating ligand access between these cavities and the distal hemepocket.     

Thin filament length regulation in striated muscle sarcomeres: pointed-end dynamics go beyond a nebulin ruler

The actin (thin) filaments in striated muscle are highly regulated and precisely specified in length to optimally overlap with the myosin (thick) filaments for efficient myofibril contraction. Here, we review and critically discuss recent evidence for how thin filament lengths are controlled in vertebrate skeletal, vertebrate cardiac, and invertebrate (arthropod) sarcomeres. Regulation of actin polymerization dynamics at the slow-growing (pointed) ends by the capping protein tropomodulin provides a unified explanation for how thin filament lengths are physiologically optimized in all three muscle types. Nebulin, a large protein thought to specify thin filament lengths in vertebrate skeletal muscle through a ruler mechanism, may not control pointed-end actin dynamics directly, but instead may stabilize a large core region of the thin filament. We suggest that this stabilizing function for nebulin modifies the lengths primarily specified by pointed-end actin dynamics to generate uniform filament lengths in vertebrate skeletal muscle. We suggest that nebulette, a small homolog of nebulin, may stabilize a correspondingly shorter core region and allow individual thin filament lengths to vary according to working sarcomere lengths in vertebrate cardiac muscle. We present a unified model for thin filament length regulation where these two mechanisms cooperate to tailor thin filament lengths for specific contractile environments in diverse muscles.     

Fast scanning calorimetry of lysozyme unfolding at scanning rates from 5 K/min to 500,000 K/min

Background

Protein denaturation is often studied using differential scanning calorimetry (DSC). However, conventional instruments are limited in the temperature scanning rate available. Fast scanning calorimetry (FSC) provides an ability to study processes at much higher rates while using extremely small sample masses [ng]. This makes it a very interesting technique for protein investigation.

SANS/SAXS study of the BSA solvation properties in aqueous urea solutions via a global fit approach

We report on the solvation properties and intermolecular interactions of a model protein (bovine serum albumine, BSA) in urea aqueous solutions, as obtained by combining small-angle neutron and X-ray scattering experiments. According to a global fit strategy, all the whole set of scattering curves are analysed by considering a unique model which includes the BSA structure, the protein-protein interactions and the thermodynamic exchange process of water/urea molecules at the protein solvent interface. As a main result, the equilibrium constant that accounts for the difference in composition between the bulk solvent and the protein solvation layer is derived. Results confirm that urea preferentially sticks to the protein surface, inducing a noticeable change in both the repulsive and the attractive interaction potentials.     

A model for water motion in crystals of lysozyme based on an incoherent quasielastic neutron-scattering study

This paper reports an incoherent quasielastic neutron scattering study of the single particle, diffusive motions of water molecules surrounding a globular protein, the hen egg-white lysozyme. For the first time such an analysis has been done on protein crystals. It can thus be directly related and compared with a recent structural study of the same sample. The measurement temperature ranged from 100 to 300 K, but focus was on the room temperature analysis. The very good agreement between the structural and dynamical studies suggested a model for the dynamics of water in triclinic crystals of lysozyme in the time range approximately 330 ps and at 300 K. Herein, the dynamics of all water molecules is affected by the presence of the protein, and the water molecules can be divided into two populations. The first mainly corresponds to the first hydration shell, in which water molecules reorient themselves fivefold to 10-fold slower than in bulk solvent, and diffuse by jumps from hydration site to hydration site. The long-range diffusion coefficient is five to sixfold less than for bulk solvent. The second group corresponds to water molecules further away from the surface of the protein, in a second incomplete hydration layer, confined between hydrated macromolecules. Within the time scale probed they undergo a translational diffusion with a self-diffusion coefficient reduced approximately 50-fold compared with bulk solvent. As protein crystals have a highly crowded arrangement close to the packing of macromolecules in cells, our conclusion can be discussed with respect to solvent behavior in intracellular media: as the mobility is highest next to the surface, it suggests that under some crowding conditions, a two-dimensional motion for the transport of metabolites can be dominant.     

Insulin modulates rat liver glucocorticoid receptor

This investigation used cytosol fraction of rat liver to examine the effects of insulin (INS) on functional properties of glucocorticoid receptor (GR). Male Wistar rats (220-250 g b.wt.) were injected with INS (50 microg/200 g b.wt, i.p.) and 18 h after INS administration used for experiments. INS-stimulated dissociation of G-R complexes was significantly increased by 133% compared to control level. However, INS treatment significantly stimulated stability of GR protein by 138% above control value. Furthermore, results show that INS stimulated activation of formed cytosol [3H] TA-R complexes by 143% in respect to control. [3H]TA-R complexes from INS treated animals could be activated and accumulated at higher rate in cell nuclei of control animals. The physiological relevance of the data was confirmed by INS-related stimulation of Tryptophan oxigenase (TO) activity. It was observed that INS stimulated TO activity while INS injected to adrenalectomized rats, exhibited less effects compared to control. The results indicate that a glucocorticoid hormone (CORT) enhances INS induced stimulation of TO activity, as evidenced by enhanced enzyme activity. Presented data suggest: that INS treatment leads to modifications of the GR protein and the nuclear components and that INS activates the rat liver CORT signaling pathway which mediates, in part, the activity of TO.