Biomedical subjects

Wataru Shinoda

Publications and source records attributed to Wataru Shinoda.

9 recordsLinked to original sources

PubMed2006

Molecular dynamics simulations of self-organized polyicosahedral Si nanowire.

A novel polyicosahedral nanowire is spontaneously formed in a series of annealing molecular dynamics simulations of liquid Si inside a nanopore of 1.36 nm in diameter. The polyicosahedral Si nanowire is stable even in a vacuum up to about 77% of the melting temperature of bulk Si. Our structural energy calculations reveal that the polyicosahedral nanowire is energetically advantageous over the pentagonal one for a wire whose diameter is less than 6.02 nm, though the latter has been recently proposed as the lowest energy wire. These results suggest the possibility of the formation of a new stable polyicosahedral Si nanowire.

Journal Article↗

PubMed2005

Calculation of heat capacities of light and heavy water by path-integral molecular dynamics.

As an application of atomistic simulation methods to heat capacities, path-integral molecular dynamics has been used to calculate the constant-volume heat capacities of light and heavy water in the gas, liquid, and solid phases. While the classical simulation based on conventional molecular dynamics has estimated the heat capacities too high, the quantum simulation based on path-integral molecular dynamics has given reasonable results based on the simple point-charge/flexible potential model. The calculated heat capacities (divided by the Boltzmann constant) in the quantum simulation are 3.1 in the vapor H2O at 300 K, 6.9 in the liquid H2O at 300 K, and 4.1 in the ice Ih H2O at 250 K, respectively, which are comparable to the experimental data of 3.04, 8.9, and 4.1, respectively. The quantum simulation also reproduces the isotope effect. The heat capacity in the liquid D2O has been calculated to be 10% higher than that of H2O, while it is 13% higher in the experiment. The results demonstrate that the path-integral simulation is a promising approach to quantitatively evaluate the heat capacities for molecular systems, taking account of quantum-mechanical vibrations as well as strongly anharmonic motions.

Journal Article↗

PubMed2005

Molecular dynamics study of the methanol effect on the membrane morphology of perfluorosulfonic ionomers.

The membranes of a perfluorosulfonic acid polymer swollen in 10-80 wt % methanol solution were investigated to elucidate the methanol effect on their morphologies, such as size of the solvent cluster, solvent location, and polymer structure, by using isothermal-isobaric molecular dynamics simulations. In higher methanol concentrations, we found less-spherical solvent aggregation and a more spread polymer structure because of the ampholytic nature of methanol. The partial radial distribution functions between solvent oxygen and fluorocarbons, which are composed of the main chain, clearly show that methanol is located closer to the polymer matrix than water. On the other hand, water is preferentially located in the vicinity of an acidic headgroup, SO(3)(-), compared with methanol, although both have similar attractive interaction energies to the acidic group. Furthermore, we discussed solvent dynamics and hydrogen bonding between sulfonic oxygen and solvent O-H groups.

Journal Article↗

PubMed2005

Molecular dynamics study of bipolar tetraether lipid membranes.

Membranes composed of bipolar tetraether lipids have been studied by a series of 25-ns molecular dynamics simulations to understand the microscopic structure and dynamics as well as membrane area elasticity. By comparing macrocyclic and acyclic tetraether and diether archaeal lipids, the effect of tail linkage of the two phytanyl-chained lipids on the membrane properties is elucidated. Tetraether lipids show smaller molecular area and lateral mobility. For the latter, calculated diffusion coefficients are indeed one order-of-magnitude smaller than that of the diether lipid. These two tetraether membranes are alike in many physical properties except for membrane area elasticity. The macrocyclic tetraether membrane shows a higher elastic area expansion modulus than its acyclic counterpart by a factor of three. Free energy profiles of a water molecule crossing the membranes show no major difference in barrier height; however, a significant difference is observed near the membrane center due to the lack of the slip-plane in tetraether membranes.

Archaea↗

PubMed2005

Quantum simulation of the heat capacity of water.

Path integral molecular dynamics simulations have been carried out to evaluate heat capacities of water in gas, liquid, and solid phases. For convergence, a 100-ps simulation run with a large path integral variable (P approximately 128) was required even by using the double-centroid virial heat capacity estimator. For all states, the quantum corrections to the heat capacities are significantly large. The calculated heat capacities for vapor and ice I(h) agreed excellently with experimental data, while that for liquid was less than the experimental value by approximately 20% due to the limit of the SPC/F potential model.

Journal Article↗

PubMed2005

Spatial confinement effect on the atomic structure of solid argon.

Molecules confined in nanopores show unusual behavior not seen in bulk systems. The present paper reports on molecular dynamics simulations of unusual freezing behavior in confined Ar. Similar to bulk Ar, liquid Ar confined in pores with a diameter D>15sigma (5.1 nm), where sigma is the diameter of the Ar atom, crystallizes when the cooling rate is lower than a critical value (Qc). We also find that the spatial confinement does not have significant influence on Qc when D>15sigma (5.1 nm). In the pore of 10sigma (3.4 nm) in diameter, on the other hand, the behavior is dramatically changed. Crystalline Ar does not appear inside the pore even when the system is cooled at a rate lower than the Qc in the bulk system by over two orders of magnitude. Instead, amorphous Ar characterized by local icosahedral configurations is formed in the pore. We further find that, even when crystalline Ar is formed outside the pore, it does not grow deeply into the pore. This supports that the amorphous Ar is actually the most stable phase in the pore. It is well known that Ar is a poor glass former. Our finding that even such an amorphous Ar is the most stable in the pore suggests that, in any system, it is possible to prepare amorphous structure selectively by using nano-molds.

Journal Article↗

PubMed2005

Molecular dynamics simulation of swollen membrane of perfluorinated ionomer.

Molecular dynamics simulations of the swollen membrane of perfluorinated ionomer, which is composed of poly(tetrafluoroethylene) backbones and perfluosulfonic pendant side chains, have been undertaken to analyze the static and dynamic properties of the water and the side chain in the membrane. The calculations were carried out for four different water contents, 5, 10, 20 and 40 wt %, at 358.15 K and 0.1 MPa. The results are summarized as follows: (1) The sulfonic acid is the unique site to which water molecules can bind, and the other sites in the pendant side chain have no bound water even at high water concentration. (2) Sulfonic acids aggregate in the short range within 4.6-7.7 A despite the electrostatic repulsion between them. In such aggregates, a water molecule bridges two sulfonic acids. (3) Pendant side chains prefer to orient perpendicular to the hydrophilic/hydrophobic interface, and long-range correlation of side chain orientations is observed at 20 and 40 wt % water uptake membranes. (4) In a low water uptake membrane, the dynamics of water is substantially restricted due to strong attractive interactions with acidic sites. In contrast, at high water content, even the water locating near the sulfonic acid is relatively mobile. The short residence time of the bound water reveals that such water can frequently exchange position with relatively free water, which locates in the center of water cluster, in highly swollen membranes.

Journal Article↗

PubMed2004

Comparative molecular dynamics study of ether- and ester-linked phospholipid bilayers.

The lipid membranes found in archaea have high bilayer stability and low permeability. The molecular structure of their constituent lipids is characterized by ether-linked, branched hydrophobic chains, whereas the conventional lipids obtained from eukaryotic or eubacterial sources have ester linked straight chains. In order to elucidate the influence of the ether linkage, instead of an ester one, on the physical properties of the lipid bilayers, we have carried out comparative 10 ns molecular dynamics simulations of diphytanyl phosphatidylcholine (ether-DPhPC) and diphytanoyl phosphatidylcholine (ester-DPhPC) bilayers in water, respectively. We analyze bilayer structures, hydration of the lipids, membrane dipole potentials, and free energy profiles of water and oxygen across the bilayers. We observe that the membrane dipole potential for the ether-DPhPC bilayer, which arises mainly from the ether linkage, is about half of that of the ester-DPhPC. The calculated free energy barrier for a water molecule in the ether-DPhPC bilayer system is slightly higher than that in the ester-DPhPC counterpart, which is in accord with experimental data.

Journal Article↗

PubMed2003

Rigid-body dynamics in the isothermal-isobaric ensemble: a test on the accuracy and computational efficiency.

We have developed a time-reversible rigid-body (rRB) molecular dynamics algorithm in the isothermal-isobaric (NPT) ensemble. The algorithm is an extension of rigid-body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non-Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state-of-the-art constraint-dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time-step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS-MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL.

Journal Article↗