Olsen BN, Schlesinger PH, Ory DS, Baker NA. Side-chain oxysterols: from cells to membranes to molecules. BBA Biomembranes, 1818, 330-336, 2012.
This review discusses the application of cellular biology, molecular biophysics, and computational simulation to understand membrane-mediated mechanisms by which oxysterols regulate cholesterol homeostasis. Side-chain oxysterols, which are produced enzymatically in vivo, are physiological regulators of cholesterol homeostasis and primarily serve as cellular signals for excess cholesterol. These oxysterols regulate cholesterol homeostasis through both transcriptional and non-transcriptional pathways; however, many molecular details of their interactions in these pathways are still not well understood. Cholesterol trafficking provides one mechanism for regulation. The current model of cholesterol trafficking regulation is based on the existence of two distinct cholesterol pools in the membrane: a low and a high availability/activity pool. It is proposed that the low availability/activity pool of cholesterol is integrated into tightly packing phospholipids and relatively inaccessible to water or cellular proteins, while the high availability cholesterol pool is more mobile in the membrane and is present in membranes where the phospholipids are not as compressed. Recent results suggest that that oxysterols may promote cholesterol egress from membranes by shifting cholesterol from the low to the high activity pools. Furthermore, molecular simulations suggest a potential mechanism for oxysterol “activation” of cholesterol through its displacement in the membrane. This review discusses these results as well as several other important interactions between oxysterols and cholesterol in cellular and model lipid membranes.
Olsen BN, Schlesinger PH, Ory DS, Baker NA. 25-Hydroxycholesterol increases the availability of cholesterol in phospholipid membranes. Biophysical Journal, 100, 948-56, 2011.
Side-chain oxysterols are enzymatically generated oxidation products of cholesterol that serve a central role in mediating cholesterol homeostasis. Recent work has shown that side-chain oxysterols, such as 25-hydroxycholesterol (25-HC), alter membrane structure in very different ways from cholesterol, suggesting a possible mechanism for how these oxysterols regulate cholesterol homeostasis. Here we extend our previous work, using molecular dynamics simulations of 25-HC and cholesterol mixtures in 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers to examine interactions between 25-HC and cholesterol in the same bilayer. When added to cholesterol-containing membranes, 25-HC causes larger changes in membrane structure than when added to cholesterol-free membranes, demonstrating interactions between the two sterols. We also find that the presence of 25-HC changes the position, orientation, and solvent accessibility of cholesterol, shifting it into the water interface and therefore its availability to external acceptors. This is consistent with experimental results showing that oxysterols can trigger cholesterol trafficking from the plasma membrane to the endoplasmic reticulum. These interactions provide a potential mechanism for 25-HC-mediated regulation of cholesterol trafficking and homeostasis through direct modulation of cholesterol availability.
Chen Z, Baker NA, Wei GW. Differential geometry based solvation model I: Eulerian formulation, J Comput Phys, 229, 8231-58, 2010.
This paper presents a differential geometry based model for the analysis and computation of the equilibrium property of solvation. Differential geometry theory of surfaces is utilized to define and construct smooth interfaces with good stability and differentiability for use in characterizing the solvent–solute boundaries and in generating continuous dielectric functions across the computational domain. A total free energy functional is constructed to couple polar and nonpolar contributions to the solvation process. Geometric measure theory is employed to rigorously convert a Lagrangian formulation of the surface energy into an Eulerian formulation so as to bring all energy terms into an equal footing. By optimizing the total free energy functional, we derive coupled generalized Poisson–Boltzmann equation (GPBE) and generalized geometric flow equation (GGFE) for the electrostatic potential and the construction of realistic solvent–solute boundaries, respectively. By solving the coupled GPBE and GGFE, we obtain the electrostatic potential, the solvent–solute boundary profile, and the smooth dielectric function, and thereby improve the accuracy and stability of implicit solvation calculations. We also design efficient second-order numerical schemes for the solution of the GPBE and GGFE. Matrix resulted from the discretization of the GPBE is accelerated with appropriate preconditioners. An alternative direct implicit (ADI) scheme is designed to improve the stability of solving the GGFE. Two iterative approaches are designed to solve the coupled system of nonlinear partial differential equations. Extensive numerical experiments are designed to validate the present theoretical model, test computational methods, and optimize numerical algorithms. Example solvation analysis of both small compounds and proteins are carried out to further demonstrate the accuracy, stability, efficiency and robustness of the present new model and numerical approaches. Comparison is given to both experimental and theoretical results in the literature.
Olsen BN, Schlesinger PH, Baker NA. Perturbations of membrane structure by cholesterol and cholesterol derivatives are determined by sterol orientation. J Am Chem Soc, 131, 4854-65, 2009.
Cholesterol is essential for proper function and regulation of eukaryotic membranes, and significant amounts of metabolic energy are dedicated to controlling cellular cholesterol levels. Oxidation products of cholesterol, the oxysterols, are enzymatically produced molecules that play a major role in mediating cholesterol homeostasis through mechanisms which have not yet been fully elucidated. Certain oxysterols are known to have direct effects on membrane permeability and structure; effects that are strikingly different from that of cholesterol. We use molecular dynamics simulations of these oxysterols in 1-palmitoyl 2-oleoyl phosphatidylcholine (POPC) bilayers to explain the structural origins for the differing effects of cholesterol and 25-hydroxycholesterol on bilayer properties. In particular, we demonstrate that the source for these differing perturbations is the much wider range of molecular orientations accessible to 25-hydroxycholesterol when compared to cholesterol. This study shows that direct membrane perturbation by side-chain oxysterols is significant, and suggests that these membrane perturbations may play a role in the oxysterol regulation of cholesterol homeostasis.
Dolinsky TJ, Czodrowski P, Li H, Nielsen JE, Jensen JH, Klebe G, Baker NA. PDB2PQR: Expanding and upgrading automated preparation of biomolecular structures for molecular simulations. Nucleic Acids Res, 35, W522-5, 2007.
Real-world observable physical and chemical characteristics are increasingly being calculated from the 3D structures of biomolecules. Methods for calculating pKa values, binding constants of ligands, changes in protein stability are readily available, but often the limiting step in computational biology is the conversion of PDB structures into formats ready for use with biomolecular simulation software. The continued sophistication and integration of biomolecular simulation methods for systems- and genome-wide studies requires a fast, robust, physically realistic and standardized protocol for preparing macromolecular structures for biophysical algorithms. As described previously, the PDB2PQR web server addresses this need for electrostatic field calculations (Dolinsky et al, NAR, 2004). Here we report the significantly-expanded PDB2PQR that includes the following features: robust standalone command line support, improved pKa estimation via the PROPKA framework, ligand parameterization via PEOE_PB charge methodology, expanded set of force fields and easily-incorporated user-defined parameters via XML input files, and improvement of atom addition and optimization code. These features are available through a new web interface (http://pdb2pqr.sourceforge.net/) which offers users a wide range of options for PDB file conversion, modification, and parameterization.
Song Y, Guallar V, Baker NA. Molecular dynamics simulations of salicylate effects on the micro- and mesoscopic properties of a dipalmitoylphosphatidylcholine bilayer. Biochemistry, 44, 13425-33, 2005.
Salicylate, an amphiphilic molecule and a popular member of non-steroidal anti-inflammatory drug family, is known to affect hearing through reduction of the electromechanical coupling in the outer hair cells of the ear. This reduction of electromotility by salicylate has been widely studied but the molecular mechanism of the phenomenon is still unknown. In this study, we investigated one aspect of salicylate’s action; namely, the perturbation of electrical and mechanical membrane properties by salicylate in the absence of cytoskeletal or membrane-bound motor proteins such as prestin. In particular, we simulated the interaction of salicylate with a dipalmitoylphosphatidylcholine (DPPC) bilayer via atomically-detailed molecular dynamics simulations to observe the effect of salicylate on the microscopic and mesoscopic properties of the bilayer. The results demonstrate that salicylate interacts with the bilayer by associating at the water-DPPC interface in a nearly perpendicular orientation and penetrating more deeply into the bilayer than either sodium or chloride. This association has several affects on the membrane properties. First, binding of salicylate to the membrane displaces chloride from the bilayer-water interface. Second, salicylate influences the electrostatic potential and dielectric properties of the bilayer, with significant changes at the water-lipid bilayer interface. Third, salicylate association results in structural changes including decreased head group area per lipid and increased lipid tail order. However, salicylate does not significantly alter the mechanical properties of the DPPC bilayer; bulk compressibility, area compressibility, and bending modulus were only perturbed by small, statistically-insignificant amounts, by the presence of salicylate. The observations from these simulations are in qualitative agreement with experimental data and support the conclusion that salicylate influences the electrical but not the mechanical properties of DPPC membranes.