Prediction of thermal expansion properties of carbon nanotubes using molecular dynamics simulations

The axial coefficients of thermal expansion (CTE) of various carbon nanotubes (CNTs), i.e., single-wall carbon nanotubes (SWCNTs), and some multi-wall carbon nanotubes (MWCNTs), were predicted using molecular dynamics (MDs) simulations. The effects of two parameters, i.e., temperature and the CNT diameter, on CTE were investigated extensively. For all SWCNTs and MWCNTs, the obtained results clearly revealed that within a wide low temperature range, their axial CTEs are negative. As the diameter of CNTs decreases, this temperature range for negative axial CTEs becomes narrow, and positive axial CTEs appear in high temperature range. It was found that the axial CTEs vary nonlinearly with the temperature, however, they decrease linearly as the CNT diameter increases. Moreover, within a wide temperature range, a set of empirical formulations was proposed for evaluating the axial CTEs of armchair and zigzag SWCNTs using the above two parameters. Finally, it was found that the absolute value of the negative axial CTE of any MWCNT is much smaller than those of its constituent SWCNTs, and the average value of the CTEs of its constituent SWCNTs. The present fundamental study is very important for understanding the thermal behaviors of CNTs in such as nanocomposite temperature sensors, or nanoelectronics devices using CNTs.

Formation of carbon nanoscrolls from graphene nanoribbons : a molecular dynamics study

Carbon nanoscrolls (CNSs) are one of the carbon-based nanomaterials similar to carbon nanotubes (CNTs) but are not widely studied in spite of their great potential applications. Their practical applications are hindered by the challenging fabrication of the CNSs. A physical approach has been proposed recently to fabricate the CNS by rolling up a monolayer graphene nanoribbon (GNR) around a CNT driven by the interaction energy between them. In this study, we perform extensive molecular dynamics (MD) simulations to investigate the various factors that impact the formation of the CNS from GNR. Our simulation results show that the formation of the CNS is sensitive to the length of the CNT and temperature. When the GNR is functionalized with hydrogen, the formation of the CNS is determined by the density and distribution of the hydrogen atoms. Graphyne, the allotrope of graphene, is inferior to graphene in the formation of the CNS due to the weaker bonds and the associated smaller atom density. The mechanism behind the rolling of GNR into CNS lies in the balance between the GNR–CNT van der Waals (vdW) interactions and the strain energy of GNR. The present work reveals new important insights and provides useful guidelines for the fabrication of the CNS.

Thermal behavior of kaolinite–urea intercalation complex and molecular dynamics simulation for urea molecule orientation

The thermal behavior of kaolinite–urea intercalation complex was investigated by thermogravimetry–differential scanning calorimetry (TG–DSC), X-ray diffraction (XRD), and fourier transform infrared spectroscopy (FTIR). In addition, the interaction mode of urea molecules intercalated into the kaolinite gallery was studied by means of molecular dynamics simulation. Three main mass losses were observed at 136 °C, in the range of 210–270 °C, and at 500 °C in the TG–DSC curves, which were, respectively, attributed to (1) melting of the surface-adsorbed urea, (2) removal of the intercalated urea, and (3) dehydroxylation of the deintercalated kaolinite. The three DSC endothermic peaks at 218, 250, and 261 °C were related to the successive removals of intercalated urea with three different distribution structures. Based on the angle between the dipole moment vector of urea and the basal surface of kaolinite, the three urea models could be described as follows: (1) Type A, the dipole moment vector is nearly parallel to the basal surface of kaolinite; (2) Type B, the dipole moment vector points to the silica tetrahedron with the angle between it and the basal surface of kaolinite ranging from 20°to 40°; and (3) Type C, the dipole moment vector is nearly perpendicular to the basal surface of kaolinite. The three distribution structures of urea molecules were validated by the results of the molecular dynamics simulation. Furthermore, the thermal behavior of the kaolinite–urea intercalation complex investigated by TG–DSC was also supported by FTIR and XRD analyses.

Fractal and complex network analyses of protein molecular dynamics

Based on protein molecular dynamics, we investigate the fractal properties of energy, pressure and volume time series using the multifractal detrended fluctuation analysis (MF-DFA) and the topological and fractal properties of their converted horizontal visibility graphs (HVGs). The energy parameters of protein dynamics we considered are bonded potential, angle potential, dihedral potential, improper potential, kinetic energy, Van der Waals potential, electrostatic potential, total energy and potential energy. The shape of the h(q)h(q) curves from MF-DFA indicates that these time series are multifractal. The numerical values of the exponent h(2)h(2) of MF-DFA show that the series of total energy and potential energy are non-stationary and anti-persistent; the other time series are stationary and persistent apart from series of pressure (with H≈0.5H≈0.5 indicating the absence of long-range correlation). The degree distributions of their converted HVGs show that these networks are exponential. The results of fractal analysis show that fractality exists in these converted HVGs. For each energy, pressure or volume parameter, it is found that the values of h(2)h(2) of MF-DFA on the time series, exponent λλ of the exponential degree distribution and fractal dimension dBdB of their converted HVGs do not change much for different proteins (indicating some universality). We also found that after taking average over all proteins, there is a linear relationship between 〈h(2)〉〈h(2)〉 (from MF-DFA on time series) and 〈dB〉〈dB〉 of the converted HVGs for different energy, pressure and volume.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.

Adsorption and Unfolding of Lysozyme at a Polarized Aqueous-Organic Liquid Interface

The adsorption of proteins at the interface between two immiscible electrolyte solutions has been found to be key to their bioelectroactivity at such interfaces. Combined with interfacial complexation of organic phase anions by cationic proteins, this adsorption process may be exploited to achieve nanomolar protein detection. In this study, replica exchange molecular dynamics simulations have been performed to elucidate for the first time the molecular mechanism of adsorption and subsequent unfolding of hen egg white lysozyme at low pH at a polarized 1,2-dichloroethane/water interface. The unfolding of lysozyme was observed to occur as soon as it reaches the organic−aqueous interface,which resulted in a number of distinct orientations at the interface. In all cases, lysozyme interacted with the organic phase through regions rich in nonpolar amino acids, such that the side chains are directed toward the organic phase, whereas charged and polar residues were oriented toward the aqueous phase. By contrast, as expected, lysozyme in neat water at low pH does not exhibit significant structural changes. These findings demonstrate the key influence of the organic phase upon adsorption of lysozyme under the influence of an electric field, which results in the unfolding of its structure.

Dynamics deflection analysis of nanoelectromechanical switches based on carbon nanotubes

This paper aims to develop an effective numerical simulation technique for the dynamic deflection analysis of nanotubes-based nanoswitches. The nanoswitch is simplified to a continuum structure, and some key material parameters are extracted from typical molecular dynamics (MD). An advanced local meshless formulation is applied to obtain the discretized dynamic equations for the numerical solution. The developed numerical technique is firstly validated by the static deflection analyses of nanoswitches, and then, the fundamental dynamic properties of nanoswitches are analyzed. A parametric comparison with the results in the literature and from experiments shows that the developed modelling approach is accurate, efficient and effective.

Molecular investigation of the mechanical properties of single actin filaments based on vibration analyses

The mechanical vibration properties of single actin filaments from 50 to 288 nm are investigated by the molecular dynamics simulation in this study. The natural frequencies obtained from the molecular simulations agree with those obtained from the analytical solution of the equivalent Euler–Bernoulli beam model. Through the convergence study of the mechanical properties with respect to the filament length, it was found that the Euler–Bernoulli beam model can only be reliably used when the single actin filament is of the order of hundreds of nanometre scale. This molecular investigation not only provides the evidence for the use of the continuum beam model in characterising the mechanical properties of single actin filaments, but also clarifies the criteria for the effective use of the Euler–Bernoulli beam model.

Nuclear instruments and methods in physics research section B : beam interactions with materials and atoms

The impact-induced deposition of Al13 clusters with icosahedral structure on Ni(0 0 1) surface was studied by molecular dynamics (MD) simulation using Finnis–Sinclair potentials. The incident kinetic energy (Ein) ranged from 0.01 to 30 eV per atom. The structural and dynamical properties of Al clusters on Ni surfaces were found to be strongly dependent on the impact energy. At much lower energy, the Al cluster deposited on the surface as a bulk molecule. However, the original icosahedral structure was transformed to the fcc-like one due to the interaction and the structure mismatch between the Al cluster and Ni surface. With increasing the impinging energy, the cluster was deformed severely when it contacted the substrate, and then broken up due to dense collision cascade. The cluster atoms spread on the surface at last. When the impact energy was higher than 11 eV, the defects, such as Al substitutions and Ni ejections, were observed. The simulation indicated that there exists an optimum energy range, which is suitable for Al epitaxial growth in layer by layer. In addition, at higher impinging energy, the atomic exchange between Al and Ni atoms will be favourable to surface alloying.

The role of the molecular footprint of EGFR in tailoring treatment decisions in NSCLC

The majority of patients with non-small-cell lung cancer (NSCLC) present with advanced disease, with targeted therapies providing some improvement in clinical outcomes. The epidermal growth factor receptor (EGFR) tyrosine kinase (TK) plays an important role in the pathogenesis of NSCLC. Tyrosine kinase inhibitors (TKIs), which target the EGFR TK domain, have proven to be an effective treatment strategy; however, patient responses to treatment vary considerably. Therefore, the identification of patients most likely to respond to treatment is essential to optimise the benefit of TKIs. Tumour-associated activating mutations in EGFR can identify patients with NSCLC who are likely to have a good response to TKIs. Nonetheless, the majority of patients relapse within a year of starting treatment. Studies of tumours at relapse have demonstrated expression of a T790M mutation in exon 20 of the EGFR TK domain in approximately 50% of cases. Although conferring resistance to reversible TKIs, these patients may remain sensitive to new-generation irreversible/panerb inhibitors. A number of techniques have been employed for genotypic assessment of tumourassociated DNA to identify EGFR mutations, each of which has advantages and disadvantages. This review presents an overview of the current methodologies used to identify such molecular markers. Recent developments in technology may make the monitoring of changes in patients' tumour genotypes easier in clinical practice, which may enable patients' treatment regimens to be tailored during the course of their disease, potentially leading to improved patient outcomes.

Accurate measurement of the molecular thickness of thin organic shells on small inorganic cores using dynamic light scattering

Dynamic light scattering (DLS) has become a primary nanoparticle characterization technique with applications from materials characterization to biological and environmental detection. With the expansion in DLS use from homogeneous spheres to more complicated nanostructures, comes a decrease in accuracy. Much research has been performed to develop different diffusion models that account for the vastly different structures but little attention has been given to the effect on the light scattering properties in relation to DLS. In this work, small (core size < 5 nm) core-shell nanoparticles were used as a case study to measure the capping thickness of a layer of dodecanethiol (DDT) on Au and ZnO nanoparticles by DLS. We find that the DDT shell has very little effect on the scattering properties of the inorganic core and hence can be ignored to a first approximation. However, this results in conventional DLS analysis overestimating the hydrodynamic size in the volume and number weighted distributions. By introducing a simple correction formula that more accurately yields hydrodynamic size distributions a more precise determination of the molecular shell thickness is obtained. With this correction, the measured thickness of the DDT shell was found to be 7.3 ± 0.3 Å, much less than the extended chain length of 16 Å. This organic layer thickness suggests that on small nanoparticles, the DDT monolayer adopts a compact disordered structure rather than an open ordered structure on both ZnO and Au nanoparticle surfaces. These observations are in agreement with published molecular dynamics results.

The molecular structure of kaolinite–potassium acetate intercalation complexes: A combined experimental and molecular dynamic simulation study

The kaolinite (Kaol) intercalated with potassium acetate (Ac) was prepared and characterized by X-ray diffraction (XRD), Fourier transform infrared spectroscopy (FTIR), and thermogravimetry. Molecular dynamic simulation was performed to investigate the structure of Kaol–Ac intercalation complex and the hydrogen bonds between Kaol and intercalated Ac andwater using INTERFACE forcefield. The acetate anions andwater arranged in a bilayer structure in the interlayer space of Kaol. The potassium cations distributed in the interlayer space and strongly coordinated with acetate anions aswell aswater rather than keyed into the ditrigonal holes of tetrahedral surface of Kaol. Strong hydrogen bonds formed between the hydrogen atoms of hydroxyl on the octahedral surface and oxygen atoms of both acetate anions and water. The acetate anions andwater also weakly bonded hydrogen to the silica tetrahedral surface through their hydrogen atoms with the oxygen atoms of silica tetrahedral surface.

Comparing charge transport predictions for a ternary electrolyte using the Maxwell–Stefan and Nernst–Planck equations

In this work, we investigate and compare the Maxwell–Stefan and Nernst–Planck equations for modeling multicomponent charge transport in liquid electrolytes. Specifically, we consider charge transport in the Li+/I−/I3−/ACN ternary electrolyte originally found in dye-sensitized solar cells. We employ molecular dynamics simulations to obtain the Maxwell–Stefan diffusivities for this electrolyte. These simulated diffusion coefficients are used in a multicomponent charge transport model based on the Maxwell– Stefan equations, and this is compared to a Nernst–Planck based model which employs binary diffusion coefficients sourced from the literature. We show that significant differences between the electrolyte concentrations at electrode interfaces, as predicted by the Maxwell–Stefan and Nernst–Planck models, can occur. We find that these differences are driven by a pressure term that appears in the Maxwell–Stefan equations. We also investigate what effects the Maxwell–Stefan diffusivities have on the simulated charge transport. By incorporating binary diffusivities found in the literature into the Maxwell–Stefan framework, we show that the simulated transient concentration profiles depend on the diffusivities; however, the simulated equilibrium profiles remain unaffected.