Search

Claudin-4 (CLDN4) is a class of transmembrane protein in the family of tight junction (TJ) proteins. Overexpression of CLDN4 is reported in the case of ovarian cancer and epithelial malignancies. The current study is focused on the identification of lead compounds for CLDN4 adopting the structure-based drug design method. The Schrodinger glide is used as a molecular docking tool for the initial docking of CLDN4 with Asinex Database by performing high throughput virtual screening, top hits were identified. Then, compounds BDF 33196188 and BDE 30874918 were identified by molecular docking based on binding energy in the active site of CLDN4. Subsequently, critical residues were identified such as Asp146 and Arg158 with the least binding energy from Extra Precision method. Further, molecular dynamics simulations of claudin-4 protein were used for the optimization of best ligands with claudin-4 in a dynamic system. Molecular docking and molecular dynamics simulations predicted critically important residues ASP146 and ARG158 involved in claudin-4 binding. The hits retrieved from screening were docked into protein by relevant procedures including HTVS, SP, and XP. Finally, two molecules were identified as potential claudin-4 inhibitors. The two ligands BDF 33196188 and BDE 30874918 are suggested as potential inhibitors for CLDN4. In summary, our computational strategy established novel leads against CLDN4 from Asinex Database and recommended as anti-cancer agents.

Inhibitors of phosphodiesterases (PDEs) have been widely studied as therapeutics for the treatment of human diseases, but improvement of inhibitor selectivity is still desirable for the enhancement of inhibitor potency. Here, we report identification of a water-containing subpocket as a PDE4-specific pocket for inhibitor binding. We designed against the pocket and synthesized two enantiomers of PDE4 inhibitor Zl-n-91. The ( S)-Zl-n-91 enantiomer showed IC 50 values of 12 and 20 nM for the catalytic domains of PDE4D2 and PDE4B2B, respectively, selectivity several thousand-fold greater than those of other PDE families, and potent neuroprotection activities. Crystal structures of the PDE4D2 catalytic domain in complex with each Zl-n-91 enantiomer revealed that ( S)-Zl-n-91 but not ( R)-Zl-n-91 formed a hydrogen bond with the bound water in the pocket, thus explaining its higher affinity. The structural superposition between the PDE families revealed that this water-containing subpocket is unique to PDE4 and thus valuable for the design of PDE4 selective inhibitors.

With the growing number of labor-intensive data in the pharmaceutical industries and public domain for protein-ligand complexes, a significant challenge is still remaining in managing and leveraging this vast information. Here, a standalone application is presented for analysis, organization, and illustration of structural data and molecular interactions for exploiting 3D-structures into simple 1D fingerprints. The utility of the approach was shown in unraveling a feasible solution for post-processing of docking results in parallel with providing fruitful analysis for users in order to investigate molecular interactions. Remarkably, all interaction possibilities including (hydrogen bond, water-bridged, electrostatic, and hydrophobic as well as π- π and cation-π interactions) are supported both in the form of fingerprints and compelling reports. These investigations are mainly considered based on right orientation, location, and geometry of the interacting pairs rather than the acquisition of the energy terms. The reasonable efficiency of our application in different models was comparable to recent methods It is clearly presented that FLIP provides a faster way to generate usable fingerprints for ligand and protein binding modes. FLIP is free for academic use and is available at: http://zistrayan.com/development/download/flip/package.zip., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Molecular recognition is a highly interdependent process. Subsite couplings within the active site of proteases are most often revealed through conditional amino acid preferences in substrate recognition. However, the potential effect of these couplings on inhibition and thus inhibitor design is largely unexplored. The present study examines the interdependency of subsites in HIV-1 protease using a focused library of protease inhibitors, to aid in future inhibitor design. Previously a series of darunavir (DRV) analogs was designed to systematically probe the S1' and S2' subsites. Co-crystal structures of these analogs with HIV-1 protease provide the ideal opportunity to probe subsite interdependency. All-atom molecular dynamics simulations starting from these structures were performed and systematically analyzed in terms of atomic fluctuations, intermolecular interactions, and water structure. These analyses reveal that the S1' subsite highly influences other subsites: the extension of the hydrophobic P1' moiety results in 1) reduced van der Waals contacts in the P2' subsite, 2) more variability in the hydrogen bond frequencies with catalytic residues and the flap water, and 3) changes in the occupancy of conserved water sites both proximal and distal to the active site. In addition, one of the monomers in this homodimeric enzyme has atomic fluctuations more highly correlated with DRV than the other monomer. These relationships intricately link the HIV-1 protease subsites and are critical to understanding molecular recognition and inhibitor binding. More broadly, the interdependency of subsite recognition within an active site requires consideration in the selection of chemical moieties in drug design; this strategy is in contrast to what is traditionally done with independent optimization of chemical moieties of an inhibitor.

We report the high-resolution (1.9 Å) crystal structure of oligomycin bound to the subunit c(10) ring of the yeast mitochondrial ATP synthase. Oligomycin binds to the surface of the c(10) ring making contact with two neighboring molecules at a position that explains the inhibitory effect on ATP synthesis. The carboxyl side chain of Glu59, which is essential for proton translocation, forms an H-bond with oligomycin via a bridging water molecule but is otherwise shielded from the aqueous environment. The remaining contacts between oligomycin and subunit c are primarily hydrophobic. The amino acid residues that form the oligomycin-binding site are 100% conserved between human and yeast but are widely different from those in bacterial homologs, thus explaining the differential sensitivity to oligomycin. Prior genetics studies suggest that the oligomycin-binding site overlaps with the binding site of other antibiotics, including those effective against Mycobacterium tuberculosis, and thereby frames a common "drug-binding site." We anticipate that this drug-binding site will serve as an effective target for new antibiotics developed by rational design.

STAT3 is a transcription factor that has been found to be constitutively activated in a number of human cancers. Dimerization of STAT3 via its SH2 domain and the subsequent translocation of the dimer to the nucleus leads to transcription of anti-apoptotic genes. Prevention of the dimerization is thus an attractive strategy for inhibiting the activity of STAT3. Phosphotyrosine-based peptidomimetic inhibitors, which mimic pTyr-Xaa-Yaa-Gln motif and have strong to weak binding affinities, have been previously investigated. It is well-known that structures of protein-inhibitor complexes are important for understanding the binding interactions and designing stronger inhibitors. Experimental structures of inhibitors bound to the SH2 domain of STAT3 are, however, unavailable. In this paper we describe a computational study that combined molecular docking and molecular dynamics to model structures of 12 peptidomimetic inhibitors bound to the SH2 domain of STAT3. A detailed analysis of the modeled structures was performed to evaluate the characteristics of the binding interactions. We also estimated the binding affinities of the inhibitors by combining MMPB/GBSA-based energies and entropic cost of binding. The estimated affinities correlate strongly with the experimentally obtained affinities. Modeling results show binding modes that are consistent with limited previous modeling studies on binding interactions involving the SH2 domain and phosphotyrosine(pTyr)-based inhibitors. We also discovered a stable novel binding mode that involves deformation of two loops of the SH2 domain that subsequently bury the C-terminal end of one of the stronger inhibitors. The novel binding mode could prove useful for developing more potent inhibitors aimed at preventing dimerization of cancer target protein STAT3.

In this study we have calculated a 3D structure of apoptin and through modeling and docking approaches, we show its interaction with Bcr-Abl oncoprotein and its downstream signaling components, following which we confirm some of the newly-found interactions by biochemical methods. Bcr-Abl oncoprotein is aberrantly expressed in chronic myelogenous leukaemia (CML). It has several distinct functional domains in addition to the Abl kinase domain. The SH3 and SH2 domains cooperatively play important roles in autoinhibiting its kinase activity. Adapter molecules such as Grb2 and CrkL interact with proline-rich region and activate multiple Bcr-Abl downstream signaling pathways that contribute to growth and survival. Therefore, the oncogenic effect of Bcr-Abl could be inhibited by the interaction of small molecules with these domains. Apoptin is a viral protein with well-documented cancer-selective cytotoxicity. Apoptin attributes such as SH2-like sequence similarity with CrkL SH2 domain, unique SH3 domain binding sequence, presence of proline-rich segments, and its nuclear affinity render the molecule capable of interaction with Bcr-Abl. Despite almost two decades of research, the mode of apoptin's action remains elusive because 3D structure of apoptin is unavailable. We performed in silico three-dimensional modeling of apoptin, molecular docking experiments between apoptin model and the known structure of Bcr-Abl, and the 3D structures of SH2 domains of CrkL and Bcr-Abl. We also biochemically validated some of the interactions that were first predicted in silico. This structure-property relationship of apoptin may help in unlocking its cancer-selective toxic properties. Moreover, such models will guide us in developing of a new class of potent apoptin-like molecules with greater selectivity and potency.

The regulatory (R) subunit of protein kinase A serves to modulate the activity of protein kinase A in a cAMP-dependent manner and exists in two distinct and structurally dissimilar, end point cAMP-bound "B" and C-subunit-bound "H"-conformations. Here we report mechanistic details of cAMP action as yet unknown through a unique approach combining x-ray crystallography with structural proteomics approaches, amide hydrogen/deuterium exchange and ion mobility mass spectrometry, applied to the study of a stereospecific cAMP phosphorothioate analog and antagonist((Rp)-cAMPS). X-ray crystallography shows cAMP-bound R-subunit in the B form but surprisingly the antagonist Rp-cAMPS-bound R-subunit crystallized in the H conformation, which was previously assumed to be induced only by C-subunit-binding. Apo R-subunit crystallized in the B form as well but amide exchange mass spectrometry showed large differences between apo, agonist and antagonist-bound states of the R-subunit. Further ion mobility reveals the apo R-subunit as an ensemble of multiple conformations with collisional cross-sectional areas spanning both the agonist and antagonist-bound states. Thus contrary to earlier studies that explained the basis for cAMP action through "induced fit" alone, we report evidence for conformational selection, where the ligand-free apo form of the R-subunit exists as an ensemble of both B and H conformations. Although cAMP preferentially binds the B conformation, Rp-cAMPS interestingly binds the H conformation. This reveals the unique importance of the equatorial oxygen of the cyclic phosphate in mediating conformational transitions from H to B forms highlighting a novel approach for rational structure-based drug design. Ideal inhibitors such as Rp-cAMPS are those that preferentially "select" inactive conformations of target proteins by satisfying all "binding" constraints alone without inducing conformational changes necessary for activation.

Small changes in ligand structure can lead to large unexpected changes in activity yet it is often not possible to rationalize these effects using empirical modeling techniques, suggesting more effective methods are required. In this study we investigate the use of high level QM methods to study the interactions found within protein-ligand complexes as improved understanding of these could help in the design of new, more active molecules. We study aspects of ligand binding in a set of protein ligand complexes containing ligand efficient, fragment-like inhibitors as these structures are often challenging to determine experimentally. To assess the reliability of our theoretical models we compare the MP2/6-31+G** QM results to the original X-ray coordinates and to QM/MM B3LYP/6-31G*//UFF results which we have previously reported. We also contrast these results with data obtained from an analysis of the distribution of comparable interactions found in (a) high resolution kinase complexes (≤ 1.8Å) from the PDB and (b) more generic, small molecule crystal structures from the CSD., (Copyright © 2010 Elsevier Inc. All rights reserved.)

The worldwide TB structural genomics initiative has identified several new drug targets for Mycobacterium tuberculosis (M. tb). Dihydrofolate reductase (DHFR) catalyzes the NADPH-dependent reduction of dihydrofolate to tetrahydrofolate that is essential for DNA synthesis. Inhibition of its activity leads to arrest of DNA synthesis and hence cell death. Thus, M. tb DHFR (mtDHFR) is an attractive novel drug target for developing anti-TB drugs. Structural comparison of mtDHFR and human DHFR (hDHFR) reveals key differences in the active sites. These differences can be exploited for the design of selective inhibitors for mtDHFR. Based on the recently determined high resolution crystal structure of mtDHFR complexed with known inhibitor methotrexate (MTX) and cofactor NADPH, a tri-peptide inhibitor has been identified using a structure-based drug design approach. Docking studies indicate that the designed tripeptide inhibitor has a high potency (K (d) = 1.78 nM) and is a selective (approximately 120 fold over hDHFR) inhibitor for mtDHFR. Hence, the tripeptide is a suitable lead compound for the development of novel anti-TB drugs.

In the quest for the rational design of selective and potent inhibitors for members of the pancreatic ribonuclease A (RNase A) family of biomedical interest, the binding of uridine 5'-phosphate (U5P) and uridine 5'-diphosphate (UDP) to RNase A have been investigated using kinetic studies and X-ray crystallography. Both nucleotides are competitive inhibitors of the enzyme, with K(i) values of 4.0 and 0.65 mM, respectively. They bind to the active site of the enzyme by anchoring two molecules connected to each other by hydrogen bonds and van der Waals interactions. While the first of the inhibitor molecules binds with its nucleobase in the pyrimidinyl-binding subsite, the second is bound at the purine-preferring subsite. The unexpected binding of a pyrimidine at the purine-binding subsite has added new important elements to the rational design approach for the discovery of new potent inhibitors of the RNase A superfamily.

Severe acute respiratory syndrome coronavirus (SARS-CoV) main protease (M(pro)), a protein required for the maturation of SARS-CoV, is vital for its life cycle, making it an attractive target for structure-based drug design of anti-SARS drugs. The structure-based virtual screening of a chemical database containing 58,855 compounds followed by the testing of potential compounds for SARS-CoV M(pro) inhibition leads to two hit compounds. The core structures of these two hits, defined by the docking study, are used for further analogue search. Twenty-one analogues derived from these two hits exhibited IC50 values below 50 microM, with the most potent one showing 0.3 microM. Furthermore, the complex structures of two potent inhibitors with SARS-CoV M(pro) were solved by X-ray crystallography. They bind to the protein in a distinct manner compared to all published SARS-CoV M(pro) complex structures. They inhibit SARS-CoV M(pro) activity via intensive H-bond network and hydrophobic interactions, without the formation of a covalent bond. Interestingly, the most potent inhibitor induces protein conformational changes, and the inhibition mechanisms, particularly the disruption of catalytic dyad (His41 and Cys145), are elaborated.

Background: The aspartic proteinase renin plays an important physiological role in the regulation of blood pressure. It catalyses the first step in the conversion of angiotensinogen to the hormone angiotensin II. In the past, potent peptide inhibitors of renin have been developed, but none of these compounds has made it to the end of clinical trials. Our primary aim was to develop novel nonpeptide inhibitors. Based on the available structural information concerning renin-substrate interactions, we synthesized inhibitors in which the peptide portion was replaced by lipophilic moieties that interact with the large hydrophobic S1/S3-binding pocket in renin., Results: Crystal structure analysis of renin-inhibitor complexes combined with computational methods were employed in the medicinal-chemistry optimisation process. Structure analysis revealed that the newly designed inhibitors bind as predicted to the S1/S3 pocket. In addition, however, these compounds interact with a hitherto unrecognised large, distinct, sub-pocket of the enzyme that extends from the S3-binding site towards the hydrophobic core of the enzyme. Binding to this S3(sp) sub-pocket was essential for high binding affinity. This unprecedented binding mode guided the drug-design process in which the mostly hydrophobic interactions within subsite S3(sp) were optimised., Conclusions: Our design approach led to compounds with high in vitro affinity and specificity for renin, favourable bioavailability and excellent oral efficacy in lowering blood pressure in primates. These renin inhibitors are therefore potential therapeutic agents for the treatment of hypertension and related cardiovascular diseases.