Essential dynamics/factor analysis for the interpretation of molecular dynamics trajectories

Subject of this work is the analysis of molecular dynamics (MD) trajectories of neurophysins I (NPI) and II (NPII) and their complexes with the neurophyseal nonapeptide hormones oxytocin (OT) and vasopresssin (VP), respectively, simulated in water. NPs serve in the neurosecretory granules as carrier proteins for the hormones before their release to the blood. The starting data consisted of two pairs of different trajectories for each of the (NPII/VP)2 and (NPI/OT)2 heterotetramers and two more trajectories for the NPII2 and NPI2 homodimers (six trajectories in total). Using essential dynamics which, to our judgement, is equivalent to factor analysis, we found that only about 10 degrees of freedom per trajectory are necessary and sufficient to describe in full the motions relevant for the function of the protein. This is consistent with these motions to explain about 90% of the total variance of the system. These principal degrees of freedom represent slow anharmonic motional modes, clearly pointing at distinguished mobility of the atoms involved in the protein's functionality.     

Interactive essential dynamics

Essential dynamics(ED) is a useful method for analyzing trajectories generated by molecular dynamics (MD), but current tools are awkward to use, limiting the usefulness of the technique. This paper describes a new interactive graphical interface for visualization of ED results, including filtering a trajectory on an arbitrary set of eigenvectors and manipulation of a structure's projection along any eigenvector.     

A hybrid method of molecular dynamics and harmonic dynamics for docking of flexible ligand to flexible receptor

We have developed a new docking method to consider receptor flexibility, a hybrid method of molecular dynamics and harmonic dynamics. The global motions of the whole receptor were approximately introduced into those of the receptor in the docking simulation as harmonic dynamics. On the other hand, the local flexibility of the side chains was also considered by conventional molecular dynamics. We confirmed that this new method can reproduce the fluctuations of the whole receptor by making a comparison of the directions and amplitudes of the global fluctuations. Then this method was applied to the docking of HIV-1 protease and its ligand. As a result, we observed a docking process where the ligand enters into the binding pocket well, which implies that this method is effective enough to reproduce a molecular complex formation.     

Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis

ABSTRACT

It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.     

Grcarma: A fully automated task‐oriented interface for the analysis of molecular dynamics trajectories

We report the availability of grcarma, a program encoding for a fully automated set of tasks aiming to simplify the analysis of molecular dynamics trajectories of biological macromolecules. It is a cross‐platform, Perl/Tk‐based front‐end to the program carma and is designed to facilitate the needs of the novice as well as those of the expert user, while at the same time maintaining a user‐friendly and intuitive design. Particular emphasis was given to the automation of several tedious tasks, such as extraction of clusters of structures based on dihedral and Cartesian principal component analysis, secondary structure analysis, calculation and display of root‐mean‐square deviation (RMSD) matrices, calculation of entropy, calculation and analysis of variance–covariance matrices, calculation of the fraction of native contacts, etc. The program is free‐open source software available immediately for download. © 2013 Wiley Periodicals, Inc.     

Pteros: Fast and easy to use open‐source C++ library for molecular analysis

An open‐source Pteros library for molecular modeling and analysis of molecular dynamics trajectories for C++ programming language is introduced. Pteros provides a number of routine analysis operations ranging from reading and writing trajectory files and geometry transformations to structural alignment and computation of nonbonded interaction energies. The library features asynchronous trajectory reading and parallel execution of several analysis routines, which greatly simplifies development of computationally intensive trajectory analysis algorithms. Pteros programming interface is very simple and intuitive while the source code is well documented and easily extendible. Pteros is available for free under open‐source Artistic License from http://sourceforge.net/projects/pteros/ . © 2012 Wiley Periodicals, Inc.     

PuCO体系的分子反应动力学研究

基于多体项展式理论方法导出的PuCO分子基态(X^7A")的分析势能函数,用准经典的Monte-Carlo轨线法对Pu(^7Fg)+CO(0,0)和O(^3Pg)+PuC(0,0)的分子反应动力学过程进行了计算。结果表明：Pu(^7Fg)与CO(0,0)碰撞易生成PuCO配合物分子,该反应是无阈能反应,反应截面σ随能量Et的升高而下降,当Et=502.1kJ.mol^-^1时,σ几乎为零。O(^3Pg)与PuC(0,0)碰撞易发生生成Pu+CO的交换反应,该反应无阈能。     

Single‐pass incremental force updates for adaptively restrained molecular dynamics

Adaptively restrained molecular dynamics (ARMD) allows users to perform more integration steps in wall‐clock time by switching on and off positional degrees of freedoms. This article presents new, single‐pass incremental force updates algorithms to efficiently simulate a system using ARMD. We assessed different algorithms for speedup measurements and implemented them in the LAMMPS MD package. We validated the single‐pass incremental force update algorithm on four different benchmarks using diverse pair potentials. The proposed algorithm allows us to perform simulation of a system faster than traditional MD in both NVE and NVT ensembles. Moreover, ARMD using the new single‐pass algorithm speeds up the convergence of observables in wall‐clock time. © 2017 Wiley Periodicals, Inc.     

New insights into the meaning and usefulness of principal component analysis of concatenated trajectories

A comparison between different conformations of a given protein, relating both structure and dynamics, can be performed in terms of combined principal component analysis (combined‐PCA). To that end, a trajectory is obtained by concatenating molecular dynamics trajectories of the individual conformations under comparison. Then, the principal components are calculated by diagonalizing the correlation matrix of the concatenated trajectory. Since the introduction of this approach in 1995 it has had a large number of applications. However, the interpretation of the eigenvectors and eigenvalues so obtained is based on intuitive foundations, because analytical expressions relating the concatenated correlation matrix with those of the individual trajectories under consideration have not been provided yet. In this article, we present such expressions for the cases of two, three, and an arbitrary number of concatenated trajectories. The formulas are simple and show what is to be expected and what is not to be expected from a combined‐PCA. Their correctness and usefulness is demonstrated by discussing some representative examples. The results can be summarized in a simple sentence: the correlation matrix of a concatenated trajectory is given by the average of the individual correlation matrices plus the correlation matrix of the individual averages. From this it follows that the combined‐PCA of trajectories belonging to different free energy basins provides information that could also be obtained by alternative and more straightforward means. © 2014 Wiley Periodicals, Inc.     

Rotamer decomposition and protein dynamics: Efficiently analyzing dihedral populations from molecular dynamics

Molecular mechanics methods have matured into powerful methods to understand the dynamics and flexibility of macromolecules and especially proteins. As multinanosecond to microsecond length molecular dynamics (MD) simulations become commonplace, advanced analysis tools are required to generate scientifically useful information from large amounts of data. Some of the key degrees of freedom to understand protein flexibility and dynamics are the amino acid residue side chain dihedral angles. In this work, we present an easily automated way to summarize and understand the relevant dihedral populations. A tremendous reduction in complexity is achieved by describing dihedral timeseries in terms of histograms decomposed into Gaussians. Using the familiar and widely studied protein lysozyme, it is demonstrated that our approach captures essential properties of protein structure and dynamics. A simple classification scheme is proposed that indicates the rotational state population for each dihedral angle of interest and allows a decision if a given side chain or peptide backbone fragment remains rigid during the course of an MD simulation, adopts a converged distribution between conformational substates or has not reached convergence yet. © 2012 Wiley Periodicals, Inc.     

Lightweight object oriented structure analysis: Tools for building tools to analyze molecular dynamics simulations

LOOS (Lightweight Object Oriented Structure‐analysis) is a C++ library designed to facilitate making novel tools for analyzing molecular dynamics simulations by abstracting out the repetitive tasks, allowing developers to focus on the scientifically relevant part of the problem. LOOS supports input using the native file formats of most common biomolecular simulation packages, including CHARMM, NAMD, Amber, Tinker, and Gromacs. A dynamic atom selection language based on the C expression syntax is included and is easily accessible to the tool‐writer. In addition, LOOS is bundled with over 140 prebuilt tools, including suites of tools for analyzing simulation convergence, three‐dimensional histograms, and elastic network models. Through modern C++ design, LOOS is both simple to develop with (requiring knowledge of only four core classes and a few utility functions) and is easily extensible. A python interface to the core classes is also provided, further facilitating tool development. © 2014 Wiley Periodicals, Inc.     

Compressing molecular dynamics trajectories: Breaking the one‐bit‐per‐sample barrier

Molecular dynamics simulations yield large amounts of trajectory data. For their durable storage and accessibility an efficient compression algorithm is paramount. State of the art domain‐specific algorithms combine quantization, Huffman encoding and occasionally domain knowledge. We propose the high resolution trajectory compression scheme (HRTC) that relies on piecewise linear functions to approximate quantized trajectories. By splitting the error budget between quantization and approximation, our approach beats the current state of the art by several orders of magnitude given the same error tolerance. It allows storing samples at far less than one bit per sample. It is simple and fast enough to be integrated into the inner simulation loop, store every time step, and become the primary representation of trajectory data. © 2016 Wiley Periodicals, Inc.     

Thermal analysis via molecular dynamics simulation

Thermal analysis by classical molecular dynamics simulations is discussed on hand of heat capacity of crystals of 9600 atoms. The differences between quantum mechanical and classical mechanical calculations are shown. Anharmonicity is proven to be an important factor. Finally, it is found that defects contribute to an increase in heat capacity before melting. The energy of conformational gauche defects within the crystal is only about 10% due to internal rotation. The other energy must be generated by cooperative strain. The conclusion is that the next generation of faster computers may permit wider use of molecular dynamics simulations in support of the interpretation of thermal analysis.Dedicated to Professor Bernhard Wunderlich on the occasion of his 65th birthday     

StreaMD: Advanced analysis of molecular dynamics using R

Gromacs is one of the most popular molecular simulation suites currently available. In this contribution we present streaMD , the first interface between Gromacs trajectory files and the statistical language R . The amount of data created due to ever increasing computational power renders fast and efficient analysis of trajectories into a challenge. Especially as standard approaches such as root‐mean square fluctuations and the like provide only limited physical insight. In our streaMD package integration of the Gromacs I/O libraries with advanced, graph‐based analysis methods as the java library Stream leads to both: improved speed and analysis depth. We benchmark our results and highlight the applicability of the package by an interesting problem in RNA design, namely the interaction of tetracycline with an aptamer. © 2018 Wiley Periodicals, Inc.     

Dynamic conformational ensembles regulate casein kinase-1 isoforms: Insights from molecular dynamics and molecular docking studies

Casein kinase-1 (CK1) isoforms actively participate in the down-regulation of canonical Wnt signaling pathway; however recent studies have shown their active roles in oncogenesis of various tissues through this pathway. Functional loss of two isoforms (CK1-α/ε) has been shown to activate the carcinogenic pathway which involves the stabilization of of cytoplasmic β-catenin. Development of anticancer therapeutics is very laborious task and depends upon the structural and conformational details of the target. This study focuses on, how the structural dynamics and conformational changes of two CK1 isoforms are synchronized in carcinogenic pathway. The conformational dynamics in kinases is the responsible for their action as has been supported by the molecular docking experiments.     

Wavelet transform analysis of ab initio molecular dynamics simulation: application to core-excitation dynamics of BF3

We propose a novel analysis method of ab initio molecular dynamics (AIMD) simulation using a continuous wavelet transform (c-WT) technique. The c-WT technique, one of the time-frequency signal analysis methods, provides a clear view of the dynamical information in time developments. Combined with the auto-correlation function of velocity by AIMD simulation, c-WT analysis enables us to well understand dynamical distribution, such as the vibrational properties following a change of electronic structure in a molecular system. As a practical application, AIMD simulation of core-excited BF(3) (B1s --> 2a(2) (')) is illustrated. AIMD simulation leads to the change of vibrational motion as well as structural deformation by core-excitation. The c-WT analysis clarifies the relationship between structural deformation and the related significant vibrational modes in core-excitation within 50 fs.     

Software news and updates. Basis-set completeness profiles in two dimensions

A two-electron basis-set completeness profile is proposed by analogy with the one-electron profile introduced by D. P. Chong (Can J Chem 1995, 73, 79). It is defined as Y(alpha, beta) = sigmam sigman (Galpha(1)Gbeta(2)/(1/r12)/ psim(1)psin(2)) (psim(1)psin(2)/r12/Galpha(1)Gp(2)) and motivated by the expression for the basis-set truncation correction that occurs in the framework of explicitly correlated methods (Galpha is a scanning Gaussian-type orbital of exponent alpha and [psim] is the orthonormalized one-electron basis under study). The two-electron basis-set profiles provide a visual assessment of the suitability of basis sets to describe electron-correlation effects. Furthermore, they provide the opportunity to assess the quality of the basis set as a whole--not only of the individual angular momentum subspaces, as is the case for the one-electron basis-set profiles. The two-electron completeness profiles of the cc-pVXZ (X = D, T, Q), aug-cc-pVTZ, cc-pCVTZ, and SVP-auxiliary basis sets for the carbon atom are presented as illustrative examples.     

Self-assembly and molecular dynamics of oligoindenofluorenes.

The optical properties, self-assembly mechanism, thermal properties, and the associated molecular dynamics of a series of stable blue-emitting oligoindenofluorenes up to the polymer were investigated by photoluminescence, wide-angle X-ray scattering, polarizing optical microscopy, differential scanning calorimetry, and dielectric spectroscopy. Oligomers of indenofluorenes possess high structural order and form smectic mesophases, as opposed to the nematic mesophases formed in the corresponding oligofluorenes. The alpha process associated with the liquid-to-glass transition displays an asymmetric broadening of relaxation times (non-Debye process), with a non-Arrhenius temperature dependence of relaxation times and a strong molecular weight dependence. The low dielectric loss and low ionic conductivity revealed the absence of ketone structural defects, which is consistent with the observed stable blue emission.