• Title/Summary/Keyword: Molecular Dynamics (MD) simulations

Search Result 123, Processing Time 0.017 seconds

Analysis of Cis- Trans Photoisomerization Mechanism of Rhodopsin Based on the Tertiary Structure of Rhodopsin

  • Yamada, Atsushi;Yamato, Takahisa;Kakitani, Toshiaki;Yamamoto, Shigeyoshi
    • Journal of Photoscience
    • /
    • v.9 no.2
    • /
    • pp.51-54
    • /
    • 2002
  • We propose a novel mechanism (Twist Sharing Mechanism) for the cis-trans photoisomerization of rhodopsin, based on the molecular dynamics (MD) simulation study. New things devised in our simulations are (1) the adoption of Mt. Fuji potentials in the excited state for twisting of the three bonds C9=C10, C11=C12 and C13=14 which are modeled using the detailed ab initio quantum chemical calculations and (2) to use the rhodopsin structure which was resolved recently by the X-ray crystallographic study. As a result, we found the followings: Due to the intramolecular steric hindrance between 20-methyl and 10-H in the retinal chromophore, the C12-C13 and C10-C11 bonds are considerably twisted counterclockwise in rhodopsin, allowing only counterclockwise rotation of the C11 =C12 in the excited state. The movement of 19-methyl in rhodopsin is blocked by the surrounding three amino acids, Thr 118, Met 207 and Tyr 268, prohibiting the rotation of C9=C10. As a result only all-trans form of the chromophore is obtainable as a photoproduct. At the 90$^{\circ}$ twisting of C11=C12 in the course of photoisomerization, twisting energies of the other bonds amount to about 20 kcal/mol. If the transition state for the thermal isomerization is assumed to be similar to this structure, the activation energy for the thermal isomerization around C11=C12'in rhodopsin is elevated by about 20 kcal/mol and the thermal isomerization rate is decelerated by 10$\^$-14/ times than that of the retinal chromophore in solution, protecting photosignal from the thermal noise.

  • PDF

A Study on the Sequential Multiscale Homogenization Method to Predict the Thermal Conductivity of Polymer Nanocomposites with Kapitza Thermal Resistance (Kapitza 열저항이 존재하는 나노복합재의 열전도 특성 예측을 위한 순차적 멀티스케일 균질화 해석기법에 관한 연구)

  • Shin, Hyunseong;Yang, Seunghwa;Yu, Suyoung;Chang, Seongmin;Cho, Maenghyo
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.25 no.4
    • /
    • pp.315-321
    • /
    • 2012
  • In this study, a sequential multiscale homogenization method to characterize the effective thermal conductivity of nano particulate polymer nanocomposites is proposed through a molecular dynamics(MD) simulations and a finite element-based homogenization method. The thermal conductivity of the nanocomposites embedding different-sized nanoparticles at a fixed volume fraction of 5.8% are obtained from MD simulations. Due to the Kapitza thermal resistance, the thermal conductivity of the nanocomposites decreases as the size of the embedded nanoparticle decreases. In order to describe the nanoparticle size effect using the homogenization method with accuracy, the Kapitza interface in which the temperature discontinuity condition appears and the effective interphase zone formed by highly densified matrix polymer are modeled as independent phases that constitutes the nanocomposites microstructure, thus, the overall nanocomposites domain is modeled as a four-phase structure consists of the nanoparticle, Kapitza interface, effective interphase, and polymer matrix. The thermal conductivity of the effective interphase is inversely predicted from the thermal conductivity of the nanocomposites through the multiscale homogenization method, then, exponentially fitted to a function of the particle radius. Using the multiscale homogenization method, the thermal conductivities of the nanocomposites at various particle radii and volume fractions are obtained, and parametric studies are conducted to examine the effect of the effective interphase on the overall thermal conductivity of the nanocomposites.

Chain Length Effect on the Configurational Properties of an n-Alkane Chain in Solution

  • Jeon, Seung-Ho;Ree, Tai-Kyue;Oh, In-Joon
    • Bulletin of the Korean Chemical Society
    • /
    • v.7 no.5
    • /
    • pp.367-371
    • /
    • 1986
  • Dynamic and equilibrium properties of n-alkane chains immersed in solvent molecules have been investigated by a molecular dynamics method. The n-alkane chain is assumed to be a chain of elements (CH$_2$) interconnected by bonds having a fixed bond length and bond angle, but each bond of the chain is allowed to execute hindered internal rotation. We studied the effect of the number of the chain elements (N$_c$ = 10, 15 and 20) on the equilibrium properties of the system, e.g., the pair correlation functions between a chain element and solvent molecules, g$_{cs}$(r), and between the chain elements, g$_{cc}$(r), and the configurational properties such as the mean-square end-to-end distance < R$^2$ >, the mean-square radius of gyration < S$^2$ >, and the eigenvalues of the moment-of-inertia tensor < S$_i^2$ > / < S$^2$ > (i = 1, 2 and 3). We also studied the dynamic properties of the system, e.g., the autocorrelation function C(A;t) where A = R$^2$(t), = S$^2$(t), or = ${\vec{V}}(t)({\vec{V}}$ = velocity of the center of mass), and the diffusion coefficient D. The g$_{cs}$(r)'s are almost equal irrespective of the change of Nc while g$_{cc}$(r) becomes larger as N$_c$ increases; The MD computed configurational properties < R$^2$2 > and < S$^2$ > were found to be a little different from the values calculated from the statistical equations of < R$^2$ > and < S$^2$ >, it may be due to the fact that our model for the MD simulations includes a long-range volume effect. From the < S$_i^2$ > / < S$^2$ >, it is found that the chain molecule has a nearly spherical shape irrespective of the variation of N$_c$. For the dynamic properties we found that the C(R$^2$;t) and C(S$^2$;t) of lower N$_c$ decay faster than those of higher N$_c$, while the C($\vec V$;t) of the center of mass in the chain is weakly dependent on the N$_c$. The center of mass diffusion coefficient D$_c$ decreases as N$_c$ increases while the end point diffusion coefficient D$_e$ is nearly equal irrespective of the change of N$_c$.