Molecular Dynamics Simulation of the Bio-Adhesion in Molecular Motors

Molecular dynamics (MD) simulation is employed to study the bio-adhesion in F1 ATP molecular motor. Histidine-peptide is widely used as linkage in micro systems because of its strong binding strength to metals. This paper focuses on the adhesion between a synthetic peptide containing 6xHis-tag (Gly-Gly-Lys-Gly-Gly-Lys-Gly-Gly-His-His-His-His-His-His) and metal substrate, which is used to define the position of the F1 ATP molecular motor on the metal substrate. It is shown that the binding strength between histidine and nickel substrate is the strongest, while that of copper is smaller and that of gold substrate is the smallest. From the result of simulation, we find that the stability of adhesion between histidine and the metal substate result of the ringed structure in histidine.

Molecular dynamics studies of brittle fracture of SiO2

By applying for molecular dynamics (MD) simulation and Griffith fracture criterion, the brittle behavior of crack extension of mode I type is investigated. The critical stress intensity factor (SIF)K-Ic(MD) of crack extension is calculated, and the evolution of atoms near crack tip is observed. It is found that K-Ic(MD) is in good agreement with the Griffith ftacture criterion K-Ic(Griffith).

Size-dependent elastic properties of Ni nanofilms by molecular dynamics simulation

Size-dependent elastic properties of Ni nanofilms are investigated by molecular dynamics ( MD) simulations with embedded atom method (EAM). The surface effects are considered by calculating the surface relaxation, surface energy, and surface stress. The Young's modulus and yield stress are obtained as functions of thickness and crystallographic orientation. It is shown that the surface relaxation has important effects on the the elastic properties at nanoscale. When the surface relaxation is outward, the Young's modulus decreases with the film thickness decreasing, and vice versa. The results also show that the yield stresses of the films increase with the films becoming thinner. With the thickness of the nanofilms decreasing, the surface effects on the elastic properties become dominant.

Molecular dynamics simulation of adhesion forces in a dipalmitoylphosphatidylcholine membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Effect Of Water On Brittle Fracture Of Sio2 By Molecular Dynamics Study

The influence of water on the brittle behavior of beta-cristobalite is studied by means of molecular dynamics (MD) simulation With the TTAM potential. Crack extension of mode 1 type is observed as the crack opening is filled LIP With water. The critical stress intensity factor K-lc(MD) is used to characterize the crack extension of MD simulation. The surface energy of SiO2 covered with layers of water is calculated at temperature of 300 K. Based oil the Griffith fracture criterion, the critical stress intensity factor K-lc(Griffith) is calculated, and it is in good agreement with that of MD simulation. (C) 2008 Elsevier B.V. All rights reserved.

Shape effects on the yield stress and deformation of silicon nanowires: A molecular dynamics simulation

The tension and compression of single-crystalline silicon nanowires (SiNWs) with different cross-sectional shapes are studied systematically using molecular dynamics simulation. The shape effects on the yield stresses are characterized. For the same surface to volume ratio, the circular cross-sectional SiNWs are stronger than the square cross-sectional ones under tensile loading, but reverse happens in compressive loading. With the atoms colored by least-squares atomic local shear strain, the deformation processes reveal that the failure modes of incipient yielding are dependent on the loading directions. The SiNWs under tensile loading slip in {111} surfaces, while the compressive loading leads the SiNWs to slip in the {110} surfaces. The present results are expected to contribute to the design of the silicon devices in nanosystems.

Splitting the fast and slow motions in molecular dynamics simulations based on the change of cold potential well bottom

Abstract. The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

Molecular dynamics research of G249 and S249 substitutions of p53 protein.

,The molecular dynamics research of the core domain of p53 protein crystal structure shows that besides the stability in biochemistry this domain also shows a high stability in molecular mechanics. Based on that work, the residue R249 was substituted with amino acids Gly and Ser respectively, and molecular dynamics researches were performed separately. The results show that these substitutions cause a relax tendency between loop2 and 3 domains, leading to an alteration of the whole conformation of p53 core domain and ruining its stability. The results visually explains the mechanism of p53 changes in immunological and biochemical reactions, which are caused by 249 residue substitutions from 3-D structure variations.

Structural analysis of human CCR2b and primate CCR2b by molecular modeling and molecular dynamics simulation

CCR2b, a chemokine receptor for MCP-1, -2, -3, -4, plays an important role in a variety of diseases involving infection, inflammation, and/or injury, as well as being a coreceptor for HIV-1 infection. Two models of human CCR2b (hCCR2b) were generated by h

Molecular dynamics simulation on the complex of the tobramycin and 16S rRNA a site

A 3. 6 ns molecular dynamics simulation was carried out on the complex system of tobramycin and 16S rRNA in order to understand the speciality recognition mechanism between tobramycin and 16S rRNA at the molecular level. The results demonstrate that two l

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Controlling electronic structures by irradiation in single-walled SiC nanotubes: a first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, the displacement threshold energy and defect configurations are determined in SiC nanotubes. The simulation results reveal that a rich variety of defect structures (vacancies, Stone-Wales defects and antisite defects) are formed with threshold energies from 11 to 64 eV. The threshold energy shows an anisotropic behavior and exhibits a dramatic decrease with decreasing tube diameter. The electronic structure can be altered by the defects formed by irradiation, which suggests that the electron irradiation may be a way to use defect engineering to tailor electronic properties of SiC nanotubes.

Molecular Dynamics Simulation on Characteristics of Decomposition of Methane Hydrate in the Presence of Hydrogen Peroxide Solution

In this work, the characteristics of the decomposition of methane hydrate Structure I (SI) in the presence of hydrogen peroxide solution is investigated using the molecular dynamics simulation. The mechanism of the transformation process from the solid hydrate to the liquid is analyzed with the effect of hydrogen peroxide (HP) solution. In addition, the effect of ethylene glycol (EG) with the same molar concentration with HP on the methane hydrate dissociation is also studied. The results illustrate that both HP and EG promote well the hydrate dissociation. The work provides the important reference value for the experimental investigation into the promotion effect of HP on the hydrate dissociation.