PEPCAT—A new tool for conformational analysis of peptides

We report a new technique for the efficient analysis and visualization of peptide and protein conformations and conformational relationships, which we have implemented in a computer program called PEPCAT. PEPCAT (an abbreviation for Peptide Conformational Analysis Tool) provides a simple, graphical, and flexible framework that allows the user to define a specific structural feature or juxtaposition of amino acids and to follow the fate of the motif during a molecular dynamics simulation. Here we describe the PEPCAT analysis of the effects of environmental and chemical modifications on conformational preferences of a regulator of hemopoiesis, namely the pentapeptide pyro‐EEDCK, and of a conformational transition in the immunosuppressant drug cyclosporin A. PEPCAT, however, can be applied to the conformational analysis of peptides and proteins in general. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 446–461, 2000     

Cluster and Molecular Mechanical Analysis of Cobalt(III) 14-Membered Macrocyclic Complexes

The configurational distribution of all the cobalt(III) complexes containing a 14-membered tetra-aza macrocyclic backbone in the Cambridge Structural Database was determined by cluster analysis, and the reason for some of the complexes adopting their configuration was established by molecular mechanics.     

CENCALC: A computational tool for conformational entropy calculations from molecular simulations

We present the CENCALC software that has been designed to estimate the conformational entropy of single molecules from extended Molecular Dynamics (MD) simulations in the gas‐phase or in solution. CENCALC uses both trajectory coordinates and topology information in order to characterize the conformational states of the molecule of interest by discretizing the time evolution of internal rotations. The implemented entropy methods are based on the mutual information expansion, which is built upon the converged probability density functions of the individual torsion angles, pairs of torsions, triads, and so on. Particularly, the correlation‐corrected multibody local approximation selects an optimum cutoff in order to retrieve the maximum amount of genuine correlation from a given MD trajectory. We illustrate these capabilities by carrying out conformational entropy calculations for a decapeptide molecule either in its unbound form or in complex with a metalloprotease enzyme. CENCALC is distributed under the GNU public license at http://sourceforge.net/projects/cencalc/ .     

IMiCMO: a new integrated ab initio multicenter molecular orbitals method for molecular dynamics calculations in solvent cluster systems

A new computational scheme integrating ab initio multicenter molecular orbitals for determining forces of individual atoms in large cluster systems is presented. This method can be used to treat systems of many molecules, such as solvents by quantum mechanics. The geometry parameters obtained for three models of water clusters by the present method are compared with those obtained by the full ab initio MO method. The results agree very well. The scheme proposed in this article also intended for use in modeling cluster systems using parallel algorithms. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1107–1112, 2001     

A new soft‐core potential function for molecular dynamics applied to the prediction of protein loop conformations

We have developed a new soft‐core potential function for the conformational search of complex systems with molecular dynamics. The potential function was designed to maintain the main equilibrium properties of the original force field, which means that the soft‐core potential gives physically realistic performance also without additional restraints, different from most of the previous soft‐core potential functions. The performance of the method was demonstrated by applying it to the problem of finding native conformations for protein loops. Short loops from neocarzinostatin and parvalbumin were used as the first test cases. The use of the new soft‐core potential function was shown to improve significantly the performance of molecular dynamics in the search of the native conformation of protein loops. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 388–397, 2000     

Solvent effect on the synthesis of clarithromycin: A molecular dynamics study

Clarithromycin (6-O-methylerythromycin A) is a 14-membered macrolide antibiotic which is active in vitro against clinically important gram-positive and gram-negative bacteria. The selectivity of the methylation of the C-6 OH group is studied on erythromycin A derivatives. To understand the effect of the solvent on the methylation process, detailed molecular dynamics (MD) simulations are performed in pure DMSO, pure THF and DMSO:THF (1:1) mixture by using the anions at the C-6, C-11 and C-12 positions of 2',4"-[O-bis(TMS)]erythromycin A 9-[O-(dimethylthexylsilyl)oxime] under the assumption that the anions are stable on the sub-nanosecond time scale. The conformations of the anions are not affected by the presence of the solvent mixture. The radial distribution functions are computed for the distribution of different solvent molecules around the 'O-' of the anions. At distances shorter than 5 A, DMSO molecules are found to cluster around the C-11 anion, whereas the anion at the C-12 position is surrounded by the THF molecules. The anion at the C-6 position is not blocked by the solvent molecules. The results are consistent with the experimental finding that the methylation yield at the latter position is increased in the presence of a DMSO:THF (1:1) solvent mixture. Thus, the effect of the solvent in enhancing the yield during the synthesis is not by changing the conformational properties of the anions, but rather by creating a suitable environment for methylation at the C-6 position.     

GridMAT‐MD: A grid‐based membrane analysis tool for use with molecular dynamics

GridMAT‐MD is a new program developed to aid in the analysis of lipid bilayers from molecular dynamics simulations. It reads a GROMACS coordinate file and generates two types of data: a two‐dimensional contour plot depicting membrane thickness, and a polygon‐based tessellation of the individual lipid headgroups. GridMAT‐MD can also account for proteins or small molecules within the headgroups of the lipids, closely approximating their occupied lateral area. The program requires no installation, is fast, and is freely available. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

A method for predicting protein conformational pathways by using molecular dynamics simulations guided by difference distance matrices

Here, an efficient method that predicts natural transition pathways between two endpoint states of an allosteric protein has been proposed. This method helps create structures that bridge these endpoints through multiple iterative and unbiased molecular dynamics simulations with explicit water. Difference distance matrices provide an approach for identifying states involving concerted slow motion. A series of structures are readily generated along the transition pathways of adenylate kinase. Predicted structures may be useful for an initial pathway to evaluate free energy landscapes via umbrella sampling and chain‐of‐states methods. © 2016 Wiley Periodicals, Inc.     

Estimation of active conformations of drugs by a new molecular superposing procedure

We have developed a new program, SUPERPOSE, to superpose two molecules based on the physicochemical properties of functional atoms within individual molecules. SUPERPOSE treats a pseudo-molecule consisting of functional atoms instead of a real molecule. Four types of physicochemical properties – hydrophobicity, presence of a hydrogen-bonding donor, presence of a hydrogen-bonding acceptor and presence of a hydrogen-bonding donor/acceptor – were supposed and a score was given to each overlap. When functional atoms with the same physicochemical properties were overlapped, points were added to the score, and when the functional atoms with different physicochemical properties were overlapped, points were subtracted. We applied SUPERPOSE to 12 pairs of 24 enzyme inhibitors and found that the best scored overlay for each inhibitor pair could successfully reproduce the superposition obtained from X-ray crystallography. Next, we applied SUPERPOSE to estimate the active conformations of the thrombin inhibitors MQPA, 4-TAPAP and NAPAP. Superpositions of conformers sampled by the high-temperature molecular dynamics calculation with respect to the three inhibitors were performed, and 13 sets of conformers having the best common overlay to the three inhibitors were selected. One among 13 sets was consistent with the superposition of the active conformations derived from the X-ray crystallography of the thrombin–inhibitor complexes.     

Molecular Dynamics Simulations of Cellulose Oligomers: Conformational Analysis

Summary: The results of classical molecular simulations of cellulose oligomers are presented here. The conformations of the chains in the high temperature melt, room temperature quenched melt and gas phase are compared with respect to various geometrical parameters including square end‐to‐end distances, glycosidic link torsion correlations, ring puckering and hydrogen bonding. The cellulose oligomer melts were relaxed at 800 K with molecular dynamics, and then cooled down in three different ways to obtain dense amorphous systems at 500 K and at room temperature. The sample resulting from the quench (step) shows too much similarity with the melt at 800 K. The two other cooling schemes (ramp, 2ramps) give very similar results for all quantities investigated. The relevance of previous single molecule calculations with respect to the dense amorphous systems is called into question. Comparisons between the chains in the dense systems and those in the gas phase reveal that, even for these relatively short stiff chains, differences exist in the preferred conformations. At high temperatures, where both systems are in equilibrium, the distribution of square end‐to‐end distances are both fairly smooth, but the gas phase clearly prefers more compact conformations. At 300 K, the differences are exacerbated as the equilibrium distribution for the gas phase shows a high proportion of folded conformers, whereas the nonequilibrium quenched systems necessarily retain the extended envelope of the higher temperature. Differences are also evident in the puckering, the rotation of the hydroxymethyl groups and the pattern of hydrogen bonds.

The probability density distribution for the square end‐to‐end distance for octaose in the gas phase (light line) and in the dense phase (dark line) at 300 K.     

Cluster analysis of molecular simulation trajectories for systems where both conformation and orientation of the sampled states are important

Clustering methods have been widely used to group together similar conformational states from molecular simulations of biomolecules in solution. For applications such as the interaction of a protein with a surface, the orientation of the protein relative to the surface is also an important clustering parameter because of its potential effect on adsorbed‐state bioactivity. This study presents cluster analysis methods that are specifically designed for systems where both molecular orientation and conformation are important, and the methods are demonstrated using test cases of adsorbed proteins for validation. Additionally, because cluster analysis can be a very subjective process, an objective procedure for identifying both the optimal number of clusters and the best clustering algorithm to be applied to analyze a given dataset is presented. The method is demonstrated for several agglomerative hierarchical clustering algorithms used in conjunction with three cluster validation techniques. © 2016 Wiley Periodicals, Inc.     

Conformational analysis of synthetic neolignans active against leishmaniasis,using the molecular mechanics method (MM2)

Conformational analysis of 20 neolignans was performed to determine the most probable conformer that may fit the receptor. The molecular mechanics method (MM2) was employed to construct conformational maps in both a vacuum and a biological environment. Boltzmann's distribution among several local minima was calculated. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18: 712–721, 1997     

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.     

Conformational analysis of sesquiterpene lactones of the germacrane type 1. 1(10)E,4E-Germacranolides

The structure of 1(10)E,4E-germacranolides was studied by a molecular mechanics method. Possible conformers, the probabilities of their occurrence, and the barriers to conformational transitions were determined.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No, 12, pp, 2890–2893, December, 1996.     

A molecular mechanics and semiempirical conformational analysis of the herbicide diuron inhibitor of photosystem II

The potential energy surface (PES) for the herbicide diuron (DCMU), a photosystem II inhibitor, has been extensively investigated using the quantum-mechanical semiempirical molecular orbital methods AM1 and PM3 and molecular mechanics method. A detailed conformational search has been carried out which revealed the occurrence of four genuine minimum energy structures. The relative stability of the conformers and rotational barriers to conformational interconversion were evaluated using distinct theoretical approaches. The results showed that thetrans form of the diuron molecule is more stable than thecis form in all methods, and so it may possibly be the biologically active isomer.     

MkVsites: A tool for creating GROMACS virtual sites parameters to increase performance in all-atom molecular dynamics simulations

The absolute performance of any all-atom molecular dynamics simulation is typically limited by the length of the individual timesteps taken when integrating the equations of motion. In the GROMACS simulation software, it has for a long time been possible to use so-called virtual sites to increase the length of the timestep, resulting in a large gain of simulation efficiency. Up until now, support for this approach has in practice been limited to the standard 20 amino acids however, shrinking the applicability domain of virtual sites. MkVsites is a set of python tools which provides a convenient way to obtain all parameters necessary to use virtual sites for virtually any molecules in a simulation. Required as input to MkVsites is the molecular topology of the molecule(s) in question, along with a specification of where to find the parent force field. As such, MkVsites can be a very valuable tool suite for anyone who is routinely using GROMACS for the simulation of molecular systems.     

A simple erosion dynamics model of molecular sputter depth profiling

A simple model which describes the essential features commonly observed in a molecular sputter depth profile is presented. General predictions of the dependence of measured molecular ion signals on the primary ion fluence are derived for the specific case where a mass spectrometric technique such as SIMS or secondary neutral mass spectrometry (SNMS) is used to analyze the momentary surface. The results are compared with recent experimental data on molecular depth profiles obtained by cluster‐ion‐initiated SIMS of organic overlayers. Copyright © 2008 John Wiley & Sons, Ltd.     

MDAnalysis: A toolkit for the analysis of molecular dynamics simulations

MDAnalysis is an object‐oriented library for structural and temporal analysis of molecular dynamics (MD) simulation trajectories and individual protein structures. It is written in the Python language with some performance‐critical code in C. It uses the powerful NumPy package to expose trajectory data as fast and efficient NumPy arrays. It has been tested on systems of millions of particles. Many common file formats of simulation packages including CHARMM, Gromacs, Amber, and NAMD and the Protein Data Bank format can be read and written. Atoms can be selected with a syntax similar to CHARMM's powerful selection commands. MDAnalysis enables both novice and experienced programmers to rapidly write their own analytical tools and access data stored in trajectories in an easily accessible manner that facilitates interactive explorative analysis. MDAnalysis has been tested on and works for most Unix‐based platforms such as Linux and Mac OS X. It is freely available under the GNU General Public License from http://mdanalysis.googlecode.com . © 2011 Wiley Periodicals, Inc. J Comput Chem 2011     

LINCS: A linear constraint solver for molecular simulations

In this article, we present a new LINear Constraint Solver (LINCS) for molecular simulations with bond constraints. The algorithm is inherently stable, as the constraints themselves are reset instead of derivatives of the constraints, thereby eliminating drift. Although the derivation of the algorithm is presented in terms of matrices, no matrix matrix multiplications are needed and only the nonzero matrix elements have to be stored, making the method useful for very large molecules. At the same accuracy, the LINCS algorithm is three to four times faster than the SHAKE algorithm. Parallelization of the algorithm is straightforward. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1463–1472, 1997