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.

Microbead dynamics on quartz crystal microbalance at elevated amplitudes

The application of the Quartz Crystal Microbalance (QCM) for biochemical sensing is well known. However, utilizing the nonlinear response of the QCM at elevated amplitudes has received sporadic attention. This study presents results for QCM-analyte interaction that provide insight into the nonlinear dynamics of the QCM with attached analyte. In particular, interactions of the QCM with polystyrene microbeads physisorbed via self-assembled monolayer (SAM) were studied through experiments and modelling. It was found that the response of the QCM coupled to these surface adsorbents is anharmonic even at low oscillation amplitudes and that the nonlinear signals from such interactions are much higher than those for bare quartz. Therefore, these signals can potentially be used as sensitive signatures of adsorbents and their kinetics on the surface. ©2009 IEEE.

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.

Transient Ventilation Dynamics Following a Change in Strength of a Point Source of Heat

We investigate the transient ventilation flow within a confined ventilated space, with high- and low-level openings, when the strength of a low-level point source of heat is changed instantaneously. The steady-flow regime in the space involves a turbulent buoyant plume, which rises from the point source to a well-mixed warm upper layer. The steady-state height of the interface between this layer and the lower layer of exterior fluid is independent of the heat flux, but the upper layer becomes progressively warmer with heat flux. New analogue laboratory experiments of the transient adjustment between steady states identify that if the heat flux is increased, the continuing plume propagates to the top of the room forming a new, warmer layer. This layer gradually deepens, and as the turbulent plume entrains fluid from the original warm layer, the original layer is gradually depleted and disappears, and a new steady state is established. In contrast, if the source buoyancy flux is decreased, the continuing plume is cooler than the original plume, so that on reaching the interface it is of intermediate density between the original warm layer and the external fluid. The plume supplies a new intermediate layer, which gradually deepens with the continuing flow. In turn, the original upper layer becomes depleted, both as a result of being vented through the upper opening of the space, but also due to some penetrative entrainment of this layer by the plume, as the plume overshoots the interface before falling back to supply the new intermediate layer. We develop quantitative models which are in good accord with our experimental data, by combining classical plume theory with models of the penetrative entrainment for the case of a decrease in heating. Typically, we find that the effect of penetrative entrainment on the density of the intruding layer is relatively weak, provided the change in source strength is sufficiently large. However, penetrative entrainment measurably increases the rate at which the depth of the draining layer decreases. We conclude with a discussion of the importance of these results for the control of naturally ventilated spaces.

Cross-relaxation dynamics on anomalous saturation processes in low pressure supersonic cw HF chemical laser amplifier

It is assumed that both translational and rotational nonequilibrium cross-relaxations play a role simultaneoulsy in low pressure supersonic cw HF chemical laser amplifier. For two-type models of gas flow medium with laminar and turbulent flow diffusion mixing, the expressions of saturated gain spectrum are derived respectively, and the numerical calculations are performed as well. The numerical results show that turbulent flow diffusion mixing model is in the best agreement with the experimental result.

Intracellular demography and the dynamics of Salmonella enterica infections.

An understanding of within-host dynamics of pathogen interactions with eukaryotic cells can shape the development of effective preventive measures and drug regimes. Such investigations have been hampered by the difficulty of identifying and observing directly, within live tissues, the multiple key variables that underlay infection processes. Fluorescence microscopy data on intracellular distributions of Salmonella enterica serovar Typhimurium (S. Typhimurium) show that, while the number of infected cells increases with time, the distribution of bacteria between cells is stationary (though highly skewed). Here, we report a simple model framework for the intensity of intracellular infection that links the quasi-stationary distribution of bacteria to bacterial and cellular demography. This enables us to reject the hypothesis that the skewed distribution is generated by intrinsic cellular heterogeneities, and to derive specific predictions on the within-cell dynamics of Salmonella division and host-cell lysis. For within-cell pathogens in general, we show that within-cell dynamics have implications across pathogen dynamics, evolution, and control, and we develop novel generic guidelines for the design of antibacterial combination therapies and the management of antibiotic resistance.

Molecular Dynamics Simulation of Thermal and Mechanical Properties of Polyimide-Carbon-Nanotube Composites

An aromatic polyimide and its mixture with randomly distributed carbon nanotubes (NTs) are simulated by using molecular dynamics, repeated energy minimization and cooling processes. The glass transition temperatures are identified through volume-temperature curves. Stress-strain curves, Young's moduli, densities and Poisson ratios are computed at different temperatures. It is demonstrated that the carbon NT reduces the softening effects of temperature on mechanical properties and increases the ability to resist deformation.

Modelling nonlinear vehicle dynamics with neural networks

The nonlinear modelling ability of neural networks has been widely recognised as an effective tool to identify and control dynamic systems, with applications including nonlinear vehicle dynamics which this paper focuses on using multi-layer perceptron networks. Existing neural network literature does not detail some of the factors which effect neural network nonlinear modelling ability. This paper investigates into and concludes on required network size, structure and initial weights, considering results for networks of converged weights. The paper also presents an online training method and an error measure representing the network's parallel modelling ability over a range of operating conditions. Copyright © 2010 Inderscience Enterprises Ltd.

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated.We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard–Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

The design of an agent for modelling supply chain network dynamics

A useful insight into managerial decision making can be found from simulation of business systems, but existing work on simulation of supply chain behaviour has largely considered non-competitive chains. Where competitive agents have been examined, they have generally had a simple structure and been used for fundamental examination of stability and equilibria rather than providing practical guidance to managers. In this paper, a new agent for the study of competitive supply chain network dynamics is proposed. The novel features of the agent include the ability to select between competing vendors, distribute orders preferentially among many customers, manage production and inventory, and determine price based on competitive behaviour. The structure of the agent is related to existing business models and sufficient details are provided to allow implementation. The agent is tested to demonstrate that it recreates the main results of the existing modelling and management literature on supply chain dynamics. A brief exploration of competitive dynamics is given to confirm that the proposed agent can respond to competition. The results demonstrate that overall profitability for a supply chain network is maximised when businesses operate collectively. It is possible for an individual business to achieve higher profits by adopting a more competitive stance, but the consequence of this is that the overall profitability of the network is reduced. The agent will be of use for a broad range of studies on the long-run effect of management decisions on their network of suppliers and customers.

Dynamics of collapsing fluidized beds and its application in the simulation of pulsed fluidized beds

A new method to study collapsing process of fluidized bed was proposed. The method is based on the analysis of the pressure variation during collapse. A model is proposed to describe the pressure variation on any location of the bed during collapse. Three kind of particles were fluidized by air and the pressure variation after shutting down the gas supply was measured by pressure transducers and simulated by the proposed model. The simulated results were in good agreement with the experimental data. The parameters of the fluidized bed, such as bubble fraction alpha(b), particle fraction in dense phase alpha(p), bubble velocity u(b) and slip velocity of particle u(i), can be estimated from the measured pressure data by using the proposed model. The model was also used to calculate the pressure variation in pulsed fluidized beds with low frequency.

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.

A Molecular Dynamics Simulation: the Effect of Finite Size on the Thermal Conductivity in a Single Crystal Silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.

A Constrained Particle Dynamics for Continuum-Particle Hybrid Method in Micro- and Nano-Fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier-Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.