A molecular dynamics investigation of buckling behaviour of hydrogenated graphene

Molecular dynamics simulations have been performed to characterise the stability behaviour of graphene nanoribbons having different hydrogen coverage, subject to a uniaxial compressive load. The temperature is set at a very low value to circumvent the contribution of thermal agitations. The results show that hydrogen coverage promotes to a rapid drop in the strain of buckling onset due to the effects of easy rotation of newly unsupported sp³ bonds. Furthermore, we have also found a critical value of the hydrogen adsorption above which the declining trend in the stability behaviour of hydrogenated graphene nanoribbons is reversed.     

Interaction mechanism of insulin with ZnO nanoparticles by replica exchange molecular dynamics simulation

The interaction of ZnO nanoparticles with biological molecules such as proteins is one of the most important and challenging problems in molecular biology. Molecular dynamics (MD) simulations are useful technique for understating the mechanism of various interactions of proteins and nanoparticles. In the present work, the interaction mechanism of insulin with ZnO nanoparticles was studied. Simulation methods including MD and replica exchange molecular dynamics (REMD) and their conditions were surveyed. According to the results obtained by REMD simulation, it was found that insulin interacts with ZnO nanoparticle surface via its polar and charged amino acids. Unfolding insulin at ZnO nanoparticle surface, the terminal parts of its chains play the main role. Due to the linkage between chain of insulin and chain of disulfide bonds, opposite directional movements of N terminal part of chain A (toward nanoparticle surface) and N termini of chain B (toward solution) make insulin unfolding. In unfolding of insulin at this condition, its helix structures convert to random coils at terminal parts chains.     

Dynamic injection behaviour of carbon monoxide nano-jets: analysis based on molecular dynamics simulation

Molecular dynamics simulation was used to characterise the dynamical injection behaviour of CO through a gold nano-injector with a Gr coating. We also varied the nozzle outlet size, system temperature, and extrusion velocity to elucidate their influence on the flow patterns, injection pressure, and flow rate of the CO nano-jets. Simulation results revealed the following important findings. (1) At 100?K, the liquefaction of a CO jet led to a wider spray angle (ø_s?=?84～96°) and allowed molecules to attach to the Gr layer, resulting in agglomeration at the orifice. (2) At 500?K and 55.824?m/s, the nebulisation of the CO nano-jet was induced, which produced a narrower spray angle (ø_s?=?47°). (3) The flow rate of CO molecules was essentially linear under the following conditions: low extrusion velocity (≤13.956?m/s), large orifice (d?=?1.5?nm), and high system temperature (≥300?K). (4) Due to the compressibility of CO molecules, the pressure inside the chamber under a high extrusion speed (≥27.912?m/s) presented a sharp increase in the middle and final extrusion stages. A delay in the pressure increase enabled the liquefaction of the extruded CO molecules, resulting in an unstable flow rate.     

Adsorption behaviors of hydrogen sulphide on Au(110) nanoslit array surfaces using molecular dynamics simulations

The adsorption behaviour of gas molecules on detector surfaces has a profound influence on the sensitivity of the detector. For this reason, this study used molecular dynamics simulation to explore the dynamic adsorption behaviour of hydrogen sulphide (H₂S) molecules on various types of Au surfaces, including a planar Au(1?1?0) structure and three types of slit array structures. The influence of system temperature, adsorbate concentration and the slit width of nanoarrays on diffusivity, average adsorption energy and static adsorption amount were systematically examined. Simulation results indicate that the self-diffusivity of the adsorbate molecules increases with temperature but decreases with adsorbate concentration. At low concentrations (~3 mol/L), each type of Au(1?1?0) surface structure shows good capacity to adsorb all H₂S molecules. With increasing concentration at 6.5 mol/L, the high concentration leads to adsorption saturation and many free H₂S molecules in the planar Au(1?1?0) structure. Moreover, desorption also begins to appear on the planar structures at a temperature of 300 K (at 6.5 mol/L). The simulation results indicate that the columnar array structures with a slit width ≥5.76 Å allow molecules to swiftly spread into the slits and provide more stable adsorption sites (i.e. with a higher adsorption energy), which can effectively address the issues of high-temperature desorption and adsorption saturation. Particularly at low temperatures (≤100 K), slit structures presented a level of static adsorption of H₂S that was 30% to 35 higher than that of planar structures.     

Exploring single-domain antibody thermostability by molecular dynamics simulation

Abstract

Single-domain antibodies also known as nanobodies are recombinant antigen-binding domains that correspond to the heavy-chain variable region of camelid antibodies. Previous experimental studies showed that the nanobodies have stable and active structures at high temperatures. In this study, the thermal stability and dynamics of nanobodies have been studied by employing molecular dynamics simulation at different temperatures. Variations in root mean square deviation, native contacts, and solvent-accessible surface area of the nanobodies during the simulation were calculated to analyze the effect of different temperatures on the overall conformation of the nanobody. Then, the thermostability mechanism of this protein was studied through calculation of dynamic cross-correlation matrix, principal component analyses, native contact analyses, and root mean square fluctuation. Our results manifest that the side chain conformation of some residues in the complementarity-determining region 3 (CDR3) and also the interaction between α-helix region of CDR3 and framework2 play a critical role to stabilize the protein at a high temperature.

Communicated by Ramaswamy H. Sarma     

Effect of chirality,length and diameter of carbon nanotubes on the adsorption of 20 amino acids: a molecular dynamics simulation study

We carried out molecular dynamics simulations to study the adsorption of all the 20 amino acids (AAs; aromatic, polar and non-polar) on the surface of chiral, zigzag and armchair single-walled carbon nanotubes. The adsorption was occurring in all systems. In the aromatic AAs, the π–π stacking and the semi-hydrogen bond formation cause a strong interaction with the carbon nanotubes (CNTs). We also investigated the chirality, length and diameter dependencies on adsorption energies. We found that all AAs have more tendency to adsorption on the chiral and zigzag CNTs over the armchair. The results show that increasing both the diameter and the length causes the enhancement of the adsorption energy. But, the effect of the length is more than of the diameter. For example, the adsorption energy of Trp on the surface of CNT (4,1), with 2 nm length, is 20.4 kcal/mol. When the length of CNT becomes twice, the adsorption energy increases by 24 ± 0.3%. But by doubling the diameter, the adsorption energy increased only by 9.8 ± 0.25%.     

Neon atoms oscillating inside carbon and boron nitride nanotubes: a fully atomistic molecular dynamics investigation

In the present work, based on extensive fully atomistic molecular dynamics simulations, we discuss the dynamics of neon atoms oscillating inside (5,5) single-walled carbon nanotubes (CNTs) and boron nitride nanotubes (BNNTs). Our results show that sustained high-frequency oscillatory regimes are possible for a large range of temperatures. Our results also show that the general features of the oscillations are quite similar to those observed in CNT and BNNT, in contrast with some speculations in previous works in the literature about the importance of broken symmetry and chirality exhibited by BNNTs.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Prediction of chlortetracycline adsorption on the Fe3O4 nanoparticle using molecular dynamics simulation

Abstract

Molecular dynamics (MD) simulation was applied to investigate the adsorption mechanism of chlortetracycline (CTC) antibiotic molecule as the aqueous pollutant on the Fe₃O₄ nanoparticle (NP). Two different NP sizes with a diameter of about 1.4?nm and 3.5?nm were selected. Initially, the stability of both NPs in water was investigated by calculating radial distribution function curves of NP atoms. Simulation results confirmed the stable crystallographic structures of both NPs. However, small NP induce greater structural stabilization. Then, CTC molecules were adsorbed on NPs surface in various pollutant concentrations. Electrostatic and hydrogen bond were the major types of interactions between CTC molecules and the adsorbent surface. CTC molecules formed a complex with NP surface from their amine side chains; while they were parallel to each other in their aromatic rings and π-π bond between two CTC molecules was formed. Diffusion rate of CTC molecules could predict the adsorption mechanism. At lower concentration of CTC, CTC molecules tend to adsorb on the NP surface. At these concentrations, the diffusion rate of CTC was high. By increasing the CTC concentration, the pollutant agglomeration was enhanced which decreased the diffusion rate. At this time, the surface of NP was saturated. In addition, the results of isotherm curves showed that CTC adsorption on small NPs could be defined with both Langmuir and Freundlich isotherm models, while Freundlich isotherm model was more appropriate for larger NPs. In conclusion, observations confirmed that MD simulation could successfully predict the behavior of CTC adsorption on the Fe₃O₄ NP surface.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics simulation of lactate dehydrogenase adsorption onto pristine and carboxylic-functionalized graphene

ABSTRACT

Lactate dehydrogenase (LDH) is a tetrameric enzyme which is composed of two subunits known as LDHA and LDHB, which are encoded by the LDHA and LDHB genes respectively. LDH catalyses the last step in anaerobic glycolysis through the reversible conversion of pyruvate to lactate via coupled oxidation of NADH cofactor. The LDHA plays an important regulatory role in anaerobic glycolysis, by catalysing the final step of the process. Therefore, it is likely that increases in the expression level of LDHA in cancer cells could facilitate the efficiency of anaerobic glycolysis. Measuring the level of serum LDHA is a key step in the diagnosis of many cancer types. In this study, the adsorption, stability, and dynamics of LDHA on the surface of pristine graphene (PG) and carboxylated graphene (COOH-Graphene) were investigated using its molecular dynamics simulation. Variations in root mean square deviation, root mean square fluctuation, solvent accessible surface area and adsorption energy of the LDHA during the simulation were calculated to analyse the effect of PG and COOH-Graphene on the overall conformation of LDHA. Results showed that the adsorption of LDHA on COOH-Graphene is mostly mediated by electrostatic interactions, whereas on the PG, both Van der Waals and π-π interactions are prominent.     

Molecular dynamics simulation,binding free energy calculation and molecular docking of human D-amino acid oxidase (DAAO) with its inhibitors

Schizophrenia is a mental illness; most affected people live in developing countries, and neither appropriate treatment nor commercial drugs are currently available. One possibility is to inhibit human-d-amino acid oxidase (h-DAAO). In this study, molecular dynamic simulations of the monomer, dimer and tetramer forms of h-DAAO complexed with the inhibitor 3-hydroxyquinolin-2(1H)-one(2) were performed. Seven residues, Leu51, Gln53, Leu215, Tyr228, Ile230, Arg283 and Gly313, were identified as essential for interacting with the inhibitor. Molecular docking of h-DAAO with pyrrole, quinoline and kojic acid derivatives, representing 69 known or potential h-DAAO inhibitors, was also performed. The results indicated that the activity of the inhibitor can be improved by modifying the compounds to have a substituent group capable of interacting with the side chain of Tyr228. Van der Waals interactions of the inhibitor with the hydrophobic pocket of h-DAAO and electrostatic interactions or H-bonds with Arg283 and Gly313 were important elements in determining the efficiency of the inhibitor. These results provide information on the interaction between h-DAAO and its inhibitors at the molecular level and can aid in the design of novel inhibitors against h-DAAO for new drug development in the treatment of schizophrenia.     

Structure and stability of the insulin dimer investigated by molecular dynamics simulation

Molecular dynamics simulation indicates that the dynamical behaviour of the insulin dimer is asymmetric. Atomic level knowledge of the interaction modes and protein conformation in the solvation state identifies dynamical structures, held by hydrogen bonds that stabilize, mainly in one monomer, the interaction between the chains. Dynamic cross-correlation analysis shows that the two insulin monomers behave asymmetrically and are almost independent. Solvation energy, calculated to evaluate the contribute of each interface residue to the dimer association pattern, well compares with the experimental association state found in protein mutants indicating that this parameter is an important factor to explain the association properties of mutated insulin dimers.     

Understanding interaction and dynamics of water molecules in the epoxy via molecular dynamics simulation

The robust structural integrity of the epoxy plays an important role in ensuring the long-term service life of its applications, which is affected by the absorbed moisture. In order to understand the mechanism of the moisture effect, the knowledge of the interaction and dynamics of the water molecules inside the epoxy is of great interest. Molecular dynamics simulation is used in this work to investigate the structure and bonding behaviour of the water molecules in the highly cross-linked epoxy network. When the moisture concentration is low, the water molecules are well dispersed in the cross-linked structure and located in the vicinity of the epoxy functional groups, which predominantly form the hydrogen bond (H-bond) with the epoxy network, resulting in the low water mobility in the epoxy. At the high concentration, the water favourably forms the large cluster due to the predominant water–water H-bond interaction, and the water molecules diffuse primarily inside the cluster, which leads to the high water mobility and the accelerated H-bond dynamics. The variation of the bonding behaviour and dynamics of the water molecules reported here could be exploited to understand the material change and predict the long-term performance of the epoxy-based products during the intended service life.     

Interfacial adhesion properties of graphene sheet on nanoscale corrugated surface: a molecular dynamics simulation study

Adhesive contacts between graphene sheet (GS) and corrugated substrates made of an ordered array of atomic pillars with variable geometries were investigated by molecular dynamics simulations. Depending on the height and interval distance of the pillars, GS can conformably coat the surface, partially adhere, or remain flat on top of the pillars. The relationship between the geometries of the pillar and the final adhesion configurations of GS was partially established. A critical adsorption energy was determined to achieve stable adsorption configuration of GS on corrugated substrates made of ordered pillar arrays. Besides the geometries of pillars, the effects of initial coating angle of GS were also considered as an important factor that affects the final adsorption configuration. We observed two interesting morphologies of GS, ‘I shape’ and ‘L shape’, which were determined by the initial coating angles.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Topological and dynamical properties of Azurin anchored to a gold substrate as investigated by molecular dynamics simulation

A classical molecular dynamics study of the electron transfer protein azurin, covalently bound to a gold substrate through its native disulphide group, is carried out at full hydration. With the aim of investigating the effects on the protein structure and dynamics as induced by the presence of an electric field, simulations are performed on neutral, positively and negatively charged substrates. A number of parameters, such as the average structure, the root mean square deviations and fluctuations, the intraprotein hydrogen bonds and solvent accessible surface of the protein, are monitored during 10 ns of run. The orientation, the height and the lateral size of the protein, with respect to the substrate are evaluated and compared with the experimental data obtained by scanning probe nanoscopies. The electron transfer properties between the copper redox center and the disulphide bridge bound to the substrate are investigated and briefly discussed.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.