Solvation of cholesterol in different solvents: a molecular dynamics simulation study

Molecular dynamics simulations of the interfacial behaviors and photo-oxidation of phytosterol under different emulsion oil content

Phytosterol, recognized for its health benefits, is predominantly extracted from plants and exhibits significantly reduced stability under varying light conditions. Their photooxidation is significantly influenced by emulsion interfaces. This study examined the mechanism of interface structure on phytosterol photooxidation with unparalleled molecular precision, utilizing molecular dynamics simulations and experimental procedures. Hydrogen bonding between the hydroxyl group at the C3 position of phytosterols and water molecules, coupled with van der Waals forces between the hydrophobic regions and the oil phase, induced phytosterol molecules to disperse toward the interface. The elevated polarity of the oil phase, specifically in tributyrin, facilitated the permeation of water molecules into the oil phase. This was achieved by diminishing the emulsion's interfacial tension, thereby fostering the development of more interface or micelles, and accelerating the photooxidation process of phytosterols. These simulations unraveled that the preponderance of phytosterol distribution is localized and oxidized at the oil-water interface.

Structural descriptions of ligand interactions to RNA quadruplexes folded from the non-coding region of pseudorabies virus

To rationalise the binding of specific ligands to RNA-quadruplex we investigated several naphthalene diimide ligands that interact with the non-coding region of Pseudorabies virus (PRV). Herein we report on the x-ray structure of the naphthalene diimide ND11 with an RNA G-quadruplex putative forming sequence from rPRV. Consistent with previously observed rPRV sequence it assembles into a bimolecular RNA G-quadruplex consisting of a pair of two tetrads stacked 3' to 5'. We observe that ND11 interacts by binding on both the externally available 5' and 3' quartets. The CUC (loop 1) is structurally altered to enhance the 5' mode of interaction. These loop residues are shifted significantly to generate a new ligand binding pocket whereas the terminal A14 residue is lifted away from the RNA G-quadruplex tetrad plane to be restacked above the bound ND11 ligand NDI core. CD analysis of this family of NDI ligands shows consistency in the spectra between the different ligands in the presence of the rPRV RNA G-quadruplex motif, reflecting a common folded topology and mode of ligand interaction. FRET melt assay confirms the strong stabilising properties of the tetrasubstituted NDI compounds and the contributions length of the substituted groups have on melt temperatures.

Sterol-activated amyloid beta fibril formation

The metabolic processes that link Alzheimer's disease (AD) to elevated cholesterol levels in the brain are not fully defined. Amyloid beta (Aβ) plaque accumulation is believed to begin decades prior to symptoms and to contribute significantly to the disease. Cholesterol and its metabolites accelerate plaque formation through as-yet-undefined mechanisms. Here, the mechanism of cholesterol (CH) and cholesterol 3-sulfate (CS) induced acceleration of Aβ42 fibril formation is examined in quantitative ligand binding, Aβ42 fibril polymerization, and molecular dynamics studies. Equilibrium and pre-steady-state binding studies reveal that monomeric Aβ42•ligand complexes form and dissociate rapidly relative to oligomerization, that the ligand/peptide stoichiometry is 1-to-1, and that the peptide is likely saturated in vivo. Analysis of Aβ42 polymerization progress curves demonstrates that ligands accelerate polymer synthesis by catalyzing the conversion of peptide monomers into dimers that nucleate the polymerization reaction. Nucleation is accelerated ∼49-fold by CH, and ∼13,000-fold by CS - a minor CH metabolite. Polymerization kinetic models predict that at presumed disease-relevant CS and CH concentrations, approximately half of the polymerization nuclei will contain CS, small oligomers of neurotoxic dimensions (∼12-mers) will contain substantial CS, and fibril-formation lag times will decrease 13-fold relative to unliganded Aβ42. Molecular dynamics models, which quantitatively predict all experimental findings, indicate that the acceleration mechanism is rooted in ligand-induced stabilization of the peptide in non-helical conformations that readily form polymerization nuclei.

Monosaccharide induced temporal delay in cholesterol self-aggregation

Self-assembly of cholesterol (CHL) is infamous for its diverse deleterious effects on human health. Clinical research over several decades indicates that a diet rich in CHL typically leads to arterial plaques, cataracts and gall stones among others. Carbohydrates like the β-glucans efficiently lower serum CHL, possibly by inhibiting CHL absorption in the digestive tract. Using molecular dynamics simulations, we explore how β-D-glucose (BGLC), the building block of β-glucans, interferes with CHL aggregation. BGLC slows down CHL diffusion and disrupts the formation of the robust hydrophobic CHL assembly. Estimation of the translational entropy of the CHL molecules shows the extent of retardation induced by BGLC. Coordination numbers obtained from the adjacency matrix and collective variable analysis of the packing of the CHL molecules in presence of BGLC show the time evolution of CHL aggregation. In presence of BGLC, small isolated CHL islands form, consolidate and disintegrate over time as compared to the blank CHL system. The predominance of smaller CHL clusters is an effect of the significant retardation of the translational motion of CHL molecules induced by BGLC.Communicated by Ramaswamy H. Sarma.

Direct Measurement of Intermolecular Mechanical Force for Nonspecific Interactions between Small Molecules

Mechanical force can evaluate intramolecular interactions in macromolecules. Because of the rapid motion of small molecules, it is extremely challenging to measure mechanical forces of nonspecific intermolecular interactions. Here, we used optical tweezers to directly examine the intermolecular mechanical force (IMMF) of nonspecific interactions between two cholesterols. We found that IMMFs of dimeric cholesterol complexes were dependent on the orientation of the interaction. The surprisingly high IMMF in cholesterol dimers (∼30 pN) is comparable to the mechanical stability of DNA secondary structures. Using Hess-like cycles, we quantified that changes in free energy of solubilizing cholesterol (ΔGsolubility) by β-cyclodextrin (βCD) and methylated βCD (Me-βCD) were as low as -16 and -27 kcal/mol, respectively. Compared to the ΔGsolubility of cholesterols in water (5.1 kcal/mol), these values indicated that cyclodextrins can easily solubilize cholesterols. Our results demonstrated that the IMMF can serve as a generic and multipurpose variable to dissect nonspecific intermolecular interactions among small molecules into orientational components.

On the Consistency of the Exfoliation Free Energy of Graphenes by Molecular Simulations

Monolayer graphene is now produced at significant yields, by liquid phase exfoliation of graphites in solvents. This has increased the interest in molecular simulation studies to give new insights in the field. We use decoupling simulations to compute the exfoliation free energy of graphenes in a liquid environment. Starting from a bilayer graphene configuration, we decouple the Van der Waals interactions of a graphene monolayer in the presence of saline water. Then, we introduce the monolayer back into water by coupling its interactions with water molecules and ions. A different approach to compute the graphene exfoliation free energy is to use umbrella sampling. We apply umbrella sampling after pulling the graphene monolayer on the shear direction up to a distance from a bilayer. We show that the decoupling and umbrella methods give highly consistent free energy results for three bilayer graphene samples with different size. This strongly suggests that the systems in both methods remain closely in equilibrium as we move between the states before and after the exfoliation. Therefore, the amount of nonequilibrium work needed to peel the two layers apart is minimized efficiently.

Photophysical and Structural Characterization of Intrinsically Fluorescent Sterol Aggregates

Self-association of cholesterol into aggregates and crystals is a hallmark of developing atherosclerosis. Intrinsically fluorescent sterols, such as dehydroergosterol (DHE), can be used to study sterol aggregation by fluorescence spectroscopy and microscopy, but a thorough understanding of DHE's photophysical and structural properties in the aggregated state is missing. Here, we show that DHE forms submicron fluorescent aggregates when evaporated from an ethanol solution. Using atomic force microscopy, we find that DHE, like cholesterol, forms compact oblate-shape aggregates of <100 nm in diameter. DHE's fluorescence is lowered in the aggregate compared to the monomeric form, and characteristic spectral changes accompany the aggregation process. Electronic structure calculations of DHE dimers in water indicate that Frenkel-type exciton coupling contributes to the lowered DHE fluorescence in the aggregates. Using molecular dynamics (MD) simulations, we show that DHE forms compact aggregates on the nanosecond scale and with strong intermolecular attraction, in which a broad range of orientations, and therefore electronic couplings, will take place. Tight packing of DHE in aggregates also lowers the apparent absorption cross section, further reducing the molecular brightness of the aggregates. Our results pave the way for systematic solubility studies of intrinsically fluorescent analogues of biologically relevant sterols.

Combined Use of Structure Analysis, Studies of Molecular Association in Solution, and Molecular Modelling to Understand the Different Propensities of Dihydroxybenzoic Acids to Form Solid Phases

The arrangement of hydroxyl groups in the benzene ring has a significant effect on the propensity of dihydroxybenzoic acids (diOHBAs) to form different solid phases when crystallized from solution. All six diOHBAs were categorized into distinctive groups according to the solid phases obtained when crystallized from selected solvents. A combined study using crystal structure and molecule electrostatic potential surface analysis, as well as an exploration of molecular association in solution using spectroscopic methods and molecular dynamics simulations were used to determine the possible mechanism of how the location of the phenolic hydroxyl groups affect the diversity of solid phases formed by the diOHBAs. The crystal structure analysis showed that classical carboxylic acid homodimers and ring-like hydrogen bond motifs consisting of six diOHBA molecules are prominently present in almost all analyzed crystal structures. Both experimental spectroscopic investigations and molecular dynamics simulations indicated that the extent of intramolecular bonding between carboxyl and hydroxyl groups in solution has the most significant impact on the solid phases formed by the diOHBAs. Additionally, the extent of hydrogen bonding with solvent molecules and the mean lifetime of solute-solvent associates formed by diOHBAs and 2-propanol were also investigated.

Interaction of lipid-free apolipoprotein A-I with cholesterol revealed by molecular modeling

We report the modeling of the interaction of differently self-associated lipid-free apoA-I with cholesterol monomer and tail-to-tail (TT) or face-to-face (FF) cholesterol dimer. Cholesterol dimerization is exploited to reconcile the existing experimental data on cholesterol binding to apoA-I with extremely low critical micelle concentration of cholesterol. Two crystal structures of 1-43 N-truncated apolipoprotein Δ(1-43)A-I tetramer (PDB ID: 1AV1, structure B), 185-243 C-truncated apolipoprotein Δ(185-243)A-I dimer (PDB ID: 3R2P, structure M) were analyzed. Cholesterol monomers bind to multiple binding sites in apoA-I monomer, dimer and tetramer with low, moderate and high energy (-10 to -28 kJ/mol with Schrödinger package), still insufficient to overcome the thermodynamic restriction by cholesterol micellization (-52.8 kJ/mol). The binding sites partially coincide with the putative cholesterol-binding motifs. However, apoA-I monomer and dimer existing in structure B, that contain nonoverlapping and non-interacting pairs of binding sites with high affinity for TT and FF cholesterol dimers, can bind in common 14 cholesterol molecules that correspond to existing values. ApoA-I monomer and dimer in structure M can bind in common 6 cholesterol molecules. The values of respective total energy of cholesterol binding up to 64.5 and 67.0 kJ/mol for both B and M structures exceed the free energy of cholesterol micellization. We hypothesize that cholesterol dimers may simultaneously interact with extracellular monomer and dimer of lipid-free apoA-I, that accumulate at acid pH in atheroma. The thermodynamically allowed apolipoprotein-cholesterol interaction outside the macrophage may represent a new mechanism of cholesterol transport by apoA-I from atheroma, in addition to ABCA1-mediated cholesterol efflux.