QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by “GAMESS”, and the rest atoms are treated as MM part calculated by “TINKER”. The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with theQMpart with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(1 0 0) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the imidazole rings are attached to the substrate more tightly than other bases in this peptide.

Forced Dissociation of Selectin-ligand Complexes Using Steered Molecular Dynamics Simulation

Selectin-ligand interactions are crucial to such biological processes as inflammatory cascade or tumor metastasis. How transient formation and dissociation of selectin-ligand bonds in blood flow are coupled to molecular conformation at atomic level, however, has not been well understood. In this study, steered molecular dynamics (SMD) simulations were used to elucidate the intramolecular and intermolecular conformational evolutions involved in forced dissociation of three selectin-ligand systems: the construct consisting of P-selectin lectin (Lec) and epidermal growth factor (EGF)-like domains (P-LE) interacting with synthesized sulfoglycopeptide or SGP-3, P-LE with sialyl Lewis X (sLeX), and E-LE with sLeX. SMD simulations were based on newly built-up force field parameters including carbohydrate units and sulfated tyrosine(s) using an analogy approach. The simulations demonstrated that the complex dissociation was coupled to the molecular extension. While the intramolecular unraveling in P-LESGP-3 system mainly resulted from the destroy of the two anti-parallel sheets of EGF domain and the breakage of hydrogen-bond cluster at the Lec-EGF interface, the intermolecular dissociation was mainly determined by separation of fucose (FUC) from Ca2+ ion in all three systems. Conformational changes during forced dissociations depended on pulling velocities and forces, as well as on how the force was applied. This work provides an insight into better understanding of conformational changes and adhesive functionality of selectin-ligand interactions under external forces.

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by "GAMESS", and the rest atoms are treated as MM part calculated by "TINKER". The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with the QM part with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(100) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the in-tidazole rings are attached to the substrate more tightly than other bases in this peptide.

Accuracy analysis of parallel plate interferometer for angular displacement measurement

The measurement accuracy of a parallel-plate interferometer for angular displacement measurement is analyzed. The measurement accuracy of angular displacement is not only related to the accuracy of phase extraction, but also related to initial incident angle, refraction index and thickness of plane-parallel plate as well as wavelength's stability of laser diode, etc. Theoretical analysis and computer simulation show that the measurement error of the angular displacement bears a minimum value when choosing an optimal initial incident angle in a large range. These analytical results serve as a guide in practical measurement. In this interferometer, reducing the refraction index or increasing the thickness of the parallel plate can improve the measurement accuracy; and the relative error of the phase measurement is 3.0 x 10(-4) corresponding to 1 degrees C temperature variation. Based on these theoretical and experimental results, the measurement accuracy of the parallel-plate interferometer is up to an order of 10(-8) rad. (c) 2007 Elsevier Ltd. All rights reserved.

Simulation of an apodizer's effect for high-density optical storage

The effect of an apodizer with two parallel taper refractive surfaces is theoretically investigated for high-density optical storage. The apodizer may modulate an incident Gaussian beam into an annular beam. Simulation shows that with the increasing inner radius of the modulated beam, the focal spot shrinks obviously. The depolarization effect gets strong simultaneously, which induces the circular symmetry loss of the focal spot. In this process, pattern density of the orthogonal and longitudinal diffractive fields increases remarkably.

Monte Carlo simulation on symmetric ABA/AB copolymer blends in confined thin films

The effects of blend composition on morphology, order-disorder transition (ODT), and chain conformation of symmetric ABA/AB copolymer blends confined between two neutral hard walls have been investigated by lattice Monte Carlo simulation. Only lamellar structure is observed in all the simulation morphologies under thermodynamic equilibrium state, which is supported by theoretical prediction. When the composition of AB diblock copolymer (phi) increases, both lamellar spacing and the corresponding ODT temperature increase, which can be attributed to the variation of conformation distribution of the diblock and the triblock copolymer chains. In addition, both diblock and triblock copolymer, chains with bridge conformation extend dramatically in the direction parallel to the surface when the system is in ordered state. Finally, the copolymer chain conformation depends strongly on both the blend composition and the incompatibility parameter chi N.

Monte Carlo simulation of the aggregation of rod-flexible triblock copolymers in a thin film

The aggregation of rod-flexible ABA and BAB triblock (A was rod block and repulsive with block B) copolymers in a thin film was studied as a function of varying the rigidity (eta) and the length of the rod block by Monte Carlo simulation. The rigidity of block A was defined as eta = R-c/R-max in this study. R-c, was the end-to-end distance below which the conformation of the block was not allowed, whereas R-max, was the longest end-to-end distance that the block could be. If eta = 0 the block was flexible, whereas if eta = 1 the block was a straight rod. The simulation results showed that the ABA triblock copolymer film were likely to form lamella structure with increasing the rigidity (eta) of block A. The lamellas were parallel each other and perpendicular to the film surface. However, the aggregation of BAB triblock copolymers tended to change from lamella to cylinder structure with increasing the rigidity (eta) of block A. Typical lamella and cylinder co-exist structure was obtained at eta = 0.504 for the BAB copolymer film. On the other hand, the simulation results indicated that the film changed from disorder to order, then to disorder structure with increasing the relative length of B block for both ABA and BAB copolymer films.

Energy conversion in shape memory alloy heat engine - Part II: Simulation

In this paper, a theoretical model proposed in Part I (Zhu et al., 2001a) is used to simulate the behavior of a twin crank NiTi SMA spring based heat engine, which has been experimentally studied by Iwanaga et al. (1988). The simulation results are compared favorably with the measurements. It is found that (1) output torque and heat efficiency decrease as rotation speed increase; (2) both output torque and output power increase with the increase of hot water temperature; (3) at high rotation speed, higher water temperature improves the heat efficiency. On the contrary, at low rotation speed, lower water temperature is more efficient; (4) the effects of initial spring length may not be monotonic as reported. According to the simulation, output torque, output power and heat efficiency increase with the decrease of spring length only in the low rotation speed case. At high rotation speed, the result might be on the contrary.

Molecular dynamics simulation of the role of dislocations in microcrack healing

The molecular dynamics method is used to simulate microcrack healing during heating or/and under compressive stress. A centre microcrack in Cu crystal would be sealed under compressive stress or by heating. The role of compressive stress and heating in crack healing was additive. During microcrack healing, dislocation generation and motion occurred. When there were pre-existing dislocations around the microcrack, the critical temperature or compressive stress necessary for microcrack healing would decrease, and, the higher the number of dislocations, the lower the critical temperature or compressive stress. The critical temperature necessary for microcrack healing depended upon the orientation of the crack plane. For example, the critical temperature for the crack along the (001) plane was the lowest, i.e. 770K.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Structural Failure Analysis and Numerical Simulation of Micro-Accelerometers under Impulsive Loading

Micromachined accelerometer is a kind of inertial MEMS devices, which usually operate under intensive impact loading. The reliability of micromachined accelerometers is one of the most important performance indices for their design, manufacture and commer