[en] 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₂O-H₂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 disorder 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