Solvation dynamics of a protein in the pre molten globule state

The nature of solvent molecules around proteins in native and different non-native states is crucial for understanding the protein folding problem. We have characterized two compact denatured states of glutaminyl-tRNA synthetase (GlnRS) under equilibrium conditions in the presence of a naturally occurring osmolyte, l-glutamate. The solvation dynamics of the compact denatured states and the fully unfolded state has been studied using a covalently attached probe, acrylodan, near the active site. The solvation dynamics progressively becomes faster as the protein goes from the native to the molten globule to the pre molten globule to the fully unfolded state. Anisotropy decay measurements suggest that the pre-molten-globule intermediate is more flexible than the molten globule although the secondary structure is largely similar. Dynamic light scattering studies reveal that both the compact denatured states are aggregated under the measurement conditions. The implications of solvation dynamics in aggregated compact denatured states have been discussed.     

Coupling between hydration layer dynamics and unfolding kinetics of HP-36

We have performed atomistic molecular dynamics simulations of aqueous solutions of HP-36 at 300 K in its native state, as well as at high temperatures to explore the unfolding dynamics of the protein and its correlation with the motion of water around it. On increasing the temperature a partially unfolded molten globule state is formed where the smallest alpha helix (helix 2) unfolds into a coil. It is observed that the unfolding is initiated around the residue Phe-18 which shows a sharp displacement during unfolding. We have noticed that the unfolding of the protein affects the density of water near the protein surface. Besides, the dynamics of water in the protein hydration layer has been found to be strongly correlated with the time evolution of the unfolding process. We have introduced and calculated a displacement time correlation function to monitor the change in water motion relative to the protein backbone during unfolding. We find that the unfolding of helix 2 is associated with an increase in mobility of water around it as compared to water around the other two helices. We have also explored the microscopic aspects of secondary structure specific and site specific solvation dynamics of the protein. The calculations reveal that unfolding influences the solvation dynamics of the protein molecule in a heterogeneous manner depending on the location of the polar probe residues. This seems to be in agreement with recent experimental findings.     

Dynamics of water in the hydration layer of a partially unfolded structure of the protein HP-36

Atomistic molecular dynamics simulations of the folded native structure and a partially unfolded molten globule structure of the protein villin headpiece subdomain or HP-36 have been carried out with explicit solvent to explore the effects of unfolding on the dynamical behavior of water present in the hydration layers of different segments (three alpha-helices) of the protein. The calculations revealed that the unfolding of helix-2 influences the translational and rotational motions of water present in the hydration layers of the three helices in a heterogeneous manner. It is observed that a correlation exists between the unfolding of helix-2 and the microscopic kinetics of protein-water hydrogen bonds formed by its residues. This in turn has an influence on the rigidity of the hydration layers of the helices in the unfolded structure versus that in the folded native structure. These results should provide a microscopic explanation to recent solvation dynamics experiments on folded native and unfolded structures of proteins.     

Relationship between the structure of the protein globule and bioluminescence spectra of firefly luciferase

Bioluminescence spectra for various native and mutant luciferases from fireflies and beetles were analyzed in the light of the known theoretical concepts on the influence of the microenvironment of the emitter on its emission spectra. The mechanism for the explanation of the nature of changing bioluminescence spectra for natural and artificial mutations of the amino acid residues in the protein globule of luciferases was proposed.     

Probing the effect of temperature on the backbone dynamics of the human alpha-lactalbumin molten globule

Molten globules are compact, partially folded proteins postulated to be general intermediates in protein folding. Human alpha-lactalbumin (alpha-LA) is a two-domain Ca(2+)-binding protein that partially unfolds at low pH to form a molten globule. NMR spectra of molten globules are characterized by broadened resonances due to conformational fluctuations on microsecond to millisecond time scales. These species are often studied at high temperature where NMR resonances are observed to sharpen. The effect of higher temperatures on fast time-scale backbone dynamics of molten globules has not been investigated previously. Here, 1D (15)N direct-detection and 2D indirect-detection (1)H-(15)N heteronuclear NOE experiments have been used to probe fast time-scale dynamics at low and high temperatures for three disulfide-bond variants of human alpha-LA that form molten globules. Disulfide bonds are found to have a significant effect on backbone dynamics within the beta-domain of the molten globule; within the alpha-domain, dynamics are not significantly influenced by these bonds. At 20 degrees C, backbone mobility is significantly decreased in both domains of the molten globule compared to the mobility at 40-50 degrees C. Heteronuclear NOE values determined at 20 degrees C for the alpha-domain are closely similar to those observed for native alpha-LA, indicating that the alpha-LA molten globule has even more native-like character than suggested by studies conducted at higher temperature. Our results highlight the importance of considering the temperature dependence of the molten globule ensemble when making comparisons between experimental data obtained under different conditions.     

Thermodynamic insights into the binding of ANS with the salt induced molten globule states of cytochrome c

The molten globule (MG) state or the A-state of cytochrome c (cyt c) has been induced by addition of salts sodium perchlorate (NaClO₄), sodium thiocyanate (NaSCN), and sodium sulphate (Na₂SO₄) at pH 2.0. Isothermal titration calorimetry (ITC) has been used for determining the energetics of binding of 8-anilino naphthalene sulphonate (ANS) to the salt induced A-state of cyt c, and the accompanying thermodynamic parameters have been analyzed to elucidate the nature of the interactions between ANS and the A-state of cyt c. Temperature dependent studies of the binding process reveal that the binding is not a two state process and there are more than a single type of interactions involved. Addition of a bulky salt tetraethylammonium bromide (TEAB) increases the stoichiometry of binding significantly, with a reduction in the binding affinity at a higher concentration. The results provide quantitative information on the binding of ANS to the salt induced molten globule states of cyt c. It is further inferred that the binding involves a combination of hydrophobic and electrostatic interactions.     

Site-specific hydration dynamics in the nonpolar core of a molten globule by dynamic nuclear polarization of water

Water-protein interactions play a direct role in protein folding. The chain collapse that accompanies protein folding involves extrusion of water from the nonpolar core. For many proteins, including apomyoglobin (apoMb), hydrophobic interactions drive an initial collapse to an intermediate state before folding to the final structure. However, the debate continues as to whether the core of the collapsed intermediate state is hydrated and, if so, what the dynamic nature of this water is. A key challenge is that protein hydration dynamics is significantly heterogeneous, yet suitable experimental techniques for measuring hydration dynamics with site-specificity are lacking. Here, we introduce Overhauser dynamic nuclear polarization at 0.35 T via site-specific nitroxide spin labels as a unique tool to probe internal and surface protein hydration dynamics with site-specific resolution in the molten globular, native, and unfolded protein states. The (1)H NMR signal enhancement of water carries information about the local dynamics of the solvent within ～10 ? of a spin label. EPR is used synergistically to gain insights on local polarity and mobility of the spin-labeled protein. Several buried and solvent-exposed sites of apoMb are examined, each bearing a covalently bound nitroxide spin label. We find that the nonpoloar core of the apoMb molten globule is hydrated with water bearing significant translational dynamics, only 4-6-fold slower than that of bulk water. The hydration dynamics of the native state is heterogeneous, while the acid-unfolded state bears fast-diffusing hydration water. This study provides a high-resolution glimpse at the folding-dependent nature of protein hydration dynamics.     

Molecular dynamics studies show solvation structure of type III antifreeze protein is disrupted at low pH

Antifreeze proteins are a class of biological molecules of interest in many research and industrial applications due to their highly specialized function, but there is little information of their stability and properties under varied pH derived from computational studies. To gain novel insights in this area, we conducted molecular dynamics (MD) simulations with the antifreeze protein 1KDF at varied temperatures and pH. Water solvation and H-bond formation around specific residues – ASN14, THR18 and GLN44 – involved in its antifreeze activity were extensively studied. We found that at pH1 there was a disruption in water solvation around the basal and the ice binding surfaces of the molecule. This was induced by a small change in the secondary structure propensities of some titrable residues, particularly GLU35. This change explains the experimentally observed reduction in antifreeze activity previously reported for this protein at pH1. We also found that THR18 showed extremely low H-bond formation, and that the three antifreeze residues all had very low average H-bond lifetimes. Our results confirm long-standing assumptions that these small, compact molecules can maintain their antifreeze activity in a wide range of pH, while demonstrating the mechanism that may reduce antifreeze activity at low pH. This aspect is useful when considering industrial and commercial use of antifreeze proteins subject to extreme pH environments, in particular in food industrial applications.     

The structure of lanthanide complexes with five membered sulfones and the specific benzene solvation

By means of a computer program based on the pseudo-contact equation, the structures of lanthanide complexes with five-membered sulfones have been established. The compounds are 1:1 bidentate complexes with the axis of complexation nearly parallel to the bisecting line of the O-S-O angle, and with the lanthanide position being at an average distance of 2.4 A from both oxygen atoms. The Ledal model, which assumes the specific benzene solvation, has been proposed as an adequate model for molecules with planar ring structure, and the distance between the benzene and solute along the preferential axis of complexation has been established as ~ 3.0 Å. The most probable conformation of the ring was taken into account because the accuracy of location of the lanthanide or benzene molecule in these complexes depends mainly on the correct model of the molecule assumed. Detailed analysis of the Eu(fod),-induced shifts in the PMR spectrum of NTBN sulfone confirmed the advanced structure of the complexes of five-membered ring sulfones with the lanthanide shift reagents.     

Microscopic structure and dynamics of molten Se50Te50 alloys

In this work we investigate the microscopic structure and dynamics of the molten equimolar alloy, Se(50)Te(50) using a combination of neutron and x-ray diffraction experiments, reverse Monte Carlo analysis, and first principles molecular dynamics. The range of temperatures studied covers the semiconductor/metal transition. From our results it can be seen that the latter is associated with an increase in coordination numbers and a reduced tendency to heterocoordination. In agreement with previous inelastic neutron scattering experiments, our molecular dynamics calculation predict a certain widening of the stretching vibrational modes band in connection with the increase of coordination and the presence of longer bonds in the metallic phase.     

The CN stretching band of aliphatic thiocyanate is sensitive to solvent dynamics and specific solvation

Infrared absorption spectra in the C[triple bond]N stretching frequency region were collected for methyl thiocyanate, the simplest model aliphatic thiocyanate, in several common solvents to establish the dependence of the C[triple bond]N spectral band of aliphatic thiocyanate on its local solvation environment. Systematic changes in the C[triple bond]N bandwidth indicate that it reports on fast solvation dynamics. Anomalous asymmetry and temperature dependence of the C[triple bond]N band in fluorinated alcohol solvents indicates that these solvents participate in formation of a discrete hydrogen-bonded complex with the C[triple bond]N end of methyl thiocyanate. These observations indicate that the C[triple bond]N band of thiocyanate could be an effective site-specific probe of both specific hydrogen bonding and local dynamics in more complex systems, such as peptides and proteins.     

Comparison between self-guided Langevin dynamics and molecular dynamics simulations for structure refinement of protein loop conformations

This article presents a comparative analysis of two replica‐exchange simulation methods for the structure refinement of protein loop conformations, starting from low‐resolution predictions. The methods are self‐guided Langevin dynamics (SGLD) and molecular dynamics (MD) with a Nosé–Hoover thermostat. We investigated a small dataset of 8‐ and 12‐residue loops, with the shorter loops placed initially from a coarse‐grained lattice model and the longer loops from an enumeration assembly method (the Loopy program). The CHARMM22 + CMAP force field with a generalized Born implicit solvent model (molecular‐surface parameterized GBSW2) was used to explore conformational space. We also assessed two empirical scoring methods to detect nativelike conformations from decoys: the all‐atom distance‐scaled ideal‐gas reference state (DFIRE‐AA) statistical potential and the Rosetta energy function. Among the eight‐residue loop targets, SGLD out performed MD in all cases, with a median of 0.48 Å reduction in global root‐mean‐square deviation (RMSD) of the loop backbone coordinates from the native structure. Among the more challenging 12‐residue loop targets, SGLD improved the prediction accuracy over MD by a median of 1.31 Å, representing a substantial improvement. The overall median RMSD for SGLD simulations of 12‐residue loops was 0.91 Å, yielding refinement of a median 2.70 Å from initial loop placement. Results from DFIRE‐AA and the Rosetta model applied to rescoring conformations failed to improve the overall detection calculated from the CHARMM force field. We illustrate the advantage of SGLD over the MD simulation model by presenting potential‐energy landscapes for several loop predictions. Our results demonstrate that SGLD significantly outperforms traditional MD in the generation and populating of nativelike loop conformations and that the CHARMM force field performs comparably to other empirical force fields in identifying these conformations from the resulting ensembles. Published 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

19F NMR studies of the native and denatured states of green fluorescent protein

Biosynthetic preparation and (19)F NMR experiments on uniformly 3-fluorotyrosine-labeled green fluorescent protein (GFP) are described. The (19)F NMR signals of all 10 fluorotyrosines are resolved in the protein spectrum with signals spread over 10 ppm. Each tyrosine in GFP was mutated in turn to phenylalanine. The spectra of the Tyr --> Phe mutants, in conjunction with relaxation data and results from (19)F photo-CIDNP (chemically induced dynamic nuclear polarization) experiments, yielded a full (19)F NMR assignment. Two (19)F-Tyr residues (Y92 and Y143) were found to yield pairs of signals originating from ring-flip conformers; these two residues must therefore be immobilized in the native structure and have (19)F nuclei in two magnetically distinct positions depending on the orientation of the aromatic ring. Photo-CIDNP experiments were undertaken to probe further the structure of the native and denatured states. The observed NMR signal enhancements were found to be consistent with calculations of the HOMO (highest occupied molecular orbital) accessibilities of the tyrosine residues. The photo-CIDNP spectrum of native GFP shows four peaks corresponding to the four tyrosine residues that have solvent-exposed HOMOs. In contrast, the photo-CIDNP spectra of various denatured states of GFP show only two peaks corresponding to the (19)F-labeled tyrosine side chains and the (19)F-labeled Y66 of the chromophore. These data suggest that the pH-denatured and GdnDCl-denatured states are similar in terms of the chemical environments of the tyrosine residues. Further analysis of the sign and amplitude of the photo-CIDNP effect, however, provided strong evidence that the denatured state at pH 2.9 has significantly different properties and appears to be heterogeneous, containing subensembles with significantly different rotational correlation times.     

Density dependence of the entropy and the solvation shell structure in supercritical water via molecular dynamics simulation

We perform molecular dynamics simulations of supercritical water (SCW) with a wide range of densities along a near critical isotherm using the simple point charge extended (SPC/E) pair potential in order to study the entropy and the solvation shell structure around a central water molecule. It is shown that both the translational and orientational two-particle correlation entropy terms can serve as the metrics of the translational-orientational structural orders in water and it is revealed that the translational structural order is very sensitive to the density variation in the gas-like and liquid-like region, while the orientational structural order is much more dependent upon compression in the medium-density SCW region. The comparison of the magnitudes of the full thermodynamic excess entropy and two-particle correlation entropy confirms the recent findings that the many-body terms other than two-body ones also present significant and non-neglectable contributions to the full excess entropy for the highly anomalous fluids like water. The analysis of entropy terms as a function of intermolecular distance and the orientational distribution functions as well as the three-dimensional spatial distribution functions indicate that the structural order occurs only in a much more diffused first solvation shell due to the elongated hydrogen bonds under supercritical conditions. It is revealed that no obvious second or higher neighbor shells occur in SCW, in contrast with the feature of normal liquid water that the anomalous decrease of translational order upon compression occurs mainly in the second shell.     

Insights into the interplay between electronic structure and protein dynamics: The case of ubiquitin

A quantum mechanical analysis of an experimental ensemble comprising 128 conformers of the protein ubiquitin has been carried out with the aid of LMO–SCF–COSMO calculations. The permanent dipole moment of the protein fluctuates in the range from 131 to 283 D while the energy-weighted average dipole has a magnitude of 197 D. The HOMO–LUMO energy gap of the conformational ensemble ranges from 7.389 to 8.397 eV and appears to being affected mainly by fluctuations in the HOMO energy. An inspection of the frontier orbitals of the 128 conformers indicates that their localization is not affected by the protein dynamics.     

Surface-enhanced resonance Raman spectroscopic characterization of the protein native structure

Surface-enhanced resonance Raman scattering (SERRS) spectra of biological species are often different from their resonance Raman (RR) spectra. A home-designed Raman flow system is used to determine the factors that contribute to the difference between the SERRS and RR of met-myoglobin (metMb). The results indicate that both the degree of protein-nanoparticles interaction and the laser irradiation contribute to the structural changes and are responsible for the observed differences between the SERRS and RR spectra of metMb. The prolonged adsorption of the protein molecules on the nanoparticle surface, which is the condition normally used for the conventional SERRS experiments, disturbs the heme pocket structure and facilitates the charge transfer process and the photoinduced transformation of proteins. The disruption of the heme pocket results in the loss of the distal water molecule, and the resulting SERRS spectrum of metMb shows a 5-coordinated high-spin heme. The flow system, when operated at a moderately high flow rate, can basically eliminate the factors that disturb the protein structure while maintaining a high enhancement factor. The SERRS spectrum obtained from a 1 x 10 (-7) M metMb solution using this flow system is basically identical to the RR spectrum of a 5 x 10 (-4) M metMb solution. Therefore, the Raman flow system reported here should be useful for characterizing the protein-nanoparticles interaction and the native structure of proteins using SERRS spectroscopy.     

Sensitivity of polar solvation dynamics to the secondary structures of aqueous proteins and the role of surface exposure of the probe

The structure and dynamics of water around a protein is expected to be sensitive to the details of the adjacent secondary structure of the protein. In this article, we explore this sensitivity by calculating both the orientational dynamics of the surface water molecules and the equilibrium solvation time correlation function of the polar amino acid residues in each of the three helical segments of the protein HP-36, using atomistic molecular dynamics simulations. The solvation dynamics of polar amino acid residues in helix-2 is found to be faster than that of the other two helices (the average time constant is smaller by a factor of 2), although the interfacial water molecules around helix-2 exhibit much slower orientational dynamics than that around the other two helices. A careful analysis shows that the origin of such a counterintuitive behavior lies in the dependence of the solvation time correlation function on the surface exposure of the probe-the more exposed is the probe, the faster the solvation dynamics. We discuss that these results are useful in explaining recent solvation dynamics experiments.     

Molecular dynamics simulation study on the transient response of solvation structure during the translational diffusion of solute

The transient response function of the density profile of the solvent around a solute during the translational diffusion of the solute is formulated based on the generalized Langevin formalism. The resultant theory is applied to both neat Lennard-Jones fluids and cations in liquid water, and the response functions are obtained from the analysis of the molecular dynamics simulations. In the case of the self-diffusion of Lennard-Jones fluids, the responses of the solvation structures are in harmony with conventional pictures based on the mode-coupling theory, that is, the binary collision in the low-density fluids, the backflow effect from medium to high density fluids, and the backscatter effect in the liquids near the triple point. In the case of cations in water, the qualitative behavior is strongly dependent on the size of cations. The pictures similar to simple dense liquids are obtained for the large ion and the neutral molecule, while the solvent waters within the first solvation shell of small ions show an oscillatory response in the short-time region. In particular, the oscillation is remarkably underdumped for lithium ion. The origin of the oscillation is discussed in relation to the theoretical treatment of the translational diffusion of ions in water.     

Electronic structure of the PYP chromophore in its native protein environment

We report on supermolecular ab initio calculations which clarify the role of the local amino acid environment in determining the unique electronic structure properties of the photoactive yellow protein (PYP) chromophore. The extensive ab initio calculations, at the level of the CC2 and EOM-CCSD methods, allow us to explicitly address how the interactions between the deprotonated p-coumaric thio-methyl ester (pCTM-) chromophore and the surrounding amino acids act together to create a specifically stabilized pCTM- species. Particularly noteworthy is the role of the Arg52 amino acid in stabilizing the chromophore against autoionization, and the role of the Tyr42 and Glu46 amino acids in determining the hydrogen-bonding properties that carry the dominant energetic effects.