Performance of the DFT-D method, paired with the PCM implicit solvation model, for the computation of interaction energies of solvated complexes of biological interest

Copper-zinc/chitosan complex hydrogels: Rheological, degradation and biological properties

Polymer hydrogels crosslinked by therapeutic metal ions have attracted increased interest in recent years due to their unique and versatile properties. Chitosan hydrogels are widely investigated for various biomedical applications such as tissue engineering and drug delivery. Copper and zinc ions are considered as therapeutic metal ions, that have important roles in bone regeneration. The aim of this study was to investigate the physicochemical and biological properties of bimetallic-chitosan complex hydrogels with different cupric and zinc ions content. Scanning electron microscopy (SEM) revealed changes in the morphology from the microstructure with larger, tubular pores for aerogels with higher Zn content, to the sheets-like structure with long pores for samples with higher Cu content. FTIR analysis indicated the formation of bimetallic-chitosan aerogels. The obtained X-ray diffraction patterns showed a broadening of chitosan's characteristic diffraction maximum, while characterization of physical properties showed decreased swelling ability and increased shear modulus with higher Cu content. ICP-MS results showed a negligible amount of copper and zinc ions released under physiological conditions during 24 h indicating a strong physical crosslink between metal ions and chitosan chains. Furthermore, accelerated in vitro degradation showed that hydrogels maintained good stability during four weeks of lysozyme activity. The MTT assay indicated that the cytotoxicity of Cu2+-Zn2+/chitosan complexes could be adjusted by the amount of cupric ions. All results imply that Cu2+ and Zn2+ ions act as physical crosslinkers of the polymer network. Also, results are in agreement with the prediction of density functional theory (DFT) which indicated stronger chitosan-Cu tetrahedral aqua complex interactions in comparison to the chitosan-[Zn(H2O)4]2+ interactions.

Investigating Paracetamol's Role as a Potential Treatment for Parkinson's Disease: Ab Initio Analysis of Dopamine, l-DOPA, Paracetamol, and NAPQI Interactions with Enzymes Involved in Dopamine Metabolism

Recently, it was found that paracetamol can extend the therapeutic window of l-DOPA treatment for Parkinson's disease [Golding (2019) BJPharm, 4(2), Article 619]. It has been posited that the effect could be due to paracetamol and its metabolite, NAPQI, inhibiting pain signals in the spinal column. In this work, we examine the possibility that the therapeutic effect of the paracetamol for the Parkinson's disease patient may be due to an inhibition of the enzymes that metabolize dopamine and/or l-DOPA, thus effectively extending the lifetime of the l-DOPA treatment. In this work, we use the M062X/6-311+G* level of theory to calculate the electronic binding energies (including explicit desolvation) of several ligands (paracetamol, NAPQI, dopamine, and l-DOPA) with a series of enzymes important to the production and metabolism of dopamine and compare them to calculated binding energy values for the natural substrates for those enzymes in order to predict possible inhibition. Benchmark interaction energies for a subset of the systems studied are calculated using the more accurate second-order Møller-Plesset perturbation (MP2) method in order to calibrate the accuracy of the M062X method. If we assume that the interaction energies calculated here can serve as a proxy for in vivo inhibition, then we can predict that paracetamol and NAPQI should not inhibit the natural production of dopamine and may in fact inhibit the metabolism of l-DOPA and dopamine, thus extending the length of l-DOPA treatments.

Assessing the persistence of chalcogen bonds in solution with neural network potentials

Non-covalent bonding patterns are commonly harvested as a design principle in the field of catalysis, supramolecular chemistry and functional materials to name a few. Yet, their computational description generally neglects finite temperature and environment effects, which promote competing interactions and alter their static gas-phase properties. Recently, neural network potentials (NNPs) trained on Density Functional Theory (DFT) data have become increasingly popular to simulate molecular phenomena in condensed phase with an accuracy comparable to ab initio methods. To date, most applications have centered on solid-state materials or fairly simple molecules made of a limited number of elements. Herein, we focus on the persistence and strength of chalcogen bonds involving a benzotelluradiazole in condensed phase. While the tellurium-containing heteroaromatic molecules are known to exhibit pronounced interactions with anions and lone pairs of different atoms, the relevance of competing intermolecular interactions, notably with the solvent, is complicated to monitor experimentally but also challenging to model at an accurate electronic structure level. Here, we train direct and baselined NNPs to reproduce hybrid DFT energies and forces in order to identify what are the most prevalent non-covalent interactions occurring in a solute-Cl$^-$-THF mixture. The simulations in explicit solvent highlight competition with chalcogen bonds formed with the solvent and the short-range directionality of the interaction with direct consequences for the molecular properties in the solution. The comparison with other potentials (e.g., AMOEBA, direct NNP and continuum solvent model) also demonstrates that baselined NNPs offer a reliable picture of the non-covalent interaction interplay occurring in solution.

B-DNA Structure and Stability: The Role of Nucleotide Composition and Order

We have quantum chemically analyzed the influence of nucleotide composition and sequence (that is, order) on the stability of double-stranded B-DNA triplets in aqueous solution. To this end, we have investigated the structure and bonding of all 32 possible DNA duplexes with Watson-Crick base pairing, using dispersion-corrected DFT at the BLYP-D3(BJ)/TZ2P level and COSMO for simulating aqueous solvation. We find enhanced stabilities for duplexes possessing a higher GC base pair content. Our activation strain analyses unexpectedly identify the loss of stacking interactions within individual strands as a destabilizing factor in the duplex formation, in addition to the better-known effects of partial desolvation. Furthermore, we show that the sequence-dependent differences in the interaction energy for duplexes of the same overall base pair composition result from the so-called "diagonal interactions" or "cross terms". Whether cross terms are stabilizing or destabilizing depends on the nature of the electrostatic interaction between polar functional groups in the pertinent nucleobases.

Mechanistic Insights into the Inhibition of SARS-CoV-2 Main Protease by Clovamide and Its Derivatives: In Silico Studies

The novel coronavirus SARS-CoV-2 Main Protease (Mpro) is an internally encoded enzyme that hydrolyzes the translated polyproteins at designated sites. The protease directly mediates viral replication processes; hence, a promising target for drug design. Plant-based natural products, especially polyphenols and phenolic compounds, provide the scaffold for many effective antiviral medications, and have recently been shown to be able to inhibit Mpro of SARS-CoV-2. Specifically, polyphenolic compounds found in cacao and chocolate products have been shown by recent experimental studies to have strong inhibitory effects against Mpro activities. This work aims to uncover the inhibition processes of Mpro by a natural phenolic compound found in cacao and chocolate products, clovamide. Clovamide (caffeoyl-DOPA) is a naturally occurring caffeoyl conjugate that is found in the phenolic fraction of Theobroma Cacao L. and a potent radical-scavenging antioxidant as suggested by previous studies of our group. Here, we propose inhibitory mechanisms by which clovamide may act as a Mpro inhibitor as it becomes oxidized by scavenging reactive oxygen species (ROS) in the body, or becomes oxidized as a result of enzymatic browning. We use molecular docking, annealing-based molecular dynamics, and Density Functional Theory (DFT) calculations to study the interactions between clovamide with its derivatives and Mpro catalytic and allosteric sites. Our molecular modelling studies provide mechanistic insights of clovamide inhibition of Mpro, and indicate that clovamide may be a promising candidate as a drug lead molecule for COVID-19 treatments.

Amino Acid Interactions (INTAA) web server v2.0: a single service for computation of energetics and conservation in biomolecular 3D structures

Interactions among amino acid residues are the principal contributor to the stability of the three-dimensional structure of a protein. The Amino Acid Interactions (INTAA) web server (https://bioinfo.uochb.cas.cz/INTAA/) has established itself as a unique computational resource, which enables users to calculate the contribution of individual residues in a biomolecular structure to its total energy using a molecular mechanical scoring function. In this update, we describe major additions to the web server which help solidify its position as a robust, comprehensive resource for biomolecular structure analysis. Importantly, a new continuum solvation model was introduced, allowing more accurate representation of electrostatic interactions in aqueous media. In addition, a low-overhead pipeline for the estimation of evolutionary conservation in protein chains has been added. New visualization options were introduced as well, allowing users to easily switch between and interrelate the energetic and evolutionary views of the investigated structures.

Exploring Epigenetic Marks by Analysis of Noncovalent Interactions

Epigenetic marks are modest chemical modifications on DNA and histone proteins that regulate the activation or silencing of genes through modulation of the intermolecular interactions between the DNA strands and the protein machinery. The process is complex and not always well understood. One of the systems studied in greater detail is the epigenetic mark on H3K9: lysine 9 of histone 3. The degree of methylation or acetylation of this histone is linked to silencing or activation of the corresponding gene, but it is not clear which effect each mark has in gene expression. We shed light on this particular methylation process by using density functional theory (DFT) calculations. First, we built a model consisting of a DNA double strand containing three base pairs and a sequence of three amino acids of the histone's tail. Then, we computed the modulation introduced into the intermolecular interactions by each epigenetic modification: from mono- to trimethylation and acetylation. The calculations show that whereas acetylation and trimethylation result in a reduction of the DNA-peptide interaction; non-, mono-, and dimethylation increase the intermolecular interactions. Such observations compare well with the findings reported in the literature, and highlight the correlation between the balance of intermolecular forces and biological properties, simultaneously advancing quantum-mechanical studies of large biochemical systems at molecular level through the use of DFT methods.

Effect analysis of quantum chemical descriptors and substituent characteristics on Henry's law constants of polybrominated diphenyl ethers at different temperatures

Twelve substituent descriptors, 17 quantum chemical descriptors and 1/T were selected to establish a quantitative structure-property relationship (QSPR) model of Henry's law constants for 7 polybrominated diphenyl ethers (PBDEs) at five different temperatures. Then, the lgH of 202 congeners at different temperatures were predicted. The variation rule and regulating mechanism of lgH was studied from the perspectives of both quantum chemical descriptors and substituent characteristics. The R² for modeling and testing sets of the final QSPR model are 0.977 and 0.979, respectively, thus indicating good fitness and predictive ability for Henry' law constants of PBDEs at different temperatures. The favorable hydrogen binding sites are the 5,5',6,6'-positions for high substituent congeners and the O atom of the ether bond for low substituent congeners, which affects the interaction between PBDEs and water molecules. lgH is negatively and linearly correlated with 1/T, and the variation trends of lgH with temperature are primarily regulated by individual substituent characteristics, wherein: the more substituents involved, the smaller the lgH. The significant sequence for the main effect of substituent positions is para>meta>ortho, where the ortho-positions are mainly involved in second-order interaction effect (64.01%). Having two substituents in the same ring also provides a significant effect, with 81.36% of second-order interaction effects, particularly where there is an adjacent distribution (55.02%).

A theoretical study on the characteristics of the intermolecular interactions in the active site of human androsterone sulphotransferase: DFT calculations of NQR and NMR parameters and QTAIM analysis

A theoretical study at the level of density functional theory (DFT) was performed to characterize noncovalent intermolecular interactions, especially hydrogen bond interactions, in the active site of enzyme human androsterone sulphotransferase (SULT2A1/ADT). Geometry optimization, interaction energy, (2)H, (14)N, and (17)O electric field gradient (EFG) tensors, (1)H, (13)C, (17)O, and (15)N chemical shielding (CS) tensors, Natural Bonding Orbital (NBO) analysis, and quantum theory of atoms in molecules (QTAIM) analysis of this active site were investigated. It was found that androsterone (ADT) is able to form hydrogen bonds with residues Ser80, Ile82, and His99 of the active site. The interaction energy calculations and NBO analysis revealed that the ADT molecule forms the strongest hydrogen bond with Ser80. Results revealed that ADT interacts with the other residues through electrostatic and Van der Waals interactions. Results showed that these hydrogen bonds influence on the calculated (2)H, (14)N, and (17)O quadrupole coupling constants (QCCs), as well as (1)H, (13)C, (17)O, and (15)N CS tensors. The magnitude of the QCC and CS changes at each nucleus depends directly on its amount of contribution to the hydrogen bond interaction.

Representative Amino Acid Side Chain Interactions in Proteins. A Comparison of Highly Accurate Correlated ab Initio Quantum Chemical and Empirical Potential Procedures

Interactions between amino acid side chains play a crucial role both within a folded protein and between the interacting protein molecules. Here we have selected a representative set of 24 of the 400 (20 × 20) possible interacting side chain pairs based on data from Atlas of Protein Side-Chain Interactions. For each pair, we obtained its most favorable interaction geometry from the structural data and computed the interaction energy in the gas phase using several different, commonly used, ab initio and force field methods, namely Møller-Plesset perturbation theory (MP2), density functional theory combined with symmetry-adapted perturbation theory (DFT-SAPT), density functional theory empirically augmented with an empirical dispersion term (DFT-D), and empirical potentials using the OPLS-AA/L and Amber03 force fields. All the methods were compared against a reference method taken to be the CCSD(T) level of theory extrapolated to the complete basis set limit. We found a high degree of agreement between the different methods, even though the range of binding energies obtained was extremely large. The most computationally intensive methods yielded the best results. Among the less computationally time-consuming methods, the DFT-D method as well as parm03 force field provided consistently good results when compared to the reference values. We also tested how representative the chosen geometries of the side chains were and investigated the effect on the binding energies of the dielectric constant of the surrounding medium.

DFT STUDY OF THE PROTONATION AND DEPROTONATION ENTHALPIES OF BENZOXAZOLE, 1,2-BENZISOXAZOLE AND 2,1-BENZISOXAZOLE AND IMPLICATIONS FOR THE STRUCTURES AND ENERGIES OF THEIR ADDUCTS WITH EXPLICIT WATER MOLECULES

Benzoxazole, 1,2-benzisoxazole and 2,1-benzisoxazole are biologically active molecules with potential applications in drug design. Their interaction with aqueous medium in biological systems may be simulated by considering their interaction with explicit water molecules. Such studies provide information on the structures, energies and type of interactions stabilizing the resulting geometric systems. The objective of the current study was to utilize theoretical approaches to investigate the structures, stabilization energy and binding energy of benzoxazole–water, 1,2-benzisoxazole–water and 2,1-benzisoxazole–water complexes. The calculations were performed utilizing the density functional theory (DFT)/M06-2X/6-311 ++ G(d,p) method and the DFT/ωB97XD method with both the 6-311 ++ G(d,p) and the aug-cc-pVDZ basis sets. The results suggest that the stability of the different clusters depends on interrelated factors including the rings formed by intermolecular hydrogen bonds and the proton affinity (PA) or acidity of the atoms forming the intermolecular hydrogen bonds with the water molecules. A comparison across methods indicates that the results follow similar trends with different methods.

B-DNA Structure and Stability as Function of Nucleic Acid Composition: Dispersion-Corrected DFT Study of Dinucleoside Monophosphate Single and Double Strands

We have computationally investigated the structure and stability of all 16 combinations of two out of the four natural DNA bases A, T, G and C in a di-2′-deoxyribonucleoside-monophosphate model DNA strand as well as in 10 double-strand model complexes thereof, using dispersion-corrected density functional theory (DFT-D). Optimized geometries with B-DNA conformation were obtained through the inclusion of implicit water solvent and, in the DNA models, of sodium counterions, to neutralize the negative charge of the phosphate groups. The results obtained allowed us to compare the relative stability of isomeric single and double strands. Moreover, the energy of the Watson–Crick pairing of complementary single strands to form double-helical structures was calculated. The latter furnished the following increasing stability trend of the double-helix formation energy: d(TpA)₂ <d(CpA)₂ <d(ApT)₂ <d(ApA)₂ <d(GpT)₂ <d(GpA)₂ <d(ApG)₂ <d(CpG)₂ <d(GpG)₂ <d(GpC)₂, where the energy differences between the last four dimers, d(ApG)₂, d(CpG)₂, d(GpG)₂ and d(GpC)₂, is within 4.0 kcal mol⁻¹, and the energy between the most and the least stable isomers is 13.4 kcal mol⁻¹. This trend shows that the formation energy essentially increases with the number of hydrogen bonds per base pair, that is two between A and T and three between G and C. Superimposed on this main trend are more subtle effects that depend on the order in which bases occur within a strand from the 5’- to the 3’-end.

Lattice model of oligonucleotide hybridization in solution. I. Model and thermodynamics

A coarse-grained lattice model of DNA oligonucleotides is proposed to investigate the general mechanisms by which single-stranded oligonucleotides hybridize to their complementary strands in solution. The model, based on a high-coordination cubic lattice, is simple enough to allow the direct simulation of DNA solutions, yet capturing how the fundamental thermodynamic processes are microscopically encoded in the nucleobase sequences. Physically relevant interactions are considered explicitly, such as interchain excluded volume, anisotropic base-pairing and base-stacking, and single-stranded bending rigidity. The model is studied in detail by a specially adapted Monte Carlo simulation method, based on parallel tempering and biased trials, which is designed to overcome the entropic and enthalpic barriers associated with the sampling of hybridization events of multiple single-stranded chains in solution. This methodology addresses both the configurational complexity of bringing together two complementary strands in a favorable orientation (entropic barrier) and the energetic penalty of breaking apart multiple associated bases in a double-stranded state (enthalpic barrier). For strands with sequences restricted to nonstaggering association and homogeneous pairing and stacking energies, base-pairing is found to dominate the hybridization over the translational and conformational entropy. For strands with sequence-dependent pairing corresponding to that of DNA, the complex dependence of the model's thermal stability on concentration, sequence, and degree of complementarity is shown to be qualitatively and quantitatively consistent both with experiment and with the predictions of statistical mechanical models.

Differential stabilization of adenine quartets by anions and cations

We have investigated the structures and stabilities of four different adenine quartets with alkali and halide ions in the gas phase and in water, using dispersion-corrected density functional theory at the BLYP-D/TZ2P level. First, we examine the empty quartets and how they interact with alkali cations and halide anions with formation of adenine quartet–ion complexes. Second, we examine the interaction in a stack, in which a planar adenine quartet interacts with a cation or anion in the periphery as well as in the center of the quartet. Interestingly, for the latter situation, we find that both cations and anions can stabilize a planar adenine quartet in a stack.