Optimization of a molecular mechanics force field for type-II polyoxometalates focussing on electrostatic interactions: a case study

In this study, we have focussed on type-II polyanions such as [M(7)O(24)](6-), and we have developed and validated optimized force fields that include electrostatic and van der Waals interactions. These contributions to the total steric energy are described by the nonbonded term, which encompasses all interactions between atoms that are not transmitted through the bonds. A first validation of a stochastic technique based on genetic algorithms was previously made for the optimization of force fields dedicated to type-I polyoxometalates. To describe the new nonbonded term added in the functional, a fixed-charged model was chosen. Therefore, one of the main issues was to analyze that which partial atomic charges could be reliably used to describe these interactions in such inorganic compounds. Based on several computational strategies, molecular mechanics (MM) force field parameters were optimized using different types of atomic charges. Moreover, the influence of the electrostatic and van der Waals buffering constants and 1,4-interactions scaling factors used in the force field were also tested, either being optimized as well or fixed with respect to the values of CHARMM force field. Results show that some atomic charges are not well adapted to CHARMM parameters and lead to unrealistic MM-optimized structures or a MM divergence. As a result, a new scaling factor has been optimized for Quantum Theory of Atoms in Molecules charges and charges derived from the electrostatic potential such as ChelpG. The force fields optimized can be mixed with the CHARMM force field, without changing it, to study for the first time hepta-anions interacting with organic molecules.     

Determination of van der Waals forces in monolayer films of lipids & biopolymers. Equation of state for two-dimensional films

Amphiphile molecules are characterized by the dual property arising from the interactions between the apolar [alkyl] and the polar part and the surrounding solvent, i.e., water. In assemblies which amphiphiles form in diverse systems, e.g., micelles, soap bubbles, monolayers or bilayers at interfaces, the attractive forces are attributed to the van der Waals forces. It is not easy to estimate the magnitude of van der Waals forces in some of these systems by any direct method.The magnitude of van der Waals forces in spread monolayers of lipids and biopolymers has been reported to be estimated from experimental data. The magnitude of these forces has been estimated by using an equation of state of a very general form, as delineated herein. In the current literature no such attempt has been reported in the analyses of these monolayers spread on aqueous surfaces. These analyses suggest that the predominant surface forces arise from van der Waals interactions, if the magnitude of electrostatic charge repulsions is weak. The equation-of-state as derived indicates that it is useful in providing information about the molecular interaction in monolayers, for both lipids and biopolymers.     

溶致液晶中金属纳米粒子的掺杂及其作用机制研究

在用琥珀酸二异辛酯磺酸钠(AOT)构建的具有长程有序结构的层状溶致液晶内, 用不同方式导入预制的亲油或亲水贵金属纳米粒子, 可得到纳米粒子分布在不同介观空间内的无机/有机杂合体. 依据小角X射线散射和偏光显微镜结果, 通过分析掺杂纳米粒子与液晶模板的相互作用, 对掺杂前后体系结构的变化及制得杂合体的稳定性进行了表征. 结果表明, 除考虑掺杂粒子与层状模板空间的匹配外, 体系中静电斥力、范德华引力和Helfrich涨落力之间的平衡是维持液晶结构稳定的基本条件.     

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.     

Crystallization Force—A Density Functional Theory Concept for Revealing Intermolecular Interactions and Molecular Packing in Organic Crystals

Organic molecules are prone to polymorphic formation in the solid state due to the rich diversity of functional groups that results in comparable intermolecular interactions, which can be greatly affected by the selection of solvent and other crystallization conditions. Intermolecular interactions are typically weak forces, such as van der Waals and stronger short‐range ones including hydrogen bonding, that are believed to determine the packing of organic molecules during the crystal‐growth process. A different packing of the same molecules leads to the formation of a new crystal structure. To disclose the underlying causes that drive the molecule to have various packing motifs in the solid state, an electronic concept or function within the framework of conceptual density functional theory has been developed, namely, crystallization force. The concept aims to describe the local change in electronic structure as a result of the self‐assembly process of crystallization and may likely quantify the locality of intermolecular interactions that directs the molecular packing in a crystal. To assess the applicability of the concept, 5‐methyl‐2‐[(2‐nitrophenyl)amino]‐3‐thiophenecarbonitrile, so‐called ROY, which is known to have the largest number of solved polymorphs, has been examined. Electronic calculations were conducted on the seven available crystal structures as well as on the single molecule. The electronic structures were analyzed and crystallization force values were obtained. The results indicate that the crystallization forces are able to reveal intermolecular interactions in the crystals, in particular, the close contacts that are formed between molecules. Strong correlations exist between the total crystallization force and lattice energy of a crystal structure, further suggesting the underlying connection between the crystallization force and molecular packing.     

Trajectory Analysis and Collision Efficiency during Microbubble Flotation

The hydrodynamic interaction between a rising bubble and a sedimenting particle during microbubble flotation is considered. The effects of attractive van der Waals forces and attractive or repulsive electrostatic forces are included. A mathematical model is presented which is used to perform a trajectory analysis and to calculate collision efficiencies between the bubble and particle. It is shown that collision efficiencies and the nature of the bubble-particle interactions are strongly dependent on the relative strengths of the van der Waals and electrostatic forces and on the lengthscales over which these forces act. It is demonstrated that optimal operating conditions can be suggested to achieve efficient microbubble flotation by correctly accounting for the interaction of van der Waals, electrostatic, and hydrodynamic forces. Copyright 1999 Academic Press.     

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.     

Fluoroaromatic-fluoroaromatic interactions between inhibitors bound in the crystal lattice of human carbonic anhydrase II

Intermolecular interactions of eleven different fluoroaromatic inhibitors are probed within the scaffolding of the crystal lattice of Phe-131-->Val carbonic anhydrase II. The degree and pattern of fluorine substitution on the inhibitor benzyl ring modulate its size, shape, and electronic character. In turn, these properties affect the geometry of intermolecular interactions between the fluoroaromatic rings of two different inhibitor molecules bound in the crystal lattice, as determined by X-ray crystallography. Depending on the degree and pattern of fluorine substitution, we observe a face-to-face (aromatic-aromatic) interaction, an atom-to-face (carbonyl-aromatic) interaction, or no interaction at all. These interaction geometries are analyzed with regard to van der Waals, electrostatic, and possible charge-transfer effects. For the aromatic-aromatic interactions investigated in this study, with aromatic ring quadrupoles specifically "tuned" by the degree and pattern of fluorination, the structural results suggest that London forces and charge-transfer complexation dominate over weakly polar electrostatic interactions in the association of aromatic ring pairs.     

Modeling intermolecular interactions of physisorbed organic molecules using pair potential calculations

The understanding and control of epitaxial growth of organic thin films is of crucial importance in order to optimize the performance of future electronic devices. In particular, the start of the submonolayer growth plays an important role since it often determines the structure of the first layer and subsequently of the entire molecular film. We have investigated the structure formation of 3,4,9,10-perylene-tetracarboxylic dianhydride and copper-phthalocyanine molecules on Au(111) using pair-potential calculations based on van der Waals and electrostatic intermolecular interactions. The results are compared with the fundamental lateral structures known from experiment and an excellent agreement was found for these weakly interacting systems. Furthermore, the calculations are even suitable for chemisorptive adsorption as demonstrated for copper-phthalocyanine/Cu(111), if the influence of charge transfer between substrate and molecules is known and the corresponding charge redistribution in the molecules can be estimated. The calculations are of general applicability for molecular adsorbate systems which are dominated by electrostatic and van der Waals interaction.     

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.     

Influence of chlorine substitution on the self-assembly of zinc phthalocyanine

The adsorption and ordering of zinc phthalocyanine (ZnPc) and octachloro zinc phthalocyanine (ZnPcCl(8)) on an Ag(111) surface is studied in situ by scanning tunneling microscopy under ultrahigh vacuum. Two-dimensional self-assembled supramolecular domains are observed for these two molecules. We show how substituting chlorine atoms for half of the peripheral hydrogen atoms on ZnPc influences the self-assembly mechanisms. While intermolecular interactions are dominated by van der Waals forces in ZnPc molecular networks, ZnPcCl(8) molecular packing undergoes a sequential phase evolution driven by the creation of C-Cl...H-C hydrogen bonds between adjacent molecules. At the end of this evolution, the final molecular assembly involves all possible hydrogen bonds. Our study also reveals the influence of molecule-substrate interactions through the presence of fault lines generating a stripe structure in the molecular film.     

Assembly,Thermodynamics, and Structure of a Two‐Wheeled Composite of a Dumbbell‐Shaped Molecule and Cylindrical Molecules with Different Edges

A carbonaceous dumbbell was able to spontaneously glue two tubular receptors to form a unique two‐wheeled composite through van der Waals interactions, thus forcing the wheel components into contact with each other at the edges. In the present study, two tubular receptors with enantiomeric carbon networks were assembled on the dumbbell joint, and the handedness of the receptors was discriminated, thus leading to the self‐sorting of homomeric receptors from a mixture of enantiomeric tubes. The crystal structures of the composites revealed the structural origins of the molecular recognition driven by van der Waals forces as well as the presence of a columnar array of C₁₂₀ molecules in a 1:1 composite.     

Assembly,Thermodynamics, and Structure of a Two‐Wheeled Composite of a Dumbbell‐Shaped Molecule and Cylindrical Molecules with Different Edges

A carbonaceous dumbbell was able to spontaneously glue two tubular receptors to form a unique two‐wheeled composite through van der Waals interactions, thus forcing the wheel components into contact with each other at the edges. In the present study, two tubular receptors with enantiomeric carbon networks were assembled on the dumbbell joint, and the handedness of the receptors was discriminated, thus leading to the self‐sorting of homomeric receptors from a mixture of enantiomeric tubes. The crystal structures of the composites revealed the structural origins of the molecular recognition driven by van der Waals forces as well as the presence of a columnar array of C₁₂₀ molecules in a 1:1 composite.     

Uncovering the interactions driving carotenoid binding in light-harvesting complexes

Carotenoids are essential constituents of plant light-harvesting complexes (LHCs), being involved in protein stability, light harvesting, and photoprotection. Unlike chlorophylls, whose binding to LHCs is known to require coordination of the central magnesium, carotenoid binding relies on weaker intermolecular interactions (such as hydrogen bonds and van der Waals forces), whose character is far more elusive. Here we addressed the key interactions responsible for carotenoid binding to LHCs by combining molecular dynamics simulations and polarizable quantum mechanics/molecular mechanics calculations on the major LHC, LHCII. We found that carotenoid binding is mainly stabilized by van der Waals interactions with the surrounding chlorophyll macrocycles rather than by hydrogen bonds to the protein, the latter being more labile than predicted from structural data. Furthermore, the interaction network in the binding pockets is relatively insensitive to the chemical structure of the embedded carotenoid. Our results are consistent with a number of experimental data and challenge the role played by specific interactions in the assembly of pigment-protein complexes.

Carotenoids are essential constituents of plant light-harvesting complexes. This in silico study shows that carotenoid binding is mainly driven by van der Waals interactions with the surrounding chlorophylls rather than hydrogen bonds to the protein.     

Structural changes in the molecular sheets along (hk0) planes derived from cellulose Iβ by molecular dynamics simulations

We studied the stability of molecular sheets with four cellotetraoses in an aqueous environment by molecular dynamics simulation to identify the molecular details of first structure as one of the possibilities in the course of crystallization of cellulose I. After simulation, the molecular sheets formed by van der Waals forces along the (11?0) and (110) crystal plane did not change their structures in an aqueous environment, whereas the other ones formed by hydrogen bonds along the (100) and (200) crystal plane changed into a van der Waals associated molecular sheet, similar to the former. These simulated molecular sheets formed by van der Waals forces were structurally stable in water because of their hydrophilic exterior and hydrophobic interior. Therefore, if the molecular sheet structures are formed in the real system, the sheets formed by van der Waals forces are probably the initial structure of crystallization. A close analysis indicated that these sheets could be classified into two groups in terms of the hydrogen bonding networks, camber angle, and main and side chain conformations. One group was the molecular sheets corresponding to the (110) after simulation. This sheet is probably rigid because intramolecular hydrogen bonds of the chains in the sheet are highly developed. The other group was the molecular sheets corresponding to (200), (100), and (11?0) crystal plane: the chains in these sheets seemed to be rather flexible due to their moderately developed intramolecular hydrogen bonds.     

Structural properties of liquidN,N-dimethylformamide

Monte-Carlo simulation of liquidN,N-dimethylformamide was performed. The influence or electrostatic and van der Waals interactions on regularities of the mutual molecular arrangement was investigated. The spatial structure of liquid DMF is determined by the type of molecular packing and steric factors and is close to a random closely packed system of soft spheres. The electrostatic interactions only affect the mutual orientation of the molecular Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 21–27, January, 1998.     

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.     

QM/MM-PBSA method to estimate free energies for reactions in proteins

We have developed a method to estimate free energies of reactions in proteins, called QM/MM-PBSA. It estimates the internal energy of the reactive site by quantum mechanical (QM) calculations, whereas bonded, electrostatic, and van der Waals interactions with the surrounding protein are calculated at the molecular mechanics (MM) level. The electrostatic part of the solvation energy of the reactant and the product is estimated by solving the Poisson-Boltzmann (PB) equation, and the nonpolar part of the solvation energy is estimated from the change in solvent-accessible surface area (SA). Finally, the change in entropy is estimated from the vibrational frequencies. We test this method for five proton-transfer reactions in the active sites of [Ni,Fe] hydrogenase and copper nitrite reductase. We show that QM/MM-PBSA reproduces the results of a strict QM/MM free-energy perturbation method with a mean absolute deviation (MAD) of 8-10 kJ/mol if snapshots from molecular dynamics simulations are used and 4-14 kJ/mol if a single QM/MM structure is used. This is appreciably better than the original QM/MM results or if the QM energies are supplemented with a point-charge model, a self-consistent reaction field, or a PB model of the protein and the solvent, which give MADs of 22-36 kJ/mol for the same test set.     

EGFR和4-苯胺喹唑啉类抑制剂之间相互作用模式的研究

采用分子动力学和MM/PBSA相结合的方法预测了表皮生长因子受体和4-苯胺喹 啉类抑制剂的相互作用模式。在分子动力学采样的基础上,采用MM/PBSA的方法分 别预测了四种可能结合模式下表皮生长因子受体和4-苯胺喹唑啉类抑制剂间的结合 自由能。在MM/PBSA计算中,受体和抑制剂之间的非键相互作用能采用分子力学 (MM)的方法得到;溶剂效应中极性部分对自由能的贡献通过解Possion- Boltzmanne (PB)方程的方法得到;溶液效应中非极性部分对自由能的贡献则通过 分子表面积计算(SA)的方法得到。计算表明,在四种结合模式下,表皮生长因子受 体和4-苯胺喹唑啉类抑制剂之间的结合自由能有较大的差别。在最佳的相互作用模 式中,抑制剂的苯胺部分位于活性口袋的底部,能够与受体残基的非极性侧链产生 很强的范德华和疏水相互作用。抑制剂喹唑啉环上的N(1)原子能够和Met-769上的 NH形成稳定的氢键,而抑制剂上的N(3)原子则和周围的一个水分子形成氢键。同时 ,抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏 水相互作用。最佳结合模式能够很好地解释已有抑制剂结构和活性间的关系。     

VB/MM--the validity of the underlying approximations

Two recently developed methods, VB/MM and density embedded VB/MM (DE-VB/MM), are compared, and their respective approximations are examined. The two methods combine valence-bond (VB) calculations with molecular mechanics (MM) and aim to allow VB analysis of reactions in large biological environments. Furthermore, the two methods utilize two major approximations regarding both the overlap and the reduced resonance integral between the various VB configurations. The difference between the two methods, however, is that VB/MM employs these approximations for the overall interaction of the reacting fragments with their surrounding, whereas DE-VB/MM employs the approximations only with regards to the van der Waals (VdW) interactions whereas the electrostatic interactions are calculated rigorously at the quantum level. The approximations that lay the grounds for the two methods involve the assumption that the overlap between the VB configurations and the respective reduced resonance integral are both invariant to the environment. Similar approximations are utilized in several other VB-based QM/MM methods. However, although extensively used, these approximations were never rigorously proved. Here, we exploit the development of the DE-VB/MM method to numerically examine the approximations by calculating the accurate as well as the approximated values of overlap and reduced resonance integral for systems where the environment involves only electrostatic interactions. The quality of the approximations is examined together with their effect on the absolute energies, the wave function, and the overall energetics. Three test cases are chosen, the dissociation of CH 3F and LiF and the identity S N2 reaction. It is shown that the approximations are usually good with the exception of cases where extreme changes are expected in the wave function. Furthermore, the impact of the approximations on the overall wave function and the overall energetics is found to be quite small. It is concluded that VB/MM, where the approximations are used more extensively, can serve as the first method of choice.