Hydrogen bonding, halogen bonding and lithium bonding: an atoms in molecules and natural bond orbital perspective towards conservation of total bond order, inter- and intra-molecular bonding

One hundred complexes have been investigated exhibiting D-X center dot center dot center dot A interactions, where X = H, Cl or Li and DX is the `X bond' donor and A is the acceptor. The optimized structures of all these complexes have been used to propose a generalized `Legon-Millen rule' for the angular geometry in all these interactions. A detailed Atoms in Molecules (AIM) theoretical analysis confirms an important conclusion, known in the literature: there is a strong correlation between the electron density at the X center dot center dot center dot A bond critical point (BCP) and the interaction energy for all these interactions. In addition, we show that extrapolation of the fitted line leads to the ionic bond for Li-bonding (electrostatic) while for hydrogen and chlorine bonding, it leads to the covalent bond. Further, we observe a strong correlation between the change in electron density at the D-X BCP and that at the X center dot center dot center dot A BCP, suggesting conservation of the bond order. The correlation found between penetration and electron density at BCP can be very useful for crystal structure analysis, which relies on arbitrary van der Waals radii for estimating penetration. Various criteria proposed for shared-and closed-shell interactions based on electron density topology have been tested for H/Cl/Li bonded complexes. Finally, using the natural bond orbital (NBO) analysis it is shown that the D-X bond weakens upon X bond formation, whether it is ionic (DLi) or covalent (DH/DCl) and the respective indices such as ionicity or covalent bond order decrease. Clearly, one can think of conservation of bond order that includes ionic and covalent contributions to both D-X and X center dot center dot center dot A bonds, for not only X = H/Cl/Li investigated here but also any atom involved in intermolecular bonding.

Characterization of Interactions Involving Bromine in 2,2-Dibromo2,3-dihydroinden-1-one via Experimental Charge Density Analysis

Experimental and theoretical charge density analyses on 2,2-dibromo-2,3-dihydroinden-1-one have been carried out to quantify the topological features of a short CBr....O halogen bond with nearly linear geometry (2.922 angstrom, angle CBr....O = 172.7 degrees) and to assess the strength of the interactions using the topological features of the electron density. The electrostatic potential map indicates the presence of the s-hole on bromine, while the interaction energy is comparable to that of a moderate OH....O hydrogen bond. In addition, the energetic contribution of CH.....Br interaction is demonstrated to be on par with that of the CBr....O halogen bond in stabilizing the crystal structure.

Computational Study of the Formation of Short Centrosymmetric N-center dot center dot center dot S Supramolecular Synthon and Related Weak Interactions in Crystalline 1,2,4-Triazoles

A comprehensive analysis of the crystal packing and the energetic features of a series of four biologically active molecules belonging to the family of substituted 4-(benzylideneamino)-3-(4-fluoro-3-phenoxyphenyl)-1H-1,2,4-triazole-5-(4 H)-thione derivatives have been performed based on the molecular conformation and the supramolecular packing. This involves the formation of a short centrosymmetric R-2(2)(8) NH...S supramolecular synthon in the solid state, including the presence of CH...S, CH...O, CH...N, CH...F, CH...Cl, CF...FC, CCl...ClC, and CH...pi intermolecular interactions along with pp stacking to evaluate the role of noncovalent interactions in the crystal. The presence of such synthons has a substantial contribution toward the interaction energy (-18 to -20 kcal/mol) as obtained from the PIXEL calculation, wherein the Coulombic and polarization contribution are more significant than the dispersion contribution. The geometrical characteristics of such synthons favor short distance, and the population of related molecules having these geometries is rare as has been obtained from the Cambridge Structural Database (CSD). Furthermore, their interaction energies have been compared with those present in our molecules in the solid state. The topological characteristics of the NH...S supramolecular synthon, in addition to related weak interactions, CH...N, CH...Cl, CF...FC, and CCl...ClC, have been estimated using the quantum theory of atoms in molecules (QTAIM). In addition, an analysis of the Hirshfeld surface and associated fingerprint plots of these four molecules also have provided a platform for the evaluation of the contribution of different atom...atom contacts, which contribute toward the packing of the molecules in solids.

Differential Cocrystallization Behavior of Isomeric Pyridine Carboxamides toward Antitubercular Drug Pyrazinoic Acid

Pyrazinoic acid, the active form of the antitubercular pro-drug Pyrazinamide, is an amphiprotic molecule containing carboxylic acid and pyridine groups and therefore can form both salts and cocrystals with relevant partner molecules. Cocrystallization of pyrazinoic acid with isomeric pyridine carboxamide series resulted in a dimorphic mixed-ionic complex with isonicotinamide and in eutectics with nicotinamide and picolinamide, respectively. It is observed that with alteration of the carboxamide position, steric and electrostatic compatibility issues between molecules of the combination emerge and affect intermolecular interactions and supramolecular growth, thus leading to either cocrystal or eutectic for different pyrazinoic acid-pyridine carboxamide combinations. Intermolecular interaction energy calculations have been performed to understand the role of underlying energetics on the formation of cocrystal/eutectic in different combinations. On the other hand, two molecular salts with piperazine and cytosine and a gallic acid cocrystal of the drug were obtained, and their X-ray crystal structures were also determined in this work.

The Self-Association of Graphane Is Driven by London Dispersion and Enhanced Orbital Interactions

We investigated the nature of the cohesive energy between graphane sheets via multiple CH center dot center dot center dot HC interactions, using density functional theory (DFT) including dispersion correction (Grimmes D3 approach) computations of n]graphane sigma dimers (n = 6-73). For comparison, we also evaluated the binding between graphene sheets that display prototypical pi/pi interactions. The results were analyzed using the block-localized wave function (BLW) method, which is a variant of ab initio valence bond (VB) theory. BLW interprets the intermolecular interactions in terms of frozen interaction energy (Delta E-F) composed of electrostatic and Pauli repulsion interactions, polarization (Delta E-pol), charge-transfer interaction (Delta E-CT), and dispersion effects (Delta E-disp). The BLW analysis reveals that the cohesive energy between graphane sheets is dominated by two stabilizing effects, namely intermolecular London dispersion and two-way charge transfer energy due to the sigma CH -> sigma*(HC) interactions. The shift of the electron density around the nonpolar covalent C-H bonds involved in the intermolecular interaction decreases the C-H bond lengths uniformly by 0.001 angstrom. The Delta E-CT term, which accounts for similar to 15% of the total binding energy, results in the accumulation of electron density in the interface area between two layers. This accumulated electron density thus acts as an electronic glue for the graphane layers and constitutes an important driving force in the self-association and stability of graphane under ambient conditions. Similarly, the double faced adhesive tape style of charge transfer interactions was also observed among graphene sheets in which it accounts for similar to 18% of the total binding energy. The binding energy between graphane sheets is additive and can be expressed as a sum of CH center dot center dot center dot HC interactions, or as a function of the number of C-H bonds.

Correlation equations for average deposition rate coefficients of nanoparticles in a cylindrical pore

Nanoparticle deposition behavior observed at the Darcy scale represents an average of the processes occurring at the pore scale. Hence, the effect of various pore-scale parameters on nanoparticle deposition can be understood by studying nanoparticle transport at pore scale and upscaling the results to the Darcy scale. In this work, correlation equations for the deposition rate coefficients of nanoparticles in a cylindrical pore are developed as a function of nine pore-scale parameters: the pore radius, nanoparticle radius, mean flow velocity, solution ionic strength, viscosity, temperature, solution dielectric constant, and nanoparticle and collector surface potentials. Based on dominant processes, the pore space is divided into three different regions, namely, bulk, diffusion, and potential regions. Advection-diffusion equations for nanoparticle transport are prescribed for the bulk and diffusion regions, while the interaction between the diffusion and potential regions is included as a boundary condition. This interaction is modeled as a first-order reversible kinetic adsorption. The expressions for the mass transfer rate coefficients between the diffusion and the potential regions are derived in terms of the interaction energy profile. Among other effects, we account for nanoparticle-collector interaction forces on nanoparticle deposition. The resulting equations are solved numerically for a range of values of pore-scale parameters. The nanoparticle concentration profile obtained for the cylindrical pore is averaged over a moving averaging volume within the pore in order to get the 1-D concentration field. The latter is fitted to the 1-D advection-dispersion equation with an equilibrium or kinetic adsorption model to determine the values of the average deposition rate coefficients. In this study, pore-scale simulations are performed for three values of Peclet number, Pe = 0.05, 5, and 50. We find that under unfavorable conditions, the nanoparticle deposition at pore scale is best described by an equilibrium model at low Peclet numbers (Pe = 0.05) and by a kinetic model at high Peclet numbers (Pe = 50). But, at an intermediate Pe (e.g., near Pe = 5), both equilibrium and kinetic models fit the 1-D concentration field. Correlation equations for the pore-averaged nanoparticle deposition rate coefficients under unfavorable conditions are derived by performing a multiple-linear regression analysis between the estimated deposition rate coefficients for a single pore and various pore-scale parameters. The correlation equations, which follow a power law relation with nine pore-scale parameters, are found to be consistent with the column-scale and pore-scale experimental results, and qualitatively agree with the colloid filtration theory. These equations can be incorporated into pore network models to study the effect of pore-scale parameters on nanoparticle deposition at larger length scales such as Darcy scale.

Ricin-membrane interaction: membrane penetration depth by fluorescence quenching and resonance energy transfer

The entry of the plant toxin ricin and its A- and B-subunits in model membranes in the presence as well as absence of monosialoganglioside (GM(1)) has been studied. Dioleoylphosphatidylcholine and 5-, 10-, and 12-doxyl- or 9,10-dibromophosphatidylcholines serve as quenchers of intrinsic tryptophan fluorescence of the proteins. The parallax method of Chattopadhyay and London [(1987) Biochemistry 26, 39-45] has been employed to measure the average membrane penetration depth of tryptophans of ricin and its B-chain and the actual depth of the sole Trp 211 in the A-chain. The results indicate that both of the chains as well as intact ricin penetrate the membrane deeply and the C-terminal end of the A-chain is well inside the bilayer, especially at pH 4.5. An extrinsic probe N-(iodoacetyl)-N'-(5-sulfo-1-naphthyl) ethylenediamine (I-AEDANS) has been attached to Cys 259 of the A-chain, and the kinetics of penetration has been followed by monitoring the increase in AEDANS fluorescence at 480 nm. The insertion follows first-order kinetics, and the rate constant is higher at a lower pH. The energy transfer distance analysis between Trp 211 and AEDANS points out that the conformation of the A-chain changes as it inserts into the membrane. CD studies indicate that the helicity of the proteins increases after penetration, which implies that some of the unordered structure in the native protein is converted to the ordered form during this process. Hydrophobic forces seem to be responsible for stabilizing a particular protein conformation inside the membrane.

Energy Hyperspace for Stacking Interaction in AU/AU Dinucleotide Step: Dispersion-Corrected Density Functional Theory Study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence-directed features, such as small roll values for the purine-pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near-zero slide to be most favorable for the purine-pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar-phosphate backbone with C3-endo sugars and this demands C1-C1 distance of about 5.4 angstrom along the chains. Consideration of an energy penalty term for deviation of C1-C1 distance from the mean value, to the recent DFT-D functionals, specifically B97X-D appears to predict reliable energy contour for AU/AU step. Such distance-based penalty improves energy contours for the other purine-pyrimidine sequences also. (c) 2013 Wiley Periodicals, Inc. Biopolymers 101: 107-120, 2014.

Singlet and triplet molecular interaction in excimers

A perturbative approach is developed to determine, term by term, the contributions of the various forces to the excimer potentials of the singlet and triplet excimers. The results show that the singlet excimer of naphthalene is more stable than the corresponding triplet excimer primarily due to large contributions of the exciton-resonance and the dispersion energy terms. The variation of the various terms with the conformations of the excimers suggests that the singlet and triplet excimers of naphthalene cannot have identical structure.

A study of interaction between two He+ ions in terms of the electronic forces

The potential energy curve of the He2+2 system dissociating into two He+ ions is examined in terms of the electronic force exerted on each nucleus as a function of the internuclear separation. The results are compared with the process of bond-formation in H2 from the separated atoms.

A study of interaction in the lowest singlet and triplet states of H2

The potential energy curves of the ground state and the first excited state of H2 are examined in terms of the electronic force acting on each nucleus. The results reveal the detailed course of events that occur when two hydrogen atoms with parallel and antiparallel electron spins approach one another from a large internuclear separation.

Sequence Specific Interaction between Promoter DNA and Escherichia coli RNA Polymerase: Comparative Thermodynamic Analysis with One Immobilized Partner

Sequence specific interaction between DNA and protein molecules has been a subject of active investigation for decades now. Here, we have chosen single promoter containing bacteriophage Delta D-III T7 DNA and Escherichia coli RNA polymerase and followed their recognition at the air-water interface by using the surface plasmon resonance (SPR) technique, where the movement of one of the reacting species is restricted by way of arraying them on an immobilized support. For the Langmuir monolayer studies, we used a RNA polymerase with a histidine tag attached to one of its subunits, thus making it an xcellent substrate for Ni(II) ions, while the SPR Studies were done using biotin-labeled DNA immobilized on a streptavidin-coated chip. Detailed analysis of the thermodynamic parameters as a function of concentration and temperature revealed that the interaction of RNA polymerase with T7 DNA is largely entropy driven (83 (+/- 12) kcal mol(-1)) with a positive enthalpy of 13.6 (+/- 3.6) kcal mol(-1), The free energy of reaction determined by SPR and Langmuir-Blodgett technique was -11 (+/- 2) and -15.6 kcal mol(-1), respectively. The ability of these methods to retain the specificity of the recognition process was also established.

Amino acid interaction preferences in proteins

Understanding the key factors that influence the interaction preferences of amino acids in the folding of proteins have remained a challenge. Here we present a knowledge-based approach for determining the effective interactions between amino acids based on amino acid type, their secondary structure, and the contact based environment that they find themselves in the native state structure as measured by their number of neighbors. We find that the optimal information is approximately encoded in a 60 x 60 matrix describing the 20 types of amino acids in three distinct secondary structures (helix, beta strand, and loop). We carry out a clustering scheme to understand the similarity between these interactions and to elucidate a nonredundant set. We demonstrate that the inferred energy parameters can be used for assessing the fit of a given sequence into a putative native state structure.

Excitation energy transfer from a fluorophore to single-walled carbon nanotubes

We study the process of electronic excitation energy transfer from a fluorophore to the electronic energy levels of a single-walled carbon nanotube. The matrix element for the energy transfer involves the Coulombic interaction between the transition densities on the donor and the acceptor. In the Foumlrster approach, this is approximated as the interaction between the corresponding transition dipoles. For energy transfer from a dye to a nanotube, one can use the dipole approximation for the dye, but not for the nanotube. We have therefore calculated the rate using an approach that avoids the dipole approximation for the nanotube. We find that for the metallic nanotubes, the rate has an exponential dependence if the energy that is to be transferred, h is less than a threshold and a d(-5) dependence otherwise. The threshold is the minimum energy required for a transition other than the k(i,perpendicular to)=0 and l=0 transition. Our numerical evaluation of the rate of energy transfer from the dye pyrene to a (5,5) carbon nanotube, which is metallic leads to a distance of similar to 165 A degrees up to which energy transfer is appreciable. For the case of transfer to semiconducting carbon nanotubes, apart from the process of transfer to the electronic energy levels within the one electron picture, we also consider the possibility of energy transfer to the lowest possible excitonic state. Transfer to semiconducting carbon nanotubes is possible only if>=epsilon(g)-epsilon(b). The long range behavior of the rate of transfer has been found to have a d(-5) dependence if h >=epsilon(g). But, when the emission energy of the fluorophore is in the range epsilon(g)>h >=epsilon(g)-epsilon(b), the rate has an exponential dependence on the distance. For the case of transfer from pyrene to the semiconducting (6,4) carbon nanotube, energy transfer is found to be appreciable up to a distance of similar to 175 A degrees.