Research

Membrane Structure of Chlorosulfolipids

  1. J. Phys. Chem. Lett., 2021 (accepted)


Chlorosulfolipids (CSLs) are major components of flagellar membranes in sea algae. Unlike typical biological lipids that have a hydrophilic head group and hydrophobic hydrocarbon tails, CSLs contain hydrophilic sulfate and chloride groups in the hydrocarbon tail region; this has eluded the prediction of the CSL membrane structure since 1960. In this study, we combine coarse-grained (CG) and atomistic molecular dynamics (MD) simulations to gain significant insights into the membrane structure of Danicalipin A, which is one of the typical CSLs. It is observed from the CG MD that Danicalipin A lipids form a stable monolayer membrane structure wherein the hydrocarbon moieties are sandwiched by hydrophilic sulfate and chloride groups in both the head and tail regions. Based on the mesoscopic membrane structure, we built the corresponding atomistic model to investigate the integrity of the CSL monolayer membrane structure. The monolayer membrane comprising bent lipids shows high thermal stability up to 313 K. The gel-liquid crystalline phase transition is observed around 300 K.

Organic Photovoltaics

  1. J. Mater. Chem. A, 7, 9698, 2019.

  2. Adv. Energy Mater., 9, 1901856, 2019.

  3. Nano Energy, 74, 104883, 2020.

  4. BKCS, 2021, in press.


Long-term stability of organic blends is a key factor for the practical use of organic solar cells (OSCs) in commercial fields. Here, we report the strategic incorporation of non-fullerene small molecules in polymer:fullerene blends to obtain ternary OSCs with improved efficiency and extended lifetimes. Non- fullerene small molecules employed in the polymer:fullerene blend successfully increased the photon- to-current conversion process as an efficient charge cascade acceptor. A combination of theoretical simulations and experimental measurements revealed that aggregation of meta-stable fullerene molecules was significantly alleviated in the ternary blend, thereby preventing an unintentional increase in the threshold for charge transfer during operation. Thus, the ternary OSCs could exhibit highly extended lifetimes with improved morphological stability and better resistance to performance decay under harsh real operational conditions compared to their binary counterparts. Combined with its high efficiency and improved device lifetimes, the high tolerance to the ternary blend thickness offers promise for commercially acceptable ternary OSCs fabricated by a printing process.

Molecular Modeling of Neurofilaments

  1. J. Mol. Biol. 391, 648, 2009.

  2. J. Mol. Biol. 405, 1101, 2011.

  3. J. Neurol. Sci. 307, 132-138, 2011.

  4. J. Neurosci. 32, 6209-6219, 2012.

  5. J. Biol. Phys. 39, 343-362, 2013.

  6. J. Chem. Phys. 138, 015103, 2013.

  7. J. Neurochem. 146, 631-641, 2018


Neurofilaments (NFs) are essential cytoskeletal filaments that impart mechanical integrity to nerve cells. They are assembled from three distinct molecular mass proteins that bind to each other to form a 10-nm-diameter filamentous rod with sidearm extensions. The sidearms are considered to play a critical role in modulating interfilament spacing and axonal caliber. However, the precise mechanism by which NF protrusions regulate axonal diameter remains to be well understood. In particular, the role played by individual NF protrusions in specifying interfilament distances is yet to be established. To gain insight into the role of individual proteins, we investigated the structural organization of NF architecture under different phosphorylation conditions. To this end, a physically motivated sequence- based coarse-grain model of NF brush has been developed based on the three-dimensional architecture of NFs. The model incorporates the charge distribution of sidearms, including charges from the phosphorylation sites corresponding to Lys-Ser-Pro repeat motifs. The model also incorporates the proper grafting of the real NF sidearms based on the stoichiometry of the three subunits. The equilibrium structure of the NF brush is then investigated under different phosphorylation conditions. The phosphoryla- tion of NF modifies the structural organization of sidearms. Upon phosphorylation, a dramatic change involving a transformation from a compact conformation to an extended conformation is found in the heavy NF (NF-H) protein. However, in spite of extensive phosphorylation sites present in the NF-H subunit, the tails of the medium NF subunit are found to be more extended than the NF-H sidearms. This supports the notion that medium NF protrusions are critical in regulating NF spacings and, hence, axonal caliber.

Effects of Nanoparticles on Biological Membranes

  1. J. Phys. Chem. B 110, 5073, 2006.

  2. Biophys. J. 95, 4102, 2008.

  3. BKCS 31, 3195, 2010.

  4. J. Phys. Chem. B 118, 6792, 2014.

  5. BKCS 37, 1076-1085, 2016.

  6. BKCS 39, 516-523, 2018.


We have performed molecular dynamics simulations of atomistic models of C60 clusters and a DMPC bilayer membrane to study the static and dynamic effects of carbon nanoparticles on biological membranes. All four C60-membrane systems were investigated representing dilute and concentrated solutions of C60 residing either inside or outside the membrane. The concentrated C60 molecules in water phase start forming an aggregated cluster. Due to its heavy mass, the cluster tends to adhere on the surface of the bilayer membrane, hindering both translational and rotational diffusion of individual C60. On the other hand, once C60 molecules accumulate inside the membrane, they are well dispersed in the central region of the bilayer membrane. Because of the homogeneous dispersion inside the membrane, each leaflet is pulled away from the center, making the bilayer membrane thicker. This thickening of the membrane provides more room for both translational and rotational motions of C60 inside the membrane compared to that in the water region. As a result, the dynamics of C60 inside the membrane becomes faster with increasing its concentration.

Polymer Dynamics in Porous Medium

  1. J. Phys. Chem. B 110, 5073, 2006.

  2. Biophys. J. 95, 4102, 2008.

  3. BKCS 31, 3195, 2010.

  4. J. Phys. Chem. B 118, 6792, 2014.

  5. BKCS 37, 1076-1085, 2016.

  6. BKCS 39, 516-523, 2018.


The dynamic behavior of hard chains in disordered materials composed of fixed hard spheres is studied using discontinuous molecular dynamics simulations. The matrix induces entanglements in the chain fluid, i.e., for high matrix densities the diffusion coefficient D scales with the chain length N as D ~ N-2. At high matrix densities the rotational relaxation time becomes very large but the translational diffusion is not affected significantly; i.e., the chains display a dynamic heterogeneity reminiscent of probe diffusion in supercooled liquids and glasses. We show that this is because some chains are trapped, and move via a hopping mechanism. There are no signatures of this dynamic heterogeneity in the matrix static structure, however, which is identical to that of a hard-sphere liquid.

DNA Conformation in Nanochannels

  1. Lap on a Chip 11, 1721, 2011.

  2. J. Chem. Phys. 136, 095101, 2012.

  3. ACS Macro Lett. 3, 926-930, 2014.


We have performed canonical ensemble Monte Carlo simulations of a primitive DNA model to study the conformation of 2.56 ∼ 21.8 μm long DNA molecules confined in nanochannels at various ionic concentrations with the comparison of our previous experimental findings. In the model, the DNA molecule is represented as a chain of charged hard spheres connected by fixed bond length and the nanochannels as planar hard walls. System potentials consist of explicit electrostatic potential along with short-ranged hard-sphere and angle potentials. Our primitive model system provides valuable insight into the DNA conformation, which cannot be easily obtained from experiments or theories. First, the visualization and statistical analysis of DNA molecules in various channel dimensions and ionic strengths verified the formation of locally coiled structures such as backfolding or hairpin and their significance even in highly stretched states. Although the folding events mostly occur within the region of ∼ 0.5 μm from both chain ends, significant portion of the events still take place in the middle region. Second, our study also showed that two controlling factors such as channel dimension and ionic strength widely used in stretching DNA molecules have different influence on the local DNA structure. Ionic strength changes local correlation between neighboring monomers by controlling the strength of electrostatic interaction (and thus the persistence length of DNA), which leads to more coiled local conformation. On the other hand, channel dimension controls the overall stretch by applying the geometric constraint to the non-local DNA conformation instead of directly affecting local correlation. Third, the molecular weight dependence of DNA stretch was observed especially in low stretch regime, which is mainly due to the fact that low stretch modes observed in short DNA molecules are not readily accessible to much longer DNA molecules, resulting in the increase in the stretch of longer DNA molecules.