Study on the β Peak of χ″ for a Two-Dimensional Supercooled Fluid State Via Molecular Dynamics Simulation

Abstract

We have computed the generalized susceptibility (imaginary part of χ) of two-dimensional supercooled fluid states through molecular dynamics simulations with different system sizes. It is found that when the temperature of the system is higher than T_c (critical temperature) there is no significant difference between the large (10,000-particle) and small (100-particle) systems in the χ obtained, but for much lower temperatures the spectra of the χ″ obtained reveal a remarkable system-size dependent behavior such as the β peak becomes much broader for the large system. This phenomenon is physically understood in terms of the strongly correlated motion of atoms extending over a wide range of space of the system. We will discuss in this work about the relaxations on two-dimensional supercooled fluids and their microscopic origins as well as their system-size dependence.     

Molecular dynamics simulation of proteins under high pressure: Structure,function and thermodynamics

BackgroundMolecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.Scope of reviewFirst, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.Major conclusionsRecent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.General significanceMD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.     

Combined pharmacophore modeling, 3D-QSAR and docking studies to identify novel HDAC inhibitors using drug repurposing

Abstract

Histone deacetylases (HDACs), a critical family of epigenetic enzymes, has emerged as a promising target for antitumor drugs. Here, we describe our protocol of virtual screening in identification of novel potential HDAC inhibitors through pharmacophore modeling, 3D-QSAR, molecular docking and molecular dynamics (MD) simulation. Considering the limitation of current virtual screening works, drug repurposing strategy was applied to discover druggable HDAC inhibitor. The ligand-based pharmacophore and 3D-QSAR models were established, and their reliability was validated by different methods. Then, the DrugBank database was screened, followed by molecular docking. MD simulation (100?ns) was performed to further study the stability of ligand binding modes. Finally, results indicated the hit DB03889 with high in silico inhibitory potency was suitable for further experimental analysis.

Communicated by Ramaswamy H. Sarma     

Structural Transformations of Ice at High Pressures Via Molecular Dynamics Simulations

Abstract

A simple classical model is used for the study of the structural transformations of ice under high pressures, such as ice VIII to VII and X, via classical molecular dynamics (MD) simulation. In the present MD simulation, pair potentials of a simple form between pair of atoms and a thee-body potential representing the H-O-H angle dependence, originally developed by Kawamura et al., were used. Starting with a stable ice VIII at low pressure and low temperature, we have carried out two different MD runs, one with increasing pressure keeping the temperature constant (simulation I) and the other with increasing temperature under constant pressure (simulation II). From these MD simulations we have obtained the structural transformations from ice VIII to VII for both simulations; the former was finally transformed into ice X for the simulation I. The present results are compatible with recent experiments on high pressure ices.     

Induced fit docking,free energy calculation and molecular dynamics studies on Mycobacterium tuberculosis alanine racemase inhibitor

Alar, a Pyridoxal 5′-phosphate (PLP)-dependent bacterial enzyme is responsible for the racemisation of L-alanine into D-alanine which is essential for the peptidoglycan biosynthesis in both Gram-positive and Gram-negative bacteria. In the present study, we performed induced fit docking, binding free energy calculation and molecular dynamics simulation to elucidate the Alar inhibition potential of 1,2,4-thiadiazolidine-3,5-dione-based inhibitor 1. The inhibitor binds to the hydrophobic groove of Alar and the binding was found to be stable throughout 20-ns MD simulation. Induced fit docking result showed that Lys42, Tyr46, Tyr175 and Tyr364 residues are primarily responsible for the stabilisation of inhibitor–protein complex. Further, high negative van der Waals binding free energy value of –38.88 kcal/mol, indicated it as the main driving force for the inhibitor binding. Based on the information obtained from this study, we designed few molecules as potent Alar inhibitor. In order to gain structural insight and to validate the stability of complex, we performed 20-ns MD simulation of the designed molecule D1. Results obtained from this study can be used for the design of M. tuberculosis Alar potent inhibitors lacking affinity for the co-factor PLP.     

Ab Initio Molecular Dynamics for Simple Liquid Metals Based on the Hypernetted-Chain Approximation

Abstract

We present an ab initio molecular dynamics (MD) method for simple liquid metals based on the quantal hypernetted-chain (QHNC) theory derived from exact expressions for radial distribution functions (RDF's) of the electron-ion model for liquid metals. In our method based on the QHNC equations, the classical MD is performed repeatedly to determine a self-consistent effective interionic potential, which depends on the ion-ion RDF of the system. This resultant effective ionic potential is obtained to be consistent with the density distribution of a pseudoatom and the electron-ion RDF, as well as the ion-ion RDF and the ion-ion bridge function, which are determined exactly as a result of the repeated MD simulation. We have applied this QHNC-MD method for Li, Na, K, Rb, and Cs near the melting temperature using upto 16,000 particles for the MD simulation. It is found that the convergence of the effective interionic potential is fast enough for practical applications; typically two MD runs are enough for convergence of the effective ionic potential within accuracy of 3 to 4 digits. Furthermore the resultant static structure factor is in excellent agreement with experimental data of X-ray and/or neutron scatering.     

Discovery of potential sclerostin inhibitors from plants with loop2 region of sclerostin inhibition by interacting with residues outside Pro-Asn-Ala-Ile-Gly motif

Abstract

Sclerostin, an antagonist of the Wnt/β-catenin signaling pathway, was discovered as a potential therapeutic target for stimulating bone formation in osteoporosis. In this study, molecular docking was employed to predict the binding of 29 herbal compounds, which were reported as bone formation stimulators, to the loop2 region of sclerostin. Then, the 50 ns molecular dynamics (MD) simulation of the complexes between sclerostin and the top 10 hits obtained from molecular docking were carried out. Root mean square deviations (RMSDs) analysis of MD trajectories pointed out that all ligands-complexes remain stable throughout the duration of MD simulations. In addition, the molecular mechanics/generalized born surface area (MM/GBSA) binding free energy and energy decomposition analyses were determined. The results here suggested that baicalin is the most promising inhibitor of sclerostin. Interestingly, baicalin binds to sclerostin via the hydrophobic interaction with the amino acid residues on loop2 region but outside the Pro-Asn-Ala-Ile-Gly (PNAIG) motif, particularly the Arg-Gly-Lys-Trp-Trp-Arg (RGKWWR) motif. This finding could be a novel strategy for developing new sclerostin inhibitors in the future.

Communicated by Ramaswamy H. Sarma     

Gibbs-Ensemble Molecular Dynamics: Liquid-Gas Equilibria for Lennard-Jones Spheres and n-Hexane

Abstract

We present a novel method to simulate phase equilibria in atomic and molecular systems. The method is a Molecular Dynamics version of the Gibbs-Ensemble Monte Carlo technique, which has been developed some years ago for the direct simulation of phase equilibria in fluid systems. The idea is to have two separate simulation boxes, which can exchange particles (or molecules) in a thermodynamically consistent fashion. Here we pres the derivation of the generalized equations of motion and discuss the relation of the resulting trajectory averages to the relevant ensemble. We test this Gibbs-Ensemble Molecular Dynamics algorithm by applying it to an atomic and a molecular system, i.e. to the liquid-gas coexistence in a Lennard-Jones fluid and in n-hexane. In both cases our results are in good accord with previous mean field and Gibbs-Ensemble Monte Carlo results as well as with the experimental data in the case of hexane. We also show that our Gibbs-Ensemble Molecular Dynamics algorithm like other Molecular Dynamics techniques can be used to study the dynamics of the system. Self-diffusion coefficients calculated with this method are in agreement with the result of conventional constant temperature Molecular Dynamics.     

Low-temperature molecular dynamics simulations of horse heart cytochrome c and comparison with inelastic neutron scattering data

Molecular dynamics (MD) simulation combined with inelastic neutron scattering can provide information about the thermal dynamics of proteins, especially the low-frequency vibrational modes responsible for large movement of some parts of protein molecules. We performed several 30-ns MD simulations of cytochrome c (Cyt c) in a water box for temperatures ranging from 110 to 300 K and compared the results with those from experimental inelastic neutron scattering. The low-frequency vibrational modes were obtained via dynamic structure factors, S(Q, ω), obtained both from inelastic neutron scattering experiments and calculated from MD simulations for Cyt c in the same range of temperatures. The well known thermal transition in structural movements of Cyt c is clearly seen in MD simulations; it is, however, confined to unstructured fragments of loops Ω₁ and Ω₂; movement of structured loop Ω₃ and both helical ends of the protein is resistant to thermal disturbance. Calculated and experimental S(Q, ω) plots are in qualitative agreement for low temperatures whereas above 200 K a boson peak vanishes from the calculated plots. This may be a result of loss of crystal structure by the protein–water system compared with the protein crystal.     

A force field for naproxen

A united-atom potential model for naproxen suitable for molecular dynamics (MD) simulation has been developed. The charge distribution is approximated by point charges obtained from ab initio calculations using the CHELPG method. Also the intramolecular interactions such as bond and angle vibration, and the torsion potential are obtained from ab initio calculations. The dispersive interaction contribution is taken from the literature. By MD simulation using a naproxen film in slap geometry, the temperature dependence of the density, surface tension and self-diffusion coefficient as well as the melting temperature for the developed potential model are obtained.     

Modeling the interaction between anti-cancer drug penicillamine and pristine and functionalized carbon nanotubes for medical applications: density functional theory investigation and a molecular dynamics simulation

Abstract

The present study focuses on the prediction and investigation of binding properties of penicillamine with pure (5,5) single-walled carbon nanotube (SWCNT) and functionalized SWCNT (f-SWCNT) through the B3LYP and M06-2X functionals using the 6-31G** basis set. The electronic and structural properties, adsorption energy and frontier molecular orbitals of various configurations are examined. Our theoretical results indicated that the interaction of the nanotubes with penicillamine molecule is weak so that the drug adsorption process is typically physisorption. Also, results of theoretical calculations show that the adsorption of the drug molecule on f-SWCNT is stronger with shorter intermolecular distances in comparison to pure SWCNT. The natural bond orbital (NBO) analysis of studied systems demonstrates that the charge is transferred from penicillamine molecule to the nanotubes. Furthermore, molecular dynamics (MD) simulation is employed to evaluate the dynamic and diffusion behavior of drug molecule on SWCNT and f-SWCNT. Energy results show that drug molecule spontaneously moves toward the carriers, and the van der Waals energy contributions in drug adsorption are more than electrostatic interaction. The obtained results from MD simulation confirm that the functionalization of SWCNT leads to increase in the solubility of the carrier in aqueous solution.

Communicated by Ramaswamy H. Sarma     

Accuracy and precision of NMR relaxation experiments and MD simulations for characterizing protein dynamics

Model-free parameters obtained from nuclear magnetic resonance (NMR) relaxation experiments and molecular dynamics (MD) simulations commonly are used to describe the intramolecular dynamical properties of proteins. To assess the relative accuracy and precision of experimental and simulated model-free parameters, three independent data sets derived from backbone ¹⁵N NMR relaxation experiments and two independent data sets derived from MD simulations of Escherichia coli ribonuclease HI are compared. The widths of the distributions of the differences between the order parameters for pairs of NMR data sets are congruent with the uncertainties derived from statistical analyses of individual data sets; thus, current protocols for analyzing NMR data encapsulate random uncertainties appropriately. Large differences in order parameters for certain residues are attributed to systematic differences between samples for intralaboratory comparisons and unknown, possibly magnetic field-dependent, experimental effects for interlaboratory comparisons. The widths of distributions of the differences between the order parameters for two NMR sets are similar to widths of distributions for an NMR and an MD set or for two MD sets. The linear correlations between the order parameters for an MD set and an NMR set are within the range of correlations observed between pairs of NMR sets. These comparisons suggest that the NMR and MD generalized order parameters for the backbone amide N—H bond vectors are of comparable accuracy for residues exhibiting motions on a fast time scale (<100 ps). Large discrepancies between NMR and MD order parameters for certain residues are attributed to the occurrence of “rare” motional events over the simulation trajectories, the disruption of an element of secondary structure in one of the simulations, and lack of consensus among the experimental data sets. Consequently, (easily detectable) severe distortions of local protein structure and infrequent motional events in MD simulations appear to be the most serious artifacts affecting the accuracy and precision, respectively, of MD order parameters relative to NMR values. In addition, MD order parameters for motions on a fast (<100 ps) timescale are more precisely determined than their NMR counterparts, thereby permitting more detailed dynamic characterization of biologically important residues by MD simulation than is sometimes possible by experimental methods. Proteins 28:481–493, 1997. © 1997 Wiley-Liss, Inc.     

An assessment of optimal time scale of conformational resampling for parallel cascade selection molecular dynamics

Parallel cascade selection molecular dynamics (PaCS-MD) has been proposed as a conformational sampling method for enhancing structural transitions from a given reactant to a product by repeating cycles of short-time MD simulations. In the present paper, we assessed how the time scale of a short-time MD simulation affected the computational efficiency by changing the simulation length. In conclusion, ps-order (t_ps) PaCS-MD simulations showed a higher computational efficiency as a total simulation time over the cycles than ns-order (t_ns) PaCS-MD simulations, indicating that t_ps might be suitable for generating structural transitions efficiently.     

In silico studies on p21-activated kinase 4 inhibitors: comprehensive application of 3D-QSAR analysis,molecular docking,molecular dynamics simulations,and MM-GBSA calculation

Abstract

P21-activated kinase 4 (PAK4) is a serine/threonine protein kinase, which is associated with many cancer diseases, and thus being considered as a potential drug target. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations were performed to explore the structure-activity relationship of a series of pyrropyrazole PAK4 inhibitors. The statistical parameters of comparative molecular field analysis (CoMFA, Q ² = 0.837, R ² = 0.990, and R ² _pred = 0.967) and comparative molecular similarity indices analysis (CoMSIA, Q ² = 0.720, R ² = 0.972, and R ² _pred = 0.946) were obtained from 3D-QSAR model, which exhibited good predictive ability and significant statistical reliability. The binding mode of PAK4 with its inhibitors was obtained through molecular docking study, which indicated that the residues of GLU396, LEU398, LYS350, and ASP458 were important for activity. Molecular mechanics generalized born surface area (MM-GBSA) method was performed to calculate the binding free energy, which indicated that the coulomb, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. Furthermore, through 100?ns MD simulations, we obtained the key amino acid residues and the types of interactions they participated in. Based on the constructed 3D-QSAR model, some novel pyrropyrazole derivatives targeting PAK4 were designed with improved predicted activities. Pharmacokinetic and toxicity predictions of the designed PAK4 inhibitors were obtained by the pkCSM, indicating these compounds had better absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. Above research provided a valuable insight for developing novel and effective pyrropyrazole compounds targeting PAK4.     

Cross-Scale Numerical Simulations Using Discrete Particle Models

Abstract

We propose a concept for a homogenous computational model in carrying out cross-scale numerical experiments on liquids. The model employs the particle paradigm and comprises three types of simulation techniques: molecular dynamics (MD), dissipative particle dynamics (DPD) and smoothed particle hydrodynamics (SPH). With respect to the definition of the collision operator, this model may work in different hierarchical spatial and time scales as: MD in the atomistic scale, DPD in the mesoscale and SPH in the macroscale. The optimal computational efficiency of the three types of cross-scale experiments are estimated in dependence on: the system size N-where N is the number of particles-and the number of processors P employed for computer simulation. For the three-hierarchical-stage, as embodied in the MD-DPD-SPH model, the efficiency is proportional to N ^8/7 but its dependence on P is different for each of the three types of cross-scale experiments. The problem of matching the different scales is discussed.     

Electronic Circular Dichroism of [16]Helicene With Simplified TD‐DFT: Beyond the Single Structure Approach

The electronic circular dichroism (ECD) spectrum of the recently synthesized [16]helicene and a derivative comprising two triisopropylsilyloxy protection groups was computed by means of the very efficient simplified time‐dependent density functional theory (sTD‐DFT) approach. Different from many previous ECD studies of helicenes, nonequilibrium structure effects were accounted for by computing ECD spectra on "snapshots" obtained from a molecular dynamics (MD) simulation including solvent molecules. The trajectories are based on a molecule specific classical potential as obtained from the recently developed quantum chemically derived force field (QMDFF) scheme. The reduced computational cost in the MD simulation due to the use of the QMDFF (compared to ab‐initio MD) as well as the sTD‐DFT approach make realistic spectral simulations feasible for these compounds that comprise more than 100 atoms. While the ECD spectra of [16]helicene and its derivative computed vertically on the respective gas phase, equilibrium geometries show noticeable differences, these are “washed” out when nonequilibrium structures are taken into account. The computed spectra with two recommended density functionals (ωB97X and BHLYP) and extended basis sets compare very well with the experimental one. In addition we provide an estimate for the missing absolute intensities of the latter. The approach presented here could also be used in future studies to capture nonequilibrium effects, but also to systematically average ECD spectra over different conformations in more flexible molecules. Chirality 28:365–369, 2016. © 2016 Wiley Periodicals, Inc.     

Homology modeling and examination of the effect of the D92E mutation on the H5N1 nonstructural protein NS1 effector domain

Virulent H5N1 strains of influenza virus often harbor a D92E point mutation in the nonstructural protein NS1. This crucial mutation has been correlated with increased virulence and/or cytokine resistance, but the structural implications of such a change are still unclear. Furthermore, NS1 protein could also be a potential target for the development of novel antiviral agents against H5N1 strains. Therefore, a reasonable 3D model of H5N1 NS1 is important for the understanding of the molecular basis of increased virulence and the design of novel antiviral agents. Based on the crystal structure of a non-H5N1 NS1 protein, a model of H5N1 NS1 was developed by homology modeling, molecular mechanics and molecular dynamics simulations. It was found that the D92E mutation could result in weakened interactions of the carboxylate side chain with other phosphorylated residues, thereby activating phosphorylation of NS1. Figure Superposition of snapshots picked from the two molecular dynamic (MD) trajectories: a H5N1 NS1 homology model and b non-H5N1 NS1 crystal structure after 0 (green ribbon), 5 (blue ribbon) and 10 ns (pink ribbon) MD simulation     

Molecular-Dynamics Simulations for the Density Autocorrelation Function in a Supercooled Fluid Phase

Abstract

We have re-calculated the self part of the density autocorrelation function F_s(k, t) (incoherent scattering function) for the binary soft-sphere fluid with a much longer molecular-dynamics (MD) simulation than our previous MD calculations, and with a larger system size (N = 4000) to a longer time window as well as to study a system-size dependence, if it exists. The full density autocorrelation function F(k, t) was also computed. It is found that all F(k, t)'s that we have computed in this work can be fitted over a wide range of time steps (at least over three figures of the decay) by a Williams-Watts stretched exponential function F_s(k, t) = A exp [— (t/t ₀)^β], where A, β and t ₀ are adjustable parameters. Other significant dynamical behaviours were also presented in mean square displacements and non-Gaussian parameters for highly supercooled fluids with N = 4000. The present results are compatible to our previous computations with N = 500, but a significant size dependence is suggested.     

Molecular Dynamics Study of BaB2O4 in Crystalline and Molten States

Abstract

Structural aspects of BaB₂O₄ liquids have been investigated by the molecular dynamics simulation including the determination on the parameters of the interatomic potential applicable to BaB₂O₄ in both crystalline and molten states. The structure and physical properties of BaB₂O₄ crystals were successfully reproduced by the MD simulation for both α and β phases. The simulated interference function in the liquid state was also in good agreement with the experimental one. Several interesting features on the relaxation phenomena just after melting were reproduced by the simulation that the structure factors of simulated liquid maintain the characteristic features of the original crystal structure for more than 40ps after melting, and the variation of the number of rings formed by B-O bondings was found to increase after melting.     

On the ability of molecular dynamics force fields to recapitulate NMR derived protein side chain order parameters

Molecular dynamics (MD) simulations have become a central tool for investigating various biophysical questions with atomistic detail. While many different proxies are used to qualify MD force fields, most are based on largely structural parameters such as the root mean square deviation from experimental coordinates or nuclear magnetic resonance (NMR) chemical shifts and residual dipolar couplings. NMR derived Lipari–Szabo squared generalized order parameter (O²) values of amide N? H bond vectors of the polypeptide chain were also often employed for refinement and validation. However, with a few exceptions, side chain methyl symmetry axis order parameters have not been incorporated into experimental reference sets. Using a test set of five diverse proteins, the performance of several force fields implemented in the NAMDD simulation package was examined. It was found that simulations employing explicit water implemented using the TIP3 model generally performed significantly better than those using implicit water in reproducing experimental methyl symmetry axis O² values. Overall the CHARMM27 force field performs nominally better than two implementations of the Amber force field. It appeared that recent quantum mechanics modifications to side chain torsional angles of leucine and isoleucine in the Amber force field have significantly hindered proper motional modeling for these residues. There remained significant room for improvement as even the best correlations of experimental and simulated methyl group Lipari–Szabo generalized order parameters fall below an R² of 0.8.