We report experimental and theoretical studies on water and protein dynamics following photoexcitation of apomyoglobin. Using site-directed mutation and with femtosecond resolution, we experimentally observed relaxation dynamics with a biphasic distribution of time scales, 5 and 87 ps, around the site Trp7. Theoretical studies using both linear response and direct nonequilibrium molecular dynamics (MD) calculations reproduced the biphasic behavior. Further constrained MD simulations with either frozen protein or frozen water revealed the molecular mechanism of slow hydration processes and elucidated the role of protein fluctuations. Observation of slow water dynamics in MD simulations requires protein flexibility, regardless of whether the slow Stokes shift component results from the water or protein contribution. The initial dynamics in a few picoseconds represents fast local motions such as reorientations and translations of hydrating water molecules, followed by slow relaxation involving strongly coupled water-protein motions. We observed a transition from one isomeric protein configuration to another after 10 ns during our 30 ns ground-state simulation. For one isomer, the surface hydration energy dominates the slow component of the total relaxation energy. For the other isomer, the slow component is dominated by protein interactions with the chromophore. In both cases, coupled water-protein motion is shown to be necessary for observation of the slow dynamics. Such biologically important water-protein motions occur on tens of picoseconds. One significant discrepancy exists between theory and experiment, the large inertial relaxation predicted by simulations but clearly absent in experiment. Further improvements required in the theoretical model are discussed.

Molecular dynamics simulations are used to calculate the time dependent Stokes shift following photoexcitation of Trp-7 (W7) in myoglobin. In agreement with experiment, a long time (~60ps) component is observed. Since the long time Stokes shift component is absent when we repeat the calculation with protein frozen at the instant of photoexcitation, we firmly establish that protein flexibility is required to observe slow Stokes shift dynamics in this case. A transition between substates near the middle of a 30ns ground state trajectory gave us an opportunity to compare solvation dynamics in two different environments. While some of the superficial features are different, we find that the underlying dynamics are shared by the two isomers. It is necessary to look beyond a decomposition of the Stokes shift into protein and water contributions, and probe the underlying dynamics of protein side groups, backbone, and water dynamics to obtain a full picture of the relaxation process. We analyze water residence times, diffusion, and reorientation dynamics in the hydration layer. We find slow components in each of these quantities, and critically examine their origin and how they affect the observed Stokes shift.

Melittin, an amphipathic peptide from honeybee venom, consists of 26 amino acid residues and adopts different conformations from a random coil, to an alpha-helix, and to a self-assembled tetramer under certain aqueous environments. We report here our systematic studies of the hydration dynamics in these conformations using single intrinsic tryptophan (W19) as a molecular probe. With femtosecond resolution, we observed the solvation dynamics occurring in 0.62 and 14.7 ps in a random-coiled primary structure. The former represents bulklike water motion, and the latter reflects surface-type hydration dynamics of proteins. As a comparison, a model tripeptide (KWK) was also studied. At a membrane-water interface, melittin folds into a secondary alpha-helical structure, and the interfacial water motion was found to take as long as 114 ps, indicating a well-ordered water structure along the membrane surface. In high-salt aqueous solution, the dielectric screening and ionic solvation promote the hydrophobic core collapse in melittin aggregation and facilitate the tetramer formation. This self-assembled tertiary structure is also stabilized by the strong hydrophilic interactions of charged C-terminal residues and associated ions with water molecules in the two assembled regions. The hydration dynamics was observed to occur in 87 ps, significantly slower than typical water relaxation at protein surfaces but similar to water motion at membrane interfaces. Thus, the observed time scale of approximately 100 ps probably implies appropriate water mobility for mediating the formation of high-order structures of melittin in an alpha-helix and a self-assembled tetramer. These results elucidate the critical role of hydration dynamics in peptide conformational transitions and protein structural stability and integrity.

It is well known that there are some similarities among various naturally occurring amino acids. Thus, the complexity in protein systems could be reduced by sorting these amino acids with similarities into groups and then protein sequences can be simplified by reduced alphabets. This paper discusses how to group similar amino acids and whether there is a minimal amino acid alphabet by which proteins can be folded. Various reduced alphabets are obtained by reserving the maximal information for the simplified protein sequence compared with the parent sequence using global sequence alignment. With these reduced alphabets and simplified similarity matrices, we achieve recognition of the protein fold based on the similarity score of the sequence alignment. The coverage in dataset SCOP40 for various levels of reduction on the amino acid types is obtained, which is the number of homologous pairs detected by program BLAST to the number marked by SCOP40. For the reduced alphabets containing 10 types of amino acids, the ability to detect distantly related folds remains almost at the same level as that by the alphabet of 20 types of amino acids, which implies that 10 types of amino acids may be the degree of freedom for characterizing the complexity in proteins.