Three entropic classes of side chain in a globular protein.

The relationship between the NMR methyl group axial order parameter and the side chain conformational entropy is investigated in inhibitor-bound and apo human HIV protease using molecular dynamics simulation. Three distinct entropic classes of methyl-bearing side chains, determined by the topological distance of the methyl group from the protein backbone (i.e.

Elastic and conformational softness of a globular protein.

Flexibility, or softness, is crucial for protein function and consists of a conformational component, involving jumps between potential wells, and an elastic component, involving fluctuations within the wells. Combining molecular dynamics simulation with incoherent neutron scattering and light scattering measurements on green fluorescent protein, we reveal a relationship between the intrawell fluctuations and elastic moduli of the protein. This finding leads to a simple means of experimentally separating the conformational from the elastic atomic displacements.

Coupled flexibility change in cytochrome P450cam substrate binding determined by neutron scattering, NMR, and molecular dynamics simulation.

Neutron scattering and nuclear magnetic resonance relaxation experiments are combined with molecular dynamics (MD) simulations in a novel, to our knowledge, approach to investigate the change in internal dynamics on substrate (camphor) binding to a protein (cytochrome P450cam). The MD simulations agree well with both the neutron scattering, which furnishes information on global flexibility, and the nuclear magnetic resonance data, which provides residue-specific order parameters. Decreased fluctuations are seen in the camphor-bound form using all three techniques, dominated by changes in specific regions of the protein.

Temperature-dependent dynamical transitions of different classes of amino acid residue in a globular protein.

The temperature dependences of the nanosecond dynamics of different chemical classes of amino acid residue have been analyzed by combining elastic incoherent neutron scattering experiments with molecular dynamics simulations on cytochrome P450cam. At T = 100-160 K, anharmonic motion in hydrophobic and aromatic residues is activated, whereas hydrophilic residue motions are suppressed because of hydrogen-bonding interactions. In contrast, at T = 180-220 K, water-activated jumps of hydrophilic side chains, which are strongly coupled to the relaxation rates of the hydrogen bonds they form with hydration water, become apparent.

Surface hydration amplifies single-well protein atom diffusion propagating into the macromolecular core.

The effect of surface hydration water on internal protein motion is of fundamental interest in molecular biophysics. Here, by decomposing the picosecond to nanosecond atomic motion in molecular dynamics simulations of lysozyme at different hydration levels into three components--localized single-well diffusion, methyl group rotation, and nonmethyl jumps--we show that the effect of surface hydration is mainly to increase the volume of the localized single-well diffusion. These diffusive motions are coupled in such a way that the hydration effect propagates from the protein surface into the dry core.

REACH coarse-grained simulation of a cellulose fiber.

A molecular level understanding of the structure, dynamics and mechanics of cellulose fibers can aid in understanding the recalcitrance of biomass to hydrolysis in cellulosic biofuel production. Here, a residue-scale REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-grained force field was derived from all-atom molecular dynamics (MD) simulations of the crystalline Iβ cellulose fibril. REACH maps the atomistic covariance matrix onto coarse-grained elastic force constants.