The quantitative characterization of residue contributions to protein-protein binding across extensive flexible interfaces poses a significant challenge for biophysical computations. It is attributable to the inherent imperfections in the experimental structures themselves, as well as to the lack of reliable computational tools for the evaluation of all types of noncovalent interactions. This study leverages recent advancements in semiempirical quantum-mechanical and implicit solvent approaches embodied in the PM6-D3H4S/COSMO2 method for the development of a hierarchical computational protocols encompassing molecular dynamics, fragmentation, and virtual glycine scan techniques for the investigation of flexible protein-protein interactions. As a model, the binding of insulin to its receptor is selected, a complex and dynamic process that has been extensively studied experimentally. The interaction energies calculated at the PM6-D3H4S/COSMO2 level in ten molecular dynamics snapshots did not correlate with molecular mechanics/generalized Born interaction energies because only the former method is able to describe nonadditive effects. This became evident by the examination of the energetics in small-model dimers featuring all the present types of noncovalent interactions with respect to DFT-D3 calculations. The virtual glycine scan has identified 15 hotspot residues on insulin and 15 on the insulin receptor, and their contributions have been quantified using PM6-D3H4S/COSMO2. The accuracy and credibility of the approach are further supported by the fact that all the insulin hotspots have previously been detected by biochemical and structural methods. The modular nature of the protocol has enabled the formulation of several variants, each tailored to specific accuracy and efficiency requirements. The developed computational strategy is firmly rooted in general biophysical chemistry and is thus offered as a general tool for the quantification of interactions across relevant flexible protein-protein interfaces.

• MeSH
• Insulin metabolism   chemistry   MeSH
• Protein Conformation MeSH
• Humans MeSH
• Receptor, Insulin * chemistry   metabolism   MeSH
• Molecular Dynamics Simulation * MeSH
• Thermodynamics MeSH
• Protein Binding MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Insulin MeSH
• Receptor, Insulin * MeSH

Previously studied complexes with protonic and hydridic hydrogen bonds exhibit significant similarities. The present study provides a detailed investigation of the structure, stabilization, electronic properties, and spectral characteristics of protonic and hydridic hydrogen bonds using low-temperature infrared (IR) spectroscopy and computational methods. Complexes of pentafluorobenzene with ammonia (C₆F₅H⋯NH₃) and triethylgermane with trifluoroiodomethane (Et₃GeH⋯ICF₃) were analyzed using both experimental and computational tools. Additionally, 30 complexes with protonic hydrogen bonds and 30 complexes with hydridic hydrogen bonds were studied computationally. Our findings reveal that, despite the opposite atomic charges on the hydrogens in these hydrogen bonds, and consequently the opposite directions of electron transfer in protonic and hydridic hydrogen bonds, their spectral manifestations - specifically, the red shifts in the X-H stretching frequency and the increase in intensity - are remarkably similar. The study also discusses the limitations of the current IUPAC definition of hydrogen bonding in covering both types of H-bonds and suggests a way to overcome these limitations.

• Publication type
• Journal Article MeSH

Accurate estimates of intermolecular interaction energy, ΔE, are crucial for modeling the properties of organic electronic materials and many other systems. For a diverse set of 50 dimers comprising up to 50 atoms (Set50-50, with 7 of its members being models of single-stacking junctions), benchmark ΔE data were compiled. They were obtained by the focal-point strategy, which involves computations using the canonical variant of the coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] performed while applying a large basis set, along with extrapolations of the respective energy components to the complete basis set (CBS) limit. The resulting ΔE data were used to gauge the performance for the Set50-50 of several density-functional theory (DFT)-based approaches, and of one of the localized variants of the CCSD(T) method. This evaluation revealed that (1) the proposed "silver standard" approach, which employs the localized CCSD(T) method and CBS extrapolations, can be expected to provide accuracy better than two kJ/mol for absolute values of ΔE, and (2) from among the DFT techniques, computationally by far the cheapest approach (termed "ωB97X-3c/vDZP" by its authors) performed remarkably well. These findings are directly applicable in cost-effective yet reliable searches of the potential energy surfaces of noncovalent complexes.

• Keywords
• CCSD(T), DFT, interaction energy, noncovalent interactions, supramolecular junctions,
• MeSH
• Dimerization MeSH
• Electronics * MeSH
• Physical Phenomena MeSH
• Polymers MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Polymers MeSH

There have been attempts, both experimental and based on density-functional theory (DFT) modeling, at understanding the factors that govern the electronic conductance behavior of single-stacking junctions formed by pi-conjugated materials in nanogaps. Here, a reliable description of relevant stacked configurations of some thiophene-cored systems is provided by means of high-level quantum chemical approaches. The minimal structures of these configurations, which are found using the dispersion-corrected DFT approach, are employed in calculations that apply the coupled cluster method with singles, doubles and perturbative triples [CCSD(T)] and extrapolations to the complete basis set (CBS) limit in order to reliably quantify the strength of intermolecular binding, while their physical origin is investigated using the DFT-based symmetry-adapted perturbation theory (SAPT) of intermolecular interactions. In particular, for symmetrized S-Tn dimers (where "S" and "T" denote a thiomethyl-containing anchor group and a thiophene segment comprising "n" units, respectively), the CCSD(T)/CBS interaction energies are found to increase linearly with n ≤ 6, and significant conformational differences between the flanking 2-thiophene group in S-T1 and S-T2 are described by the CCSD(T)/CBS and SAPT/CBS computations. These results are put into the context of previous work on charge transport properties of S-Tn and other types of supramolecular junctions.

• Keywords
• CCSD(T), DFT, intermolecular stacking, oligothiophenes, supramolecular junctions,
• MeSH
• Electronics * MeSH
• Polymers * MeSH
• Thiophenes MeSH
• Virion MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Polymers * MeSH
• Thiophenes MeSH

Submolecular charge distribution significantly affects the physical-chemical properties of molecules and their mutual interaction. One example is the presence of a π-electron-deficient cavity in halogen-substituted polyaromatic hydrocarbon compounds, the so-called π-holes, the existence of which was predicted theoretically, but the direct experimental observation is still missing. Here we present the resolution of the π-hole on a single molecule using the Kelvin probe force microscopy, which supports the theoretical prediction of its existence. In addition, experimental measurements supported by theoretical calculations show the importance of π-holes in the process of adsorption of molecules on solid-state surfaces. This study expands our understanding of the π-hole systems and, at the same time, opens up possibilities for studying the influence of submolecular charge distribution on the chemical properties of molecules and their mutual interaction.

• Publication type
• Journal Article MeSH

It has been recognized that the C-H⋯O structural motif can be present in destabilizing as well as highly stabilizing intermolecular environments. Thus, it should be of interest to describe the strength of the C-H⋯O hydrogen bond for constant structural factors so that this intrinsic strength can be quantified and compared to other types of interactions. This description is provided here for C2h-symmetric dimers of acrylic acid by means of the calculations that employ the coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] together with an extrapolation to the complete basis set (CBS) limit. Dimers featuring the C-H⋯O and O-H⋯O hydrogens bonds are carefully investigated in a wide range of intermolecular separations by the CCSD(T)/CBS approach, and also by the symmetry-adapted perturbation theory (SAPT) method, which is based on the density-functional theory (DFT) treatment of monomers. While the nature of these two types of hydrogen bonding is very similar according to the SAPT-DFT/CBS calculations and on the basis of a comparison of the intermolecular potential curves, the intrinsic strength of the C-H⋯O interaction is found to be about a quarter of its O-H⋯O counterpart that is less than one might anticipate.

• Keywords
• CCSD(T), DFT, SAPT, hydrogen bonding, noncovalent interactions,
• MeSH
• Polymers * MeSH
• Density Functional Theory MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Polymers * MeSH

We have performed a hierarchical ab initio benchmark and DFT performance study of D2 Ch•••A- chalcogen bonds (Ch = S, Se; D, A = F, Cl). The ab initio benchmark study is based on a series of ZORA-relativistic quantum chemical methods [HF, MP2, CCSD, CCSD(T)], and all-electron relativistically contracted variants of Karlsruhe basis sets (ZORA-def2-SVP, ZORA-def2-TZVPP, ZORA-def2-QZVPP) with and without diffuse functions. The highest-level ZORA-CCSD(T)/ma-ZORA-def2-QZVPP counterpoise-corrected complexation energies (ΔECPC ) are converged within 1.1-3.4 kcal mol-1 and 1.5-3.1 kcal mol-1 with respect to the method and basis set, respectively. Next, we used the ZORA-CCSD(T)/ma-ZORA-def2-QZVPP (ΔECPC ) as reference data for analyzing the performance of 13 different ZORA-relativistic DFT approaches in combination with the Slater-type QZ4P basis set. We find that the three-best performing functionals are M06-2X, B3LYP, and M06, with mean absolute errors (MAE) of 4.1, 4.2, and 4.3 kcal mol-1 , respectively. The MAE for BLYP-D3(BJ) and PBE amount to 8.5 and 9.3 kcal mol-1 , respectively.

• Keywords
• benchmark study, chalcogen bonds, coupled-cluster, density functional calculations, noncovalent interactions,
• Publication type
• Journal Article MeSH

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

• MeSH
• DNA chemistry   MeSH
• Catalysis MeSH
• Nucleic Acid Conformation * MeSH
• Computer Simulation MeSH
• RNA chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Review MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• DNA MeSH
• RNA MeSH

General and reliable description of structures and energetics in protein-ligand (PL) binding using the docking/scoring methodology has until now been elusive. We address this urgent deficiency of scoring functions (SFs) by the systematic development of corrected semiempirical quantum mechanical (SQM) methods, which correctly describe all types of noncovalent interactions and are fast enough to treat systems of thousands of atoms. Two most accurate SQM methods, PM6-D3H4X and SCC-DFTB3-D3H4X, are coupled with the conductor-like screening model (COSMO) implicit solvation model in so-called "SQM/COSMO" SFs and have shown unique recognition of native ligand poses in cognate docking in four challenging PL systems, including metalloprotein. Here, we apply the two SQM/COSMO SFs to 17 diverse PL complexes and compare their performance with four widely used classical SFs (Glide XP, AutoDock4, AutoDock Vina, and UCSF Dock). We observe superior performance of the SQM/COSMO SFs and identify challenging systems. This method, due to its generality, comparability across the chemical space, and lack of need for any system-specific parameters, gives promise of becoming, after comprehensive large-scale testing in the near future, a useful computational tool in structure-based drug design and serving as a reference method for the development of other SFs.

• Publication type
• Journal Article MeSH

Understanding strength and nature of noncovalent binding to surfaces imposes significant challenge both for computations and experiments. We explored the adsorption of five small nonpolar organic molecules (acetone, acetonitrile, dichloromethane, ethanol, ethyl acetate) to fluorographene and fluorographite using inverse gas chromatography and theoretical calculations, providing new insights into the strength and nature of adsorption of small organic molecules on these surfaces. The measured adsorption enthalpies on fluorographite range from -7 to -13 kcal/mol and are by 1-2 kcal/mol lower than those measured on graphene/graphite, which indicates higher affinity of organic adsorbates to fluorographene than to graphene. The dispersion-corrected functionals performed well, and the nonlocal vdW DFT functionals (particularly optB86b-vdW) achieved the best agreement with the experimental data. Computations show that the adsorption enthalpies are controlled by the interaction energy, which is dominated by London dispersion forces (∼70%). The calculations also show that bonding to structural features, like edges and steps, as well as defects does not significantly increase the adsorption enthalpies, which explains a low sensitivity of measured adsorption enthalpies to coverage. The adopted Langmuir model for fitting experimental data enabled determination of adsorption entropies. The adsorption on the fluorographene/fluorographite surface resulted in an entropy loss equal to approximately 40% of the gas phase entropy.

• Publication type
• Journal Article MeSH