Picosecond dynamics of ligand interconversion in the primary docking site of heme proteins

We have used femtosecond IR spectroscopy to probe interconversion dynamics of ligand in the primary docking site of heme proteins under physiological conditions. The docking site, fashioned with highly conserved amino acid residues, modulates ligand-binding activity by mediating the passage of ligand to and from the active binding site. Ligands in two states of the docking site interconvert on the picosecond time scale, and the rates are about 4 times slower in hemoglobin than that in myoglobin. The accurate interconversion rates on the time scale readily accessible by MD simulations can be used to refine computer simulations, which could in turn provide a detailed mechanistic picture of ligand binding in heme proteins.     

Pharmacophore model-aided virtual screening combined with comparative molecular docking and molecular dynamics for identification of marine natural products as SARS-CoV-2 papain-like protease inhibitors

Targeting SARS-CoV-2 papain-like protease using inhibitors is a suitable approach for inhibition of virus replication and dysregulation of host anti-viral immunity. Engaging all five binding sites far from the catalytic site of PLpro is essential for developing a potent inhibitor. We developed and validated a structure-based pharmacophore model with 9 features of a potent PLpro inhibitor. The pharmacophore model-aided virtual screening of the comprehensive marine natural product database predicted 66 initial hits. This hit library was downsized by filtration through a molecular weight filter of ≤ 500 g/mol. The 50 resultant hits were screened by comparative molecular docking using AutoDock and AutoDock Vina. Comparative molecular docking enables benchmarking docking and relieves the disparities in the search and scoring functions of docking engines. Both docking engines retrieved 3 same compounds at different positions in the top 1 % rank, hence consensus scoring was applied, through which CMNPD28766, aspergillipeptide F emerged as the best PLpro inhibitor. Aspergillipeptide F topped the 50-hit library with a pharmacophore-fit score of 75.916. Favorable binding interactions were predicted between aspergillipeptide F and PLpro similar to the native ligand XR8-24. Aspergillipeptide F was able to engage all the 5 binding sites including the newly discovered BL2 groove, site V. Molecular dynamics for quantification of Cα-atom movements of PLpro after ligand binding indicated that it exhibits highly correlated domain movements contributing to the low free energy of binding and a stable conformation. Thus, aspergillipeptide F is a promising candidate for pharmaceutical and clinical development as a potent SARS-CoV-2 PLpro inhibitor.     

Conformational and docking studies of acyl homoserine lactones as a robust method to investigate bioactive conformations

A method aiming at investigating possible bioactive conformations of acyl homoserine lactone (AHL) quorum sensing (QS) modulators is established. The method relies on the exhaustive conformational analysis of AHLs by varying torsion angles around the amide group then on the selection of the closest conformation to those known from co-crystallized XRD data of AHL-receptor complexes. These latter are then docked as rigid ligand within the receptor binding site, leading to interactions with binding site residues which are highly consistent as compared with the data arising from XRD studies. The method is first validated using AHLs for which XRD data of their complexes with their cognate receptor are available, then extended to examples for which the binding mode is still unknown.Three compounds were used to validate the method: hexanoyl homoserine lactone (HHL) as an example of autoinducer, 3-oxo-butanoyl homoserine lactone (OBHL), as a representative model of 3-oxo-AHLs, and 4-(4-chlorophenoxy)butanoyl homoserine lactone (CPOBHL) as an example of a QS inhibitor. The conformational analysis of these three compounds to their cognate protein (TraR, SdiA, LasR and CviR) provides the data which enable the next rigid docking step. Further rigid docking of the closest conformations compared to the known bioactive ones within the binding sites allows to recover the expected binding mode with high precision (atomic RMSD < 2 Å). This “conformational analysis/torsion angle filter/rigid ligand docking” method was then used for investigating three non-natural AHL-type QS inhibitors without known co-crystallized XRD structures, namely was 2-hexenoyl homoserine lactone (HenHL), 3-oxo-4-phenylbutanoyl homoserine lactone (OPBHL) and 3-(4-bromophenyl)propanoyl homoserine lactone (BPPHL). Results provide insights into their possible binding mode by identifying specific interactions with some key residues within the receptor binding site, allowing discussion of their biological activity.     

HIV蛋白酶抑制剂——利迪链菌素的分子对接研究

用分子对接的方法, 对利迪链菌素的抗HIV蛋白酶活性进行了研究. 为了更准确地反映利迪链菌素分子与酶蛋白结合的情况, 充分考虑受体活性部位的柔性, 采用了FlexX(初步对接)和Flexidock(精确对接)分两步将配体与受体进行对接. 在初步对接中, 设计了不同的受体活性部位来考察是否有结合水分子参与抑制剂与酶的结合. 对一种作用方式已知的非肽类HIV蛋白酶抑制剂Aha006进行的对接研究显示, 分子模拟的结果与实际情况吻合得较好, 证明了本文所采用的方法的可靠性. 利迪链菌素与蛋白酶活性部位的对接结果显示, 配体分子与受体之间的结合没有结合水分子的参与, 两者通过5对氢键作用结合成为稳定的复合物. 利迪链菌素占据结合腔, 覆盖了蛋白酶的活性三联体Asp25-Thr26-Gly27, 从而起到抑制其生物活性的作用.     

Binding affinity of substituted ureido‐benzenesulfonamide ligands to the carbonic anhydrase receptor: A theoretical study of enzyme inhibition

The binding properties of a series of benzenesulfonamide inhibitors (4‐substituted‐ureido‐benzenesulfonamides, UBSAs) of human carbonic anhydrase II (hCA II) enzyme with active site residues have been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) model. To account for the important docking interactions between the UBSAs ligand and hCA II enzyme, a molecular docking program AutoDock Vina is used. The molecular docking results obtained by AutoDock Vina revealed that the docked conformer has root mean square deviation value less than 1.50 Å compared to X‐ray crystal structures. The inhibitory activity of UBSA ligands against hCA II is found to be in good agreement with the experimental results. The thermodynamic parameters for inhibitor binding show that hydrogen bonding, hydrophilic, and hydrophobic interactions play a major role in explaining the diverse inhibitory range of these derivatives. Additionally, natural bond orbital analysis is performed to characterize the ligand–metal charge transfer stability. The insights gained from this study have great potential to design new hCA‐II inhibitor, 4‐[3‐(1‐p‐Tolyl‐4‐trifluoromethyl‐1H‐pyrazol‐3‐yl)‐ureido]‐benzenesulfonamide, which belongs to the family of UBSA inhibitors and shows similar type of inhibitor potency with hCA II. This work also reveals that a QM/MM model and molecular docking method are computationally feasible and accurate for studying substrate–protein inhibition. © 2013 Wiley Periodicals, Inc.     

Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

Computational analysis of protein–ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions, such as in docking programs, all rely on equations that sum each type of protein–ligand interactions in order to predict the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions constitute one of the most important part of total interactions between macromolecules. Unlike dispersion forces, they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this study, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and darunavir) were analyzed by using charge densities from the transferable aspherical-atom University at Buffalo Databank (UBDB). Moreover, we analyzed the electrostatic interaction energy for an ensemble of structures, using molecular dynamic simulations to highlight the main features of electrostatic interactions important for binding affinity.     

Profiling the structural determinants for the selectivity of representative factor-Xa and thrombin inhibitors using combined ligand-based and structure-based approaches

The current study deciphers the combined ligand- and structure-based computational insights to profile structural determinants for the selectivity of representative diverse classes of FXa-selective and thrombin-selective as well as dual FXa-thrombin high affinity inhibitors. The thrombin-exclusive insertion 60-loop (D-pocket) was observed to be one of the most notable recognition sites for the known thrombin-selective inhibitors. Based on the topological comparison of four common active-site pockets (S1-S4) of FXa and thrombin, the greater structural disparity was observed in the S4-pocket, which was more symmetrical (U-shaped) in FXa as compared to thrombin mainly due to the presence of L99 and I174 residues in latter in place of Y99 and F174 respectively in former protease. The S2 pocket forming partial roof at the entry of 12 ? deep S1-pocket, with two extended β-sheets running antiparallel to each other by undergoing U-turn (～180?), has two conserved glycine residues forming H-bonds with the bound ligand for governing ligand binding affinity. The docking, scoring, and binding pose comparison of the representative high-affinity and selective inhibitors into the active sites of FXa and thrombin revealed critical residues (S214, Y99, W60D) mediating selectivity through direct- and long-range electrostatic interactions. Interestingly, most of the thrombin-selective inhibitors attained S-shaped conformation in thrombin, while FXa-selective inhibitors attained L-shaped conformations in FXa. The role of residue at 99th position of FXa and thrombin toward governing protease selectivity was further substantiated using molecular dynamics simulations on the wild-type and mutated Y99L FXa bound to thrombin-selective inhibitor 2. Furthermore, predictive CoMFA (FXa q2 = 0.814; thrombin q2 = 0.667) and CoMSIA (FXa q2 = 0.807; thrombin q2 = 0.624) models were developed and validated (FXa r2(test) = 0.823; thrombin r(2)(test) = 0.816) to feature molecular determinants of ligand binding affinity using the docking-based conformational alignments (DBCA) of 141 (88(train)+53(test)) and 39 (27(train)+11(test)) nonamidine class of potent FXa (0.004 ≤ K(i) (nM) ≤ 4700) and thrombin (0.001 ≤ K(i) (nM) ≤ 940) inhibitors, respectively. Interestingly, the ligand-based insights well corroborated with the structure-based insights in terms of the role of steric, electrostatic, and hydrophobic parameters for governing the selectivity for the two proteases. The new computational insights presented in this study are expected to be valuable for understanding and designing potent and selective antithrombotic agents.     

Characterization and crystal structure of a high-affinity pentavalent receptor-binding inhibitor for cholera toxin and E. coli heat-labile enterotoxin

Multivalent ligand design constitutes an attractive avenue to the inhibition of receptor recognition and other biological events mediated by oligomeric proteins with multiple binding sites. One example is the design of multivalent receptor blockers targeting members of the AB(5) bacterial toxin family. We report here the synthesis and characterization of a pentavalent inhibitor for cholera toxin and Escherichia coli heat-labile enterotoxin. This inhibitor is an advance over the symmetric pentacyclen-derived inhibitor described in our earlier work in that it presents five copies of m-nitrophenyl-alpha-D-galactoside (MNPG) rather than five copies of beta-D-galactose. The approximately 100-fold higher single-site affinity of MNPG for the toxin receptor binding site relative to galactose is found to yield a proportionate increase in the affinity and IC50 measured for the respective pentavalent constructs. We show by dynamic light scattering that inhibition of receptor binding by the pentavalent inhibitor is due to 1:1 inhibitor:toxin association rather than to inhibitor-mediated aggregation. This 1:1 association is in complete agreement with a 1.46 A resolution crystal structure of the pentavalent inhibitor:toxin complex, which shows that the favorable single-site binding interactions of MNPG are retained by the five arms of the 5256 Da pentavalent MNPG-based inhibitor and that the initial segment of the linking groups interacts with the surface of the toxin B pentamer.     

Validation of the AmpC β-lactamase binding site and identification of inhibitors with novel scaffolds

AmpC β-lactamase confers resistance to β-lactam antibiotics in multiple Gram-negative bacteria. Therefore, identification of non-β-lactam compounds that inhibit the enzyme is considered crucial to the development of novel antibacterial therapies. Given the highly solvent-exposed active site, it is important to study the induced-fit movements and water-mediated interactions to improve docking accuracy and virtual screening enrichments in structure-based design of new AmpC inhibitors. Here, we tested multiple models of the AmpC binding site to investigate the importance of conserved water molecules and binding site plasticity on molecular docking. The results indicate that at least one conserved water molecule greatly improves the binding pose predictions and virtual screening enrichments of known noncovalent AmpC inhibitors. The best model was tested prospectively in the virtual screening of about 6 million commercially available compounds. Sixty-one chemically diverse top-scoring compounds were experimentally tested, which led to the identification of seven previously unknown inhibitors. These findings validate the essential features of the AmpC binding site for molecular recognition and are useful for further optimization of identified inhibitors.     

Binding of alkaloids into the S1 specificity pocket of α-chymotrypsin: evidence from induced circular dichroism spectra

Non-covalent binding of planar aromatic molecules into the S1 specificity pocket of the serine protease α-chymotrypsin (αCHT) can be detected by measuring induced circular dichroism (CD) spectroscopic signals. Utilizing this phenomenon, αCHT association of proflavine (PRF), the well known serine protease inhibitor has been investigated together with plant-derived compounds including isoquinoline, pyridocarbazole and indoloquinoline alkaloids, of which αCHT binding has never been reported. Non-degenerate exciton coupling between π-π* transitions of the ligand molecules and two tryptophan residues (Trp172 and Trp215) near to the binding site is proposed to be responsible for the induced CD activity. The association constants calculated from CD titration data indicated strong αCHT association of sanguninarine, ellipticine, desmethyl-isocryptolepine and isoneocryptolepine (K(a) ≈ 10(5) M(-1)) while berberine, coptisine and chelerythrine bind to the enzyme with lower, PRF-like affinity (K(a) ≈ 10(4) M(-1)). PRF-trypsin and ellipticine-trypsin binding interactions have also been demonstrated. The binding of the alkaloids into the S1 pocket of αCHT has been confirmed by CD competition experiments. Molecular docking calculations showed the inclusion of PRF as well as the alkaloid molecules in the S1 cavity where they are stabilized by hydrophobic and H-bonding interactions. These novel nonpeptidic scaffolds can be used for developing selective inhibitors of serine proteases having chymotrypsin-like folds. Furthermore, the results provide a novel, CD spectroscopic based approach for probing the ligand binding of αCHT and related proteases.     

Pharmacophore anchor models of ATAT1 to discover potential inhibitors and lead optimization

Post-translation modification of microtubules is associated with many diseases like cancer. Alpha Tubulin Acetyltransferase 1 (ATAT1) is a major enzyme that acetylates ‘Lys-40’ in alpha-tubulin on the luminal side of microtubules and is a drug target that lacks inhibitors. Here, we developed pharmacophore anchor models of ATAT1 which were constructed statistically using thousands of docked compounds, for drug design and investigating binding mechanisms. Our models infer the compound moiety preferences with the physico-chemical properties for the ATAT1 binding site. The results from the pharmacophore anchor models show the three main sub-pockets, including S1 acetyl site, S2 adenine site, and S3 diphosphate site with anchors, where conserved moieties interact with respective sub-pocket residues in each site and help in guiding inhibitor discovery. We validated these key anchors by analyzing 162 homologous protein sequences (>99 species) and over 10 structures with various bound ligands and mutations. Our results were consistent with previous works also providing new interesting insights. Our models applied in virtual screening predicted several ATAT1 potential inhibitors. We believe that our model is useful for future inhibitor discovery and for guiding lead optimization.     

Design of multi-binding-site inhibitors, ligand efficiency, and consensus screening of avian influenza H5N1 wild-type neuraminidase and of the oseltamivir-resistant H274Y variant

The binding sites of wild-type avian influenza A H5N1 neuraminidase, as well as those of the Tamiflu (oseltamivir)-resistant H274Y variant, were explored computationally to design inhibitors that target simultaneously several adjacent binding sites of the open conformation of the virus protein. The compounds with the best computed free energies of binding, in agreement by two docking methods, consensus scoring, and ligand efficiency values, suggest that mimicking a polysaccharide, beta-lactam, and other structures, including known drugs, could be routes for multibinding site inhibitor design. This new virtual screening method based on consensus scoring and ligand efficiency indices is introduced, which allows the combination of pharmacodynamic and pharmacokinetic properties into unique measures.     

Hydration in drug design. 3. Conserved water molecules at the ligand-binding sites of homologous proteins

Summary Water molecules are known to play an important rôle in mediating protein-ligand interactions. If water molecules are conserved at the ligand-binding sites of homologous proteins, such a finding may suggest the structural importance of water molecules in ligand binding. Structurally conserved water molecules change the conventional definition of binding sites by changing the shape and complementarity of these sites. Such conserved water molecules can be important for site-directed ligand/drug design. Therefore, five different sets of homologous protein/protein-ligand complexes have been examined to identify the conserved water molecules at the ligand-binding sites. Our analysis reveals that there are as many as 16 conserved water molecules at the FAD binding site of glutathione reductase between the crystal structures obtained from human and E. coli. In the remaining four sets of high-resolution crystal structures, 2–4 water molecules have been found to be conserved at the ligand-binding sites. The majority of these conserved water molecules are either bound in deep grooves at the protein-ligand interface or completely buried in cavities between the protein and the ligand. All these water molecules, conserved between the protein/protein-ligand complexes from different species, have identical or similar apolar and polar interactions in a given set. The site residues interacting with the conserved water molecules at the ligand-binding sites have been found to be highly conserved among proteins from different species; they are more conserved compared to the other site residues interacting with the ligand. These water molecules, in general, make multiple polar contacts with protein-site residues.     

Structural Insights into Novel 15-Prostaglandin Dehydrogenase Inhibitors

We discovered SW033291 in a high throughput chemical screen aimed at identifying 15-prostaglandin dehydrogenase (15-PGDH) modulators. The compound exhibited inhibitory activity in in vitro biochemical and cell-based assays of 15-PGDH activity. We subsequently demonstrated that this compound, and several analogs thereof, are effective in in vivo mouse models of bone marrow transplant, colitis, and liver regeneration, where increased levels of PGE2 positively potentiate tissue regeneration. To better understand the binding of SW033291, we carried out docking studies for both the substrate, PGE2, and an inhibitor, SW033291, to 15-PGDH. Our models suggest similarities in the ways that PGE2 and SW033291 interact with key residues in the 15-PGDH-NAD+ complex. We carried out molecular dynamics simulations (MD) of SW033291 bound to this complex, in order to understand the dynamics of the binding interactions for this compound. The butyl side chain (including the sulfoxide) of SW033291 participates in crucial binding interactions that are similar to those observed for the C15-OH and the C16-C20 alkyl chain of PGE2. In addition, interactions with residues Ser138, Tyr151, and Gln148 play key roles in orienting and stabilizing SW033291 in the binding site and lead to enantioselectivity for the R-enantiomer. Finally, we compare the binding mode of (R)-S(O)-SW033291 with the binding interactions of published 15-PGDH inhibitors.     

FAST-NMR: functional annotation screening technology using NMR spectroscopy

An abundance of protein structures emerging from structural genomics and the Protein Structure Initiative (PSI) are not amenable to ready functional assignment because of a lack of sequence and structural homology to proteins of known function. We describe a high-throughput NMR methodology (FAST-NMR) to annotate the biological function of novel proteins through the structural and sequence analysis of protein-ligand interactions. This is based on basic tenets of biochemistry where proteins with similar functions will have similar active sites and exhibit similar ligand binding interactions, despite global differences in sequence and structure. Protein-ligand interactions are determined through a tiered NMR screen using a library composed of compounds with known biological activity. A rapid co-structure is determined by combining the experimental identification of the ligand binding site from NMR chemical shift perturbations with the protein-ligand docking program AutoDock. Our CPASS (Comparison of Protein Active Site Structures) software and database are then used to compare this active site with proteins of known function. The methodology is demonstrated using unannotated protein SAV1430 from Staphylococcus aureus.     

Identification of novel human USP2 inhibitor and its putative role in treatment of COVID-19 by inhibiting SARS-CoV-2 papain-like (PLpro) protease

Human ubiquitin carboxyl-terminal hydrolase-2 (USP2) inhibitors, such as thiopurine analogs, have been reported to inhibit SARS-CoV papain-like proteases (PLpro). The PLpro have significant functional implications in the innate immune response during SARS-CoV-2 infection and considered an important antiviral target. Both proteases share strikingly similar USP fold with right-handed thumb–palm–fingers structural scaffold and conserved catalytic triad Cys-His-Asp/Asn. In this urgency situation of COVID-19 outbreak, there is a lack of in-vitro facilities readily available to test SARS-CoV-2 inhibitors in whole-cell assays. Therefore, we adopted an alternate route to identify potential USP2 inhibitor through integrated in-silico efforts. After an extensive virtual screening protocol, the best compounds were selected and tested. The compound Z93 showed significant IC₅₀ value against Jurkat (9.67 μM) and MOTL-4 cells (11.8 μM). The binding mode of Z93 was extensively analyzed through molecular docking, followed by MD simulations, and molecular interactions were compared with SARS-CoV-2. The relative binding poses of Z93 fitted well in the binding site of both proteases and showed consensus π-π stacking and H-bond interactions with histidine and aspartate/asparagine residues of the catalytic triad. These results led us to speculate that compound Z93 might be the first potential chemical lead against SARS-CoV-2 PLpro, which warrants in-vitro evaluations.     

Potent HIV-1 protease inhibitors incorporating meso-bicyclic urethanes as P2-ligands: structure-based design, synthesis, biological evaluation and protein-ligand X-ray studies

Recently, we designed a series of novel HIV-1 protease inhibitors incorporating a stereochemically defined bicyclic fused cyclopentyl (Cp-THF) urethane as the high affinity P2-ligand. Inhibitor with this P2-ligand has shown very impressive potency against multi-drug-resistant clinical isolates. Based upon the -bound HIV-1 protease X-ray structure, we have now designed and synthesized a number of meso-bicyclic ligands which can conceivably interact similarly to the Cp-THF ligand. The design of meso-ligands is quite attractive as they do not contain any stereocenters. Inhibitors incorporating urethanes of bicyclic-1,3-dioxolane and bicyclic-1,4-dioxane have shown potent enzyme inhibitory and antiviral activities. Inhibitor (K(i) = 0.11 nM; IC(50) = 3.8 nM) displayed very potent antiviral activity in this series. While inhibitor showed comparable enzyme inhibitory activity (K(i) = 0.18 nM) its antiviral activity (IC(50) = 170 nM) was significantly weaker than inhibitor . Inhibitor maintained an antiviral potency against a series of multi-drug resistant clinical isolates comparable to amprenavir. A protein-ligand X-ray structure of -bound HIV-1 protease revealed a number of key hydrogen bonding interactions at the S2-subsite. We have created an active model of inhibitor based upon this X-ray structure.