Hydroxyl and water molecule orientations in trypsin: Comparison to molecular dynamics structures

A comparison is presented of experimentally observed hydroxyl and water hydrogens in trypsin determined from neutron density maps with the results of a 140ps molecular dynamics (MD) simulation. Experimental determination of hydrogen and deuterium atom positions in molecules as large as proteins is a unique capability of neutron diffraction. The comparison addresses the degree to which a standard force-field approach can adequately describe the local electrostatic and van der Waals forces that determine the orientations of these hydrogens. Neutron densities, derived from 2.1{Angstrom} D{sub 2}O-H{sub 2}O difference Fourier maps, provide a database of 27 well-ordered hydroxyl hydrogens. Most of the simulated hydroxyl orientations are within a standard deviation of the experimentally-observed positions, including several examples in which both the simulation and the neutron density indicate that a hydroxyl group is shifted from a {open_quote}standard{close_quote} rotamer. For the most highly ordered water molecules, the hydrogen distributions calculated from the trajectory were in good agreement with neutron density; simulated water molecules that displayed multiple hydrogen bonding networks had correspondingly broadened neutron density profiles. This comparison was facilitated by development of a method to construct a pseudo 2{Angstrom} density map based on the hydrogen atom distributions from the simulation. The degree of disordermore » of internal water molecules is shown to result primarily from the electrostatic environment surrounding that water molecule as opposed to the cavity size available to the molecule. A method is presented for comparing the discrete observations sampled in a dynamics trajectory with the time- averaged data obtained from X-ray or neutron diffraction studies. This method is particularly useful for statically-disordered water molecules, in which the average location assigned from a trajectory may represent a site of relatively low occupancy.« less