A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

Stress-induced phase transformation and pseudo-elastic/pseudo-plastic recovery in intermetallic Ni-Al nanowires

Extensive molecular dynamics (MD) simulations have been performed in a B2-NiAl nanowire using an embedded atom method (EAM) potential. We show a stress induced B2 -> body-centered-tetragonal (BCT) phase transformation and a novel temperature and cross-section dependent pseudo-elastic/pseudo-plastic recovery from such an unstable BCT phase with a recoverable strain of similar to 30% as compared to 5-8% in polycrystalline materials. Such a temperature and cross-section dependent pseudo-elastic/pseudo-plastic strain recovery can be useful in various interesting applications of shape memory and strain sensing in nanoscale devices. Effects of size, temperature, and strain rate on the structural and mechanical properties have also been analyzed in detail. For a given size of the nanowire the yield stress of both the B2 and the BCT phases is found to decrease with increasing temperature, whereas for a given temperature and strain rate the yield stress of both the B2 and the BCT phase is found to increase with increase in the cross-sectional dimensions of the nanowire. A constant elastic modulus of similar to 80 GPa of the B2 phase is observed in the temperature range of 200-500 K for nanowires of cross-sectional dimensions in the range of 17.22-28.712 angstrom, whereas the elastic modulus of the BCT phase shows a decreasing trend with an increase in the temperature.

The designed environment and how it affects brain morphology and mental health

Background The environment is inextricably related to mental health. Recent research replicates findings of a significant, linear correlation between a childhood exposure to the urban environment and psychosis. Related studies also correlate the urban environment and aberrant brain morphologies. These findings challenge common beliefs that the mind and brain remain neutral in the face of worldly experience. Aim There is a signature within these neurological findings that suggests that specific features of design cause and trigger mental illness. The objective in this article is to work backward from the molecular dynamics to identify features of the designed environment that may either trigger mental illness or protect against it. Method This review analyzes the discrete functions putatively assigned to the affected brain areas and a neurotransmitter called dopamine, which is the primary target of most antipsychotic medications. The intention is to establish what the correlations mean in functional terms, and more specifically, how this relates to the phenomenology of urban experience. In doing so, environmental mental illness risk factors are identified. Conclusions Having established these relationships, the review makes practical recommendations for those in public health who wish to use the environment itself as a tool to improve the mental health of a community through design.

Experimental and computational studies on lipoprotein lipolysis at acidic pH

Atherosclerosis is a disease of the arteries; its characteristic features include chronic inflammation, extra- and intracellular lipid accumulation, extracellular matrix remodeling, and an increase in extracellular matrix volume. The underlying mechanisms in the pathogenesis of advanced atherosclerotic plaques, that involve local acidity of the extracellular fluid, are still incompletely understood. In this thesis project, my co-workers and I studied the different mechanisms by which local extracellular acidity could promote accumulation of the atherogenic apolipoprotein B-100 (apoB-100)-containing plasma lipoprotein particles in the inner layer of the arterial wall, the intima. We found that lipolysis of atherogenic apoB-100-containing plasma lipoprotein particles (LDL, IDL, and sVLDL) by the secretory phospholipase A2 group V (sPLA2-V) enzyme, was increased at acidic pH. Also, the binding of apoB-100-containing plasma lipoprotein particles to human aortic proteoglycans was dramatically enhanced at acidic pH. Additionally, lipolysis by sPLA2-V enzyme further increased this binding. Using proteoglycan-affinity chromatography, we found that sVLDL lipoprotein particles consist of populations, differing in their affinities toward proteoglycans. These populations also contained different amounts of apolipoprotein E (apoE) and apolipoprotein C-III (apoC-III); the amounts of apoC-III and apoE per particle were highest in the population with the lowest affinity toward proteoglycans. Since PLA2-modification of LDL particles has been shown to change their aggregation behavior, we also studied the effect of acidic pH on the monolayer structure covering lipoprotein particles after PLA2-induced hydrolysis. Using molecular dynamics simulations, we found that, in acidity, the monolayer is more tightly packed laterally; moreover, its spontaneous curvature is negative, suggesting that acidity may promote lipoprotein particles fusion. In addition to extracellular lipid accumulation, the apoB-100-containing plasma lipoprotein particles can be taken up by inflammatory cells, namely macrophages. Using radiolabeled lipoprotein particles and cell cultures, we showed that sPLA2-V-modification of LDL, IDL, and sVLDL lipoproteins particles, at neutral or acidic pH, increased their uptake by human monocyte-derived macrophages.

Graphene with patterned fluorination: Morphology modulation and implications

Recent years have witnessed a large volume of works on the modification of graphene; however, an understanding of the associated morphology or mechanical properties changes is still lacking, which is vital for its engineering implementation. By taking the C4F fluorination as an example, we find that the morphology of both graphene sheet (GS) and graphene nanoribbon (GNR) can be effectively tailored by fluorination patterning via molecular dynamics simulations. The fluorine atom produces out-of-plane forces which trigger several intriguing morphology changes to monolayer graphene, including zigzag, folded, ruffle, nanoscroll, and chain structures. Notably, for multilayer GNR, the delamination and climbing phenomena of the surface layer are observed. Further studies show that the fluorination pattern can also be utilized to modulate the mechanical properties of graphene, e.g., about 40% increase of the effective yield strain is observed for the examined GNR with fluorination patterns. This study not only demonstrates the significant impacts on the morphology of graphene from fluorination but also suggests an effective avenue to tailor the morphology and thus mechanical properties of GS and GNR.

Jump Reorientation of Water Molecules Confined in Narrow Carbon Nanotubes

We used molecular dynamics (MD) simulations to study the reorientational dynamics of water molecules confined inside narrow carbon nanotubes immersed in a bath of water. Our simulations show that the confined water molecules exhibit bistability in their reorientational relaxation, which proceeds by angular jumps between the two stable states. The angular jump of a water molecule in the bulk involves the breaking of a hydrogen bond with one of its neighbors and the formation of a hydrogen bond with a different neighbor. In contrast, the angular jump of a confined water molecule corresponds to an interchange of the two hydrogen atoms that can form a hydrogen bond with the same neighbor. The free energy barrier between these two states is a few k(B)T. The analytic solution of a simplified two-state jump model that qualitatively explains the reorientational behavior observed in simulations is also presented.

Effects of Isothermal and Adiabatic Thermal Loadings on Size and Strain Rate Dependence of Copper Nanowire

In the present paper, the size and strain rate effects on ultra-thin < 100 >/{100} Cu nanowires at an initial temperature of 10 K have been discussed. Extensive molecular dynamics (MD) simulations have been performed using Embedded atom method (EAM) to investigate the structural behaviours and properties under high strain rate. Velocity-Verlet algorithm has been used to solve the equation of motions. Two different thermal loading cases have been considered: (i) Isothermal loading, in which Nose-Hoover thermostat is used to maintain the constant system temperature, and (ii) Adiabatic loading, i.e., without any thermostat. Five different wire cross-sections were considered ranging from 0.723 x 0.723 nm(2) to 2.169 x 2.169 nm(2) The strain rates used in the present study were 1 x 10(9) s(-1), 1 x 10(8) s(-1), and 1 x 10(7) s(-1). The effect of strain rate on the mechanical properties of copper nanowires was analysed, which shows that elastic properties are independent of thermal loading for a given strain rate and cross-sectional dimension of nanowire. It showed a decreasing yield stress and yield strain with decreasing strain rate for a given cross- section. Also, a decreasing yield stress and increasing yield strain were observed for a given strain rate with increasing cross-sectional area. Elastic modulus was found to be similar to 100 GPa, which was independent of processing temperature, strain rate, and size for a given initial temperature. Reorientation of < 100 >/{100} square cross-sectional copper nanowire into a series of stable ultra-thin Pentagon copper nanobridge structures with dia of similar to 1 nm at 10 K was observed under high strain rate tensile loading. The effect of isothermal and adiabatic loading on the formation of such pentagonal nanobridge structure has been discussed.

The Influence of Bilayer Composition on the Gel to Liquid Crystalline Transition

We report molecular dynamics simulations of bilayers using a united atom model with explicit solvent molecules. The bilayer consists of the single tail cationic surfactant behenyl trimethyl ammonium chloride (BTMAC) with stearyl alcohol (SA) as the cosurfactant. We study the gel to liquid crystalline transitions in the bilayer by varying the amount of water at fixed BTMAC to SA ratio as well as by varying the BTMAC to SA ratio at fixed water content. The bilayer is found to exist in the tilted, Lβ′ phase at low temperatures, and for the compositions investigated in this study, the Lβ′ to Lα melting transition occurred in the temperature range 330−338 K. For the highest BTMAC to SA composition (2:3 molar ratio), a diffuse headgroup−water interface is observed at lower temperatures, and an increase in the d-spacing occurs prior to the melting transition. This pretransition swelling is accompanied by a sharpening in the water density variation across the headgroup region of the bilayer. Signatures of this swelling effect which can be observed in the alkane density distributions, area per headgroup, and membrane thickness are attributed to the hydrophobic effect. At a fixed bilayer composition, the transition temperature (>338 K) from the Lβ′ to Lα transition obtained for the high water content bilayer (80 wt %) is similar to that obtained with low water content (54.3 wt %), confirming that the melting transition at these water contents is dominated by chain melting.

A computational study of Compton profiles of water - methanol solutions

The molecular level structure of mixtures of water and alcohols is very complicated and has been under intense research in the recent past. Both experimental and computational methods have been used in the studies. One method for studying the intra- and intermolecular bindings in the mixtures is the use of the so called difference Compton profiles, which are a way to obtain information about changes in the electron wave functions. In the process of Compton scattering a photon scatters inelastically from an electron. The Compton profile that is obtained from the electron wave functions is directly proportional to the probability of photon scattering at a given energy to a given solid angle. In this work we develop a method to compute Compton profiles numerically for mixtures of liquids. In order to obtain the electronic wave functions necessary to calculate the Compton profiles we need some statistical information about atomic coordinates. Acquiring this using ab-initio molecular dynamics is beyond our computational capabilities and therefore we use classical molecular dynamics to model the movement of atoms in the mixture. We discuss the validity of the chosen method in view of the results obtained from the simulations. There are some difficulties in using classical molecular dynamics for the quantum mechanical calculations, but these can possibly be overcome by parameter tuning. According to the calculations clear differences can be seen in the Compton profiles of different mixtures. This prediction needs to be tested in experiments in order to find out whether the approximations made are valid.

Reorientation of Water Inside Carbon Nanorings by Large Angular Jumps

Molecular dynamics simulations of the orientational dynamics of water molecules confined inside narrow carbon nanorings reveal that reorientational relaxation is mediated by large amplitude angular jumps. The distribution of waiting time between jumps peaks at about 60 fs, and has a slowly decaying exponential tail with a timescale of about 440 fs. These time scales are much faster than the mean waiting time between jumps of the water molecules in bulk.

Higher-resolution structure of the human insulin receptor ectodomain: Multi-modal inclusion of the insert domain

Insulin receptor (IR) signaling is critical to controlling nutrient uptake and metabolism. However, only a low-resolution (3.8 Å) structure currently exists for the IR ectodomain, with some segments ill-defined or unmodeled due to disorder. Here, we revise this structure using new diffraction data to 3.3 Å resolution that allow improved modeling of the N-linked glycans, the first and third fibronectin type III domains, and the insert domain. A novel haptic interactive molecular dynamics strategy was used to aid fitting to low-resolution electron density maps. The resulting model provides a foundation for investigation of structural transitions in IR upon ligand binding.

A phosphorylation-induced turn defines the Alzheimer's disease AT8 antibody epitope on the tau protein

Post mortem biochemical staging of Alzheimer’s disease is currently based on immunochemical analysis of brain slices with the AT8 antibody. The epitope of AT8 is described around the pSer202/pThr205 region of the hyperphosphorylated form of the neuronal protein tau. In this study, NMR spectroscopy was used to precisely map the AT8 epitope on phosphorylated tau, and derive its defining structural features by a combination of NMR analyses and molecular dynamics. A particular turn conformation is stabilized by a hydrogen bond of the phosphorylated Thr205 residue to the amide proton of Gly207, and is further stabilized by the two Arg residues opposing the pSer202/pThr205.

Free energy calculations of glycosaminoglycan - Protein interactions

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5-6 of PECAM-1 and a region of high affinity in domains 2-3, consistent with experimental data and ligand-protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.