Nuclear quantum effects and hydrogen bonding in liquids

We have employed ab initio path integral molecular dynamics simulations to investigate the role of nuclear quantum effects on the strength of hydrogen bonds in liquid hydrogen fluoride. Nuclear quantum effects are shown to be responsible for a stronger hydrogen bond and an enhanced dipole-dipole interaction, which lead, in turn, to a shortening of the H...F intrachain distance. The simulation results are analyzed in terms of the electronic density shifts with respect to a purely classical treatment of the nuclei. The observed enhanced hydrogen-bond interaction, which arises from a coupling of intra- and intermolecular effects, should be a general phenomenon occurring in all hydrogen-bonded systems.     

Nuclear quantum and H/D isotope effects on three-centered bonding diborane: Path integral molecular dynamics simulations

Nuclear quantum and H/D isotope effects of bridging and terminal hydrogen atoms of diborane (B₂H₆) molecules were systematically studied by classical ab initio molecular dynamics (CLMD) and ab initio path integral molecular dynamics (PIMD) simulations with BHandHLYP/6-31++G** level of theory at room temperature (298.15 K). Calculated results clearly show that H/D isotope effect appears in the distribution of hydrogen (deuterium) of B₂H₆ (B₂D₆). Geometry of B₂H₆ also plays a significant role in the nuclear quantum effect proved by PIMD simulations, but slightly deviated from its equilibrium structure when simulated via CLMD simulation. The bond lengths between boron atoms R (B1 … B2) and the bridging hydrogen atoms R_HH (H_B1 … H_B2) of the B₂H₆ molecule obtained from PIMD simulations are slightly longer than those of the deuterated form of the diborane (B₂D₆) molecule. The principal component analysis (PCA) was also employed to distinguish the important modes of bridging hydrogen as related to the nuclear quantum and H/D isotope effects. The highest level of contribution obtained from PCA of PIMD simulations is bending, while various mixed vibrations with less contribution were also found. Therefore, the nuclear quantum and H/D isotope effects need to be taken into account for a better understanding of diborane geometry.     

Two-center molecular repulsion integrals over slater functions

For the general two-electron two-center integral over Slater functions, use of the Neumann expansion for the electron-electron interaction term yields the standard auxiliary functions. These are expanded and integrated explicitly by two independent methods. The resulting simple analytic formula for the total integral is completely general, requiring only the Slater function quantum numbers and exponents and the internuclear separation. Hence all two-electron hydrid, coulomb, exchange, and one-center integrals are considered. The efficiency of calculation of this expression is compared with those of other methods, indicating an order of magnitude improvement in speed over recursion for the exchange integral.     

Exponential repulsion improves structural predictability of molecular docking

Molecular docking is a powerful tool for theoretical prediction of the preferred conformation and orientation of small molecules within protein active sites. The obtained poses can be used for estimation of binding energies, which indicate the inhibition effect of designed inhibitors, and therefore might be used for in silico drug design. However, the evaluation of ligand binding affinity critically depends on successful prediction of the native binding mode. Contemporary docking methods are often based on scoring functions derived from molecular mechanical potentials. In such potentials, nonbonded interactions are typically represented by electrostatic interactions between atom‐centered partial charges and standard 6–12 Lennard–Jones potential. Here, we present implementation and testing of a scoring function based on more physically justified exponential repulsion instead of the standard Lennard–Jones potential. We found that this scoring function significantly improved prediction of the native binding modes in proteins bearing narrow active sites such as serine proteases and kinases. © 2016 Wiley Periodicals, Inc.     

Hydrogen bonding effects on molecular structure: crystallographic studies of barbiturates

X-ray crystal structure determinations show systematic perturbations of the order of 0.01 Å in the oxopyrimidine rings of twelve barbiturate molecules. They appear to be related to the mode of intermolecular NHβO hydrogen bonding.     

Nuclear quantum effects induce metallization of dense solid molecular hydrogen

We present an accurate computational study of the electronic structure and lattice dynamics of solid molecular hydrogen at high pressure. The band‐gap energies of the , Pc, and structures at pressures of 250, 300, and 350 GPa are calculated using the diffusion quantum Monte Carlo (DMC) method. The atomic configurations are obtained from ab initio path‐integral molecular dynamics (PIMD) simulations at 300 K and 300 GPa to investigate the impact of zero‐point energy and temperature‐induced motion of the protons including anharmonic effects. We find that finite temperature and nuclear quantum effects reduce the band‐gaps substantially, leading to metallization of the and Pc phases via band overlap; the effect on the band‐gap of the structure is less pronounced. Our combined DMC‐PIMD simulations predict that there are no excitonic or quasiparticle energy gaps for the and Pc phases at 300 GPa and 300 K. Our results also indicate a strong correlation between the band‐gap energy and vibron modes. This strong coupling induces a band‐gap reduction of more than 2.46 eV in high‐pressure solid molecular hydrogen. Comparing our DMC‐PIMD with experimental results available, we conclude that none of the structures proposed is a good candidate for phases III and IV of solid hydrogen. © 2017 Wiley Periodicals, Inc.     

Hydrogen bonding in clusters and molecular crystals

Ab initio calculations on the Hartree–Fock level are used to study the effects of cooperativity in ternary complexes and in infinite chains on the example of hydrogen fluoride. SCF energy partitioning demonstrates that polarization forces dominate three-body energies in these aggregates. The crystal orbital method is applied to investigate the structural differences between dimers and molecular crystals.     

Nuclear quantum effects at aqueous metal interfaces captured by molecular dynamics simulations

In this review, we summarize the recent development in modeling nuclear quantum effects at aqueous metal interfaces. First, we review the nuclear quantum effects on the water-metal interface at ultrahigh vacuum. Then, we illustrate the nuclear quantum effects at the potential of zero charge conditions. At last, we give some outlook for the perspective work in modeling the nuclear quantum effects at electrochemical interfaces and some practical simulation strategies.     

Hydrogen bonding control of molecular self-assembly

In this short review we describe approaches to the design and construction of synthetic molecules that mimic the process of self organization that is at the heart of biological complexity. Multi-subunit enzymes, viruses, and higher order DNA structures are formed by the non-covalent association of many smaller components. This self-assembly is controlled by the nature, number and orientation of interacting groups on the surface of the subunits. The central problem lies in overcoming the unfavorable entropy of multi-subunit association by significant enthalpic contribution from the binding of complementary regions on the subunits. We will place particular emphasis on the design of synthetic molecules that use hydrogen bonding interactions to control the formation of aggregates of well-defined structure.     

Fluorescence quenching in molecular recognition by hydrogen bonding

Steady-state fluorescence and single photon timing have been used to study the effect of the presence of hydrogen bonding on the intermolecular quenching of pyrene covalently linked to a guanine-like receptor I by an aliphatic amine (N,N-dimethylpropylamine) covalently linked to cytosine derivative II. By comparing the fluorescence quenching of I by II with that of 10methylpyrene (1-MP) by triethylamine (TEA), as a model system in which no hydrogen bonding can occur, one could possibly analyze the effect of the hydrogen bonding between receptor and substrate as a quenching as it leads to a higher local concentration of donor and acceptor. While the quenching of I by II was observed with an apparent rate constant k_q of (1.78 ± 0.10) × 10⁹ M⁻¹ s⁻¹ and (8.72 ± 0.42) × 10⁸ M⁻¹ s⁻¹ in toluene and acetonitrile, respectively, no quenching could be observed in methanol. Upon excitation of 1-MP, no quenching by II could be detected in the same concentration range as used in the quenching of I. Quenching of I and of 1-MP by TEA (⩾ 10⁻² M) in toluene leads to exciplex formation with maxima centred at 540 and 514 nm, respectively. The rate constants of exciplex formation and dissociation of I with TEA were analyzed using a global compartmental analysis. The following values were obtained for the rate constants: k₀₁ = (9.70 ± 0.01) × 10⁶ s⁻¹, k₂₁ = (1.12 ± 0.003) × 10⁹ M⁻¹ s⁻¹, k₀₂ = (5.24 ± 0.01) × 10⁷ s⁻¹ and k₁₂ = (7.74 ± 0.08) × 10⁶ s⁻¹. Quenching of I by TEA in the presence of III, a hydrogen-bonding system without an alkyl amine substituent, leads to exciplex formation centred at 538 nm. The rate constant values for the exciplex formation and dissociation of I with TEA in the presence of III were: k₀₁ = (9.32 ± 0.08) × 10⁶ s⁻¹, k₂₁ = (9.32 ± 0.003) × 10⁸ M⁻¹ s⁻¹, k₀₂ = (6.16 ± 0.03) × 10⁷ s⁻¹ and k₁₂ = (21.90 ± 0.3) × 10⁶ s⁻¹. The apparent rate constants k_q for this system was (7.26 ± 0.56) × 10⁶ M⁻¹ s⁻¹. The observed decrease in the rate of exciplex formation of I with TEA in the presence of III could suggest that the guanine-like moiety in I forms hydrogen bonds with the cytosine-like moiety and this could decrease the electron affinity of I. The rate constant of exciplex dissociation increased, indicating that the exciplex is less stable in the presence of III. Because of the single exponential decay of I in the presence and absence of II and of the agreement between steady-state and transient fluorescence measurements, the information available for quantitative analysis of the association between I and II is limited.     

Poisson equation for molecular exchange,hybrid and coulomb electron repulsion integrals

The various multicenter exchange, hybrid and Coulomb electron repulsion integrals that occur in molecular quantum mechanics are shown to satisfy a Poisson equation in which an overlap integral plays the role of a source distribution function. Two-, three-and four-center exchange integrals arise from four-center source functions; two- and three-center hybrid integrals arise from three-center distributions; and one- and two-center Coulomb integrals have two-center sources.     

Nuclear electric dipole shielding in molecular ions

A formula, which includes the effects of finite nuclear masses, is derived for the force on a nucleus in a (possibly molecular) ion in a spatially uniform (possibly time dependent) electric field.     

Hydrogen bonding in molecular recognition by HIV-1 protease

Molecular recognition of the cleavage sites of the substrates by HIV-1 protease is analyzed in terms of hydrogen bonding. Crystal structures of an inactive enzyme complexed with six different substrates were used as reference structures. Applying molecular mechanics calculations it can be shown that the interaction energies between the real substrate and the enzyme are larger than with other peptides. From the analysis, it can be concluded that water molecules are essential in the recognition process. Moreover, the hydrogen bonds between the protease and various substrates are characterized in detail.     

Hydrogen bonding and molecular association in methylated 3-hydroxypyridines

Association in derivatives of 3-hydroxypyridine has been examined via the chemical shifts for OH in PMR for equimolar mixtures of the derivatives of pyridine and phenol in CCl₄ and CH₂Cl₂. It is found that spatial screening of the pyridine N and phenol OH by methyl groups has an anomalous effect on the strength of the complex. A possible explanation is that the ionic form NH . . . O is stabilized on account of increased basicity of the N and acidity of the H in the methyl derivatives of pyridine and phenol. Domination by the intermolecular hydrogen bond between OH and N is the distinctive feature of association in 3-hydroxypyridines relative to association in phenol. The temperature dependence of _OH gives the enthalpy and entropy of association for 2, 4, 6-trimethyl-3-hydroxypyridine and of complex formation for various model systems. IR measurements confirm the conclusions.See [1] for the previous communication.We are indebted to G. G. Dvoryantseva for making the IR measurements.     

Elastic effects behind cooperative bonding in beta-sheets

We present extensive density functional theory calculations of the bonding between strands in beta-sheets. We identify a significant cooperative effect whereby the interaction increases in strength with the number of strands. We show that the effect is related to a coupling between interstrand bonding and intrastrand elastic properties. It is found that a direct consequence of this coupling is that the pitch of beta-sheets should contract with increasing number of strands, and we show that the effect can be observed directly in experimental data from the Protein Data Bank.