Non-covalent interactions in receptor-ligand complexes. A study based on the electron charge density

Exploring protein-ligand binding affinity prediction with electron density-based geometric deep learning

Rational structure-based drug design relies on accurate predictions of protein-ligand binding affinity from structural molecular information. Although deep learning-based methods for predicting binding affinity have shown promise in computational drug design, certain approaches have faced criticism for their potential to inadequately capture the fundamental physical interactions between ligands and their macromolecular targets or for being susceptible to dataset biases. Herein, we propose to include bond-critical points based on the electron density of a protein-ligand complex as a fundamental physical representation of protein-ligand interactions. Employing a geometric deep learning model, we explore the usefulness of these bond-critical points to predict absolute binding affinities of protein-ligand complexes, benchmark model performance against existing methods, and provide a critical analysis of this new approach. The models achieved root-mean-squared errors of 1.4-1.8 log units on the PDBbind dataset, and 1.0-1.7 log units on the PDE10A dataset, not indicating significant advantages over benchmark methods, and thus rendering the utility of electron density for deep learning models context-dependent. The relationship between intermolecular electron density and corresponding binding affinity was analyzed, and Pearson correlation coefficients r > 0.7 were obtained for several macromolecular targets.

How procyanidin C1 sticks to collagen: The role of proline rings

Molecular interactions between proteins and polyphenols are responsible for many natural phenomena like colloidal turbidity, astringency, denaturation of enzymes and leather tanning. Although these phenomena are well known, there are open questions about the specific interactions involved in the complexation process. In this work, Molecular Dynamic (MD) simulations and the topology of the electron density analysis were used to study the interactions between the flavonoid procyanidin C1 and a collagen fragment solvated in water. Root mean square deviation; root mean square fluctuation and hydrogen bonds occupancy were examined after 50 ns. The interactions were also analyzed by means of the quantum theory of atoms in molecules. Our results show that the main interactions are hydrogen bonds between -OH groups of the polyphenol and CO groups of the peptide bond. Stacking interactions between proline rings and phenol rings, that is CH⋯π hydrogen bonds, also stabilize the dynamic structure of the complex.

Design, synthesis, biological evaluation and molecular modelling of substituted pyrrolo[2,1-a]isoquinolinone derivatives: discovery of potent inhibitors of AChE and BChE

Study of the molecular interactions in L–R complexes of acetyl- and butyryl-cholinesterase using MD/QTAIM calculations for designing new potent cholinesterase inhibitors.

Combined MD/QTAIM techniques to evaluate ligand-receptor interactions. Scope and limitations

In medicinal chemistry, it is extremely important to evaluate, as accurately as possible, the molecular interactions involved in the formation of different ligand-receptor (L-R) complexes. Evaluating the different molecular interactions by quantum mechanics calculations is not a simple task, since formation of an L-R complex is a dynamic process. In this case, the use of combined techniques of molecular dynamics (MD) and quantum calculations is one the best possible approaches. In this work we report a comparative study using combined MD and QTAIM (Quantum Theory of Atoms In Molecules) calculations for five biological systems with different levels of structural complexity. We have studied Acetylcholinesterase (AChE), D2 Dopamine Receptor (D2DR), beta Secretase (BACE1), Dihydrofolate Reductase (DHFR) and Sphingosine Kinase 1 (SphK1). In these molecular targets, we have analyzed different ligands with diverse structural characteristics. The inhibitory activities of most of them have been previously measured in our laboratory. Our results indicate that QTAIM calculations can be extremely useful for in silico studies. It is possible to obtain very accurate information about the strength of the molecular interactions that stabilize the formation of the different L-R complexes. Better correlations can be obtained between theoretical and experimental data by using QTAIM calculations, allowing us to discriminate among ligands with similar affinities. QTAIM analysis gives fairly accurate information for weak interactions which are not well described by MD simulations. QTAIM study also allowed us to evaluate and determine which parts of the ligand need to be modified in order to increase its interactions with the molecular target. In this study we have discussed the importance of combined MD/QTAIM calculations for this type of simulations, showing their scopes and limitations.

Arylaminopropanone Derivatives as Potential Cholinesterase Inhibitors: Synthesis, Docking Study and Biological Evaluation

Neurodegenerative diseases in which the decrease of the acetylcholine is observed are growing worldwide. In the present study, a series of new arylaminopropanone derivatives with N-phenylcarbamate moiety (1–16) were prepared as potential acetylcholinesterase and butyrylcholinesterase inhibitors. In vitro enzyme assays were performed; the results are expressed as a percentage of inhibition and the IC₅₀ values. The inhibitory activities were compared with reference drugs galantamine and rivastigmine showing piperidine derivatives (1–3) as the most potent. A possible mechanism of action for these compounds was determined from a molecular modelling study by using combined techniques of docking, molecular dynamics simulations and quantum mechanics calculations.

A quantum chemical approach representing a new perspective concerning agonist and antagonist drugs in the context of schizophrenia and Parkinson's disease

Schizophrenia and Parkinson's disease can be controlled with dopamine antagonists and agonists. In order to improve the understanding of the reaction mechanism of these drugs, in this investigation we present a quantum chemical study of 20 antagonists and 10 agonists. Electron donor acceptor capacity and global hardness are analyzed using Density Functional Theory calculations. Following this theoretical approach, we provide new insights into the intrinsic response of these chemical species. In summary, antagonists generally prove to be better electron acceptors and worse electron donors than dopamine, whereas agonists present an electron donor-acceptor capacity similar to that of dopamine. The chemical hardness is a descriptor that captures the resistance of a chemical compound to change its number of electrons. Within this model, harder molecules are less polarizable and more stable systems. Our results show that the global hardness is similar for dopamine and agonists whilst antagonists present smaller values. Following the Hard and Soft Acid and Bases principle, it is possible to conclude that dopamine and agonists are hard bases while antagonists are soft acids, and this can be related to their activity. From the electronic point of view, we have evolved a new perspective for the classification of agonist and antagonist, which may help to analyze future results of chemical interactions triggered by these drugs.

Combining Charge Density Analysis with Machine Learning Tools To Investigate the Cruzain Inhibition Mechanism

Trypanosoma cruzi, a flagellate protozoan parasite, is responsible for Chagas disease. The parasite major cysteine protease, cruzain (Cz), plays a vital role at every stage of its life cycle and the active-site region of the enzyme, similar to those of other members of the papain superfamily, is well characterized. Taking advantage of structural information available in public databases about Cz bound to known covalent inhibitors, along with their corresponding activity annotations, in this work, we performed a deep analysis of the molecular interactions at the Cz binding cleft, in order to investigate the enzyme inhibition mechanism. Our toolbox for performing this study consisted of the charge density topological analysis of the complexes to extract the molecular interactions and machine learning classification models to relate the interactions with biological activity. More precisely, such a combination was useful for the classification of molecular interactions as "active-like" or "inactive-like" according to whether they are prevalent in the most active or less active complexes, respectively. Further analysis of interactions with the help of unsupervised learning tools also allowed the understanding of how these interactions come into play together to trigger the enzyme into a particular conformational state. Most active inhibitors induce some conformational changes within the enzyme that lead to an overall better fit of the inhibitor into the binding cleft. Curiously, some of these conformational changes can be considered as a hallmark of the substrate recognition event, which means that most active inhibitors are likely recognized by the enzyme as if they were its own substrate so that the catalytic machinery is arranged as if it is about to break the substrate scissile bond. Overall, these results contribute to a better understanding of the enzyme inhibition mechanism. Moreover, the information about main interactions extracted through this work is already being used in our lab to guide docking solutions in ongoing prospective virtual screening campaigns to search for novel noncovalent cruzain inhibitors.

Searching new structural scaffolds for BRAF inhibitors. An integrative study using theoretical and experimental techniques

The identification of the V600E activating mutation in the protein kinase BRAF in around 50% of melanoma patients has driven the development of highly potent small inhibitors (BRAFi) of the mutated protein. To date, Dabrafenib and Vemurafenib, two specific BRAFi, have been clinically approved for the treatment of metastatic melanoma. Unfortunately, after the initial response, tumors become resistant and patients develop a progressive and lethal disease, making imperative the development of new therapeutic options. The main objective of this work was to find new BRAF inhibitors with different structural scaffolds than those of the known inhibitors. Our study was carried out in different stages; in the first step we performed a virtual screening that allowed us to identify potential new inhibitors. In the second step, we synthesized and tested the inhibitory activity of the novel compounds founded. Finally, we conducted a molecular modelling study that allowed us to understand interactions at the molecular level that stabilize the formation of the different molecular complexes. Our theoretical and experimental study allowed the identification of four new structural scaffolds, which could be used as starting structures for the design and development of new inhibitors of BRAF. Our experimental data indicate that the most active compounds reduced significantly ERK½ phosphorylation, a measure of BRAF inhibition, and cell viability. Thus, from our theoretical and experimental results, we propose new substituted hydroxynaphthalenecarboxamides, N-(hetero)aryl-piperazinylhydroxyalkylphenylcarbamates, substituted piperazinylethanols and substituted piperazinylpropandiols as initial structures for the development of new inhibitors for BRAF. Moreover, by performing QTAIM analysis, we are able to describe in detail the molecular interactions that stabilize the different Ligand-Receptor complexes. Such analysis indicates which portion of the different molecules must be changed in order to obtain an increase in the binding affinity of these new ligands.

Synthesis, Analysis, Cholinesterase-Inhibiting Activity and Molecular Modelling Studies of 3-(Dialkylamino)-2-hydroxypropyl 4-[(Alkoxy-carbonyl)amino]benzoates and Their Quaternary Ammonium Salts

Tertiary amines 3-(dialkylamino)-2-hydroxypropyl 4-[(alkoxycarbonyl)amino]benzoates and their quaternary ammonium salts were synthesized. The final step of synthesis of quaternary ammonium salts was carried out by microwave-assisted synthesis. Software-calculated data provided the background needed to compare fifteen new resulting compounds by their physicochemical properties. The acid dissociation constant (pK_a) and lipophilicity index (log P) of tertiary amines were determined; while quaternary ammonium salts were characterized by software-calculated lipophilicity index and surface tension. Biological evaluation aimed at testing acetylcholinesterase and butyrylcholinesterase-inhibiting activity of synthesized compounds. A possible mechanism of action of these compounds was determined by molecular modelling study using combined techniques of docking; molecular dynamics simulations and quantum mechanics calculations.

An integrative study to identify novel scaffolds for sphingosine kinase 1 inhibitors

Sphingosine kinase 1 (SphK1), the enzyme that produces the bioactive sphingolipid metabolite, sphingosine-1-phosphate, is a promising new molecular target for therapeutic intervention in cancer and inflammatory diseases. In view of its importance, the main objective of this work was to find new and more potent inhibitors for this enzyme possessing different structural scaffolds than those of the known inhibitors. Our theoretical and experimental study has allowed us to identify two new structural scaffolds (three new compounds), which could be used as starting structures for the design and then the development of new inhibitors of SphK1. Our study was carried out in different steps: virtual screening, synthesis, bioassays and molecular modelling. From our results, we propose a new dihydrobenzo[b]pyrimido[5,4-f]azepine and two alkyl{3-/4-[1-hydroxy-2-(4-arylpiperazin-1-yl)ethyl]phenyl}carbamates as initial structures for the development of new inhibitors. In addition, our molecular modelling study using QTAIM calculations, allowed us to describe in detail the molecular interactions that stabilize the different Ligand-Receptor complexes. Such analyses indicate that the cationic head of the different compounds must be refined in order to obtain an increase in the binding affinity of these ligands.

Direct Determination of Site-Specific Noncovalent Interaction Strengths of Proteins from NMR-Derived Fast Side Chain Motional Parameters

A novel approach to accurately determine residue-specific noncovalent interaction strengths (ξ) of proteins from NMR-measured fast side chain motional parameters (O_axis²) is presented. By probing the environmental sensitivity of side chain conformational energy surfaces of individual residues of a diverse set of proteins, the microscopic connections between ξ, O_axis², conformational entropy (S_conf), conformational barriers, and rotamer stabilities established here are found to be universal among proteins. The results reveal that side chain flexibility and conformational entropy of each residue decrease with increasing ξ and that for each residue type there exists a critical range of ξ, determined primarily by the mean side chain conformational barriers, within which flexibility of any residue can be reversibly tuned from highly flexible (with O_axis² ∼ 0) to highly restricted (with O_axis² ∼ 1) by increasing ξ by ∼3 kcal/mol. Beyond this critical range of ξ, both side chain flexibility and conformational entropy are insensitive to ξ. The interrelationships between conformational dynamics, conformational entropy, and noncovalent interactions of protein side chains established here open up new avenues to probe perturbation-induced (for example, ligand-binding, temperature, pressure) changes in fast side chain dynamics and thermodynamics of proteins by comparing their conformational energy surfaces in the native and perturbed states.

Theoretical investigation of gas-phase molecular complex formation between 2-hydroxy thiophenol and a water molecule

A systematic study of the interaction energies and hydrogen bonding interaction of a gas-phase molecular complex between 2-hydroxy thiophenol and a water molecule.

Halogen bonding in biological context: a computational study of D2 dopamine receptor

In this work, Halogen Bond (X‐bond) interactions formed by halogenated ligands (LX) at the Dopamine Receptor D2 (DRD2) binding pocket were studied by Molecular Dynamics (MD) and charge density analysis. The X‐bonds were contrasted with the Hydrogen Bond (H‐bond) interactions established by hydroxylated analogs (LOH, where X was replaced by OH). The ligands for this study were extracted from a dataset of compounds deposited in ZINC database that were active in binding assays to DRD2. This dataset was subjected to the filtering rules by employing cheminformatics tools to find the LX/LOH pairs that were then submitted to MD simulations. A homology model of DRD2 was employed for the simulations because no crystal structure is yet available for the receptor. To mimic the positive cap (σ‐hole) on the halogen atom, a massless, positive charged extra‐point was introduced in the force field. An analysis of the charge density (QTAIM) was performed on reduced models of simulated complexes to explain their binding differences. Results show that the halogen atom tends to form X‐bond with protein backbone oxygen atom. Two out of the four halogenated ligands studied form a specific X‐bond with the carbonyl oxygen of Ser193. This specific X‐bond decreases the inherent propensity of transmembrane 5 to unfolding. These results suggest a possible role of the X‐bond as a protein secondary structure modulator because of the ability of the halogen to interact with the protein backbone. Copyright © 2016 John Wiley & Sons, Ltd.

Using Density Functional Theory To Study Neutral and Ionized Stacked Thymine Dimers

Stacking interactions in thymine dimers are studied with density functional theory. According to our calculations, six dimers of comparable stability can be prepared at low temperature, but dimerization is impossible at room temperature due to the large entropy contribution that accompanies it. Analysis of vibrational anharmonic coupling terms shows that each of the dimers exhibits distinct vibrational dynamics. Properties of electron density in the intermolecular region are used to analyze neutral stacked species and their ionized forms. Bond paths and critical points in the intermolecular region are identified, but a simple relationship between binding energy and total electron density in the intermolecular critical points could not be found due to an uneven electron distribution in the binding region. The reduced density gradient was confirmed to be a useful tool for analysis of weak stacking interactions. Those interactions also affect vertical and adiabatic ionization energies, which are computed to be slightly lower for the dimers compared to the monomer.

Evolution of the hydrogen-bonding motif in the melamine-cyanuric acid co-crystal: a topological study

The melamine (M)/cyanuric acid (CA) supramolecular system is perhaps one of the most exploited in the field of self-assembly because of the high complementarity of the components. However, it is necessary to investigate further the factors involved in the assembly process. In this study, we analyzed a set of 13 M n /CA m clusters (with n , m = 1, 2, 3), taken from crystallographic data, to characterize the nature of the hydrogen bonds involved in the self-assembly of these components as well as to provide greater understanding of the phenomenon. The calculations were performed at the B3LYP/6-311++G(d,p) and ω-B97XD (single point) levels of theory, and the interactions were analyzed within the framework of the quantum theory of atoms in molecules and by means of molecular electrostatic potential maps. Our results show that the stablest structure is the rosette-type motif and the aggregation mechanism is governed by a combination of cooperative and anticooperative effects. Our topological results explain the polymorphism in the self-assembly of coadsorbed monolayers of M and CA. Graphical abstract The aggregation steps of the melamine-cyanuric co-crystal is driven by a hydrogen-bonded network which is governed by a complex combination of cooperative and anticooperative effects.

Dopaminergic isoquinolines with hexahydrocyclopenta[ij]-isoquinolines as D2-like selective ligands

Dopamine receptors (DR) ligands are potential drug candidates for treating neurological disorders including schizophrenia or Parkinson's disease. Three series of isoquinolines: (E)-1-styryl-1,2,3,4-tetrahydroisoquinolines (series 1), 7-phenyl-1,2,3,7,8,8a-hexahydrocyclopenta[ij]-IQs (HCPIQs) (series 2) and (E)-1-(prop-1-en-1-yl)-1,2,3,4- tetrahydroisoquinolines (series 3), were prepared to determine their affinity for both D1 and D2-like DR. The effect of different substituents on the nitrogen atom (methyl or allyl), the dioxygenated function (methoxyl or catechol), the substituent at the β-position of the THIQ skeleton, and the presence or absence of the cyclopentane motif, were studied. We observed that the most active compounds in the three series (2c, 2e, 3a, 3c, 3e, 5c and 5e) possessed a high affinity for D2-like DR and these remarkable features: a catechol group in the IQ-ring and the N-substitution (methyl or allyl). The series showed the following trend to D2-RD affinity: HCPIQs > 1-styryl > 1-propenyl. Therefore, the substituent at the β-position of the THIQ and the cyclopentane ring also modulated this affinity. Among these dopaminergic isoquinolines, HCPIQs stood out for unexpected selectivity to D2-DR since the Ki D1/D2 ratio reached values of 2465, 1010 and 382 for compounds 3a, 3c and 3e, respectively. None of the most active THIQs in D2 DR displayed relevant cytotoxicity in human neutrophils and HUVEC. Finally, and in agreement with the experimental data, molecular modeling studies on DRs of the most characteristic ligands of the three series revealed stronger molecular interactions with D2 DR than with D1 DR, which further supports to the encountered enhanced selectivity to D2 DR.