Molecular Dynamics Simulations of End-to-End Contact Formation in Hydrocarbon Chains in Water and Aqueous Urea Solution

Abstract

We probe the urea-denaturation mechanism using molecular dynamics simulations of an elementary "folding" event, namely, the formation of end-to-end contact in the linear hydrocarbon chain (HC) CH(3)(CH(2))(18)CH(3). Electrostatic effects are examined using a model HC in which one end of the chain is positively charged (+0.2e) and the other contains a negative charge (-0.2e). For these systems multiple transitions between "folded" (conformations in which the chain ends are in contact) and "unfolded" (end-to-end contact is broken) can be observed during 4 ns molecular dynamics simulations. In water and 6 M aqueous urea solution HC and the charged HC fluctuate between collapsed globular conformations and a set of expanded structures. The collapsed conformation adopted by the HC in water is slightly destablized in 6 M urea. In contrast, the end-to-end contact is disrupted in the charged HC only in aqueous urea solution. Despite the presence of a large hydrophobic patch, on length scales on the order of approximately 8-10 A "denaturation" (transition to the expanded unfolded state) occurs by a direct interaction of urea with charges on the chain ends. The contiguous patch of hydrophobic moieties leads to "mild dewetting", which becomes more pronounced in the charged HC in 6 M aqueous urea solution. Our simulations establish that the urea denaturation mechanism is most likely electrostatic in origin.