Quantum chemical molecular dynamics simulation of carbon nanotube–graphene fusion

Abstract

Here we report a quantum mechanical molecular dynamics (QM/MD) study of a fusion process of an open-ended carbon nanotube on a graphene hole, which results in the formation of a so-called pillared graphene structure – a three-dimensional nanomaterial consisting entirely of sp²-carbons. The self-consistent-charge density-functional tight-binding potential was adopted in this study. Two different sizes of graphene holes with 12 or 24 central carbon atoms removed from a graphene flake, and a (6,6) carbon nanotube with a compatible diameter were adopted. Formations of 6–7–6/5–8–5 defect structures were found on the fusion border between tube and graphene hole. The 6–7–6 structure was found to bear less curvature-induced strain energy and therefore to be more stable and much easier to form than the 5–8–5 structure.     

Deformation mechanisms in nanotwinned γ-TiAl by molecular dynamics simulation

ABSTRACT

In this work, the plastic deformation mechanisms and fracture toughness of nanotwinned γ-TiAl with different twin boundary (TB) spacing are investigated by using molecular dynamics simulation. The simulation results reveal that there are pronounced shifts in the mechanical behaviour of nanotwinned γ-TiAl when the TB spacing is 3.50, 4.20 and 4.90?nm. In addition, the variation of the dislocation density with strain at these three TB spacing illustrates that a smaller TB spacing induces a higher dislocation density. Different TB spacing has an influence on the dislocation behaviour. The dislocation pile-up, dislocation–dislocation, dislocation–twin and twin–twin reactions, hierarchical twins including their generation and density, step formation, dislocation emission from steps and TB migration are the main plastic deformation mechanisms. The results also show that TB migration, twinning formation and interaction of crack and TB dominate the deformation mechanism of nanotwinned γ-TiAl with crack. The generation of hierarchical twins, lower distance between crack surface plane and twin plane, dislocation–twin, twin–twin interaction and crack deflection increase the fracture toughness of nanotwinned γ-TiAl.     

Understanding LuxS-based quorum sensing and its inhibition – molecular dynamics simulation study

Molecular dynamics simulations were successfully applied to LuxS protein and its protein–ligand complex using the newly developed force field parameters for the iron containing active site. To the best of our knowledge, this was the first attempt to develop force field parameters for the iron containing active site of the LuxS protein. From the simulations, catalytically important amino acid residues were identified which were found to stabilise the ligand. Two residues Glu57 and Asp77 were involved in polar interactions while the protein region in the range between amino acid residue 125 and 131 were predicted to facilitate the entry of ligand to the active site. Other residues like Arg65, Asp77, Ile78 and Ser79 were also recognised as ligand stabilising factors deduced from the simulation. These results were also found to be in good agreement with earlier studies and thus demonstrated the successful application of MD simulations to the LuxS protein. Moreover, the simulation data were expected to be considered for the development of rational approaches in order to identify new LuxS-based quorum sensing antagonists for the treatment of pathologies caused by resistant bacteria.     

Correlated ab initio molecular dynamics simulations of the acetone–carbon dioxide complex: implications for solubility in supercritical CO2

Ab initio molecular dynamics simulations of the acetone–CO₂ complex (MP2/6-31G(d) level) were performed to investigate the effect of dynamics at finite temperature on the weak electron donor–acceptor intermolecular interactions. In addition, we carried out a study of the free energy of formation of the complex by means of umbrella sampling technique at the MP2 level with a perturbative CCSD(T) correction. The potential of mean force was obtained along a reaction coordinate describing the acetone–CO₂ interaction. The results obtained here support some hypothesis that we already explored in past works using static electronic calculations. In particular, when interacting with a molecule having a carbonyl function, carbon dioxide displays both Lewis acid and Lewis base behaviour. This property can be exploited to design molecular systems that are easily solubilised in supercritical CO₂.     

Cα torsion angles as a flexible criterion to extract secrets from a molecular dynamics simulation

A multipolar, polarizable electrostatic method for future use in a novel force field is described. Quantum Chemical Topology (QCT) is used to partition the electron density of a chemical system into atoms, then the machine learning method Kriging is used to build models that relate the multipole moments of the atoms to the positions of their surrounding nuclei. The pilot system serine is used to study both the influence of the level of theory and the set of data generator methods used. The latter consists of: (i) sampling of protein structures deposited in the Protein Data Bank (PDB), or (ii) normal mode distortion along either (a) Cartesian coordinates, or (b) redundant internal coordinates. Wavefunctions for the sampled geometries were obtained at the HF/6-31G(d,p), B3LYP/apc-1, and MP2/cc-pVDZ levels of theory, prior to calculation of the atomic multipole moments by volume integration. The average absolute error (over an independent test set of conformations) in the total atom-atom electrostatic interaction energy of serine, using Kriging models built with the three data generator methods is 11.3 kJ mol^-1 (PDB), 8.2 kJ mol^-1 (Cartesian distortion), and 10.1 kJ mol^-1 (redundant internal distortion) at the HF/6-31G(d,p) level. At the B3LYP/apc-1 level, the respective errors are 7.7 kJ mol^-1, 6.7 kJ mol^-1, and 4.9 kJmol^-1, while at the MP2/cc-pVDZ level they are 6.5 kJ mol^-1, 5.3 kJ mol^-1, and 4.0 kJmol^-1. The ranges of geometries generated by the redundant internal coordinate distortion and by extraction from the PDB are much wider than the range generated by Cartesian distortion. The atomic multipole moment and electrostatic interaction energy predictions for the B3LYP/apc-1 and MP2/cc-pVDZ levels are similar, and both are better than the corresponding predictions at the HF/6-31G(d,p) level.     

Study on properties of liquid ammonia via molecular dynamics simulation based on ABEEMσπ polarisable force field

Abstract

The molecular dynamics simulation has been performed to investigate the charge distribution, structural and dynamical properties of liquid ammonia at 273 K using a polarisable force field of the atom-bond electronegativity equalisation method (ABEEMσπ). One ammonia molecule in this model has eight charge sites, one N atomic site, three H atomic sites, three N–H bond sites and one lone-pair electron site. ABEEMσπ model can present the quantitative site charges of molecular ammonias in liquid and their changing in response to their surroundings. The radial distribution functions and dynamical properties are in fair agreement with the available experimental data. The first peak of g_NN(r) appears at N–N distance of ~3.50 ± 0.05 Å where most hydrogen bonds are formed. The average coordination number of the first shell is 13.0 ± 0.1 among which a central ammonia molecule intimately connects 3 ~ 4 ammonia molecules by hydrogen bonds. The power spectrum shows the vibrations of hydrogen bonds. For a reference, a simple estimation of the average hydrogen bonding energy in liquid ammonia is 6.5 ± 0.1 kcal/mol larger than 3.8 ± 0.3 kcal/mol in dimer ammonia. Our simulation results provide more detailed information about liquid ammonia.     

Characterization of molecular recognition of Phosphoinositide-3-kinase α inhibitor through molecular dynamics simulation

Phosphatidylinositol 3-kinase α (PI3Kα) is a promising target for anticancer drug discovery due to its overactivation in tumor cells. To systematically investigate the interactions between PI3Kα and PIK75 which is the most selective PI3Kα inhibitor reported to date, molecular docking, molecular dynamics simulation, and ensuing energetic analysis were utilized. The binding free energy between PI3Kα and PIK75 is −10.04 kcal•mol⁻¹ using MMPBSA method, while −13.88 kcal•mol⁻¹ using MMGBSA method, which is beneficial for the binding. The van der Waals/hydrophobic and electrostatic interactions play critical roles for the binding. The binding mode of PIK75 for PI3Kα is predicted. The conserved hydrophobic adenine region of PI3Kα made up of Ile800, Ile848, Val850, Val851, Met922, Phe930, and Ile932 accommodates the flat 6-bromine imidazo[1,2-a]pyridine ring of PIK75. The 2-methyl-5-nitrophenyl group of PIK75 extends to the P-loop region, and has four hydrogen-bond arms with the backbone and side chain of Ser773 and Ser774. And the distinct conformation of the P-loop induced by PIK75 is speculated to be responsible for the selectivity profile of PIK75. The predicted binding mode of PIK75 for PI3Kα presented in this study may help design high affinity and selective compounds to target PI3Kα.     

Structure of the antimicrobial β-hairpin peptide protegrin-1 in a DLPC lipid bilayer investigated by molecular dynamics simulation

All atom molecular dynamics simulations of the 18-residue β-hairpin antimicrobial peptide protegrin-1 (PG-1, RGGRLCYCRRRFCVCVGR-NH₂) in a fully hydrated dilauroylphosphatidylcholine (DLPC) lipid bilayer have been implemented. The goal of the reported work is to investigate the structure of the peptide in a membrane environment (previously solved only in solution [R.L. Fahrner, T. Dieckmann, S.S.L. Harwig, R.I. Lehrer, D. Eisenberg, J. Feigon, Solution structure of protegrin-1, a broad-spectrum antimicrobial peptide from porcine leukocytes. Chemistry and Biology, 3 (1996) 543-550]), and to delineate specific peptide-membrane interactions which are responsible for the peptide's membrane binding properties. A novel, previously unknown, “kick” shaped conformation of the peptide was detected, where a bend at the C-terminal β-strand of the peptide caused the peptide backbone at residues 16-18 to extend perpendicular to the β-hairpin plane. This bend was driven by a highly persistent hydrogen-bond between the polar peptide side-chain of TYR7 and the unshielded backbone carbonyl oxygen atom of GLY17. The H-bond formation relieves the unfavorable free energy of insertion of polar groups into the hydrophobic membrane core. PG-1 was anchored to the membrane by strong electrostatic binding of the protonated N-terminus of the peptide to the lipid head group phosphate anions. The orientation of the peptide in the membrane, and its influence on bilayer structural and dynamic properties are in excellent agreement with solid state NMR measurements [S. Yamaguchi, T. Hong, A. Waring, R.I. Lehrer, M. Hong, Solid-State NMR Investigations of Peptide-Lipid Interaction and Orientation of a b-Sheet Antimicrobial Peptide, Protegrin, Biochemistry, 41 (2002) 9852-9862]. Importantly, two simulations which started from different initial orientations of the peptide converged to the same final equilibrium orientation of the peptide relative to the bilayer. The kick-shaped conformation was observed only in one of the two simulations.     

Investigation of the ‘double cross’ splitting mechanism of single-crystal diamond under nanoindentation via molecular dynamics simulation

Elucidating the mechanical response of diamond is a difficult task due to its ultrahard nature. Here, we applied a molecular dynamics (MD) method to investigate the mechanical response of single-crystal diamond under nanoindentation. There was no obvious “pop in” phenomenon on the load–depth curve, and the elastic modulus deduced from the curve was 1128 GPa, which was similar to the value obtained from experimental measurements. Results from computed tomography (CT) and the coordination number showed that the distribution of the mismatched C atoms around the deformation zone took the form of a ‘double cross.’ The atoms around the indenter tip could be divided into two zones, a translation zone and a lattice distortion zone, based on their movements. Subsequent first-principles calculations revealed that the C-atom displacement barrier varied significantly with direction, which resulted in shear stress between the two zones and the formation of the double-cross splitting.

Open image in new window

Graphical Abstract The displacement of the atoms around the indenter tip

     

IB–LBM simulation of the haemocyte dynamics in a stenotic capillary

To study the behaviour of a haemocyte when crossing a stenotic capillary, the immersed boundary–lattice Boltzmann method was used to establish a quantitative analysis model. The haemocyte was assumed to be spherical and to have an elastic cell membrane, which can be driven by blood flow to adopt a highly deformable character. In the stenotic capillary, the spherical blood cell was stressed both by the flow and the wall dimension, and the cell shape was forced to be stretched to cross the stenosis. Our simulation investigated the haemocyte crossing process in detail. The velocity and pressure were anatomised to obtain information on how blood flows through a capillary and to estimate the degree of cell damage caused by excessive pressure. Quantitative velocity analysis results demonstrated that a large haemocyte crossing a small stenosis would have a noticeable effect on blood flow, while quantitative pressure distribution analysis results indicated that the crossing process would produce a special pressure distribution in the cell interior and to some extent a sudden change between the cell interior and the surrounding plasma.     

Role of invariant water molecules in retaining the active site geometry of β-lactamase: a molecular dynamics simulation study

Water molecules play a critical role in stabilising the three-dimensional architecture, dynamics and function of biological macromolecules. Comparative analysis of structurally similar proteins has shown that there are water molecules conserved in the same relative positions and make similar hydrogen bonds with proteins in all crystal structures. These invariant water molecules are essential for the maintenance of the native structure of proteins. The present study explores the role of invariant water molecules to maintain the active site geometry of β-lactamase enzyme. Thirteen crystal structures of class-A β-lactamase from Staphylococcus aureus have been used in this study. Molecular dynamics simulations of the protein structures were performed in hydrated as well as in dehydrated conditions. The analysis showed that significant changes occur in the active site geometry due to dehydration. These changes can be attributed to the removal of water molecules at the active site.     

Whales in the carbon cycle: can recovery remove carbon dioxide?

The great whales (baleen and sperm whales), through their massive size and wide distribution, influence ecosystem and carbon dynamics. Whales directly store carbon in their biomass and contribute to carbon export through sinking carcasses. Whale excreta may stimulate phytoplankton growth and capture atmospheric CO₂; such indirect pathways represent the greatest potential for whale-carbon sequestration but are poorly understood. We quantify the carbon values of whales while recognizing the numerous ecosystem, cultural, and moral motivations to protect them. We also propose a framework to quantify the economic value of whale carbon as populations change over time. Finally, we suggest research to address key unknowns (e.g., bioavailability of whale-derived nutrients to phytoplankton, species- and region-specific variability in whale carbon contributions).     

Molecular dynamics simulation of α-melanocyte stimulating hormone in a water-membrane model interface

The conformation of the tridecapeptide α-melanocyte stimulating hormone in the presence of a double water-membrane interface was studied by molecular dynamics simulation, using the computational package THOR. In this program the solvent is represented by a continuous medium with dielectric constant ɛ, and the interface between different media is simulated by a surface of discontinuity of the dielectric constant. The electrostatic image method was used to write down the terms, added to the force field, that describe the polarisation effects induced in the interface by the atomic charges. The program was further improved by the introduction of a second surface, parallel to the first one, to mimic the membrane. A conformational search using the software Prelude was employed to find an initial geometry for the peptide in water. The molecular dynamics simulation performed during 10 ns showed that the peptide structure is flexible in water, without stabilisation of any preferential conformation. In the presence of the model membrane, the peptide moved to the medium representing the interior of the membrane. Inside the low dielectric constant medium, the structure of the peptide showed a turn in the central sequence of amino acids and a packed conformation remained stabilised during more than 7.0 ns of simulation. Received: 27 November 1998 / Revised version: 11 March 1999 / Accepted: 8 April 1999     

Effects of ground water and harvest intensity on alkaline grassland ecosystem dynamics – a simulation study

A model for the alkaline grassland ecosystems, MAGE, was applied to plant communities dominated by three species. Field observations on two communities dominated respectively by Puccinellia tenuiflora and Suaeda corniculata were used to parameterize the model for multiple species interaction. The model behaves reasonably in following the seasonal variations of water content, soluble sodium cation and calcium cation in surface soil, as well as biomass of the plant communities.Simulations were run to investigate the effects of ground water quality, ground water table depth, maximum non-capillary porosity in surface soil and harvest intensity, on ecosystem dynamics. The results indicated that ground water sodium concentration and ground water table depth had primary control on soil alkalization and vegetation status. The improvement of soil conditions by vegetation is limited to an extent with moderate ground water depth and sodium concentration. Non-capillary pores are critical for vegetation to affect the soil alkalization/de-alkalization process, but the effect of non-capillary pores tends to saturate when maximum non-capillary porosity is greater than 0.1.     

What contributions to protein side-chain dynamics are probed by NMR experiments? A molecular dynamics simulation analysis

Molecular dynamics simulations of the structurally homologous proteins TNfn3 and FNfn10 have been used to investigate the contributions to side-chain dynamics measured by NMR relaxation experiments. The results reproduce the variation in core side-chain dynamics observed by NMR and highlight the relevance of anharmonic motion and transitions between local minima for explaining NMR side-chain order parameters. A method is described for calculating converged order parameters by use of replica exchange molecular dynamics in conjunction with an implicit solvent model. These simulations allow the influence of various factors, such as the flexibility of side-chains and their free volume, on the mobility to be tested by perturbing the system. Deletion mutations are found to have the largest effect on the more densely packed FNfn10. Some counterintuitive effects are seen, such as an increase in order parameters close to deletion mutation sites, but these can be rationalized in terms of direct interactions with the modified side-chains. A statistical analysis of published order parameters supports the conclusions drawn from the simulations.     

Molecular mechanism of transdermal co-delivery of interferon-alpha protein with gold nanoparticle – a molecular dynamics study

The transdermal route provides numerous advantages over conventional drug delivery routes. However, passive delivery of large molecules such as proteins through the skin is challenging due to its barrier function. Therefore, to design a successful formulation, molecular interaction of these proteins with constituent molecules present in the skin responsible for its barrier function, is necessary. In this study, we have shown through extensive computer simulations that the therapeutic protein, interferon alpha (INF), can be co-delivered through the skin using the gold nanoparticle. We carried out both steered (umbrella sampling) and unrestrained coarse-grained molecular dynamics simulation to show the molecular mechanism of absorption/permeation of protein on/through skin layer in the absence/presence of gold nanoparticle. According to the steered simulations, when INF was taken alone, the free energy minimum was observed at the head group of the skin layer, whereas, when co-delivered with AuNP, it was observed in the interior of the bilayer. Unrestrained simulations have also shown that INF was adsorbed on the skin lipid bilayer head group, while in presence of AuNP, it first complexed with the AuNP and then breached the barrier. The MD simulations thus established the transdermal delivery as a possible pathway for delivery of INF protein.     

Highlighting a π–π interaction: a protein modeling and molecular dynamics simulation study on Anopheles gambiae glutathione S-transferase 1-2

Cytosolic insect theta class glutathione S-transferases (GSTs) have not been studied completely and their physiological roles are unknown. A detailed understanding of Anopheles gambiae GST (Aggst1-2) requires an accurate structure, which has not yet been determined. A high quality model structure of Aggst1-2 was constructed using homology modeling and the ligand–protein complex was obtained by the docking method. Molecular dynamics (MD) simulations were carried out to study conformational changes and to calculate binding free energy. The results of MD simulation indicate that Aggst1-2 undergoes small conformational changes after ligands dock to the protein, which facilitate the catalytic reaction. An essential hydrogen bond was found between the sulfur atom of glutathione (GSH) and the hydrogen atom of hydroxyl group in Ser9, which was in good agreement with experimental data. A π–π interaction between Phe204 and CDNB ligand was also found. This interaction seems to be important in stabilization of the ligand. Further study of binding free energy decomposition revealed a van der Waals interaction between two ligands that may play a key role in nucleophilic addition reaction. This work will be a good starting point for further determination of the biological role of cytosolic insect theta class GSTs and will aid the design of structure-based inhibitors.