Cellobiohydrolases processively hydrolyze glycosidic linkages in individual polymer chains of cellulose microfibrils, and typically exhibit specificity for either the reducing or nonreducing end of cellulose. Here, we conduct molecular dynamics simulations and free energy calculations to examine the initial binding of a cellulose chain into the catalytic tunnel of the reducing-end-specific Family 7 cellobiohydrolase (Cel7A) from Hypocrea jecorina. In unrestrained simulations, the cellulose diffuses into the tunnel from the -7 to the -5 positions, and the associated free energy profiles exhibit no barriers for initial processivity. The comparison of the free energy profiles for different cellulose chain orientations show a thermodynamic preference for the reducing end, suggesting that the preferential initial binding may affect the directional specificity of the enzyme by impeding nonproductive (nonreducing end) binding. Finally, the Trp-40 at the tunnel entrance is shown with free energy calculations to have a significant effect on initial chain complexation in Cel7A.

No abstract

No abstract

a b s t r a c tVarious structural models for amyloid b fibrils have been derived from a variety of experimental techniques. However, these models cannot differentiate between the relative position of the two arms of the b hairpin called the stagger. Amyloid fibrils of various hierarchical levels form lefthanded helices composed of b sheets. However it is unclear if positive, negative and zero staggers all form the macroscopic left-handed helices. To address this issue we have conducted extensive molecular dynamics simulations of amyloid b sheets of various staggers and shown that only negative staggers lead to the experimentally observed left-handed helices while positive staggers generate the incorrect right-handed helices. This result suggests that the negative staggers are physiologically relevant structure of the amyloid b fibrils.

We have performed 20 ns of fully atomistic molecular dynamics simulations of Hen Egg-White Lysozyme in 0, 10, 20, 30, and 100% by weight of glycerol in water to better understand the microscopic physics behind the bioprotection offered by glycerol to naturally occuring biological systems. The solvent exposure of protein surface residues changes when glycerol is introduced. The dynamic behavior of the protein, as quantified by the incoherent intermediate scattering function, shows a nonmonotonic dependence on glycerol content. The fluctuations of the protein residues with respect to each other were found to be similar in all water-containing solvents, but different from the pure glycerol case. The increase in the number of protein-glycerol hydrogen bonds in glycerol-water binary mixtures explains the slowing down of protein dynamics as the glycerol content increases. We also explored the dynamic behavior of the hydration layer. We show that the short length scale dynamics of this layer are insensitive to glycerol concentration. However, the long length scale behavior shows a significant dependence on glycerol content. We also provide insights into the behavior of bound and mobile water molecules.