Quantum Mechanical Assessment of Protein–Ligand Hydrogen Bond Strength Patterns: Insights from Semiempirical Tight-Binding and Local Vibrational Mode Theory

Abstract

Hydrogen bonds (HB)s are the most abundant motifs in biological systems. They play a key role in determining protein–ligand binding affinity and selectivity. We designed two pharmaceutically beneficial HB databases, database A including ca. 12,000 protein–ligand complexes with ca. 22,000 HBs and their geometries, and database B including ca. 400 protein–ligand complexes with ca. 2200 HBs, their geometries, and bond strengths determined via our local vibrational mode analysis. We identified seven major HB patterns, which can be utilized as a de novo QSAR model to predict the binding affinity for a specific protein–ligand complex. Glycine was reported as the most abundant amino acid residue in both donor and acceptor profiles, and N–H⋯O was the most frequent HB type found in database A. HBs were preferred to be in the linear range, and linear HBs were identified as the strongest. HBs with HB angles in the range of 100–110, typically forming intramolecular five-membered ring structures, showed good hydrophobic properties and membrane permeability. Utilizing database B, we found a generalized Badger’s relationship for more than 2200 protein–ligand HBs. In addition, the strength and occurrence maps between each amino acid residue and ligand functional groups open an attractive possibility for a novel drug-design approach and for determining drug selectivity and affinity, and they can also serve as an important tool for the hit-to-lead process.

The Discovery of Selective Protein Arginine Methyltransferase 5 Inhibitors in the Management of β-Thalassemia through Computational Methods

“…Investigating the molecular interactions between a ligand-bound protein complex is crucial in assessing complex stability during docking and MDS studies. Stability is influenced mainly by the formation of hydrogen bonds and hydrophobic interactions [43][44][45][46],…”

Virtual Screening Using a Ligand-based Pharmacophore Model from Ashitaba (Angelica keiskei K.) Isolates and Molecular Docking to Obtained New Candidates as -Glucosidase Inhibitors

“…33 Electrostatic interactions, such as ionic bonds, play a key role in determining protein-ligand binding affinity and selectivity. 36 Van der Waals interactions, such as π/π interactions, also play a critical role in ligand-receptor interactions. 35 Acarbose, a known -glucosidase inhibitor, exhibits a binding energy of -1.25 kcal/mol in its interaction with specific residues within the target protein.…”

Chemical bonding in Uranium‐based materials: A local vibrational mode case study of Cs 2UO 2Cl 4 and UCl 4 crystals

“…29 Since its inception, LVM theory has found successful applications in characterizing a wide range of systems, providing valuable insights into covalent bonding as well as non-covalent interactions involving hydrogen, chalcogen, pnicogen, and unusual chemical bonding phenomena. 21 LVM analysis has also led to a new way of analyzing vibrational spectra, enabling novel investigations such as the examination of pK a probes, 30 gaining insights into vibrational Stark effect probes, 31 analyzing characteristic vibrational coupling in nucleobases and Watson-Crick base pairs of DNA, 32 assessing bond strength in biological systems through QM/MM, 33 evaluating underlying properties of lanthanide compounds, 34 or studying protein-ligand hydrogen bond strength patterns, 35 just to name a few.…”