Investigation and optimisation of the drying of reduced glutathione by steered molecular dynamics simulation

The drying of reduced glutathione from a series of aqueous–ethanol binary solutions at 300 K (below human body temperature) and 330 K (above human body temperature) was investigated in detail by steered molecular simulation and an umbrella sampling method with the Gromacs software package and Gromos96(53a6) united atomic force field. The results show that electrostatic interactions between glutathione and solvent represent the main resistance to drying. When the aqueous solution was gradually changed to pure ethanol, the energy of electrostatic interaction between glutathione and solvent molecules increased by 445.088 kJ/mol, and the drying potential of mean force (PMF) free energy also fell by 253.040 kJ/mol. However, an increase in temperature from 300 to 330 K in the aqueous solution only results in an increase of 23.013 kJ/mol in electrostatic interaction energy and a decrease of 34.956 kJ/mol in drying PMF free energy. Furthermore, we show that hydrogen bonding is the major form of electrostatic interaction involved, and directly affects the drying of glutathione. Therefore, choosing water-miscible solvents that minimise hydrogen-bond formation with glutathione will enhance its drying rate, and this is likely to be more efficient than increasing the temperature of the process. Thus, a power-saving technology can be used to produce the high bioactivity medicines.     

Molecular dynamics simulation of d-Benzedrine transmitting through molecular channels within D3R

Dex-Benzedrine (known as d-Benzedrine or SAT) acts in dopamine receptors of central nerve cell system. In clinic, SAT is used to treat a variety of diseases; meanwhile, it has dependence and addiction. In order to investigate the pharmacology and addiction mechanisms of SAT as a medicine, in this paper, we have studied the structure of D₃R complex protein with SAT, and based on which, using potential mean force with umbrella samplings and the simulated phospholipid bilayer membrane (or POPC bilayer membrane), the molecular dynamics simulation was performed to obtain free energy changes upon the trajectories for SAT moving along the molecular channels within D₃R. The free energy change for SAT transmitting toward the outside of cell along the functional molecular channel within D₃R is 83.5 kJ mol^?1. The change of free energy for SAT to permeate into the POPC bilayer membrane along the protective molecular channel within D₃R is 87.7 kJ mol^?1. Our previous work gave that the free energy for Levo-Benzedrine (RAT) transmitting toward the outside of cell along the functional molecular channel within D₃R is 91.4 kJ mol^?1, while it is 117.7 kJ mol^?1 for RAT to permeate into the POPC bilayer membrane along the protective molecular channel within D₃R. The values of free energy suggest that SAT relatively prefers likely to pass through the functional molecular channel within D₃R for increasing the release of dopamine molecules resulting in a variety of functional effects for SAT. The obtained results show that the pharmacology and addiction mechanisms of SAT as a drug are closely related to the molecular dynamics and mechanism for SAT transmitting along molecular channels within D₃R.     

Computational revelation of binding mechanisms of inhibitors to endocellular protein tyrosine phosphatase 1B using molecular dynamics simulations

Endocellular protein tyrosine phosphatase 1B (PTP1B) is one of the most promising target for designing and developing drugs to cure type-II diabetes and obesity. Molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy methods were applied to study binding differences of three inhibitors (ID: 901, 941, and 968) to PTP1B, the calculated results show that the inhibitor 901 has the strongest binding ability to PTP1B among the current inhibitors. Principal component (PC) analysis was also carried out to investigate the conformational change of PTP1B, and the results indicate that the associations of inhibitors with PTP1B generate a significant effect on the motion of the WPD-loop. Free energy decomposition method was applied to study the contributions of individual residues to inhibitor bindings, it is found that three inhibitors can generate hydrogen bonding interactions and hydrophobic interactions with different residues of PTP1B, which provide important forces for associations of inhibitors with PTP1B. This research is expected to give a meaningfully theoretical guidance to design and develop of effective drugs curing type-II diabetes and obesity.     

Study on mechanical properties of polyethylene with chain branching in atomic scale by molecular dynamics simulation

The effects of length and content of chain branching on the mechanical properties of polyethylene (PE) in atomic scale were examined by molecular dynamics (MD) simulations. Methyl-, ethyl- and butyl-groups were adopted as branched chains to distribute along PE backbones. Plastic flow deformation was captured by providing a uniaxial tensile loading at a given strain rate, which shows the characteristic of rate dependence. Current results are in reasonable agreements with existing experimental data. The statistical results show that the longer length of chain branching induces lower equilibrium density and higher yield strength of branched PE. In addition, higher content of chain branching brings higher equilibrium density and lower yield strength of branched PE. It is assumed that the distribution of dihedral angles influences the deformation of PE definitely. The non-bond interactions contribute to the load-bearing capacity of PE largely. Branched PE shows big differences on mechanical behaviours comparing with the linear one. Chain branching distribution also greatly affects the performance of PE, which needs a further discussion.     

The mechanism of restructuring of surfactant monolayer on mica surface in aqueous solution: molecular dynamics simulation

Molecular dynamics simulation was employed to investigate the restructuring process of CTAB monolayer at mica/water interface. The reversing process of CTAB monolayer was exploited by diffusion of water molecules, reversing of CTAB molecules with time evolution and restructuring of the surfactant monolayer. The results showed that bromide ions around surfactant head groups diffused into bulk water readily due to the electrostatic repulsion caused by negatively charged mica surface. Meanwhile, because of the electrostatic attraction between water molecules and mica surface, part of water molecules can penetrate the surfactant monolayer to form water channel which bridges bulk water and mica surface. The monolayer structure was disturbed by diffusion of bromide ions and formation of water channel. Few of the head groups of surfactants tended to reverse and enter into aqueous solution. The number of reversed surfactant molecules increased with time evolution. The monolayer restructured into bilayer structure gradually. Finally, a cylindrical aggregate was obtained.     

The molecular simulation of the miscibility,mechanical properties and physical cross-linking behavior of the poly(vinyl alcohol)/poly(acrylic acid) composited membranes

The miscibility and mechanical properties of poly vinyl alcohol (PVA) and poly acrylic acid (PAA)-composited membranes were studied with molecular simulation. The Flory–Huggins parameters (δ) were calculated to prove the good miscibility of PVA and PAA. The radial distribution functions of hydroxyl and carboxyl atoms and the average number of H-bonds were observed to indicate the degree of physical cross-linking between PVA and PAA. The influences of intermolecular physical cross-linking on the glass transition temperature and mechanical properties were estimated. The results revealed that the PVA/PAA membrane with a composition of 2:3 has the best plastic properties, which exhibits a good application value. All of the simulated results showed good agreement with the experimental data. It indicates that the method presented in this work has a promising application prospect.     

Evaluating the suitability of RNA intervention mechanism exerted by some flavonoid molecules against dengue virus MTase RNA capping site: a molecular docking,molecular dynamics simulation,and binding free energy study

Abstract

Dengue virus (DENV) is one of the most dangerous mosquito-borne human pathogens known to the mankind. Currently, no vaccines or standard therapy is avaliable to treate DENV infection. This makes the drug development against DENV more significant and challenging. The MTase domain of DENV RNA RdRp NS5 is a promising drug target, because this domain hosts the RNA capping process of DENV RNA to escape from human immune system. In the present study, we have analysed the RNA intervention mechanism exerted by flavoniod molecules against NS5 MTase RNA capping site by using molecular docking, molecular dynamics simulation and the binding free energy calculations. The results from the docking analysis confirmed that the RNA intervention mecanism is exerted by the quercetagetin (QGN) molecule with all necessary intermolecular interactions and high binding affinity. Notably, QGN forms strong hydrogen bonding interactions with Asn18, Leu20 and Ser150 residues and π???π stacking interaction with Phe25 residue. The apo and QGN bound NS5 MTase and QGN-NS5 MTase complex were used for MD simulation. The results of MD simulation reveal that the RMSD and RMSF values of QGN-MTase complex have increased on comparing the apo protein due to the effect of ligand binding. The binding free energy calulation includes prediction of total binding free energy of ligand-protein complex and per-residue free energy decomposition. The QGN binding to NS5 MTase affects it’s native motion, this result is found from Principal component analysis.

Communicated by Ramaswamy H. Sarma     

Interaction investigations of HipA binding to HipB dimer and HipB dimer + DNA complex: a molecular dynamics simulation study

We carried out molecular dynamics simulations and free energy calculations for a series of ternary and diplex models for the HipA protein, HipB dimer, and DNA molecule to address the mechanism of HipA sequestration and the binding order of events from apo HipB/HipA to 2HipA + HipB dimer + DNA complex. The results revealed that the combination of DNA with the HipB dimer is energetically favorable for the combination of HipB dimer with HipA protein. The binding of DNA to HipB dimer induces a long‐range allosteric communication from the HipB₂‐DNA interface to the HipA–HipB₂ interface, which involves the closeness of α1 helices of HipB dimer to HipA protein and formations of extra hydrogen bonds in the HipA–HipB₂ interface through the extension of α2/3 helices in the HipB dimer. These simulated results suggested that the DNA molecule, as a regulative media, modulates the HipB dimer conformation, consequently increasing the interactions of HipB dimer with the HipA proteins, which explains the mechanism of HipA sequestration reported by the previous experiment. Simultaneously, these simulations also explored that the thermodynamic binding order in a simulated physiological environment, that is, the HipB dimer first bind to DNA to form HipB dimer + DNA complex, then capturing strongly the HipA proteins to form a ternary complex, 2HipA + HipB dimer + DNA, for sequestrating HipA in the nucleoid. Copyright © 2013 John Wiley & Sons, Ltd.     

Comparative molecular dynamics simulation analysis of G20 and C92 mutations in c-di-GMP I riboswitch and the wild type with docked c-di-GMP ligand

Riboswitch, a bacterial regulatory RNA consists of an aptamer (specific ligand binding unit) and an expression platform (gene expression modulation unit), which act as a potential drug target as it regulates critical genes. Therefore, it is of interest to glean information on the binding of c-di-GMP ligand to mutated conserved G20 and C92 residues of cyclic diguanosine monophosphate I (c-di-GMP I) riboswitch using molecular dynamics simulation. The result shows that the binding energy of wild/native type riboswitch-ligand complex (3IRW) is lower than the mutant complexes suggesting that the binding affinity for c-di-GMP ligand decreases in case of mutant riboswitches. The hydrogen bonding interactions analysis also showed a high number of hydrogen bonds formation in the wild type riboswitch-ligand complex as compared to the mutant complexes illustrating stronger interaction of ligand to wild type riboswitch than the mutants. The simulation result shows that the mutations affected riboswitch-ligand interactions. The residues G14, G21, C46, A47, and U92 were identified as the key residues which contributed effectively to the binding of c-di-GMP I riboswitch with the natural ligand.     

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.     

Screening and molecular dynamics simulation of compounds inhibiting MurB enzyme of drug-resistant Mycobacterium tuberculosis: An in-silico approach

Mycobacterium tuberculosis (MTB) is becoming more and more resistant to drugs and it is a common problem, making current antimicrobials ineffective and highlighting the need for new TB drugs. One of the promising targets for treating MTB is MurB enzymes. This study aimed to identify potential inhibitors of MurB enzymes in M. tuberculosis, as drug resistance among MTB is a significant problem. Attempts are being made to conduct a virtual screening of 30,417 compounds, and thirty-two compounds were chosen for further analysis based on their binding conformations. The selected compounds were assessed for their drug-likeness, pharmacokinetics, and physiochemical characteristics, and seven compounds with binding energy lower than flavin (FAD) were identified. Further, molecular dynamics simulation analysis of these seven compounds found that four of them, namely DB12983, DB15688, ZINC084726167, and ZINC254071113 formed stable complexes with the MurB binding site, exhibiting promising inhibitory activity. These compounds have not been mentioned in any other study, indicating their novelty. The study suggests that these four compounds could be promising candidates for treating MTB, but their effectiveness needs to be validated through in vitro and in vivo experiments. Overall, the findings of this study provide new insight into potential drug targets and candidates for combating drug-resistant MTB.     

3D-QSAR-aided design of potent c-Met inhibitors using molecular dynamics simulation and binding free energy calculation

Abstract

Mesenchymal-epithelial transition factor (c-Met) is a member of receptor tyrosine kinase. It involves in various cellular signaling pathways which includes proliferation, motility, migration, and invasion. Over-expression of c-Met has been reported in various cancers. Hence, it is an ideal therapeutic target for cancer. The main objective of the study is to identify crucial residues involved in the inhibition of c-Met kinase and to design a series of potent imidazo [4,5-b] pyrazine derivatives as c-Met inhibitors. Docking was used to identify important active site residues involved in the inhibition of c-Met kinase which was further validated by 100 ns of molecular dynamics simulation and free energy calculation using molecular mechanics generalized born surface area. Furthermore, binding energy decomposition identified that residues Tyr1230, Met1211, Asp1222, Tyr1159, Met1160, Val1092, Ala1108, and Leu1157 contributed favorably to the binding stability of compound 32. Receptor-guided Comparative Molecular Field Analysis (CoMFA) (q² = 0.751, NOC = 6, r² = 0.933) and Comparative Molecular Similarity Indices Analysis (COMSIA) (q² = 0.744, NOC = 6, r² = 0.950) models were generated based on the docked conformation of the most active compound 32. The robustness of these models was tested using various validation techniques and found to be predictive. The results of CoMFA and CoMSIA contour maps exposed the regions favorable to enhance the activity. Based on this information, 27 novel c-Met inhibitors were designed. These designed compounds exhibited potent activity than the most active compound of the existing dataset.

Communicated by Ramaswamy H. Sarma     

Ligand orientation of human neuroglobin obtained from solution NMR and molecular dynamics simulation as compared with X-ray crystallography

Neuroglobin, a new member of hemoprotein family, can reversibly bind oxygen and take part in many biological processes such as enzymatic reaction, signal transduction and the mitochondria function. Different from myoglobin and hemoglobin, it has a hexacoordinated heme environment, with histidyl imidazole of proximal His⁹⁶(F8) and distal His⁶⁴(E7) directly bound to the metal ion. In the present work, solution ¹H NMR spectroscopy was employed to investigate the electronic structure of heme center of wild-type met-human neuroglobin. The resonances of heme protons and key residues in the heme pocket were assigned. Two heme orientations resulting from a 180° rotation about the α-γ-meso axis with a population ratio about 2:1 were observed. Then the ¹H NMR chemical shifts of the ferriheme methyl groups were used to predict orientations of the axial ligand. The obtained axial ligand plane angle φ is consistent with that from the molecular dynamics simulation but not with those from the crystal data. Compared with mouse neuroglobin, the obtained average ligand orientation of human neuroglobin reflects the changeability of heme environment for the Ngb family.     

A computational analysis of pyrazole-based inhibitors binding to Hsp90 using molecular dynamics simulation and the MM-GBSA method

Owing to the key role of heat-shock protein 90 (Hsp90) in the evolution, development and disease pathogenesis of cancer, it has been an important target for anti-cancer chemotherapy over the years. A five-nanosecond molecular dynamics simulation combined with the calculation of the binding free energy was carried out to investigate the binding mechanisms of three Hsp90 inhibitors 4BH, 2E1 and 2D9 to Hsp90. The binding free energy of each complex was computed using the molecular mechanics–generalised Born surface area method. Detailed binding free energies between each inhibitor and residues of Hsp90 were calculated using a per-residue basis decomposition method. The detailed inhibitor–residue interaction provides insights into binding mechanisms and in-depth understanding of the structure–affinity relationship. This study suggests that van der Waals energy is primarily responsible for driving the binding of the inhibitors to Hsp90, and the three inhibitors bind to Hsp90 in a similar binding mode. However, a substituent in 2D9 leads to higher binding free energy than the other two inhibitors. These data may assist in designing new potent drugs to combat cancer.     

Crystallographic structure versus homology model: a case study of molecular dynamics simulation of human and zebrafish histone deacetylase 10

Abstract

Histone deacetylase (HDAC) 10 has been implicated in the pathology of various cancers and neurodegenerative disorders, making the discovery of novel inhibitors of the isoform an important endeavor. However, the unavailability of crystallographic structure of human HDAC10 (hHDAC10) hinders structure-based drug design effort. Previously, we reported the homology modeled structure of human HDAC10 built using the crystallographic structure of Danio rerio (zebrafish) HDAC10 (zHDAC10) (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) as a template. Here, in continuation with our study, both hHDAC10 and zHDAC10, and their respective complexes with trichostatin A (TSA), quisinostat, and the native ligand (in 5TD7), 7-[(3-aminopropyl)amino]-1,1,1-trifluoroheptane-2,2-diol (PDB ID; FKS) were submitted to 100?ns-long unrestrained molecular dynamics (MD) simulations. Comparative analyses of the MD trajectories revealed that zHDAC10 and its complexes displayed higher stability than hHDAC10 and its corresponding complexes over time. Nonetheless, docking of active and inactive set molecules revealed that more reliable conformations of hHDAC10 could be obtained at an extended time period. This study may shed more light on the reliability of hHDAC10 modeled structure for use in selective inhibitor design.

Communicated by Ramaswamy H. Sarma     

Identification of Plasmodium falciparum apicoplast-targeted tRNA-guanine transglycosylase and its potential inhibitors using comparative genomics,molecular modelling,docking and simulation studies

tRNA modifications play an important role in the proper folding of tRNA and thereby determine its functionality as an adaptor molecule. Notwithstanding the centrality of this basic process in translation, a major gap in the genomics of Plasmodium falciparum is unambiguous identification of enzymes catalysing the various tRNA modifications. In this study, tRNA-modifying enzymes of P. falciparum were annotated using homology-based approach. Based on the presence of these identified enzymes, the modifications were compared with those of prokaryotic and eukaryotic organisms. Through sequence comparison and phylogenetic analysis, we have identified P. falciparum apicoplast tRNA-guanine 34 transglycosylase (TGT, EC: 2.4.2.29), which shows evidence of its prokaryotic origin. The docking analysis of the modelled TGT structures revealed that binding of quinazolinone derivatives is more favourable with P. falciparum apicoplast TGT as compared to human TGT. Molecular dynamic simulation and molecular mechanics/generalized Born surface area analysis of the complex confirmed the greater binding affinity of the ligand in the binding pocket of P. falciparum TGT protein. Further, evolutionary patterning analysis identified the amino acids of P. falciparum apicoplast TGT that are under purifying selection pressure and hence can be good inhibitor-targeting sites. Based on these computational studies, we suggest that P. falciparum apicoplast tRNA-guanine 34 transglycosylase can be a promising drug target.     

Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild–type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the ?uctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.     

Sensitivity of the aggregation behaviour of asphaltenes to molecular weight and structure using molecular dynamics

Asphaltenes are heavy crude oil compounds, defined as soluble in toluene and precipitating in alkanes. To understand the relation between asphaltene structure and aggregation, we perform equilibrium molecular dynamics with Large-scale Atomic and Molecular Massively Parallel Software (LAMMPS), using the atomistic force field PCFF+ in the MedeA® environment. The following three molecular models are considered: the continental model (1350 g/mol) that has a large polyaromatic core and long alkyl chains, the island model (780 g/mol) that has a smaller polyaromatic unit and shorter chains and the archipelago model (1350 g/mol) that has three polyaromatic nuclei bridged with alkyl chains. The aggregation in a given solvent is monitored by visualising solvent-free configurations over 15 ns trajectories at 350 K. Nanoaggregates are characterised by stacked polyaromatic units separated by 0.33–0.4 nm. Irreversible aggregation is found with the continental model in both solvents. Aggregation of the island model is significant in n-heptane and low in toluene. The archipelago model does not aggregate significantly. Our results confirm that the island model is a reasonable average model of asphaltenes [Headen TF, Boek ES, Skipper NT. Energy Fuels 2009;23:1220–1229]. The open structure of nanoaggregates and the limited number of stacked molecules are also in agreement with previous interpretations of experimental data [Fenistein D. et al. Langmuir 1998;14:1013–1020].