The Role of Molecular Dynamics Simulations for the Study of Slow Dynamics

Abstract

A two step strategy is proposed to study dynamical properties of a physical system much slower than the time scales accessible by molecular dynamics simulations. The strategy is applied to investigate the slow dynamics of supercooled liquids.     

A Combined 2D-NMR and Molecular Dynamics Analysis of the Structure of the Actinomycin D: d(ATGCAT)2 Complex

Abstract

We present a comparative analysis of an NMR experiment and molecular and harmonic dynamics simulations of an actinomycin D: d(ATGCAT)₂ complex. A comparison of NOE measurements and 1/R⁶ weighted proton-proton distances confirm the general correctness of the Actinomycin D-DNA model proposed by Sobell. There are, however, some substantial differences between the proton-proton distances inferred from the NOE results and the molecular and harmonic dynamics simulations. The remaining discrepancies could either come from contributions of other conformations to the average properties of the complex or from uncertainties in the NMR distance analysis. An analysis of the molecular dynamics helix properties, sugar puckers, hydrogen bonding, rms fluctuations and torsional properties are qualitatively consistent with those from previous simulations, but the presence of an intercalated drug leads to some new structural and dynamical features.     

Molecular Dynamics as a Mathematical Mapping. II. Partial Derivatives in the Microcanonical Ensemble

Abstract

Performing molecular dynamics in a fully continuous and differentiable framework can be viewed as a deterministic mathematical mapping between, on one side, the force field parameters that describe the potential energy interactions and input macroscopic conditions, and, on the other, the calculated corresponding macroscopic properties of the bulk molecular system.

Within this framework, it is possible to apply standard methods of variational calculus for the computation of the partial derivatives of the molecular dynamics mapping based on the integration of either the adjoint equations or the sensitivity equations of the classical Newtonian equations of motion. We present procedures for these computations in the standard microcanonical (N, V, E) ensemble, and compare the computational efficiency of the two approaches. The general formulations developed are applied to the specific example of bulk ethane fluid.

With these procedures in place, it is now possible to compute the partial derivatives of any property determined by molecular dynamics with respect to any input property and any potential parameter. Moreover, these derivatives are computed to essentially the same level of numerical accuracy as the output properties themselves.     

Structural and energetic insight into the isoform-selective inhibitors of tumour marker Hsp90 against Grp94

The heat shock protein 90 (Hsp90) represents a new avenue for cancer therapy. A novel benzolactam inhibitor, compound 31, showed a great selectivity for Hsp90α and Hsp90β against Grp94. However, the structural basis for the great selectivity of compound 31 for Hsp90α/β versus Grp94 remains poorly understood. In this study, we carried out molecular docking, molecular dynamics simulations and binding free energy calculations (MM-GBSA) to address the isoform selective property. Molecular docking studies indicated the different binding modes of the Hsp90 and Grp94 with compound 31. The MM-GBSA calculations revealed that the hydrophobic interactions between compound 31 and proteins contributed the most to the binding affinity and the Grp94/compound 31 complex could result in a less energy-favourable complex compared with the Hsp90α/compound 31 and the Hsp90β/compound 31 complexes. This may render the weak binding of compound 31 to the Grp94. This study may be helpful for the future design of novel and selective Hsp90 inhibitors.     

Two- and Higher Point Correlation Functions in Proteins

Abstract

A method is presented for a more efficient sampling of the configurational space of proteins as compared to conventional sampling techniques such as molecular dynamics. The method is based on the large conformational changes in proteins revealed by the “essential dynamics” analysis. A form of constrained dynamics is performed, forcing the system to move along some of the essential coordinates. This results in a broader sampling of the essential subspace than in a comparable conventional molecular dynamics simulation without constraints. The new sampling method (essential dynamics sampling) was applied to the histidine-containing phosphocarrier protein HPr. The results indicate that the essential dynamics sampling method produces physically allowed structures, as estimated by the evaluation of many geometrical properties. In addition, a study of the motions in the essential subspace reveals a diffusion-like behavior.     

On the Approximation of Solvent Effects on the Conformation and Dynamics of Cyclosporin A by Stochastic Dynamics Simulation Techniques

Abstract

The molecular simulation technique of stochastic dynamics (SD) is tested by application to the immunosuppressive drug cyclosporin A (CPA). Two stochastic dynamics simulations are performed, one (SD_CCl4 ) with atomic friction coefficients proportional to the viscosity of the nonpolar solvent CCl₄, and one (SD_H2O) with atomic friction coefficients corresponding to an aqueous solution. The atomic friction coefficients are also taken proportional to an approximate expression for the atomic accessible surface area. The properties of both stochastic dynamics simulations are compared to those of two full molecular dynamics (MD) simulations of cyclosporin A, one in a box with 591 CCl₄ molecules, and one in a box with 632 H₂O molecules.

The properties of cyclosporin A as found in the molecular dynamics simulation in CCl₄ are well reproduced by the SD_CCl4 simulation. This indicates that the neglect of a mean force reresenting the average solvent effects on the solute is justified in the case of nonpolar solvents. For polar solvents, like water, this mean force may not be neglected. The SD_H2O simulation of cyclosporin A clearly fails to reproduce the amount of hydrogen bonding found in the molecular dynamics stimulation of cyclosporin A in water.

A comparison with a molecular dynamics simulation of cyclosporin A in vacuo shows that both the SD_CCl4 and the SD_H2O simulation come closer to the properties of the molecular dynamics simulations in CCl₄ and in H₂O than a molecular dynamics simulation in vacuo.     

Generation and study of a relatively large amorphous silica surface in the liquid phase

Through the classical and ab-initio molecular dynamics computer simulations, the aim of this study was to generate and analyse a surface-type structure in the liquid phase at a temperature of 3400?K. First, a crystalline structure was used (β-cristobalite) with 216 SiO₂ molecules (648 atoms) and it was melted. Next, the resulting structure was properly cut, with a total of 546 particles remaining, and this system was put into a geometric arrangement denominated ‘sandwich’, where it was equilibrated. The prior processes were implemented through the classical molecular dynamics computer simulations using the effective classical potential of Feuston [B. P. Feuston and S. H. Garofalini, J. Chem. Phys. 89, 5818 (1988)]. Finally, the electronic charge distribution and topological and bonding defects were analysed in the first five superior layers of this system. This latter was carried out through the ab-initio molecular dynamics computer simulations, where the pseudopotentials of N. Trouiller and J. L. Martins [N. Trouiller, J. L. Martins, Phys. Rev. B 43, 1993 (1991)] were used. The results show different aspects between them, the interaction between the bonding and topological defects, and the existence of two-membered rings, nonbridging oxygens (NBOs) and three-coordinated silicon at the outermost surface.     

Computer Simulation of Liquid Methanol II. System Size Effects

Abstract

A series of molecular dynamics simulations of liquid methanol has been carried out on a supernode transputer array. Four system sizes from 125 to 512 molecules have been considered, in order to study the effect of system size on the calculated structural, orientational and dynamic properties. The dielectric constant and the dielectric relaxation time are compared with experimental data.     

Reaction Field Effects on the Simulated Properties of Liquid Water

Abstract

The effects of including a reaction field contribution on the structure and dynamics of liquid water have been investigated using molecular dynamics simulations. Reaction field effects are determined for two models of water, the simple point charge (SPC) model and the extended simple point charge (SPC/E) model, and at two temperatures (277 K and 300 K). Inclusion of the reaction field leads to a reduced system density, an increase in translational diffusion, which is model dependent, an increase in internal energy, and an increase in rotational diffusion rates, in addition to the large (known) changes in the dielectric properties of liquid water. It is concluded that continued use of the reaction field technique should involve a reparameterization of the water model and not merely a merging with the original model parameters.     

基于处方挖掘与分子动力学模拟筛选严重急性呼吸综合征冠状病毒2潜在抑制剂分子的研究

目的 从中药筛选具有潜在抑制严重急性呼吸综合征冠状病毒2 （SARS-CoV-2） 活性的成分，进一步从原子水平揭 示其抑制SARS-CoV-2 表面刺突蛋白（S 蛋白） 受体结合域（RBD） 与血管紧张素转化酶2 （ACE2） 结合的内在机制。 方法 检索新型冠状病毒（简称“新冠肺炎”） 治疗中药处方，构建“新冠肺炎中药候选活性成分数据库”。用具有ACE2 抑制活性的小分子化合物构建HipHop药效团模型，并对“新冠肺炎中药候选活性成分数据库”中活性成分筛选。采用分子 对接和分子动力学模拟方法研究候选活性成分与ACE2 的结合方式及其对SARS-CoV-2 S 蛋白与ACE2 识别的影响。 结果 本文通过中药处方挖掘和分子动力学模拟，从143 个新冠肺炎治疗中药处方中筛选出10 种可与SARS-CoV-2 S 蛋白/ 人源ACE2 识别位点结合的中药成分。其中，枇杷叶主要活性成分23-trans-p-coumaryhormentic acid 与ACE2 具有最高的亲和 力，且23-trans-p-coumaryhormentic acid 的结合可有效阻断SARS-CoV-2 S蛋白与宿主细胞ACE2 的结合。结论 本文通过虚 拟筛选发现了SARS-CoV-2 潜在抑制剂分子23-trans-p-coumaryhormentic acid，同时从原子水平预测了其抑制SARS-CoV-2 S 蛋白与ACE2 结合的内在机制，这将为SARS-CoV-2 特异性抗病毒药物的研发提供理论依据。     

Some Recent Developments in Computational Chemistry

Abstract

Some recent developments in the use of computational methods to predict the properties of condensed phases are discussed; the use of Gibbs ensemble Monte Carlo to predict the phase equilibria of bulk phases, the use of molecular dynamics to elucidate Atomic Force Microscopy experiments on organic films, and the use of combined Monte Carlo/molecular dynamics techniques to enable the direct prediction of particle fluxes along pressure gradients in model microporous materials.     

Molecular aspects of the interaction between amphotericin B and a phospholipid bilayer: molecular dynamics studies

Amphotericin B (AmB) is a polyene macrolide antibiotic used to treat systemic fungal infections. The molecular mechanism of AmB action is still only partly characterized. AmB interacts with cell-membrane components and forms membrane channels that eventually lead to cell death. The interaction between AmB and the membrane surface can be regarded as the first (presumably crucial) step on the way to channel formation. In this study molecular dynamics simulations were performed for an AmB–lipid bilayer model in order to characterize the molecular aspects of AmB–membrane interactions. The system studied contained a box of 200 dimyristoylphosphatidylcholine (DMPC) molecules, a single AmB molecule placed on the surface of the lipid bilayer and 8,065 water molecules. Two molecular dynamics simulations (NVT ensemble), each lasting 1 ns, were performed for the model studied. Two different programs, CHARMM and NAMD2, were used in order to test simulation conditions. The analysis of MD trajectories brought interesting information concerning interactions between polar groups of AmB and both DMPC and water molecules. Our studies show that AmB preferentially took a vertical position, perpendicular to the membrane surface, with no propensity to enter the membrane. Our finding may suggest that a single AmB molecule entering the membrane is very unlikely.Figure The figure presents the whole structure of the system simulated—starting point. AmB is presented as a space-filling model, DMPC molecules—green sticks, water molecules—red sticks     

Molecular modeling and biophysical analysis of the c-MYC NHE-III<Subscript>1</Subscript> silencer element

G-Quadruplex and i-Motif-forming sequences in the promoter regions of several oncogenes show promise as targets for the regulation of oncogenes. In this study, molecular models were created for the c-MYC NHE-III₁ (nuclease hypersensitivity element III₁) from two 39-base complementary sequences. The NHE modeled here consists of single folded conformers of the polypurine intramolecular G-Quadruplex and the polypyrimidine intramolecular i-Motif structures, flanked by short duplex DNA sequences. The G-Quadruplex was based on published NMR structural data for the c-MYC 1:2:1 loop isomer. The i-Motif structure is theoretical (with five cytosine–cytosine pairs), where the central intercalated cytosine core interactions are based on NMR structural data obtained for a tetramolecular [d(A₂C₄)₄] model i-Motif. The loop structures are in silico predictions of the c-MYC i-motif loops. The porphyrin meso-tetra(N-methyl-4-pyridyl)porphine (TMPyP4), as well as the ortho and meta analogs TMPyP2 and TMPyP3, were docked to six different locations in the complete c-MYC NHE. Comparisons are made for drug binding to the NHE and the isolated G-Quadruplex and i-Motif structures. NHE models both with and without bound cationic porphyrin were simulated for 100 ps using molecular dynamics techniques, and the non-bonded interaction energies between the DNA and porphyrins calculated for all of the docking interactions. Figure Molecular models of the average structure of the final 20 ps of the molecular dynamics simulation of the c-MYC NHE-III₁ (nuclease hypersensitivity element III₁) “silencer” element. The G-Quadruplex structure is at the top-center, and the i-Motif is at the bottom-center of each picture. a “Rotation #1” of the G-Quadruplex, with the T15 loop at the top and rear and the G19/A20 loop at the top and front of the picture. b “Rotation #2” of the G-Quadruplex, with the T15 loop at the top and front of the image, and the G19/A20 loop at the front and adjacent to the G-Quadruplex/i-Motif interface     

Structure prediction and network analysis of chitinases from the Cape sundew,Drosera capensis

BackgroundCarnivorous plants possess diverse sets of enzymes with novel functionalities applicable to biotechnology, proteomics, and bioanalytical research. Chitinases constitute an important class of such enzymes, with future applications including human-safe antifungal agents and pesticides. Here, we compare chitinases from the genome of the carnivorous plant Drosera capensis to those from related carnivorous plants and model organisms.MethodsUsing comparative modeling, in silico maturation, and molecular dynamics simulation, we produce models of the mature enzymes in aqueous solution. We utilize network analytic techniques to identify similarities and differences in chitinase topology.ResultsHere, we report molecular models and functional predictions from protein structure networks for eleven new chitinases from D. capensis, including a novel class IV chitinase with two active domains. This architecture has previously been observed in microorganisms but not in plants. We use a combination of comparative and de novo structure prediction followed by molecular dynamics simulation to produce models of the mature forms of these proteins in aqueous solution. Protein structure network analysis of these and other plant chitinases reveal characteristic features of the two major chitinase families.General significanceThis work demonstrates how computational techniques can facilitate quickly moving from raw sequence data to refined structural models and comparative analysis, and to select promising candidates for subsequent biochemical characterization. This capability is increasingly important given the large and growing body of data from high-throughput genome sequencing, which makes experimental characterization of every target impractical.     

Structural investigations into the binding mode of novel neolignans Cmp10 and Cmp19 microtubule stabilizers by in silico molecular docking,molecular dynamics,and binding free energy calculations

Microtubule stabilizers provide an important mode of treatment via mitotic cell arrest of cancer cells. Recently, we reported two novel neolignans derivatives Cmp10 and Cmp19 showing anticancer activity and working as microtubule stabilizers at micromolar concentrations. In this study, we have explored the binding site, mode of binding, and stabilization by two novel microtubule stabilizers Cmp10 and Cmp19 using in silico molecular docking, molecular dynamics (MD) simulation, and binding free energy calculations. Molecular docking studies were performed to explore the β-tubulin binding site of Cmp10 and Cmp19. Further, MD simulations were used to probe the β-tubulin stabilization mechanism by Cmp10 and Cmp19. Binding affinity was also compared for Cmp10 and Cmp19 using binding free energy calculations. Our docking results revealed that both the compounds bind at Ptxl binding site in β-tubulin. MD simulation studies showed that Cmp10 and Cmp19 binding stabilizes M-loop (Phe272-Val288) residues of β-tubulin and prevent its dynamics, leading to a better packing between α and β subunits from adjacent tubulin dimers. In addition, His229, Ser280 and Gln281, and Arg278, Thr276, and Ser232 were found to be the key amino acid residues forming H-bonds with Cmp10 and Cmp19, respectively. Consequently, binding free energy calculations indicated that Cmp10 (?113.655 kJ/mol) had better binding compared to Cmp19 (?95.216 kJ/mol). This study provides useful insight for better understanding of the binding mechanism of Cmp10 and Cmp19 and will be helpful in designing novel microtubule stabilizers.     

The Computational Studies on the Conformation-Stabilizing Factors of the Left-Handed Z-DNA

Abstract

The conformation-stabilizing factors of the left-handed Z-DNA was estimated by the molecular dynamics calculations using the various metal cations and the polyamines. We concluded the conformation-stabilizing effect of metal cation for Z-DNA was more significant than that of the polyamine by theoretically.     

Understanding the dual mechanism of bioactive peptides targeting the enzymes involved in Renin Angiotensin System (RAS): An in-silico approach

Abstract

Understanding the dual inhibition mechanism of food derivative peptides targeting the enzymes (Renin and Angiotensin Converting enzyme) in the Renin Angiotensin System. Two peptides RALP and WYT were reported to possess antihypertensive activity targeting both renin and ACE, and we have used molecular docking and molecular dynamics simulation, in order to understand the underlying mechanism. The selected peptides (RALP and WYT) from the series of peptides reported were docked to renin and ACE and two binding modes were selected based on the binding energy, interaction pattern and clusters of docking simulation. The enzyme-peptide complexes for renin and ACE (Renin/RALP_1,2; ACE/RALP_1,2; Renin/WYT_1,2 and ACE/WYT_1,2) were subjected to molecular dynamics simulation. Our results identified that the peptides inhibiting renin, tends to move out of the binding pockets (S1’ S2’) which is critical for potent binding and occupies the less important pockets (S4 and S3). This could possibly be the reason for its low potency. Whereas, the same peptides targeting ACE, tends to be intact in the pocket because of the metal ion coordination and there is an ample room to improve on its efficacy. Our results further pave way for the biochemist, medicinal chemist to design dual peptides targeting the RAS effectively.

Communicated by Ramaswamy H. Sarma