A detailed model for the structure and dynamics of the interior of the lipid bilayer in the liquid crystal phase is presented. The model includes two classes of motion: (i) the internal dynamics of the chains, determined from Brownian dynamics simulations with a continuous version of the Marcelja mean-field potential, and (ii) noncollective reorientation (axial rotation and wobble) of the entire molecule, introduced by a cone model. The basic unit of the model is a single lipid chain with field parameters adjusted to fit the 2H order parameters and the frequency-dependent 13C NMR T1 relaxation times of dipalmitoyl phosphatidylchollne bilayers. The chain configurations obtained from the trajectory are used to construct a representation of the bilayer. The resulting lipid assembly is consistent with NMR, neutron diffraction, surface area, and density data. It indicates that a high degree of chain disorder and entanglement exists in biological membranes.Lipid bilayers have two distinct phases in the physiological temperature range. The lower temperature Lp phase is highly ordered, with the hydrocarbon tails of the lipids primarily in the extended all-trans state (1, 2). At higher temperature, the bilayer enters the partially ordered liquid crystal (L,,) phase observed for biologically active membranes (1, 2). On going from Lp3 to La, the volume and fluidity of the bilayer increase, the thickness decreases, and the individual lipid undergoes a wide range of motions including: gauche-trans isomerizations, axial rotation, collective and noncollective tilting, out-of-plane deformations, and lateral diffusion (3)(4)(5)(6)(7)(8). The precise time scales, amplitudes, and degree of coupling of these motions remain uncertain both experimentally and theoretically. As a result, our knowledge of the underlying chain distribution is limited, and our understanding of the manner in which the bilayer adjusts to the presence of proteins and other membrane components (9) or to the high local curvature characteristic of membrane fusion (10) is incomplete.Much of the recent theoretical work on lipid structure has been based on either simplified approaches, such as the lattice (11) and mean-field (12-14) approximations, or on computationally intensive molecular dynamics (15, 16) and Monte Carlo (17) simulation methods. In lattice theories, the configurations for a system composed of a set of chain molecules, each consisting of connected segments, are generated on a lattice of specified symmetry. However, for large systems such as lipid bilayers, computer limitations require either a relatively coarse lattice or a short chain length. Consequently, the chains are idealized, and it is difficult to use the results of these calculations to interpret details of the real chain distributions. Very detailed descriptions of the bilayer come from molecular dynamics simulations with all-atom potentials (15). Unfortunately, many significant lipid motions are slow relative to the molecular dynamics subnanosecond time scale. The trans-gauche ...