Architectures and algorithms for millisecond-scale molecular dynamics simulations of proteins

MotivationThe ligand/protein binding interaction is typically investigated by docking and molecular dynamics (MD) simulations. In particular, docking-based virtual screening (VS) is used to select the best ligands from database of thousands of compounds, while MD calculations assess the energy stability of the ligand/protein binding complexes. Considering the broad use of these techniques, it is of great demand to have one single software that allows a combined and fast analysis of VS and MD results. With this in mind, we have developed the Drug Discovery Tool (DDT) that is an intuitive graphics user interface able to provide structural data and physico-chemical information on the ligand/protein interaction.ResultsDDT is designed as a plugin for the Visual Molecular Dynamics (VMD) software and is able to manage a large number of ligand/protein complexes obtained from AutoDock4 (AD4) docking calculations and MD simulations. DDT delivers four main outcomes: i) ligands ranking based on an energy score; ii) ligand ranking based on a ligands’ conformation cluster analysis; iii) identification of the aminoacids forming the most occurrent interactions with the ligands; iv) plot of the ligands’ center-of-mass coordinates in the Cartesian space. The flexibility of the software allows saving the best ligand/protein complexes using a number of user-defined options.Availability and implementationDDT_site_1 (alternative DDT_site_2); the DDT tutorial movie is available here.Supplementary informationSupplementary data are available at Bioinformatics online.

Inhibition of dipeptidyl peptidase-4 (DPP4) activity has emerged as a promising therapeutic approach for the treatment of type 2 diabetes mellitus (T2DM). Bioinformatics-driven approaches have emerged as crucial tools in drug discovery. Molecular docking and molecular dynamics (MD) simulations are effective tools in drug discovery, as they reduce the time and cost associated with experimental screening. In this study, we employed structure-assisted in-silico methods, including molecular docking and MD simulations, to identify SRT2183, a small molecule that may potentially inhibit the activity of DPP4 enzyme. The interaction between the small molecule "SRT2183" and DPP4 exhibited a binding affinity of −9.9 Kcal/Mol, leading to the formation of hydrogen bonds with the amino acid residues MET348, SER376, and THR351 of DPP4. The MD simulations over a period of 100 ns indicated stable protein-ligand interactions, with no significant conformational rearrangements observed within the simulated timeframe. In conclusion, our results suggest that the small molecule SRT2183 may have the potential to inhibit the DPP4 enzyme and pave the way for the therapeutics of T2DM.

The term molecular mechanics (MM) refers to the use of simple potential-energy functions (e.g., harmonic oscillator or Coulombic potentials) to model molecular systems. Molecular mechanics approaches are widely applied in molecular structure refinement, molecular dynamics (MD) simulations, Monte Carlo (MC) simulations, and ligand-docking simulations. Typically, molecular mechanics models consist of spherical atoms connected by springs which represent bonds. Internal forces experienced in the model structure are described using simple mathematical functions. For example, Hooke’s law is commonly used to describe bonded interactions, and the nonbonded atoms might be treated as inelastic hard spheres or may interact according to a Lennard-Jones potential. Using these simple models, a molecular dynamics simulation numerically solves Newton’s equations of motion, thus allowing structural fluctuations to be observed with respect to time. Dynamic simulation methods are widely used to obtain information on the time evolution of conformations of proteins and other biological macromolecules1–4 and also kinetic and thermodynamic information. Simulations can provide fine detail concerning the motions of individual particles as a function of time. They can be utilized to quantify the properties of a system at a precision and on a time scale that is otherwise inaccessible, and simulation is, therefore, a valuable tool in extending our understanding of model systems. Theoretical consideration of a system additionally allows one to investigate the specific contributions to a property through “computational alchemy”,5 that is, modifying the simulation in a way that is nonphysical but nonetheless allows a model’s characteristics to be probed. One particular example is the artificial conversion of the energy function from that representing one system to that of another during a simulation. This is an important technique in free-energy calculations.6 Thus, molecular dynamics simulations, along with a range of complementary computational approaches, have become valuable tools for investigating the basis of protein structure and function. This review offers an outline of the origin of molecular dynamics simulation for protein systems and how it has developed into a robust and trusted tool. This review then covers more recent advances in theory and an illustrative selection of practical studies in which it played a central role. The range of studies in which MD has played a considerable or pivotal role is immense, and this review cannot do justice to them; MD simulations of biomedical importance were recently reviewed.4 Particular emphasis will be placed on the study of dynamic aspects of protein recognition, an area where molecular dynamics has scope to provide broad and far-ranging insights. This review concludes with a brief discussion of the future potential offered to advancement of the biological and biochemical sciences and the remaining issues that must be overcome to allow the full extent of this potential to be realized. 1.1. Historical Background MD methods were originally conceived within the theoretical physics community during the 1950s. In 1957, Alder and Wainwright7 performed the earliest MD simulation using the so-called hard-sphere model, in which the atoms interacted only through perfect collisions. Rahman8 subsequently applied a smooth, continuous potential to mimic real atomic interactions. During the 1970s, as computers became more widespread, MD simulations were developed for more complex systems, culminating in 1976 with the first simulation of a protein9,10 using an empirical energy function constructed using physics-based first-principles assumptions. MD simulations are now widely and routinely applied and especially popular in the fields of materials science11,12 and biophysics. As will be discussed later in this review, a variety of experimental conditions may be simulated with modern theories and algorithms. The initial simulations only considered single molecules in vacuo. Over time, more realistic or at least biologically relevant simulations could be performed. This trend is continuing today. The initial protein MD simulation, of the small bovine pancreatic trypsin inhibitor (BPTI), covered only 9.2 ps of simulation time. Modern simulations routinely have so-called equilibration periods much longer than that, and production simulations of tens of nanoseconds are routine, with the first microsecond MD simulation being reported in 1998.13 In addition, the original BPTI simulation included only about 500 atoms rather than the 104-106 atoms that are common today. While much of this advancement results from an immense increase in availability of computing power, major theoretical and methodological developments also contribute significantly. The number of publications regarding MD theory and application of MD to biological systems is growing at an extraordinary pace. A single review cannot do justice to the recent applications of MD. Using data from ISI Web of Science, the authors estimate that during 2005 at least 800 articles will be published that discuss molecular dynamics and proteins. The historical counts are shown in Figure 1. Open in a separate window Figure 1 Articles matching ISI Web of Science query “TS=(protein) AND TS=(molecular dynamics)”.

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

We present a continuum approach for efficient and accurate calculation of reaction field forces and energies in classical molecular-dynamics (MD) simulations of proteins in water. The derivation proceeds in two steps. First, we reformulate the electrostatics of an arbitrarily shaped molecular system, which contains partially charged atoms and is embedded in a dielectric continuum representing the water. A so-called fuzzy partition is used to exactly decompose the system into partial atomic volumes. The reaction field is expressed by means of dipole densities localized at the atoms. Since these densities cannot be calculated analytically for general systems, we introduce and carefully analyze a set of approximations in a second step. These approximations allow us to represent the dipole densities by simple dipoles localized at the atoms. We derive a system of linear equations for these dipoles, which can be solved numerically by iteration. After determining the two free parameters of our approximate method we check its quality by comparisons (i) with an analytical solution, which is available for a perfectly spherical system, (ii) with forces obtained from a MD simulation of a soluble protein in water, and (iii) with reaction field energies of small molecules calculated by a finite difference method.

Molecular dynamics simulations of intrinsically disordered protein regions enable biophysical interpretation of variant effect predictors.

Molecular Dynamics (MD) simulation of biomolecules provides important insights into conformational changes and dynamic behavior, revealing critical information about folding and interactions with other molecules. The collection of simulations stored in computers across the world holds immense potential to serve as training data for future Machine Learning models that will transform the prediction of structure, dynamics, drug interactions, and more. Ideally, there should exist an open access repository that enables scientists to submit and store their MD simulations of proteins and protein-drug interactions, and to find, retrieve, analyze, and visualize simulations produced by others. However, despite the ubiquity of MD simulation in structural biology, no such repository exists; as a result, simulations are instead stored in scattered locations without uniform metadata or access protocols. Here, we introduce MDRepo, a robust infrastructure that provides a relatively simple process for standardized community contribution of simulations, activates common downstream analyses on stored data, and enables search, retrieval, and visualization of contributed data. MDRepo is built on top of the open-source CyVerse research cyber-infrastructure, and is capable of storing petabytes of simulations, while providing high bandwidth upload and download capabilities and laying a foundation for cloud-based access to its stored data.

Computer simulation of materials using parallel architectures

Selectivity and transport in aquaporins from molecular simulation studies.

Introduction: Milk kefir, a fermenting milk made with kefir grains, has shown potential in promoting apoptosis, regulating the cell cycle, and reducing tumor growth in breast cancer cells. This study aimed to investigate the stability and potential of milk kefir metabolites as inhibitors of breast cancer growth by interacting with the estrogen receptor alpha (ER-α), a key protein involved in breast cancer cell proliferation. We used computational methods, specifically molecular docking simulations with AutoDock and molecular dynamics (MD) simulations with Gromacs, to analyze how these metabolites bind to ER-α. Methods: A combination of molecular docking and MD simulations was used to explore how metabolites derived from milk kefir interact with ER-α, a crucial target in breast cancer therapy. The methodology included multiple stages: preparation of target proteins, preparation and screening of the metabolites, geometry optimization, molecular docking, and MD simulations. Results: The molecular docking simulations of 43 metabolites revealed three promising candidates: 2-Methyl (S35), benzeneethanol (S42), and 2,6-dimethyl-4-heptanone (S54), with binding affinities (ΔG) of -5.08, -5.06, and -4.90 kcal/mol, respectively. MD simulations further showed that the selected metabolites stabilized the ER-α-metabolite complex, with the 2,6-dimethyl-4-heptanone (S54) metabolite demonstrating the most negative total MM-GBSA energy value (ΔG = -22.98 kcal/mol), indicating a strong and stable binding interaction. Conclusion: 2,6-Dimethyl-4-heptanone, a metabolite from milk kefir, showed promising potential as a candidate for further development as a breast cancer treatment, offering a novel alternative to conventional therapies.

Recent advances have made all-atom molecular dynamics (MD) a powerful tool to sample the conformational energy landscape. There are still however three major challenges in the application of MD to biological systems: accuracy of force field, time scale, and the analysis of simulation trajectories. Significant progress in addressing the first two challenges has been made and extensively reviewed previously. This Account focuses on strategies of analyzing simulation data of biomolecules that also covers ways to properly design simulations and validate simulation results. In particular, we examine an approach named comparative perturbed-ensembles analysis, which we developed to efficiently detect dynamics in protein MD simulations that can be linked to biological functions. In our recent studies, we implemented this approach to understand allosteric regulations in several disease-associated human proteins. The central task of a comparative perturbed-ensembles analysis is to compare two or more conformational ensembles of a system generated by MD simulations under distinct perturbation conditions. Perturbations can be different sequence variations, ligand-binding conditions, and other physical/chemical modifications of the system. Each simulation is long enough (e.g., microsecond-long) to ensure sufficient sampling of the local substate. Then, sophisticated bioinformatic and statistical tools are applied to extract function-related information from the simulation data, including principal component analysis, residue-residue contact analysis, difference contact network analysis (dCNA) based on the graph theory, and statistical analysis of side-chain conformations. Computational findings are further validated with experimental data. By comparing distinct conformational ensembles, functional micro- to millisecond dynamics can be inferred. In contrast, such a time scale is difficult to reach in a single simulation; even when reached for a single condition of a system, it is elusive as to what dynamical motions are related to functions without, for example, comparing free and substrate-bound proteins at the minimum. We illustrate our approach with three examples. First, we discuss using the approach to identify allosteric pathways in cyclophilin A (CypA), a member of a ubiquitous class of peptidyl-prolyl cis-trans isomerase enzymes. By comparing side-chain torsion-angle distributions of CypA in wild-type and mutant forms, we identified three pathways: two are consistent with recent nuclear magnetic resonance experiments, whereas the third is a novel pathway. Second, we show how the approach enables a dynamical-evolution analysis of the human cyclophilin family. In the analysis, both conserved and divergent conformational dynamics across three cyclophilin isoforms (CypA, CypD, and CypE) were summarized. The conserved dynamics led to the discovery of allosteric networks resembling those found in CypA. A residue wise determinant underlying the unique dynamics in CypD was also detected and validated with additional mutational MD simulations. In the third example, we applied the approach to elucidate a peptide sequence-dependent allosteric mechanism in human Pin 1, a phosphorylation-dependent peptidyl-prolyl isomerase. We finally present our outlook of future directions. Especially, we envisage how the approach could help open a new avenue in drug discovery.

Molecular simulations restrained to single or multiple templates are commonly used in protein-structure modeling. However, the restraints introduce additional barriers, thus impairing the ergodicity of simulations, which can affect the quality of the resulting models. In this work, the effect of restraint types and simulation schemes on ergodicity and model quality was investigated by performing template-restrained canonical molecular dynamics (MD), multiplexed replica-exchange molecular dynamics, and Hamiltonian replica exchange molecular dynamics (HREMD) simulations with the coarse-grained UNRES force field on nine selected proteins, with pseudo-harmonic log-Gaussian (unbounded) or Lorentzian (bounded) restraint functions. The best ergodicity was exhibited by HREMD. It has been found that non-ergodicity does not affect model quality if good templates are used to generate restraints. However, when poor-quality restraints not covering the entire protein are used, the improved ergodicity of HREMD can lead to significantly improved protein models. © 2017 Wiley Periodicals, Inc.

We combine molecular dynamics (MD) and Monte Carlo (MC) simulations to study the properties of gas molecules confined between two hard walls of a microchannel or nanochannel. The coupling between MD and MC simulations is introduced by performing MD near the boundaries for accuracy and MC in the bulk because of the low computational cost. We characterize the influence of different densities and molecule sizes on the equilibrium properties of the gas in the microchannel. The effect of the particle size on the simulation results is very small in the case of a dilute gas and increases with the density. The hybrid MD-MC simulation method is validated by comparing the results for density and temperature profiles with those of pure MD and pure MC simulations. These results compare well for pure MD and pure MC, as well as hybrid MD-MC, both in the bulk and near the boundaries, when hard-sphere interactions are used. When Lennard-Jones potentials are used to accurately model the interactions between the gas and wall molecules instead, the results of pure MD simulations differ significantly from the pure MC simulations near the boundaries, but the results of the hybrid method compare well with the pure MD results near the wall, and with the pure MC and pure MD results in the middle of the channel. The hybrid method also very accurately simulates the interface between the MD and MC simulation domains. Comparisons between MD, MC, and hybrid MD-MC computational costs are outlined. The speedup when using 50% of the domain for MD simulations and 50% for MC simulations is very small compared to pure MD simulations times, but this speedup increases drastically for more realistic situations where the region near the wall is small compared to the bulk region.

HighlightsAccording to the design of the experiment (DoE), multivariate analysis models were used to optimize the critical process parameters combined with multi-objective optimization.Based on the optimized operating conditions, the MILT-HD method not only enhances the extraction efficiency from Amomi fructus but also reduces energy demands and CO2 emissions.Based on the density functional theoretical (DFT) and molecular dynamics (MD) simulations, the mechanisms for ionic liquids (ILs) to improve the extraction efficiency of essential oil was comprehensively revealed.In this paper, Amomi fructus (Latin) was used to explore the mechanism of ionic liquids (ILs) in improving the extraction efficiency of essential oils. Microwave assisted ionic liquid treatment followed by a hydro-distillation (MILT-HD) process for isolating Amomi fructus essential oil was optimized by multi-objective optimization. Under optimum operating conditions, the IL-assisted extraction method not only enhances extraction efficiency but also reduces energy demands and CO2 emissions. Since the hydrogen bond structure network of cellulose in the cell wall is an important reason for hindering diffusion of essential oils, the mechanism of ILs was explored by density functional theoretical (DFT) and molecular dynamics (MD) simulations. According to DFT calculations, ILs can facilitate the cleavage of cellulose chains and have strong non-covalent interactions with cellulose. Based on the MD simulations, the degree of destruction of the cellulose hydrogen bond structure was explored. According to the DFT and MD simulations, the ILs can significantly destroy cellulose structure, thereby promoting essential oil release from the plant. These results were confirmed by scanning electron microscopy (SEM) and Fourier transform infrared spectroscopy (FTIR). This work is conducive to better understand the MILT-HD process for isolating essential oil and comprehensively understand the mechanism of ILs in the extraction process.

Prediction of nonlocal elasticity parameters using high-throughput molecular dynamics simulations and machine learning