Halogen Bonds in Ligand-Protein Systems: Molecular Orbital Theory for Drug Design

Can Duvelisib and Eganelisib work for both cancer and COVID-19? Molecular-level insights from MD simulations and enhanced samplings

SARS-CoV-2 has caused severe illness and anxiety worldwide, evolving into more dreadful variants capable of evading the host's immunity. Cytokine storms, led by PI3Kγ, are common in cancer and SARS-CoV-2. Naturally, there is a yearning to see whether any drugs could alleviate cytokine storms for both. Upon investigation, we identified two anticancer drugs, Duvelisib and Eganelisib, that could also work against SARS-CoV-2. This report is the first to decipher their synergic therapeutic effectiveness against COVID-19 and cancer with molecular insights from atomistic simulations. In addition to PI3Kγ, these drugs exhibit specificity for the main protease among all SARS-CoV-2 targets, with significant negative binding free energies and small time-dependent conformational changes of the complexes. Complexation makes active sites and secondary structures highly mechanically stiff, with barely any deformation. Replica simulations estimated large pulling forces in enhanced sampling to dissociate the drugs from Mpro's active site. Furthermore, the radial distribution function (RDF) demonstrated that the therapeutic molecules were closest to the His41 and Cys145 catalytic dyad residues. Finally, analyses implied Duvelisib and Eganelisib as promising dual-purposed anti-COVID and anticancer drugs, potentially targeting Mpro and PI3Kγ to stop virus replication and cytokine storms concomitantly. We also distinguished hotspot residues imparting significant interactions.

Situating the phosphonated calixarene-cytochrome C association by molecular dynamics simulations

Protein-calixarenes binding plays an increasingly central role in many applications, spanning from molecular recognition to drug delivery strategies and protein inhibition. These ligands obey a specific bio-supramolecular chemistry, which can be revealed by computational approaches, such as molecular dynamics simulations. In this paper, we rely on all-atom, explicit-solvent molecular dynamics simulations to capture the electrostatically driven association of a phosphonated calix-[4]-arene with cytochome-C, which critically relies on surface-exposed paired lysines. Beyond two binding sites identified in direct agreement with the x-ray structure, the association has a larger structural impact on the protein dynamics. Then, our simulations allow a direct comparison to analogous calixarenes, namely, sulfonato, similarly reported as "molecular glue." Our work can contribute to a robust in silico predictive tool to assess binding sites for any given protein of interest for crystallization, with the specificity of a macromolecular cage whose endo/exo orientation plays a role in the binding.

Relation between Halogen Bond Strength and IR and NMR Spectroscopic Markers

The relationship between the strength of a halogen bond (XB) and various IR and NMR spectroscopic quantities is assessed through DFT calculations. Three different Lewis acids place a Br or I atom on a phenyl ring; each is paired with a collection of N and O bases of varying electron donor power. The weakest of the XBs display a C-X bond contraction coupled with a blue shift in the associated frequency, whereas the reverse trends occur for the stronger bonds. The best correlations with the XB interaction energy are observed with the NMR shielding of the C atom directly bonded to X and the coupling constants involving the C-X bond and the C-H/F bond that lies ortho to the X substituent, but these correlations are not accurate enough for the quantitative assessment of energy. These correlations tend to improve as the Lewis acid becomes more potent, which makes for a wider range of XB strengths.

Identification of defactinib derivatives targeting focal adhesion kinase using ensemble docking, molecular dynamics simulations and binding free energy calculations

Focal adhesion kinase (FAK) belongs to the nonreceptor tyrosine kinases, which selectively phosphorylate tyrosine residues on substrate proteins. FAK is associated with bladder, esophageal, gastric, neck, breast, ovarian and lung cancers. Thus, FAK has been considered as a potential target for tumor treatment. Currently, there are six adenosine triphosphate (ATP)-competitive FAK inhibitors tested in clinical trials but no approved inhibitors targeting FAK. Defactinib (VS-6063) is a second-generation FAK inhibitor with an IC50 of 0.6 nM. The binding model of VS-6063 with FAK may provide a reference model for developing new antitumor FAK-targeting drugs. In this study, the VS-6063/FAK binding model was constructed using ensemble docking and molecular dynamics simulations. Furthermore, the molecular mechanics/generalized Born (GB) surface area (MM/GBSA) method was employed to estimate the binding free energy between VS-6063 and FAK. The key residues involved in VS-6063/FAK binding were also determined using per-residue energy decomposition analysis. Based on the binding model, VS-6063 could be separated into seven regions to enhance its binding affinity with FAK. Meanwhile, 60 novel defactinib-based compounds were designed and verified using ensemble docking. Overall, the present study improves our understanding of the binding mechanism of human FAK with VS-6063 and provides new insights into future drug designs targeting FAK.Communicated by Ramaswamy H. Sarma.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Tetrel bonds involving a CF3 group participate in protein-drug recognition: a combined crystallographic and computational study

In this study, the ability of CF₃ groups to bind to the electron-rich side chains and backbone groups of proteins has been investigated by combining a Protein Data Bank (PDB) survey and ab initio quantum mechanics calculations. More precisely, an inspection of the PDB involving organic ligands containing a CF₃ group and electron-rich atoms (A = N, O and S) in the vicinity revealed 419 X-ray structures exhibiting CF₃⋯A tetrel bonds (TtBs). In a posterior stage, those hits that exhibited the most relevant features in terms of directionality and intermolecular distance were selected for theoretical calculations at the RI-MP2/def2-TZVPD level of theory. Also, Hammett's regression plots of several TtB complexes involving meta- and para-substituted benzene derivatives were computed to shed light on the substituent effects. Moreover, the TtBs were characterized through several state-of-the-art computational techniques, such as the Quantum Theory of Atoms in Molecules (QTAIM) and Noncovalent Interactions plot (NCIplot) methodologies. We believe that the results gathered from our study will be useful for rational drug design and biological communities as well as for further expanding the role of this interaction to biomedical applications.

Regium-π Bonds Involving Nucleobases: Theoretical Study and Biological Implications

In this study, we provide crystallographic (Protein Data Bank (PDB) inspection) and theoretical (RI-MP2/def2-TZVP//PBE0-D3/def2-SVP level of theory) evidence of the involvement of nucleobases in Regium-π bonds (RgBs). This noncovalent interaction involves an electrophilic site located on an element of group 11 (Cu, Ag, and Au) and an electron-rich species (lone pair, LP donor, or π-system). Concretely, an initial PDB search revealed several examples where RgBs were undertaken involving DNA bases and Cu(II), Ag(I), and Au(I/III) ions. While coordination positions (mainly at the N atoms of the base) are well known, the noncovalent binding force between these counterparts has been scarcely studied in the literature. In this regard, computational models shed light on the strength and directionality properties of the interaction, which was also further characterized from a charge-density perspective using Bader's "atoms in molecules" (AIM) theory, noncovalent interaction plot (NCIplot) visual index, and natural bonding orbital (NBO) analyses. As far as our knowledge extends, this is the first time that RgBs in metal-DNA complexes are systematically analyzed, and we believe the results might be useful for scientists working in the field of nucleic acid engineering and chemical biology as well as to increase the visibility of the interaction among the biological community.

Intermolecular Covalent Interactions: Nature and Directionality

Quantum chemical methods were employed to analyze the nature and the origin of the directionality of pnictogen (PnB), chalcogen (ChB), and halogen bonds (XB) in archetypal F_m Z⋅⋅⋅F^- complexes (Z=Pn, Ch, X), using relativistic density functional theory (DFT) at ZORA-M06/QZ4P. Quantitative Kohn-Sham MO and energy decomposition analyses (EDA) show that all these intermolecular interactions have in common that covalence, that is, HOMO-LUMO interactions, provide a crucial contribution to the bond energy, besides electrostatic attraction. Strikingly, all these bonds are directional (i.e., F-Z⋅⋅⋅F^- is approximately linear) despite, and not because of, the electrostatic interactions which, in fact, favor bending. This constitutes a breakdown of the σ-hole model. It was shown how the σ-hole model fails by neglecting both, the essential physics behind the electrostatic interaction and that behind the directionality of electron-rich intermolecular interactions. Our findings are general and extend to the neutral, weaker ClI⋅⋅⋅NH₃ , HClTe⋅⋅⋅NH₃ , and H₂ ClSb⋅⋅⋅NH₃ complexes.

Understanding the impact of halogen functional group (Br, Cl, F, OH) in amprenavir ligand of the HIV protease

We focused our attention towards the most dreadful disease that threatens the mankind of 20th century - Acquired immunodeficiency syndrome (AIDS), caused through the human immunodeficiency virus (HIV) and a sexually transmitted infection (STI). In this study, our foremost interest was to identify the potency and stability of HIV ligand- Amprenavir (APV) and its modelled functional group (Br, Cl, F, CF3, CH3, NH2) ligands through halogen and hydrogen bond contact, which will have a clear portrait on the structure activity of protein ligand interactions. This will assist chemist in synthesizing novel APV ligands, which are expected to inhibit the activity of HIV-1 protease enzyme. The binding strength of Amprenavir ligand with interacting hinge region amino acid side chains: Isoleucine (ILE 147, 150, 184), Valine (VAL 82), Alanine (ALA 28), Aspartic acid (25, 30, 125, 130) and Glycine (GLY 127, 149) were understood through interaction energy calculations at HF, B3LYP, M052X, MP2 level of theories for different basis set (6-311 G**, LANL2DZ). The present work will reveal an understandable picture about the halogen and hydrogen bond interaction that grip the contact of ligand and amino acids in the hinge region. Overall the Halogen atom (Br, Cl, F) functional groups improved the binding strength of APV in HIV protease; which provide a new novel path for the functional group preference on the ligand that enclose perfectly with the amino acid in the hinge region.Communicated by Ramaswamy H. Sarma.

Charge Assisted S/Se Chalcogen Bonds in SAM Riboswitches: A Combined PDB and ab Initio Study

In this study, we provide experimental (Protein Data Bank (PDB) inspection) and theoretical (RI-MP2/def2-TZVP level of theory) evidence of the involvement of charge assisted chalcogen bonding (ChB) interactions in the recognition and folding mechanisms of S-adenosylmethionine (SAM) riboswitches. Concretely, an initial PDB search revealed several examples where ChBs between S-adenosyl methionine (SAM)/adenosyl selenomethionine (EEM) molecules and uracil (U) bases belonging to RNA take place. While these interactions are usually described as a merely Coulombic attraction between the positively charged S/Se group and RNA, theoretical calculations indicated that the σ holes of S and Se are involved. Moreover, computational models shed light on the strength and directionality properties of the interaction, which was also further characterized from a charge-density perspective using Bader’s “Atoms in Molecules” (AIM) theory, Non-Covalent Interaction plot (NCIplot) visual index, and Natural Bonding Orbital (NBO) analyses. As far as our knowledge extends, this is the first time that ChBs in SAM–RNA complexes have been systematically analyzed, and we believe the results might be useful for scientists working in the field of RNA engineering and chemical biology as well as to increase the visibility of the interaction among the biological community.

Spodium Bonds in Biological Systems: Expanding the Role of Zn in Protein Structure and Function

Understanding the structural and functional implications of metal ions is of pivotal significance to chemical biology. Herein, we report first time the evidence of spodium bonds (SpB's, an attractive noncovalent force involving elements from group 12 and electron-rich species) in tetrahedral Zn-binding sites. Through a combined crystallographic (PDB analysis) and computational (ab initio calculations) study, we demonstrate that Zn SpB's are abundant and might be involved in protein structure and enzyme inhibition.

Nature of Alkali- and Coinage-Metal Bonds versus Hydrogen Bonds

We have quantum chemically studied the structure and nature of alkali- and coinage-metal bonds (M-bonds) versus that of hydrogen bonds between A-M and B- in archetypal [A-M⋅⋅⋅B]- model systems (A, B=F, Cl and M=H, Li, Na, Cu, Ag, Au), using relativistic density functional theory at ZORA-BP86-D3/TZ2P. We find that coinage-metal bonds are stronger than alkali-metal bonds which are stronger than the corresponding hydrogen bonds. Our main purpose is to understand how and why the structure, stability and nature of such bonds are affected if the monovalent central atom H of hydrogen bonds is replaced by an isoelectronic alkali- or coinage-metal atom. To this end, we have analyzed the bonds between A-M and B- using the activation strain model, quantitative Kohn-Sham molecular orbital (MO) theory, energy decomposition analysis (EDA), and Voronoi deformation density (VDD) analysis of the charge distribution.

Insight into the binding mode of HIF-2 agonists through molecular dynamic simulations and biological validation

Hypoxia-inducible factor-2 (HIF-2), a heterodimeric transcriptional protein consisting of HIF-2α and aryl hydrocarbon receptor nuclear translocator (ARNT) subunits, has a broad transcriptional profile that plays a vital role in human oxygen metabolism. M1001, a HIF-2 agonist identified by high-throughput screening (HTS), is capable of altering the conformation of Tyr281 of the HIF-2α PAS-B domain and enhancing the affinity of HIF-2α and ARNT for transcriptional activation. M1002, an analog of M1001, shows improved efficacy than M1001. However, the cocrystal structure of M1001 and HIF-2 has some defects in revealing the agonist binding mode due to the relatively low resolution, while the binding mode of M1002 remained unexplored. To in-depth understand agonist binding profiles, herein, the molecular dynamic (MD) simulations was applied to construct a stable agonist-protein model, and a possible binding mode was proposed through the analysis of the binding free energy and hydrogen bonding of the simulation results. Nine compounds were then synthesized and evaluated to verify the proposed binding mode. Among them, compound 10 manifested improved agonistic activity and reduced toxicity compared to M1002. This study provides deep insight into the binding mode of such HIF-2 agonists, which would be useful for designing novel agonists for HIF-2.

The C2N surface as a highly selective sensor for the detection of nitrogen iodide from a mixture of NX3 (X = Cl, Br, I) explosives

Explosives are quite toxic and destructive; therefore, it is necessary to not only detect them but also remove them. The adsorption behavior of NX₃ analytes (NCl₃, NBr₃ and NI₃) over the microporous C₂N surface was evaluated by DFT calculations. The nature of interactions between NX₃ and C₂N was characterized by adsorption energy, NCI, QTAIM, SAPT0, NBO, EDD and FMO analysis. The interaction energies of NX₃ with C₂N are in the range of −10.85 to −16.31 kcal mol⁻¹ and follow the order of NCl₃@C₂N > NBr₃@C₂N > NI₃@C₂N, respectively. The 3D isosurfaces and 2D-RGD graph of NCI analysis qualitatively confirmed the existence of halogen bonding interactions among the studied systems. Halogen bonding was quantified by SAPT0 component energy analysis. The SAPT0 results revealed that ΔE_disp (56.75%) is the dominant contributor towards interaction energy, whereas contributions from ΔE_elst and ΔE_ind are 29.41% and 14.34%, respectively. The QTAIM analysis also confirmed the presence of halogen bonding between atoms of NX₃ and C₂N surface. EDD analysis also validated NCI, QTAIM and NBO analysis. FMO analysis revealed that the adsorption of NI₃ on the C₂N surface caused the highest change in the E_HOMO–LUMO gap (from 5.71 to 4.15 eV), and resulted in high sensitivity and selectivity of the C₂N surface towards NI₃, as compared to other analytes. It is worth mentioning that in all complexes, a significant difference in the E_HOMO–LUMO gap was seen when electronic transitions occurred from the analyte to the C₂N surface.

Explosives are quite toxic and destructive; therefore, it is necessary to not only detect them but also remove them.