Computational study on the ring distortions of 1,4-dicyanobenzene in the gas phase and in the crystal

A molecular mechanics software enhanced to perform empirical energy calculations on crystals (KESSHOU) was further developed to handle intermolecular electrostatic interactions as well. The packing of the molecules of 1,4-dicyanobenzene and 1,4-diisocyanobenzene in the crystal was studied. The role of the van der Waals and the electrostatic interactions in the balance of nonbonded atom-atom interactions is analyzed. The packing forces are dominated by van der Waals forces. The electrostatic interactions have higher stabilizing contribution for the dicyano isomer than for the diisocyano form. The dependence of the results on the size of the crystal, the molecular mechanics force field (MM2 vs MM3), and the dielectric constant are also assessed. Ab initio MP2/6–311G** geometries of the isolated molecules are in accordance with the observed benzene ring distortions determined by electron diffraction.     

A new force field (ECEPP-05) for peptides, proteins, and organic molecules

Parametrization and testing of a new all-atom force field for organic molecules and peptides with fixed bond lengths and bond angles are described. The van der Waals parameters for both the organic molecules and the peptides were taken from J. Phys. Chem. B 2003, 107, 7143 and J. Phys. Chem. B 2004, 108, 12181. First, the values of the 1-4 nonbonded and electrostatic scale factors appropriate to the new force field were determined by computing the conformational energies of six model molecules, namely, ethanol, ethylamine, propanol, propylamine, 1,2-ethanediol, and 1,3-propanediol with different values of these factors. The partial atomic charges of these molecules were obtained by fitting to the electrostatic potentials calculated with the HF/6-31G quantum-mechanical method. Two different charge models (single- and multiple-conformation-derived) were also considered. We demonstrated that the charge model has a stronger effect on the conformational energies than the 1-4 scaling. The choice of a charge model affected the conformational energies of even the smallest molecules considered, whereas the effect of the 1-4 electrostatic or nonbonded scaling was apparent only for 1,3-propanediol. The best agreement with high-level ab initio data was obtained with the multiple-conformation-derived charges and with no scaling of the 1-4 nonbonded or electrostatic interactions (scale factors of 1.0). Next, the torsional parameters of a large number of neutral and charged organic molecules, assumed to be models of the side chains of the 20 naturally occurring amino acids, were computed by fitting to rotational energy profiles obtained from ab initio MP2/6-31G calculations. The quality of the fits was high with average errors for torsional profiles of less than 0.2 kcal/mol. To derive the torsional parameters for the peptide backbone, the partial atomic charges of the 20 neutral and charged amino acids were obtained by fitting to the electrostatic potentials of terminally blocked amino acids using the HF/6-31G quantum-mechanical method. Then, the phi-psi energy maps of Ac-Ala-NMe and Ac-Gly-NMe were computed using MP2/6-31G//HF/6-31G quantum-mechanical methods. The phi-psi energy map of Ac-Ala-NMe was used for refinement of the nonbonded parameters for the backbone nitrogen and hydrogen bonded to it. Subsequently, the main-chain torsional parameters were obtained by fitting the molecular mechanics energies to the phi-psi energy maps of Ac-Ala-NMe and Ac-Gly-NMe. The transferability of the entire force field was demonstrated by reproducing the main energy minima of terminally blocked Ala3 from the literature. The performance of the force field was also evaluated by simulating crystal structures of small peptides. By comparison of simulated and experimental data, examination of the torsional-angle and atom-positional root-mean-square deviations of the energy-minimized crystal structures from the corresponding X-ray model structures demonstrated high accuracy of the force field.     

Importance of van der Waals Interactions in QM/MM Simulations

The importance of accurately treating van der Waals interactions between the quantum mechanical (QM) and molecular mechanical (MM) atoms in hybrid QM/MM simulations has been investigated systematically. First, a set of van der Waals (vdW) parameters was optimized for an approximate density functional method, the self-consistent charge-tight binding density functional (SCC-DFTB) approach, based on small hydrogen-bonding clusters. The sensitivity of condensed phase observables to the SCC-DFTB vdW parameters was then quantitatively investigated by SCC-DFTB/MM simulations of several model systems using the optimized set and two sets of extreme vdW parameters selected from the CHARMM22 forcefield. The model systems include a model FAD molecule in solution and a solvated enediolate, and the properties studied include the radial distribution functions of water molecules around the solute (model FAD and enediolate), the reduction potential of the model FAD and the potential of mean force for an intramolecular proton transfer in the enediolate. Although there are noticeable differences between parameter sets for gas-phase clusters and solvent structures around the solute, thermodynamic quantities in the condensed phase (e.g., reduction potential and potential of mean force) were found to be less sensitive to the numerical values of vdW parameters. The differences between SCC-DFTB/MM results with the three vdW parameter sets for SCC-DFTB atoms were explained in terms of the effects of the parameter set on solvation. The current study has made it clear that efforts in improving the reliability of QM/MM methods for energetical properties in the condensed phase should focus on components other than van der Waals interactions between QM and MM atoms.     

Interaction Forces between Hydrophobic and Hydrophilic Self-Assembled Monolayers

An investigation is presented of the interaction of charged self-assembled monolayers (SAMs) with a monoprotic ionizable acid functional group (-COOH) and uncharged SAMs with a methyl terminated functional group (-CH(3)). The strength of the interactions are determined using an atomic force microscope. For all electrolyte conditions investigated the interactions are not well described by a summation of van der Waals attractions and electrostatic repulsions in a manner suggesting that van der Waals attractions are screened. The repulsions are accurately described as corresponding to two surfaces of different charge interacting with surface charges that are independent of separation (i.e., the constant charge model). A small adhesion force was observed under all conditions and its magnitude increased with NaCl concentration. Copyright 2000 Academic Press.     

Amino acid analogues bind to carbon nanotube via π-π interactions: comparison of molecular mechanical and quantum mechanical calculations

Understanding the interaction between carbon nanotubes (CNTs) and biomolecules is essential to the CNT-based nanotechnology and biotechnology. Some recent experiments have suggested that the π-π stacking interactions between protein's aromatic residues and CNTs might play a key role in their binding, which raises interest in large scale modeling of protein-CNT complexes and associated π-π interactions at atomic detail. However, there is concern on the accuracy of classical fixed-charge molecular force fields due to their classical treatments and lack of polarizability. Here, we study the binding of three aromatic residue analogues (mimicking phenylalanine, tyrosine, and tryptophan) and benzene to a single-walled CNT, and compare the molecular mechanical (MM) calculations using three popular fixed-charge force fields (OPLSAA, AMBER, and CHARMM), with quantum mechanical (QM) calculations using the density-functional tight-binding method with the inclusion of dispersion correction (DFTB-D). Two typical configurations commonly found in π-π interactions are used, one with the aromatic rings parallel to the CNT surface (flat), and the other perpendicular (edge). Our calculations reveal that compared to the QM results the MM approaches can appropriately reproduce the strength of π-π interactions for both configurations, and more importantly, the energy difference between them, indicating that the various contributions to π-π interactions have been implicitly included in the van der Waals parameters of the standard MM force fields. Meanwhile, these MM models are less accurate in predicting the exact structural binding patterns (matching surface), meaning there are still rooms to be improved. In addition, we have provided a comprehensive and reliable QM picture for the π-π interactions of aromatic molecules with CNTs in gas phase, which might be used as a benchmark for future force field developments.     

Atom-atom partitioning of total (super)molecular energy: the hidden terms of classical force fields

Classical force fields describe the interaction between atoms that are bonded or nonbonded via simple potential energy expressions. Their parameters are often determined by fitting to ab initio energies and electrostatic potentials. A direct quantum chemical guide to constructing a force field would be the atom-atom partitioning of the energy of molecules and van der Waals complexes relevant to the force field. The authors used the theory of quantum chemical topology to partition the energy of five systems [H2, CO, H2O, (H2O)2, and (HF)2] in terms of kinetic, Coulomb, and exchange intra-atomic and interatomic contributions. The authors monitored the variation of these contributions with changing bond length or angle. Current force fields focus only on interatomic interaction energies and assume that these purely potential energy terms are the only ones that govern structure and dynamics in atomistic simulations. Here the authors highlight the importance of self-energy terms (kinetic and intra-atomic Coulomb and exchange).     

A modeling approach to describe the adhesion of rough, asymmetric particles to surfaces

A combined theoretical and experimental study of the adhesion of alumina particles and polystyrene latex spheres to silicon dioxide surfaces was performed. A boundary element technique was used to model electrostatic interactions between micron-scale particles and planar surfaces when the particles and surfaces were in contact. This method allows quantitative evaluation of the effects of particle geometry and surface roughness on the electrostatic interaction. The electrostatic interactions are combined with a previously developed model for van der Waals forces in particle adhesion. The combined model accounts for the effects of particle and substrate geometry, surface roughness and asperity deformation on the adhesion force. Predictions from the combined model are compared with experimental measurements made with an atomic force microscope. Measurements are made in aqueous solutions of varying ionic strength and solution pH. While van der Waals forces are generally dominant when particles are in contact with surfaces, results obtained here indicate that electrostatic interactions contribute to the overall adhesion force in certain cases. Specifically, alumina particles with complex geometries were found to adhere to surfaces due to both electrostatic and van der Waals interactions, while polystyrene latex spheres were not affected by electrostatic forces when in contact with various surfaces.     

Molecular mechanics (MM3) calculations on lithium amide compounds

The MM3 force field has been extended to deal with the lithium amide molecules that are widely used as efficient catalysts for stereoselective asymmetric synthesis. The MM3 force field parameters have been determined on the basis of the ab initio MP2/6-31G* and/or DFT (B3LYP/6-31G*, B3-PW91/6-31G*) geometry optimization calculations. To evaluate the electronic interactions specific to the lithium amides derived from the diamine molecules properly, the Lewis bonding potential term for the interaction between the lithium atom and the nonbonded adjacent electronegative atom such as nitrogen was introduced into the MM3 force field. The bond dipoles were evaluated correctly from the electronic charges on the atoms calculated by fitting to the electrostatic potential at points selected. The MM3 results on the molecular structures, conformational energies, and vibrational spectra show good agreement with those from the quantum mechanical calculations.     

Konformationen von 1,3-dioxanen nach der kraftfeldmethode. nichtbindende wechselwirkungen von sauerstoff

The treatment of oxygen lone pairs in force field calculations is discussed. The geometries of 1,3-dioxanes can be calculated well without explicit account of lone pairs. For the calculation of conformational energies in 5-alkyl-1,3-dioxanes electrostatic interactions between point charges taken from quantum mechanical calculations must be included. Addition of van der Waals interactions of lone pairs as weak potentials further improves the calculated conformational energies.     

A molecular mechanics force field for lignin

A CHARMM molecular mechanics force field for lignin is derived. Parameterization is based on reproducing quantum mechanical data of model compounds. Partial atomic charges are derived using the RESP electrostatic potential fitting method supplemented by the examination of methoxybenzene:water interactions. Dihedral parameters are optimized by fitting to critical rotational potentials and bonded parameters are obtained by optimizing vibrational frequencies and normal modes. Finally, the force field is validated by performing a molecular dynamics simulation of a crystal of a lignin fragment molecule and comparing simulation-derived structural features with experimental results. Together with the existing force field for polysaccharides, this lignin force field will enable full simulations of lignocellulose.     

Calculation of alcohol conformations by molecular mechanics

The geometries and energies of simple alcohols were calculated with a molecular mechanics force field. The force field requires the application of the charge interaction model with charges calculated by the CNDO/2 method, the importance of electrostatic interactions for the equilibrium of rotamers about the C-O bond exceeds that of van der Waals interactions. The calculated rotamer populations are discussed with regard to the value of ¹H NMR coupling constants ³JHCOH and other experimental data.     

Theoretical studies of the Wilcox molecular torsion balance. Is the edge-to-face aromatic interaction important?

[formula: see text] Molecular mechanics calculations with several different force fields and the GB/SA solvation model were carried out for the molecular torsion balance developed by Wilcox et al. to investigate the intramolecular nonbonded interactions between two aromatic rings. The preference is found to arise from a balance between intramolecular van der Waals interactions and solvation effects, with the latter favoring the unfolded conformation.     

Specific and nonspecific interaction forces between Escherichia coli and silicon nitride, determined by poisson statistical analysis

The nature of the physical interactions between Escherichia coli JM109 and a model surface (silicon nitride) was investigated in water via atomic force microscopy (AFM). AFM force measurements on bacteria can represent the combined effects of van der Waals and electrostatic forces, hydrogen bonding, steric interactions, and perhaps ligand-receptor type bonds. It can be difficult to decouple these forces into their individual components since both specific (chemical or short-range forces such as hydrogen bonding) and nonspecific (long-range colloidal) forces may be present in the overall profiles. An analysis is presented based on the application of Poisson statistics to AFM adhesion data, to decouple the specific and nonspecific interactions. Comparisons with classical DLVO theory and a modified form of a van der Waals expression for rough surfaces were made in order to help explain the nature of the interactions. The only specific forces in the system were due to hydrogen bonding, which from the Poisson analysis were found to be -0.125 nN. The nonspecific forces of 0.155 nN represent an overall repulsive interaction. These nonspecific forces are comparable to the forces calculated from DLVO theory, in which electrostatic-double layer interactions are added to van der Waals attractions calculated at the distance of closest approach, as long as the van der Waals model for "rough" spherical surfaces is used. Calculated electrostatic-double layer and van der Waals interactions summed to 0.116 nN. In contrast, if the classic (i.e., smooth) sphere-sphere model was used to predict the van der Waals forces, the sum of electrostatic and van der Waals forces was -7.11 nN, which appears to be a large overprediction. The Poisson statistical analysis of adhesion forces may be very useful in applications of bacterial adhesion, because it represents an easy way to determine the magnitude of hydrogen bonding in a given system and it allows the fundamental forces to be easily broken into their components.     

An Extension of the Consistent Valence Force Field (CVFF) with the Aim to Simulate the Structures of Vanadium Phosphorus Oxides and the Adsorption of n‐Butane and of 1‐Butene on their Crystal Planes

A specific force field of Consistent Valence Force Field type was developed with the aim to simulate the structures of catalysts of vanadium phosphorus oxide type and the reversible adsorption of organic compounds on specific crystallographic planes of such catalysts by molecular modeling. The appropriate parameters were derived for the bonded (stretching, bending, and torsional deformations) and nonbonded (attractive and repulsive van der Waals and Coulomb forces) atomic interactions for V—O and P—O bonds in typical fragments of these catalysts with the vanadium atom in the oxidation state IV. The parameters for bonded interactions were computed from Hessian matrices, supplied by the program DMol for performing Density Functional Theory, by means of a program for non‐linear regression. The DMol program was applied to energy minimize structures of known vanadium phosphorus oxides, which were compared with X‐ray structures, and to obtain their Hessian matrices as a basis for the force constants needed. Some hypothetical structural models had to be added. The van der Waals parameters were estimated by means of correlations between van der Waals radii and the repulsive parameters and between polarizabilities and the dispersive parameters from the literature. The force field obtained was applied to simulate the crystal structure of vanadyl pyrophosphate and to compute the heat of adsorption of n‐butane and of 1‐butene on its (100) plane (computer codes of company Biosym/MSI/Accelrys). The experimental crystal structure and the adsorption energies were fairly well reproduced, except that the a lattice constant proves somewhat too large.     

Molecular Mechanics Study of the Inclusion Complexes of Some 1,2,4-Oxadiazole Derivatives of 3,3′-Bis (1,2,4-Oxadiazol-5(4H)-one) with β-Cyclodextrin

Molecular mechanics calculations were employed to study the inclusion of some 1,2,4-oxadiazol derivatives in β-cyclodextrin in vacuum and in the presence of water as a solvent using MM + force field. The driving forces for complexation in both environments are dominated by nonbonded van der Waals host–guest interactions with little electrostatic contribution. Among 1,2,4-oxadiazole derivatives investigated in this work, 3,3′-bis(1,2,4-oxadiazol-5(4H)-one) (H₂OD) forms the least stable 1:1 complex and the stability increases as the chain length increases.