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.

Temperature dependence of the distribution of the first passage time: Results from discontinuous molecular dynamics simulations of an all-atom model of the second beta-hairpin fragment of protein G

More than 22 000 folding kinetic simulations were performed to study the temperature dependence of the distribution of first passage time (FPT) for the folding of an all-atom Go-like model of the second beta-hairpin fragment of protein G. We find that the mean FPT (MFPT) for folding has a U (or V)-shaped dependence on the temperature with a minimum at a characteristic optimal folding temperature T-opt*. The optimal folding temperature T-opt* is located between the thermodynamic folding transition temperature and the solidification temperature based on the Lindemann criterion for the solid. Both the T-opt* and the MFPT decrease when the energy bias gap against nonnative contacts increases. The high-order moments are nearly constant when the temperature is higher than T-opt* and start to diverge when the temperature is lower than T-opt*. The distribution of FPT is close to a log-normal-like distribution at T* greater than or equal to T-opt*. At even lower temperatures, the distribution starts to develop long power-law-like tails, indicating the non-self-averaging intermittent behavior of the folding dynamics. It is demonstrated that the distribution of FPT can also be calculated reliably from the derivative of the fraction not folded (or fraction folded), a measurable quantity by routine ensemble-averaged experimental techniques at dilute protein concentrations.

Formation of Fivefold Deformation Twins in Nanocrystalline Face-Centered-Cubic Copper Based on Molecular Dynamics Simulations

Fivefold deformation twins were reported recently to be observed in the experiment of the nanocrystalline face-centered-cubic metals and alloys. However, they were not predicted previously based on the molecular dynamics (MD) simulations and the reason was thought to be a uniaxial tension considered in the simulations. In the present investigation, through introducing pretwins in grain regions, using the MD simulations, the authors predict out the fivefold deformation twins in the grain regions of the nanocrystal grain cell, which undergoes a uniaxial tension. It is shown in their simulation results that series of Shockley partial dislocations emitted from grain boundaries provide sequential twining mechanism, which results in fivefold deformation twins. (c) 2006 American Institute of Physics.

Molecular Dynamics Study Of The Mechanical Behavior Of Nickel Nanowire: Strain Rate Effects

We present the analysis of uniaxial deformation of nickel nanowires using molecular dynamics simulations, and address the strain rate effects on mechanical responses and deformation behavior. The applied strain rate is ranging from 1 x 10(8) s(-1) to 1.4 x 10(11) s(-1). The results show that two critical strain rates, i.e., 5 x 10(9) s(-1) and 8 x 10(10) s(-1), are observed to play a pivotal role in switching between plastic deformation modes. At strain rate below 5 x 10(9) s(-1), Ni nanowire maintains its crystalline structure with neck occurring at the end of loading, and the plastic deformation is characterized by {111} slippages associated with Shockley partial dislocations and rearrangements of atoms close to necking region. At strain rate above 8x10(10) s(-1), Ni nanowire transforms from a fcc crystal into a completely amorphous state once beyond the yield point, and hereafter it deforms uniformly without obvious necking until the end of simulation. For strain rate between 5 x 10(9) s(-1) and 8 x 10(10) s(-1), only part of the nanowire exhibits amorphous state after yielding while the other part remains crystalline state. Both the {111} slippages in ordered region and homogenous deformation in amorphous region contribute to the plastic deformation. (C) 2007 Published by Elsevier B.V.

Size effects on the melting of nickel nanowires: a molecular dynamics study

The melting process of nickel nanowires are simulated by using molecular dynamics with the quantum Sutten-Chen many-body force field. The wires studied were approximately cylindrical in cross-section and periodic boundary conditions were applied along their length; the atoms were arranged initially in a face-centred cubic structure with the [0 0 1] direction parallel to the long axis of the wire. The size effects of the nanowires on the melting temperatures are investigated. We find that for the nanoscale regime, the melting temperatures of Ni nanowires are much lower than that of the bulk and are linear with the reciprocal of the diameter of the nanowire. When a nanowire is heated up above the melting temperature, the neck of the nanowire begins to arise and the diameter of neck decreases rapidly with the equilibrated running time. Finally, the breaking of nanowire arises, which leads to the formation of the spherical clusters. (C) 2004 Elsevier B.V. All rights reserved.

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).

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.

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 investigation on polarization retention of barium titanate nanofilm arising from ordered oxygen vacancy

Molecular-dynamics simulations have been carried out to investigate the electric hysteresis of barium titanate nanofilm containing oxygen vacancy ordering array parallel to the {101} crystal plane. The results obtained show a significant weakening of polarization retention from non-zero value to zero as the size of the array was reduced to a critical level, which was attributed to the formation and motion of head-to-head domain wall structure under external field loading process. By comparing with materials containing isolated oxygen vacancies, it was found that the zero retention was due to the oxygen vacancy ordering array rather than to the concentration of oxygen vacancy. Copyright (C) EPLA, 2010

Splitting the Fast and Slow Motions in Molecular Dynamics Simulations Based on the Change of Cold Potential Well Bottom

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 study on the phase diagrams of linear and branched chain molecules

The particle transfer molecular dynamics is used to study the phase equilibria of linear and branched chain molecules. The scaling of the critical temperature versus chain length is obtained and the critical densities are found to decrease with increasing chain length, which are in agreement with the results of experiment and theory. The phase diagrams of the linear and the branched chain molecules nearly overlap with each other. Moreover, the radial distribution functions of linear and branched chain molecules in gas phase are very similar, but in the liquid phase, they are different for different kinds of chains.

Docking and molecular dynamics studies toward the binding of new natural phenolic marine inhibitors and aldose reductase

Phenolic marine natural product is a kind of new potential aldose reductase inhibitors (ARIs). In order to investigate the binding mode and inhibition mechanism, molecular docking and dynamics studies were performed to explore the interactions of six phenolic inhibitors with human aldose reductase (hALR2). Considering physiological environment, all the neutral and other two ionized states of each phenolic inhibitor were adopted in the simulation. The calculations indicate that all the inhibitors are able to form stable hydrogen bonds with the hALR2 active pocket which is mainly constructed by residues TYR48, HIS110 and TRP111, and they impose the inhibition effect by occupying the active space. In all inhibitors, only La and its two ionized derivatives La_ion1 and La_ion2, in which neither of the ortho-hydrogens of 3-hydroxyl is substituted by Br, bind with hALR2 active residues using the terminal 3-hydroxyl. While, all the other inhibitors, at least one of whose ortho-sites of 3- and 6-hydroxyls are substituted by Br substituent which take much electron-withdrawing effect and steric hindrance, bind with hALR2 through the lactone group. This means that the Br substituent can effectively regulate the binding modes of phenolic inhibitors. Although the lactone bound inhibitors have relatively high RMSD values, our dynamics study shows that both binding modes are of high stability. For each inhibitor molecule, the ionization does not change its original binding mode, but it does gradually increase the binding free energy, which reveals that besides hydrogen bonds, the electrostatic effect is also important to the inhibitor–hALR2 interaction.

Forced extension of P-selectin construct using steered molecular dynamics

P-selectin, a 70-nm-long cellular adhesive molecule, possesses elastic and extensible properties when neutrophils roll over the activated endotheliam of blood vessel in inflammatory reaction. Transient formation and dissociation of P-selectin/ligand bond on applied force of blood flow induces the extension of P-selectin and relevant ligands. Steered molecular dynamics simulations were performed to stretch a single P-selectin construct consisting of a lectin (Lec) domain and an epithelial growth factor (EGF)-like domain, where P-selectin construct was forced to extend in water with pulling velocities of 0.005-0.05 nm/ps and with constant forces of 1000-2500 pN respectively. Resulting force-extension profiles exhibited a dual-peak pattern on various velocities, while both plateaus and shoulders appeared in the extension-time profiles on various forces. The force or extension profiles along stretching pathways were correlated to the conformational changes, suggesting that the structural collapses of P-selectin Lec/EGF domains were mainly attributed to the burst of hydrogen bonds within the major beta sheet of EGF domain and the disruptions of two hydrophobic cores of Lee domain. This work furthers the understanding of forced dissociation of P-selectin/ligand bond.

Md Simulation Of Dislocation Mobility During Cutting With Diamond Tip On Silicon

By means of Tersoff and Morse potentials, a three-dimensional molecular dynamics simulation is performed to study atomic force microscopy cutting on silicon monocrystal surface. The interatomic forces between the workpiece and the pin tool and the atoms of workpiece themselves are simulated. Two partial edge dislocations are introduced into workpiece Si, it is found that the motion of dislocations does not occur during the atomic force microscopy cutting processing. Simulation results show that the shear stress acting on dislocations is far below the yield strength of Si. (c) 2008 Elsevier Ltd. All rights reserved.