Substituent effects on the edge-to-face aromatic interactions

The edge-to-face interactions for either axially or facially substituted benzenes are investigated by using ab initio calculations. The predicted maximum energy difference between substituted and unsubstituted systems is approximately 0.7 kcal/mol (approximately 1.2 kcal/mol if substituents are on both axially and facially substituted positions). In the case of axially substituted aromatic systems, the electron density at the para position is an important stabilizing factor, and thus the stabilization/destabilization by substitution is highly correlated to the electrostatic energy. This results in its subsequent correlation with the polarization and charge transfer. Thus, the stabilization/destabilization by substitution is represented by the sum of electrostatic energy and induction energy. On the other hand, the facially substituted aromatic system depends on not only the electron-donating ability responsible for the electrostatic energy but also the dispersion interaction and exchange repulsion. Although the dispersion energy is the most dominating interaction in both axial and facial substitutions, it is almost canceled by the exchange repulsion in the axial substitution, whereas in the facial substitution, together with the exchange repulsion it augments the electrostatic energy. The systems with electron-accepting substituents (NO2, CN, Br, Cl, F) favor the axial substituent conformation, while those with electron-donating substituents (NH2, CH3, OH) favor the facial substituent conformation. The interactions for the T-shape complex systems of an aromatic ring with other counterpart such as H2, H2O, HCl, and HF are also studied.     

Pi-stacked interactions in explosive crystals: buffers against external mechanical stimuli

The pi-stacked interactions in some explosive crystal packing are discussed. Taking a typical pi-stacked explosive 2,4,6-trinitrobenzene-1,3,5-triamine (TATB) as a sample and using molecular simulations, we investigated the nature of the pi-stacked interactions versus the external mechanical stimuli causing possible slide and compression of explosives. As a result, between the neighbor layers in the TATB unit cell, the electrostatic attraction decreases with a little decrease of vdW attraction when its top layer slides, whereas the vdW attraction increases with a decrease of electrostatic attraction when TATB crystal is compressed along its c axis. Meanwhile, we studied the correlation between the pi-stacked structures and the impact sensitivities of explosives by means of three representatives including TATB with typical planar pi-stacked structures, 2,2-dinitroethylene-1,1-diamine (Fox-7) with wavelike pi-stacked structures, and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) without pi-stacked structure. The results showed that pi-stacked structures, particularly planar layers, can effectively buffer against external mechanical stimuli. That is, pi-stacked structures can partly convert the mechanical energy acting on them into their intermolecular interaction energy, to avoid the increase of the molecular vibration resulting in the explosive decomposition, the formation of hot spots, and the final detonation. This is another reason for the low mechanical sensitivity of pi-stacked explosives besides their stable conjugated molecular structures.     

Rationalizing Aggregate Structures with Orbital Contributions to the Exchange-Repulsion Energy**

It is shown that repulsive interactions have a crucial influence on the structure of prototypical non-covalently bonded systems. To explain this, we propose a molecular orbital-based model for the exchange-repulsion contribution to the total interaction energy. As a central result, our model shows that energetically preferred aggregate structures frequently exhibit reduced exchange repulsion, which can be deduced from the nodal structure of certain occupied orbitals. This is used to explain key features of the intermolecular potentials of the Cl₂-He, benzene-benzene, and benzene-hexafluorobenzene aggregates, which are not correctly reproduced by commonly applied electrostatic models.     

Calculation of electrostatic interaction energies in molecular dimers from atomic multipole moments obtained by different methods of electron density partitioning

Accurate and fast evaluation of electrostatic interactions in molecular systems is still one of the most challenging tasks in the rapidly advancing field of macromolecular chemistry, including molecular recognition, protein modeling and drug design. One of the most convenient and accurate approaches is based on a Buckingham-type approximation that uses the multipole moment expansion of molecular/atomic charge distributions. In the mid-1980s it was shown that the pseudoatom model commonly used in experimental X-ray charge density studies can be easily combined with the Buckingham-type approach for calculation of electrostatic interactions, plus atom-atom potentials for evaluation of the total interaction energies in molecular systems. While many such studies have been reported, little attention has been paid to the accuracy of evaluation of the purely electrostatic interactions as errors may be absorbed in the semiempirical atom-atom potentials that have to be used to account for exchange repulsion and dispersion forces. This study is aimed at the evaluation of the accuracy of the calculation of electrostatic interaction energies with the Buckingham approach. To eliminate experimental uncertainties, the atomic moments are based on theoretical single-molecule electron densities calculated at various levels of theory. The electrostatic interaction energies for a total of 11 dimers of alpha-glycine, N-acetylglycine and L-(+)-lactic acid structures calculated according to Buckingham with pseudoatom, stockholder and atoms-in-molecules moments are compared with those evaluated with the Morokuma-Ziegler energy decomposition scheme. For alpha-glycine a comparison with direct "pixel-by-pixel" integration method, recently developed Gavezzotti, is also made. It is found that the theoretical pseudoatom moments combined with the Buckingham model do predict the correct relative electrostatic interactions energies, although the absolute interaction energies are underestimated in some cases. The good agreement between electrostatic interaction energies computed with Morokuma-Ziegler partitioning, Gavezzotti's method, and the Buckingham approach with atoms-in-molecules moments demonstrates that reliable and accurate evaluation of electrostatic interactions in molecular systems of considerable complexity is now feasible.     

Scope and limitations of the SCS-MP2 method for stacking and hydrogen bonding interactions

Fluorobenzenes are pi-acceptor synthons that form pi-stacked structures in molecular crystals as well as in artificial DNAs. We investigate the competition between hydrogen bonding and pi-stacking in dimers consisting of the nucleobase mimic 2-pyridone (2PY) and all fluorobenzenes from 1-fluorobenzene to hexafluorobenzene (n-FB, with n = 1-6). We contrast the results of high level ab initio calculations with those obtained using ultraviolet (UV) and infrared (IR) laser spectroscopy of isolated and supersonically cooled dimers. The 2PY.n-FB complexes with n = 1-5 prefer double hydrogen bonding over pi-stacking, as diagnosed from the UV absorption and IR laser depletion spectra, which both show features characteristic of doubly H-bonded complexes. The 2-pyridone.hexafluorobenzene dimer is the only pi-stacked dimer, exhibiting a homogeneously broadened UV spectrum and no IR bands characteristic for H-bonded species. MP2 (second-order M?ller-Plesset perturbation theory) calculations overestimate the pi-stacked dimer binding energies by about 10 kJ/mol and disagree with the experimental observations. In contrast, the MP2 treatment of the H-bonded dimers appears to be quite accurate. Grimme's spin-component-scaled MP2 approach (SCS-MP2) is an improvement over MP2 for the pi-stacked dimers, reducing the binding energy by approximately 10 kJ/mol. When applied to explicitly correlated MP2 theory (SCS-MP2-R12 approach), agreement with the corresponding coupled-cluster binding energies [at the CCSD(T) level] is very good for the pi-stacked dimers, within +/- 1 kJ/mol for the 2PY complexes with 1-fluorobenzene, 1,2-difluorobenzene, 1,2,4,5-tetrafluorobenzene, pentafluorobenzene and hexafluorobenzene. Unfortunately, the SCS-MP2 approach also reduces the binding energy of the H-bonded species, leading to disagreement with both coupled-cluster theory and experiment. The SCS-MP2-R12 binding energies follow the SCS-MP2 binding energies closely, being about 0.5 and 0.7 kJ/mol larger for the H-bonded and pi-stacked forms, respectively, in an augmented correlation-consistent polarized valence quadruple-zeta basis. It seems that the SCS-MP2 and SCS-MP2-R12 methods cannot provide sufficient accuracy to replace the CCSD(T) method for intermolecular interactions where H-bonding and pi-stacking are competitive.     

Substituent effects in the benzene dimer are due to direct interactions of the substituents with the unsubstituted benzene

The prevailing views of substituent effects in the sandwich configuration of the benzene dimer are flawed. For example, in the polar/pi model of Cozzi and co-workers (J. Am. Chem. Soc. 1992, 114, 5729), electron-withdrawing substituents enhance binding in the benzene dimer by withdrawing electron density from the pi-cloud of the substituted ring, reducing the repulsive electrostatic interaction with the nonsubstituted benzene. Conversely, electron-donating substituents donate excess electrons into the pi-system and diminish the pi-stacking interaction. We present computed interaction energies for the sandwich configuration of the benzene dimer and 24 substituted dimers, as well as sandwich complexes of substituted benzenes with perfluorobenzene. While the computed interaction energies correlate well with sigmam values for the substituents, interaction energies for related model systems demonstrate that this trend is independent of the substituted ring. Instead, the observed trends are consistent with direct electrostatic and dispersive interactions of the substituents with the unsubstituted ring.     

Model chemistry calculations of thiophene dimer interactions: origin of pi-stacking

The intermolecular interaction energies of thiophene dimers have been calculated by using an aromatic intermolecular interaction (AIMI) model (a model chemistry for the evaluation of intermolecular interactions between aromatic molecules). The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium-size basis set. The calculated interaction energies of the parallel and perpendicular thiophene dimers are -1.71 and -3.12 kcal/mol, respectively. The substantial attractive interaction in the thiophene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge transfer but rather long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases the attraction significantly. The dispersion interaction is found to be the major source of attraction in the thiophene dimer. The calculated total interaction energy of the thiophene dimer is highly orientation dependent. Although electrostatic interaction is substantially weaker than dispersion interaction, it is highly orientation dependent, and therefore electrostatic interaction play an important role in the orientation dependence of the total interaction energy. The large attractive interaction in the perpendicular dimer is the cause of the preference for the herringbone structure in the crystals of nonsubstituted oligothiophenes (alpha-terthienyls), and the steric repulsion between the beta-substituents is the cause of the pi-stacked structure in the crystals of some beta-substituted oligothiophenes.     

Theoretical study on mesogenic core structures of nematic liquid crystalline compounds

The structures and intermolecular interaction energies of 10 dimers, included in the mesogenic core structures of typical liquid crystalline (LC) compounds, are obtained at the MP2/6-31G(d) level of theory. It is proved that the dispersion energy significantly contributes to the total interaction energy of these dimers. Even when bulky substituents are introduced into the core part, the interaction energy is still large. It is also revealed that when a long intermolecular distance is provided by a high steric repulsion originating from the linkage of two phenyl groups, the dispersion energy is significantly small. However, in this range of intermolecular distances, the electrostatic energy caused by a strong quadrupole-quadrupole attractive interaction plays a dominant role, and as a result, a rather stable dimer is formed. In all 10 dimers, the dispersion, electrostatic, and exchange-repulsion energies strongly depend on the geometrical orientation of the molecules. The calculated interaction energies of these dimers are also compared with the corresponding experimentally measured viscosities. The results suggest an explicit linear relationship between the interaction energies and viscosities.     

Electrostatic interactions in peptides. Polarisation effects due to an α-helix

A quantum mechanical/molecular mechanical study of a dodecapeptide made of 11 alanine and one asparagine residues in a helical conformation is carried out by means of the local self-consistent field/molecular mechanical and integrated molecular orbital and molecular mechanics computational schemes. The electronic properties of the asparagine side chain are analysed to extract the influence of electrostatic and induction interaction. One finds that induction may play an important role in the energetic and structural features of the systems modelled with mixed methods. The importance of performing quantum computations which explicitly take account of the electrostatic interactions is pointed out.Contribution to the Jacopo Tomasi Honorary Issue     

取代基对咔咯锰(V)-氧配合物Mn―O的成键影响

用密度泛函理论(DFT)方法对一系列A₃型咔咯锰(V)氧配合物进行了理论计算. 结果表明: 咔咯锰(V)-氧配合物中Mn―O键是由1 个σ键和2 个π键构成的叁键结构; 当咔咯中位取代基由推电子过渡到拉电子性质时, 咔咯骨架紧缩, Mn―O键缩短, 其拉曼光谱的伸缩振动峰往高波数移动; 取代基与氧原子的静电作用模式由正-负吸引转化为负-负排斥, 导致Mn―O键解离能ΔE减少, 即拉电子基团有利于增强咔咯锰(V)-氧配合物氧原子的活泼性.     

Tryptophan side chain electrostatic interactions determine edge-to-face vs parallel-displaced tryptophan side chain geometries in the designed beta-hairpin "trpzip2"

The interaction geometries of the four tryptophan (Trp) side chains in the 12-residue designed beta-hairpin trpzip2 are investigated using all-atom explicit-solvent molecular dynamics simulations. The experimentally observed edge-to-face (EtF) pairwise interaction geometries are stable on a time scale of 10 ns. However, removing the electrostatic multipoles of the Trp side chains while retaining the dipoles of the side chains' NH moieties induces a conformational change to a geometry in which three of the four side chains interact in a parallel-displaced (PD) manner. Free energy simulations of the Etf to PD conformational change reveal that, with the side chain multipole moments intact (+MP), the EtF conformation is preferred by 5.79 kcal/mol. Conversely, with only the dipole moments of the side chain NH moieties intact (-MP), the PD conformation's free energy is more favorable by 1.71 kcal/mol. In contrast to energetic similarities for Trp side chain-water electrostatic and Trp side chain-Trp side chain and Trp side chain-water van der Waals, +MP Trp side chain-Trp side chain electrostatic interactions are more favorable by 4.21 kcal/mol in the EtF conformation, while in the -MP case the EtF and PD conformations' Trp side chain-Trp side chain electrostatic energies are nearly identical. The results highlight the importance of electrostatic multipole moments in determining aromatic-aromatic interaction geometries in aqueous biomolecular systems and argue for the inclusion of this physics in simplified models used for protein-ligand docking and protein structure prediction, possibly through a truncated Coulomb term between aromatic moieties.     

Substitution effect on the geometry and electronic structure of the ferrocene

The substitution effects on the geometry and the electronic structure of the ferrocene are systematically and comparatively studied using the density functional theory. It is found that -NH(2) and -OH substituents exert different influence on the geometry from -CH(3), -SiH(3), -PH(2), and -SH substituents. The topological analysis shows that all the C-C bonds in a-g are typical opened-shell interactions while the Fe-C bonds are typical closed-shell interactions. NBO analysis indicates that the cooperated interaction of d --> pi* and feedback pi --> d + 4s enhances the Fe-ligand interaction. The energy partitioning analysis demonstrates that the substituents with the second row elements lead to stronger iron-ligand interactions than those with the third row elements. The molecular electrostatic potential predicts that the electrophiles are expected to attack preferably the N, O, P, or S atoms in Fer-NH(2), Fer-OH, Fer-PH(2), and Fer-SH, and attack the ring C atoms in Fer-SiH(3) and Fer-CH(3). In turn, the nucleophiles are supposed to interact predominantly by attacking the hydrogen atoms. The simulated theoretical excitation spectra show that the maximum absorption peaks are red-shifted when the substituents going from second row elements to the third row elements.     

Conformational sensitivity of polyether macrocycles to electrostatic potential: Partial atomic charges,molecular mechanics,and conformational prediction

Molecular recognition (whether by enzymes, the immune system, or chelating ligands) depends critically on molecular conformation. Molecular mechanics predicts energetically favorable molecular conformations by locating low energy conformations using an empirical fit of molecular potential energy as a function of internal coordinates. Molecular mechanics analysis of 18-crown-6 demonstrates that the nonbonded term (primarily the electrostatic part) is the largest contributor to the conformational energy. Nevertheless, common methods of treating the electrostatic interaction for 18-crown-6 yield inconsistent values for conformational energies partly because partial charges assigned to each atom can change with conformation due to through-space inductive effects which are not considered in most molecular mechanics programs. Similar findings from several other groups are reviewed to support our conclusions. We argue for care and caution in predicting conformational preferences of molecules with two or more highly polar atoms. We also discuss the desirability of using an empirical method of partial charge determination such as the charge equilibration algorithm of Rappé and Goddard (or a suitable generalization which includes polarization) as a method of including these effects in molecular mechanics and molecular dynamics calculations.     

壳聚糖在稀溶液中的分子构象及其影响因素

结合国内外对壳聚糖在溶液性质方面的研究进展,从壳聚糖的分子量、乙酰度和溶液的离子强度三个方面,探讨了它们对壳聚糖在稀溶液中的分子构象所产生的影响,并从分子间和分子内相互作用力(主要是范德华力、氢键、疏水作用和静电作用)的角度对这种构象转变进行了理论解释。     

Synthesis of electrochemically responsive TTF-based molecular tweezers: evidence of tight intramolecular TTF pairing in solution

We report the synthesis and conformational studies of TTF-containing molecular tweezers based on a 1,2,4,5-tetramethylbenzene scaffold. In the neutral form the tweezers are expected to adopt the closed conformation, while, upon oxidation, the open conformation should be preferred due to electrostatic repulsion between the oxidized TTF moieties. Cyclic voltammetry studies demonstrate electronic pairing with formation of mixed-valence [TTF]₂⁺ species and opening of the tweezers upon the full oxidation of the TTF groups. Variable-temperature (VT) NMR studies evidence tight intramolecular TTF pairing at low temperature. Molecular modeling studies showed clear preference for an open conformation of tweezers in a fully oxidized state.     

DLVO AND NON-DLVO SURFACE FORCES AND INTERACTIONS IN COLLOIDAL DISPERSIONS

Some aspects of DLVO and non-DLVO forces in colloidal systems are over-viewed. The influence of long range interactions on some kinetic properties of dispersions, as Brownian diffusion, is discussed. It is shown, both theoretically and experimentally, that the electrostatic repulsion increases the collective diffusivity. The film stratification and oscillatory structure forces in colloidal suspensions are considered within the framework of an uniform approach The presence of small colloidal species (e. g. micelles or polymer molecules) may lead to several maxima and minima in the disjoining pressure isotherm. The particular case of interacting emulsion droplets is examined accounting for the interfacial deformability. The droplet deformation acts as a soft repulsion but affects also the remaining contributions to the interaction energy due to changes of the droplet shape. A general procedure for calculating the inter-droplet interaction energy, as well as the equilibrium film radius and thickness in a doublet of droplets, is suggested. The energy of interaction between charged colloidal particles, due to correlations of the density fluctuations in the electric double layer is also studied. It is found that this effect may lead to attraction greater than the van der Waals contribution, especially when multivale counter ions are present.     

Structural and conformational analysis of two native procyanidin trimers

The structure and conformation of two native procyanidin trimers in water have been determined using 2D NMR and molecular mechanics. The results show the existence of four rotameric forms, one of which is predominant (60 to 80%). These four rotamers are shown to be in slow to intermediate exchange on the NMR timescale. Both trimers, whose structures vary owing to a different substitution of one carbon atom, adopt conformations in which stacking between different phenolic rings is favored.     

Structural, dynamic, and electrostatic properties of fully hydrated DMPC bilayers from molecular dynamics simulations accelerated with graphical processing units (GPUs)

We present results of molecular dynamics simulations of fully hydrated DMPC bilayers performed on graphics processing units (GPUs) using current state-of-the-art non-polarizable force fields and a local GPU-enabled molecular dynamics code named FEN ZI. We treat the conditionally convergent electrostatic interaction energy exactly using the particle mesh Ewald method (PME) for solution of Poisson's Equation for the electrostatic potential under periodic boundary conditions. We discuss elements of our implementation of the PME algorithm on GPUs as well as pertinent performance issues. We proceed to show results of simulations of extended lipid bilayer systems using our program, FEN ZI. We performed simulations of DMPC bilayer systems consisting of 17,004, 68,484, and 273,936 atoms in explicit solvent. We present bilayer structural properties (atomic number densities, electron density profiles), deuterium order parameters (S(CD)), electrostatic properties (dipole potential, water dipole moments), and orientational properties of water. Predicted properties demonstrate excellent agreement with experiment and previous all-atom molecular dynamics simulations. We observe no statistically significant differences in calculated structural or electrostatic properties for different system sizes, suggesting the small bilayer simulations (less than 100 lipid molecules) provide equivalent representation of structural and electrostatic properties associated with significantly larger systems (over 1000 lipid molecules). We stress that the three system size representations will have differences in other properties such as surface capillary wave dynamics or surface tension related effects that are not probed in the current study. The latter properties are inherently dependent on system size. This contribution suggests the suitability of applying emerging GPU technologies to studies of an important class of biological environments, that of lipid bilayers and their associated integral membrane proteins. We envision that this technology will push the boundaries of fully atomic-resolution modeling of these biological systems, thus enabling unprecedented exploration of meso-scale phenomena (mechanisms, kinetics, energetics) with atomic detail at commodity hardware prices.     

On the paucity of molecular actinide complexes with unsupported metal-metal bonds: a comparative investigation of the electronic structure and metal-metal bonding in U2X6 (X = Cl, F, OH, NH2, CH3) complexes and d-block analogues

Density functional calculations have been performed on M2X6 complexes (where M = U, W, and Mo and X = Cl, F, OH, NH2, and CH3) to investigate general aspects of their electronic structures and explore the similarities and differences in metal-metal bonding between f-block and d-block elements. A detailed analysis of the metal-metal interactions has been conducted using molecular orbital theory and energy decomposition methods. Multiple (sigma and pi) bonding is predicted for all species investigated, with predominant f-f and d-d metal orbital character, respectively, for U and W or Mo complexes. The energy decomposition analysis involves contributions from orbital interactions (mixing of occupied and unoccupied orbitals), electrostatic effects (Coulombic attraction and repulsion), and Pauli repulsion (associated with four-electron two-orbital interactions). The general results suggest that the overall metal-metal interaction is stronger in the Mo and W species, relative to the U analogues, as a consequence of a significantly less destabilizing contribution from the combined Pauli and electrostatic ("pre-relaxation") effects. Although the orbital-mixing ("post-relaxation") contribution to the total bonding energy is predicted to have a larger magnitude in the U complexes, this is not sufficiently strong to compensate for the comparatively greater destabilization that originates from the Pauli-plus-electrostatic effects. Of the pre-relaxation terms, the Pauli repulsion is comparable in analogous U and d-block compounds, contrary to the electrostatic term, which is (much) less favorable in the U systems than in the W and Mo systems. This generally weak electrostatic stabilization accounts for the large pre-relaxation destabilization in the U complexes and, ultimately, for the relative weakness of the U-U bonds. The origin of the small electrostatic term in the U compounds is traced primarily to MX(3) fragment overlap effects.     

C5H4-BR2 bending in ferrocenylboranes: a delocalized through-space interaction between iron and boron

A comparison of the molecular structures of mono-, di- and tetraborylated ferrocenes [Fc{B(R(1))(R(2))}] (R(1)/R(2)=Br/Br, Br/Fc, Br/Me, Me/Me, Me/OH, OMe/OMe), 1,1'-[fc{B(R(1))(R(2))}(2)] (R(1)/R(2)=Br/Br, Br/Me, OMe/OMe), and 1,1',3,3'-[Fe{C(5)H(3)(BMe(2))(2)}(2)] revealed the boryl substituent(s) to be bent out of the Cp ring plane towards the iron center. The corresponding dip angle alpha* decreases with decreasing Lewis acidity of the boron atom and with increasing degree of borylation at the ferrocene core. This trend is well reproduced by DFT calculations (including [FcBH(2)], not yet accessible experimentally). A Bader analysis of the electron density topology of [FcBH(2)] (alpha*=26.5 degrees ; BP86/TZVP) clearly showed that there is no direct iron-boron bonding in this compound. Instead, strongly delocalized orbital interactions have been identified that involve the boron p orbital, C(ipso) of the adjacent Cp ring, d orbitals at iron, and a through-space interaction with the second Cp ring. A second important factor is attractive electrostatic interactions, which are enhanced upon ligand bending. Cyclic voltammetric measurements on the series [FcBMe(2)], 1,1'-[fc(BMe(2))(2)], and 1,1',3,3'-[Fe{C(5)H(3)(BMe(2))(2)}(2)] indicate a substantial anodic shift in the oxidation potential of the central iron atom upon introduction of BMe(2) substituents. Addition of 4-dimethylaminopyridine (DMAP) does not just counterbalance this effect, but leads to a cathodic shift of the Fe(II)/Fe(III) redox transition far beyond the half-wave potential of parent ferrocene. In the Mossbauer spectra, a continuous decrease in the quadrupole splitting (QS) is observed upon going from parent ferrocene to [FcBMe(2)], to 1,1'-[fc(BMe(2))(2)], and to 1,1',3,3'-[Fe{C(5)H(3)(BMe(2))(2)}(2)]. In contrast, no significant differences are found between the QS values of ferrocene, [Fc(BMe(2)-DMAP)], and 1,1'-[fc(BMe(2)-DMAP)(2)].