Hydrogen Bonding and Molecular Surface Shape Complementarity as a Basis for Protein Docking

doi:10.1006/jmbi.1996.0634

Journal of Molecular Biology

Volume 264, Issue 1, 22 November 1996, Pages 199-210

https://doi.org/10.1006/jmbi.1996.0634 Get rights and content

Abstract

A geometric docking algorithm based upon correlation analysis for quantification of geometric complementarity between protein molecular surfaces in close interfacial contact has been developed by a detailed optimization of the conformational search of the algorithm. In order to reduce the entire conformation space search required by the method a physico-chemical pre-filter of conformation space has been developed based upon thea prioriassumption that two or more intermolecular hydrogen bonds are intrinsic to the mechanism of binding within protein complexes. Donor sites are defined spatially and directionally by the positions of explicitly calculated donor hydrogen atoms, and the vector space within a defined range about the donor atom-hydrogen atom bond vector. Acceptor sites are represented spatially and directionally by the van der Waals molecular surface points having normal vectors within a predefined range of vector space about the acceptor atom covalent bond vector(s). Geometric conditions necessary for the simultaneous hydrogen bonding interaction between both sites of functionally congruent hydrogen bonding site pairs, located on the individual proteins, are then tested on the basis of a transformation invariant parameterization of the site pair spatial and directional properties. Sterically acceptable conformations defined by interaction of functionally, spatially, and directionally compatable site pairs are then refined to a maximum contact of complementary contact surfaces using the simplex method for the angular search and correlation techniques for the translational search.

The utility of the spatial and directional properties of hydrogen bonding donor and acceptor sites for the identification of candidate docking conformations is demonstrated by the reliable preliminary reduction of conformation space, the improved geometric ranking of the minimum RMS conformations of some complexes and the overall reduction of CPU time obtained.

References (0)

Cited by (97)

Synthesis of novel Bis-1,2,4-Triazolo[3,4-b][1,3,4]Thiadiazines from natural camphoric acid as potential anti-candidal agents
2024, Results in Chemistry
Candida species have long been attributed to various diseases like candidiasis and systemic diseases and exacerbate the symptoms of immunocompromised patients. Candida species have enzymes that could function as drug targets to decrease their pathogenicity and eradicate the fungi. This research aimed to investigate the potency of new bis-triazolothiadiazine derivatives contained in inhibiting important enzymes of C. albicans as an example, through molecular docking simulation. Thus, a novel series of bis-triazolo[3,4-b][1,3,4]thiadiazines were designed and prepared via the reaction of the most versatile, hitherto unreported 5,5′-(1,2,2-trimethylcyclopentane-1,3-diyl)bis(4-amino-2,4-dihydro-3H-1,2,4-triazole-3-thione) with the appropriate hydrazonoyl halides and phenacyl bromides. Various spectroscopic techniques were used to identify the structures of the synthesized derivatives. The synthesized derivatives were tested against different species of Candida spp. The most effective compound was 6f followed by 6b, 6d, and 6e, where the inhibition zone ranged from 45 mm to 38 mm. By using molecular docking, which highlighted the important interactions with the amino acid residues Lys57, Leu77, Glu116, Gly114, Phe36, Thr58, and Glu32 at the point of binding, it was possible to determine the binding interactions of the produced derivatives to the fluconazole target fungi. The binding interaction energy was discovered to be −6.494 kcal/mol for the fungi candida albicans (PDB ID: 1IA2). The derivatives 6f and 6d, which demonstrated the highest efficacy, displayed significant conserved interactions with the amino acid residues at the binding site of the fluconazole fungi, in conjunction with the PDB co-crystal ligand 1IA2. The study’s results also revealed that the dG scores of the novel bis-triazolo-thiadiazines 6b-f. The study shows good docking scores with acceptable binding interactions. Finally, molecular docking studies revealed the lowest binding activity of derivatives 6e and 6c with the target fungi.
Probing antiviral drugs against SARS-CoV-2 through virus-drug association prediction based on the KATZ method
2020, Genomics
Citation Excerpt :
Remdesivir, oseltamivir, and zanamivir were predicted to be the top 3 small molecules associated with SARS-CoV-2. Molecular docking is a theoretical simulation method to compute the binding mode and affinity through electrostatic interactions [16], hydrogen bond interactions [17], Van der Waals forces interactions [18] and hydrophobic interactions [19] between molecules (for example, ligands and receptors). Energy matching is the basis of the stable combination between molecules.
It is urgent to find an effective antiviral drug against SARS-CoV-2. In this study, 96 virus-drug associations (VDAs) from 12 viruses including SARS-CoV-2 and similar viruses and 78 small molecules are selected. Complete genomic sequence similarity of viruses and chemical structure similarity of drugs are then computed. A KATZ-based VDA prediction method (VDA-KATZ) is developed to infer possible drugs associated with SARS-CoV-2. VDA-KATZ obtained the best AUCs of 0.8803 when the walking length is 2. The predicted top 3 antiviral drugs against SARS-CoV-2 are remdesivir, oseltamivir, and zanamivir. Molecular docking is conducted between the predicted top 10 drugs and the virus spike protein/human ACE2. The results showed that the above 3 chemical agents have higher molecular binding energies with ACE2. For the first time, we found that zidovudine may be effective clues of treatment of COVID-19. We hope that our predicted drugs could help to prevent the spreading of COVID.
Conceptual design of hybrid PCSK9 lead inhibitors against coronary artery disease
2019, Biocatalysis and Agricultural Biotechnology
Citation Excerpt :
From Table 4, among 48 hit molecules, 28 compounds fail to meet some of the criteria such as log P, MW and hepatotoxicity profile. In particular, the berberine moieties reduces the hydrogen-bonding capacity and molecular flexibility and increases the lipophilicity of the molecule compared to the other molecules (Meyer et al., 1996). The human intestinal absorption (HIA) and BBB penetration model of ADMET are defined by 95% and 99% confidence ellipses in the ADMET_PSA_2D plot (Fig. 13) (Egan and Lauri, 2002).
The enzyme, Proprotein Convertase Subtilisin Kexin type-9 (PCSK 9) inhibition serves as a significant target to prevent the risk of Coronary Artery Disease. Hence, this work is focused on the exploration of novel PCSK 9 lead inhibitors through In silico approach. In abidance to the study, BREED algorithm has been taken to design novel hybridised ligands. Subsequently, various in silico screening studies were performed in Biovia Discovery Studio 2017. Novel ligand generation protocol generated 595 hybridised ligands from the 51 parent scaffolds. These ligands were mapped and screened by ludi interactions of the active site PCSK 9 and were reduced to 245 hit compounds. The docking studies showed that the berberine skeleton contained molecule possessed CDOCKER energy (−13.5582 Kcal/mol) and formed strong interactions with Arg⁹³, Arg⁹⁷ and His¹³⁷ residues. The developed pharmacophore model revealed that hydrogen bond acceptor and negative ionizable charges are necessary for the inhibition of PCSK 9. The mapping of pharmacophore model confirmed that 96 molecules were as best fit. Finally, toxicity and dynamic simulation studies confirmed that safe and stable complex formation between lead molecules and PCSK 9 receptor. The results concluded that the 15 lead molecules are proved as potent and safe PCSK9 inhibitors.
Aspirin-induced attenuation of adipogenic differentiation of bone marrow mesenchymal stem cells is accompanied by the disturbed epigenetic modification
2018, International Journal of Biochemistry and Cell Biology
Aspirin has positive effects on bone marrow mesenchymal stem cells (BMSCs) osteogenic differentiation. However, researchers did not give much thought to its effect on BMSCs adipogenic differentiation. Here, we analyzed the effect of aspirin on the BMSCs adipogenic differentiation. To detect whether the effect of aspirin on the adipogenic differentiation of BMSCs is associated with the disturbed epigenetic modification, the expression of histone deacetylases (HDACs), activity of HDACs and HAT, global histone H3 acetylation and H3k9 acetylation alterations were investigated. Moreover, to further explore and understand the binding mode between aspirin and HDACs, an attempt was made to identify the interaction between aspirin and the HDACs with the aid of in silico docking study. The results showed that aspirin could induce inhibition of BMSCs adipogenesis. The level of HDAC activity, global histone H3 acetylation, and H3k9 acetylation were all down regulated during adipogenic differentiation, and aspirin can reverse these decreases. Furthermore, the HDAC isoforms have different expression patterns in those progresses. The expression of HDAC9 was increased in a does-dependent manner when aspirin was introduced during BMSCs adipogenic differentiation. Docking study showed that high affinity of HDAC9 to aspirin was existed, suggesting that HDAC9 may has an important role in the process of aspirin-induced suppression of adipogenesis. Further studies are needed to define the intricate mechanisms of the HDAC isoforms, and all of these enable us to understand aspirin and its efficacy of inhibition of adipogenic differentiation and pave the way to aspirin clinical using for the tissue regenerating.
Designing a polytope for use in a broad-spectrum dengue virus vaccine
2018, Journal of Taibah University Medical Sciences
Dengue virus surface proteins are often used in the development of vaccines that protect against dengue virus infection. However, the surface proteins on the four serotypes of dengue virus display high variation, which increases the difficulty of developing a vaccine that can protect against all viral strains. In this study, a polytope that is recognized by broadly neutralizing antibodies (bnAbs) was designed using conserved epitopes from the four serotypes.
We constructed a polytope using four conserved dengue virus epitopes such that two aligned epitopes were separated from the other two epitopes by a histidyl-tRNA synthetase spacer. The epitopes were selected based on our previous docking studies. We then performed molecular docking of the polytope with the four bnAbs.
The polytope bound precisely to the four bnAbs—B7, C8, A11, and C10. Moreover, the polytope had a higher affinity for the bnAbs compared to the DENV2 antigen. The polytope and A11 antibody complex had the lowest binding energy relative to complexes between the polytope and the other three antibodies assessed. The highest total number of hydrogen bonds was found in the polytope and B7 antibody complex. The hydrogen bond length in all the complexes ranged from 2.07 to 3.03 Å, implying that hydrogen bonds stabilized the complexes.
The developed polytope interacted with four different bnAbs that recognize the four serotypes of dengue virus. The results of this study suggest that this polytope warrants further development for use in a broad-spectrum vaccine against dengue virus.
Study on the interaction of artificial and natural food colorants with human serum albumin: A computational point of view
2015, Computational Biology and Chemistry
Citation Excerpt :
To score docking refined poses for each protein–ligand system the formation of hydrogen bonds was evaluated along the refinement MD trajectories, as defined by collective variable H. The value of hydrogen bonds’ spatial and directional properties has shown to be of major interest in analyzing docking conformations (Meyer et al., 1996). Hydrogen bonds are a key contributor to the high specificity of macromolecule interactions, having distributions in good agreement with energy landscapes obtained through electronic structure calculations (Morozov et al., 2004).
Due to the high amount of artificial food colorants present in infants’ diets, their adverse effects have been of major concern among the literature. Artificial food colorants have been suggested to affect children's behavior, being hyperactivity the most common disorder. In this study we compare binding affinities of a group of artificial colorants (sunset yellow, quinoline yellow, carmoisine, allura red and tartrazine) and their natural industrial equivalents (carminic acid, curcumin, peonidin-3-glucoside, cyanidin-3-glucoside) to human serum albumin (HSA) by a docking approach and further refinement through atomistic molecular dynamics simulations. Due to the protein–ligand conformational interface complexity, we used collective variable driven molecular dynamics to refine docking predictions and to score them according to a hydrogen-bond criterion. With this protocol, we were able to rank ligand affinities to HSA and to compare between the studied natural and artificial food additives. Our results show that the five artificial colorants studied bind better to HSA than their equivalent natural options, in terms of their H-bonding network, supporting the hypothesis of their potential risk to human health.

View all citing articles on Scopus

^f1: Corresponding author

View full text

Journal of Molecular Biology

Regular articleHydrogen Bonding and Molecular Surface Shape Complementarity as a Basis for Protein Docking

Abstract

Regular article
Hydrogen Bonding and Molecular Surface Shape Complementarity as a Basis for Protein Docking