Biomedical subjects

Rakwoo Chang

Publications and source records attributed to Rakwoo Chang.

5 recordsLinked to original sources

PubMed2006

Insights into the effect of combustion-generated carbon nanoparticles on biological membranes: a computer simulation study.

Classical molecular dynamics simulations of atomistic models of combustion-generated carbon nanoparticles and lipid bilayers have been performed to explore their possible structural, dynamical, and thermodynamic effects on biological membranes. The DREIDING generic force field is used for the carbonaceous nanoparticles of different morphologies, as produced from combustion sources, and the united atom model was employed for the dimyristoylphosphatidylcholine (DMPC) bilayer. It is observed that particle shape and structure have significant effects on solvation, mobility, adsorption, and permeation behavior of the particles. While combustion-generated carbon nanoparticles with an aspect ratio close to unity prefer to stay near the membrane center, precursors with other shapes mostly reside within the hydrocarbon tail region of the membrane. Carbon nanoparticles are not trapped in a local region even inside the membranes but move freely with a speed depending on their molecular weight. The adsorption of the particles on the surface of the biological membrane is comparable to thermal fluctuations because the weak segregation effect by water molecules is the main driving force to the adsorption behavior. The bigger the precursors are, the stronger they are bound to the membrane surface. The presence of combustion-generated nanoparticles inside the membrane perturbs local lipid density by pushing the neighboring lipid molecules away from the nanoparticles. This, coupled with thermal fluctuations, can induce an instantaneous membrane pore to allow water protrusion. From the umbrella sampling method, the potential of mean force for the permeation of carbona nanoparticles into the bilayer was also obtained. Surprisingly, elongated particles have a free energy barrier an order of magnitude smaller compared with more round ones. In addition, the round carbon nanoparticles showed strong hysteresis due to the local trapping of water molecules. Although the carbon soot precursors studied in this work are not the well-known carbon nanoparticles such as fullerenes or carbon nanotubes, the qualitative features of this study may be applicable to them as well.

Carbon↗

PubMed2006

Dynamics of chain molecules in disordered materials.

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 is approximately 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.

Journal Article↗

PubMed2005

Multiscale coupling of mesoscopic- and atomistic-level lipid bilayer simulations.

A multiscale method is presented to bridge between the atomistic and mesoscopic membrane systems. The atomistic model in this case is the united atom dimyristoylphosphatidylcholine membrane system, although the method is completely general. Atomistic molecular dynamics provides the expansion modulus which is used to parametrize a mesoscopic elastic membrane model. The resulting elastic membrane model, including explicit mesoscopic solvent, shows appropriate static and dynamic undulation behaviors. Large membranes of approximately 100 nm in length can then be easily simulated using the mesoscopic membrane system. The critical feedback from the mesoscopic system back down to the atomistic-scale system is accomplished by bridging the stress (or surface tension) of a small region in the mesoscopic membrane to the corresponding atomistic membrane system. Because of long length-scale modes of membranes such as undulation and buckling, the local tension responds differently from the frame tension, when subjected to external perturbations. The effect of these membrane modes is shown for the stress response of a local membrane region and therefore the atomistic membrane system. In addition, certain equilibrium static and dynamic properties of stand-alone and multiscale coupled systems are presented for several different membrane sizes. Although static properties such as two-dimensional pair-correlation function and order parameters show no noticeable discrepancy for the different systems, lipid self-diffusion and the rotational relaxation of lipid dipoles have a strong dependence on the membrane size (or long-wavelength membrane motions), which is properly modeled by the present multiscale method.

Algorithms↗

PubMed2004

Diffusion of hard sphere fluids in disordered media: a molecular dynamics simulation study.

Molecular dynamic simulations are reported for the static and dynamic properties of hard sphere fluids in matrices (or media) composed of quenched hard spheres. The effect of fluid and matrix density, matrix structure, and fluid to matrix sphere size ratio on the static and dynamic properties is studied using discontinuous molecular dynamics. The matrix density has a stronger effect on the self-diffusion coefficient than the fluid density, especially at high matrix densities where the geometric constraints due to the quenched spheres are significant. When the ratio of the size of the fluid spheres to that of the matrix spheres is equal to or greater than one, the diffusion increases as the fluid density is increased, at constant total volume fraction. This trend is however reversed if the ratio is smaller than one. Different methods of generating the matrix have a very strong effect on the dynamic properties even though the static correlations are similar. An analysis of the single-chain structure factor of the particle trajectories shows a change in the particle diffusive behavior at different time scales, suggestive of a hopping mechanism, although normal diffusion is recovered at long times. At high matrix densities, there is considerable heterogeneity in the diffusion of the fluid particles. The simulations demonstrate that the correlations in the matrix play a significant role on the diffusion of fluid spheres. For example, the diffusion constant in matrices constructed by different methods can be an order of magnitude different even though the pair correlation functions are almost identical.

Journal Article↗

PubMed2002

A Monte Carlo study of the self-assembly of bacteriorhodopsin.

Bacteriorhodopsin (BR) and specific lipid molecules self-assemble into a quasi two-dimensional lattice structure known as the purple membrane (PM). In the PM, BR molecules exist in a trimeric form with lipid molecules present in the space enclosed by each trimeric unit and in the inter-trimer space. These trimeric units, which have a roughly circular cross-section, are arranged in hexagonal patterns with long-ranged crystalline order. In this work, we investigate the self-assembly of BR in the PM via Monte Carlo simulations of a two-dimensional model of the membrane and proteins. The protein molecules are modeled as 120 degrees sectors of a circle and the lipid molecules enter into the model through effective protein-protein interactions. The sectors cannot overlap with each other, and in addition to this excluded volume interaction there are site-site attractive interactions between specific points of the proteins to mimic interactions between helices on the proteins and lipid-induced interactions. At low values of the attractive well depth, the proteins are found in the monomeric form at all concentrations. At moderate and high values of the attractive well depth, trimers are formed as the concentration increases, and with a further increase in concentration the trimers organize into a hexagonal lattice. The interactions between the proteins and those induced by the intra-trimer lipids play an equally important role in the formation of trimers and the lattice. The lipids in the inter-trimer space cause the trimers to orient in a specific direction in the hexagonal crystal lattice.

Algorithms↗