Molecular Dynamics Simulation of the Aggregation Patterns in Aqueous Solutions of Bile Salts at Physiological Conditions

Abstract

Classical molecular dynamics simulations are employed to monitor the aggregation behavior of six bile salts (nonconjugated and glycine- and taurine-conjugated sodium cholate and sodium deoxycholate) with concentration of 10 mM in aqueous solution in the presence of 120 mM NaCl. There are 150 ns trajectories generated to characterize the systems. The largest stable aggregates are analyzed to determine their shape, size, and stabilizing forces. It is found that the aggregation is a hierarchical process and that its kinetics depends both on the number of hydroxyl groups in the steroid part of the molecules and on the type of conjugation. The micelles of all salts are similar in shape-deformed spheres or ellipsoids, which are stabilized by hydrophobic forces, acting between the steroid rings. The differences in the aggregation kinetics of the various conjugates are rationalized by the affinity for hydrogen bond formation for the glycine-modified salts or by the longer time needed to achieve optimum packing for the tauro derivatives. Evidence is provided for the hypothesis from the literature that the entirely hydrophobic core of all aggregates and the enhanced dynamics of the molecules therein should be among the prerequisites for their pronounced solubilization capacity for hydrophobic substances in vivo.