Classification of auxin plant hormones by interaction property similarity indices

Although auxins were the first type of plant hormone to be identified, little is known about the molecular mechanism of this important class of plant hormones. We present a classification of a set of about 50 compounds with measured auxin activities, according to their interaction properties. Four classes of compounds were defined: strongly active, weakly active with weak antiauxin behaviour, inactive and inhibitory. All compounds were modeled in two low-energy conformations, P and T, so as to obtain the best match to the planar and tilted conformations, respectively, of indole 3-acetic acid. Each set of conformers was superimposed separately using several different alignment schemes. Molecular interaction energy fields were computed for each molecule with five different chemical probes and then compared by computing similarity indices. Similarity analysis showed that the classes are on average distinguishable, with better differentiation achieved for the T conformers than the P conformers. This indicates that the T conformation might be the active one. Further, a screening was developed which could distinguish compounds with auxin activity from inactive compounds and most antiauxins using the T conformers. The classifications rationalize ambiguities in activity data found in the literature and should be of value in predicting the activities of new plant growth substances and herbicides.     

Putting molecular similarity into context: asymmetric indices for field-based similarity measures

Some of the most widely used indices in molecular similarity searching are intrinsically symmetric in nature, meaning that each molecule under comparison contributes equally to the similarity index. For example, the Carbó and the Hodgkin–Richards similarity indices are respectively, related to the geometric and arithmetic averages of the molecular self-similarities. This work introduces the asymmetric forms of an entire family of field-based molecular similarity indices. By incorporating a weighted contribution of each molecule into the similarity index, the newly obtained asymmetric forms allow for measuring and modulating the similarity of one molecule in the context of another and thus have the potential of alleviating the size dependency often observed in chemical similarity searching     

Cassandra: An open source Monte Carlo package for molecular simulation

Cassandra is an open source atomistic Monte Carlo software package that is effective in simulating the thermodynamic properties of fluids and solids. The different features and algorithms used in Cassandra are described, along with implementation details and theoretical underpinnings to various methods used. Benchmark and example calculations are shown, and information on how users can obtain the package and contribute to it are provided. © 2017 Wiley Periodicals, Inc.     

Monte Carlo method based QSAR modelling of natural lipase inhibitors using hybrid optimal descriptors

Obesity is one of the most provoking health burdens in the developed countries. One of the strategies to prevent obesity is the inhibition of pancreatic lipase enzyme. The aim of this study was to build QSAR models for natural lipase inhibitors by using the Monte Carlo method. The molecular structures were represented by the simplified molecular input line entry system (SMILES) notation and molecular graphs. Three sets – training, calibration and test set of three splits – were examined and validated. Statistical quality of all the described models was very good. The best QSAR model showed the following statistical parameters: r² = 0.864 and Q² = 0.836 for the test set and r² = 0.824 and Q² = 0.819 for the validation set. Structural attributes for increasing and decreasing the activity (expressed as pIC₅₀) were also defined. Using defined structural attributes, the design of new potential lipase inhibitors is also presented. Additionally, a molecular docking study was performed for the determination of binding modes of designed molecules.     

Conformational behavior of antineoplastic peptides Dolastatin 10 and Dolastatin 15 from Monte Carlo and molecular dynamics simulations

In this work, we report the results of a combined Monte Carlo (MC) and molecular dynamics (MD) approach applied to the conformational study of natural potent mitosis inhibitor Dolastatin 15. A GRID and CLogP analysis has been performed with the aim to provide some hints to the interpretation of the different behavior of Dolastatin 15 and Dolastatin 10, which has been examined in a former study. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Wavelet-accelerated Monte Carlo sampling of polymer chains

We introduce a new method for coarse-graining polymer chains, based on the wavelet transform, a multiresolution data analysis technique. This method, which assigns a cluster of particles to a coarse-grained bead located at the center of mass of the cluster, reduces the complexity of the problem significantly by dividing the simulation into several stages, each with a small fraction of the number of beads in the overall chain. At each stage, we compute the distributions of coarse-grained internal coordinates as well as potential functions required for subsequent simulation stages. We show that, with this wavelet-accelerated Monte Carlo method, coarse-grained Gaussian and self-avoiding random walks can reproduce results obtained from atomistic simulations to a high degree of accuracy in orders of magnitude less time. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 897–910, 2005     

A scalable parallel Monte Carlo method for free energy simulations of molecular systems

We present a method of parallelizing flat histogram Monte Carlo simulations, which give the free energy of a molecular system as an output. In the serial version, a constant probability distribution, as a function of any system parameter, is calculated by updating an external potential that is added to the system Hamiltonian. This external potential is related to the free energy. In the parallel implementation, the simulation is distributed on to different processors. With regular intervals the modifying potential is summed over all processors and distributed back to every processor, thus spreading the information of which parts of parameter space have been explored. This implementation is shown to decrease the execution time linearly with added number of processors.     

MCDOCK: A Monte Carlo simulation approach to the molecular docking problem

Prediction of the binding mode of a ligand (a drug molecule) to its macromolecular receptor, or molecular docking, is an important problem in rational drug design. We have developed a new docking method in which a non-conventional Monte Carlo (MC) simulation technique is employed. A computer program, MCDOCK, was developed to carry out the molecular docking operation automatically. The current version of the MCDOCK program (version 1.0) allows for the full flexibility of ligands in the docking calculations. The scoring function used in MCDOCK is the sum of the interaction energy between the ligand and its receptor, and the conformational energy of the ligand. To validate the MCDOCK method, 19 small ligands, the binding modes of which had been determined experimentally using X-ray diffraction, were docked into their receptor binding sites. To produce statistically significant results, 20 MCDOCK runs were performed for each protein–ligand complex. It was found that a significant percentage of these MCDOCK runs converge to the experimentally observed binding mode. The root-mean-square (rms) of all non-hydrogen atoms of the ligand between the predicted and experimental binding modes ranges from 0.25 to 1.84 Å for these 19 cases. The computational time for each run on an SGI Indigo2/R10000 varies from less than 1 min to 15 min, depending upon the size and the flexibility of the ligands. Thus MCDOCK may be used to predict the precise binding mode of ligands in lead optimization and to discover novel lead compounds through structure-based database searching.     

Application of combined Monte Carlo and molecular dynamics method to simulation of dipalmitoyl phosphatidylcholine lipid bilayer

We describe a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer. The procedure consists of alternating molecular dynamics trajectory calculations in a constant surface tension and temperature ensemble with configurational bias Monte Carlo moves to different regions of the configuration space of the bilayer, in a constant volume and temperature ensemble. The procedure is described in detail and is applied to a bilayer of 100 molecules of dipalmitoyl phosphatidylcholine (DPPC) and 3205 water molecules. We find that the hybrid simulation procedure enhances the equilibration of the bilayer as measured by the convergence of the area per molecule and the segmental order parameters, as compared with a simulation using only molecular dynamics (MD). Progress toward equilibration is almost three times as fast in CPU time, compared with a purely MD simulation. Equilibration is complete, as judged by the lack of energy drift in three separate 200-ps runs of continuous MD started from different initial states. Results of the simulation are presented and compared with experimental data and with other recent simulations of DPPC. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1153–1164, 1999     

流体混合物超额性质的分子模拟研究

提出了逐级取样分子模拟方法。通过逐个改变粒子的类别实现在一次模拟中对混合物整个组成范围的取样,并提出适当的处理方法可以同时获取混合物在整个组成范围内的多种超额热力学性质。用本文方法在零压、115.8K下对氩-氪体系作了模拟计算,与McDonald采用常规NPTMonteCarlo模拟结果的比较表明,本文的逐级取样模拟方法是可靠而有效的。     

Monte Carlo study of electron correlation functions for small molecules

A study is made of electron-electron correlation functions for use in trial wave functions for small molecules. New forms are proposed that have only a few variational parameters, and these parameters have physical meanings that are easily discerned. Total energies for H₂, LiH and Li₂ computed using these correlation functions are presented, and comparison is made with previous forms, including the Jastrow-Pade form often used in Monte Carlo studies. We further treat the possibility that correlation depends not only on the separation of a pair of electrons but also on the location of the electron pair relative to the nuclei — indicative of a density-dependent or many body correlation effect. Our results indicate that such a many-body correlation effect is weakly present.This work was supported by the Director, Office of Energy Research, Office of Basic Energy Sciences, Chemical Sciences Division of the U.S. Department of Energy under Contract No. DE-AC03-76SF00098     

Comparative molecular dynamics and Monte Carlo study of statistical properties for coarse-grained heteropolymers

Employing a simple hydrophobic-polar heteropolymer model, we compare thermodynamic quantities obtained from Andersen and Nosé-Hoover molecular dynamics as well as replica-exchange Monte Carlo methods. We find qualitative correspondence in the results, but serious quantitative differences using the Nosé-Hoover chain thermostat. For analyzing the deviations, we study different parameterizations of the Nosé-Hoover chain thermostat. Autocorrelations from molecular dynamics and Metropolis Monte Carlo runs are also investigated.     

Rovibrationally averaged properties of H2 using Monte Carlo methods

Using variational Monte Carlo and a simple explicitly correlated wave function, we have computed 18 molecular properties of the hydrogen molecule (X¹∑) at 24 internuclear distances. These properties have been combined with rapidly convergent rovibrational wave functions to produce rovibrationally averaged properties for several of the lowest rotational and vibrational levels of this system. Our results are in very good agreement with previous values found in the literature. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Monte Carlo simulation of the local ordering of water molecules. I. Binary spatial correlations

The pair distribution functions (DFs) of the oxygen and hydrogen atoms over space around one of the water molecules were calculated during Monte Carlo simulation (Metropolis procedure for the NVT ensemble at normal density and 300 K). An analysis of the isosurfaces of the DF_OO and DF_OH constructed around the selected water molecule for several fixed values of local density allowed us to obtain detailed information about the most probable localization of water molecules in the second and third coordination spheres.     

Quantum Monte Carlo calculations of bond dissociation energies for some nitro and amino molecules

Bond dissociation energies (BDEs) for some nitro or amino contained prototypical molecules in energetic materials are computed by fixed‐node diffusion quantum Monte Carlo method. The nodes are determined from a Slater determinant calculated within density functional theory at the B3LYP/6‐311G** level. The possible errors, the nodal error, and the cancellation of nodal errors in calculating BDE are discussed, and the accuracy is compared with other available ab initio computations and experimental results. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Monte Carlo simulations of a liquid crystal copolymer in the solid state

The condensed phase of the alternating copolyester of p-hydroxybenzoic acid (HBA) and 2-hydroxy-6-naphthoic acid (HNA) is investigated by studying the room temperature packing arrangement of the copolymer chains. A molecular modeling methodology is employed with a Monte Carlo sampling of the configurational phase space. Realistic poly(HBA-alt-HNA) polymer chains are represented by an explicit atom representation of the HBA/HNA dimers. States are sampled from the NVT ensemble using a sampling scheme consisting of (1) valence and torsional variations, (2) rigid body rotations of the chain about the chain axis, and (3) rigid body translations of the chain. The effect of chain packing on the conformation of chains, as well as the relative intra- and intermolecular orientations of aromatic rings, is investigated. Correlation of chain positioning along the chain axis is dominated by aromatic rings maintaining a center-to-center plane of registry. These layers of aromatic units pack with a preference for edge-to-face orientations in a herringbone-type pattern and have an intermolecular ring angle between the pairs of aromatic rings in the unit cell that is ca. 68°. The aromatic rings, on average, are rotated 38° out from the b–c plane. The phenylene rings of these copolyesters are less restricted in their relative orientation in comparison to the naphthalene rings. Intramolecular orientational probability density distributions indicate a preference for staggering the successive aromatic rings along the chain, with a staggering angle of ca. 66°. © 1998 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 36: 727–741, 1998