PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “molecular docking”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48Linked to original sources

Probing the cons and pros of lectin-induced immunomodulation: case studies for the mistletoe lectin and galectin-1.

When imagining to monitor animal cells through a microscope with resolution at the molecular level, a salient attribute of their surfaces will be the abundance of glycan chains. They present galactosides at their termini widely extending like tentacles into the extracellular space. Their spatial accessibility and their potential for structural variability endow especially these glycan parts with capacity to act as docking points for molecular sensors (sugar receptors such as lectins). Binding and ligand clustering account for transmission of post-binding signals into the cell interior. The range of triggered activities has turned plant lectins into popular tools in cell biology and immunology. Potential for clinical application has been investigated rigorously only in recent years. As documented in vitro and in vivo for the galactoside-specific mistletoe lectin, its apparent immunomodulatory capacity reflected in upregulation of production of proinflammatory cytokines will not necessarily be clinically favorable but a double-edged sword. In fact, lectin application has been shown to stimulate tumor growth in cell lines, histocultures of human tumors and in two animal models using chemical carcinogenesis or tumor transplantation. When testing immunological effects of the endogenous lectin galectin-1, protection against disorders mediated by activated T cells came up for consideration. Elimination of these cells via CD7-dependent induction of apoptosis, and a shift to the Th2 response by the galectin, are factors to ameliorate disease states. This result encourages further efforts with other galectins. Functional redundancy, synergism, diversity or antagonism among galectins are being explored to understand the actual role of this class of endogenous lectins in inflammation. Regardless of the results of further preclinical testing for galectin-1, these two case studies break new ground in our understanding how glycans as ligands for lectins convey reactivity to immune cells, with impact on the course of a tumor or autoimmune disease.

Adjuvants, Immunologic↗

Docking studies on monoamine oxidase-B inhibitors: estimation of inhibition constants (K(i)) of a series of experimentally tested compounds.

Monoamine oxidase (EC1.4.3.4; MAO) is a mitochondrial outer membrane flavoenzyme that catalyzes the oxidation of biogenic amines. It has two distinct isozymic forms designated MAO-A and MAO-B, each displaying different substrate and inhibitor specificities. They are the well-known targets for antidepressant and neuroprotective drugs. Elucidation of the X-ray crystallographic structure of MAO-B has opened the way for molecular modeling studies. A series of experimentally tested (1-10) model compounds has been docked computationally to the active site of the MAO-B enzyme. The AutoDock 3.0.5 program was employed to perform automated molecular docking. The free energies of binding (DeltaG) and inhibition constants (K(i)) of the docked compounds were calculated by the Lamarckian Genetic Algorithm (LGA) of AutoDock 3.0.5. Excellent to good correlations between the calculated and experimental K(i) values were obtained.

Crystallography, X-Ray↗

Development of the force field parameters for phosphoimidazole and phosphohistidine.

Phosphorylation of histidine-containing proteins is a key step in the mechanism of many phosphate transfer enzymes (kinases, phosphatases) and is the first stage in a wide variety of signal transduction cascades in bacteria, yeast, higher plants, and mammals. Studies of structural and dynamical aspects of such enzymes in the phosphorylated intermediate states are important for understanding the intimate molecular mechanisms of their functioning. Such information may be obtained via molecular dynamics and/or docking simulations, but in this case appropriate force field parameters for phosphohistidine should be explicitly defined. In the present article we describe development of the GROMOS96 force field parameters for phosphoimidazole molecule--a realistic model of the phosphohistidine side chain. The parameterization is based on the results of ab initio quantum chemical calculations with subsequent refinement and testing using molecular mechanics and molecular dynamics simulations. The set of force constants and equilibrium geometry is employed to derive force field for the phosphohistidine moiety. Resulting parameters and topology are incorporated into the molecular modeling package GROMACS and used in molecular dynamics simulations of a phosphohistidine-containing protein in explicit solvent.

Algorithms↗

Effects of variable docking conditions and scoring functions on corresponding protein-aligned comparative molecular field analysis models constructed from diverse human protein tyrosine phosphatase 1B inhibitors.

The effects of variable docking conditions and scoring functions on corresponding protein-aligned comparative molecular field analysis (CoMFA) models have been assessed. To this end, a group of diverse inhibitors were docked into the active site of human protein tyrosine phosphatase 1B (h-PTP 1B). The docked structures were utilized to construct corresponding protein-aligned CoMFA models by employing probe-based (H+, OH, CH3) energy grids and genetic partial least squares (G/PLS) statistical analysis. A total of 48 different docking configurations were evaluated, of which some succeeded in producing self-consistent and predictive CoMFA models. However, the best CoMFA model coincided with docking the un-ionized ligands into the hydrated form of the binding site via the PLP1 scoring function and restricted docking settings (r2(LOO) = 0.647, r2(PRESS) against 27 test compounds = 0.617). Interestingly, the most significant CoMFA models were orthogonal and corresponded to significantly different docked conformers/poses. To utilize the predictive potentials of the best CoMFA models collectively, it was decided to combine them in a single quantitative structure-activity relationship (QSAR) model. The combination model illustrated excellent statistical properties (r2(LOO) = 0.890, r2(PRESS) against 27 test compounds = 0.750).

Binding Sites↗

Structure-based design of amidinophenylurea-derivatives for factor VIIa inhibition.

The amidinophenylurea scaffold was earlier shown to provide an excellent template for the synthesis of novel and potent inhibitors of the blood coagulation factor VIIa. In this contribution we describe the structure-based design of potent ligands guided by X-ray crystallography, molecular modeling and docking studies. The design and synthetic efforts were directed towards novel modifications to explore the protease binding region close to the S4 subsite.

Binding Sites↗

Mutagenesis of specific amino acids converts carnitine acetyltransferase into carnitine palmitoyltransferase.

Carnitine acyltransferases catalyze the exchange of acyl groups between carnitine and CoA. The members of the family can be classified on the basis of their acyl-CoA selectivity. Carnitine acetyltransferases (CrATs) are very active toward short-chain acyl-CoAs but not toward medium- or long-chain acyl-CoAs. Previously, we identified an amino acid residue (Met(564) in rat CrAT) that was critical to fatty acyl-chain-length specificity. M564G-mutated CrAT behaved as if its natural substrates were medium-chain acyl-CoAs, similar to that of carnitine octanoyltransferase (COT). To extend the specificity of rat CrAT to other substrates, we have performed new mutations. Using in silico molecular modeling procedures, we have now identified a second putative amino acid involved in acyl-CoA specificity (Asp(356) in rat CrAT). The double CrAT mutant D356A/M564G showed 6-fold higher activity toward palmitoyl-CoA than that of the single CrAT mutant M564G and a new activity toward stearoyl-CoA. We show that by performing two amino acid replacements a CrAT can be converted into a pseudo carnitine palmitoyltransferase (CPT) in terms of substrate specificity. To change CrAT specificity from carnitine to choline, we also prepared a mutant CrAT that incorporates four amino acid substitutions (A106M/T465V/T467N/R518N). The quadruple mutant shifted the catalytic discrimination between l-carnitine and choline in favor of the latter substrate and showed a 9-fold increase in catalytic efficiency toward choline compared with that of the wild-type. Molecular in silico docking supports kinetic data for the positioning of substrates in the catalytic site of CrAT mutants.

Amino Acid Sequence↗

Structural variations of piritrexim, a lipophilic inhibitor of human dihydrofolate reductase: synthesis, antitumor activity and molecular modeling investigations.

Piritrexim (PTX) (1), a lipophilic inhibitor of the human dihydrofolate reductase, has been evaluated as an anticancer agent. The synthesis of four structural variations (2-5) of PTX is reported. The PTX analogues 2-5 were obtained by reaction of suitable C3-building blocks with pyrimidine-2,4,6-triamine (14) or with cyanacetamide (7) and guanidine (10). The evaluation of 2-4 for antitumor activity against a panel of 60 human cancer cell lines showed inhibitory effects on the growth of the cell lines. These data are supported by molecular modeling and docking studies, which show that compounds 2-4 share the same binding mode within the DHFR active site. Moreover, the estimated ligand binding energies are in good agreement with the experimental activity data.

Antineoplastic Agents↗

Inhibitors of urokinase reduce size of prostate cancer xenografts in severe combined immunodeficient mice.

Proteolytic enzymes are required to mediate tumor cell invasion and metastasis. The urokinase plasminogen activator (uPA) is commonly overexpressed by many human cancers. Therefore, uPA is a logical target to inhibit cancer invasion and metastasis. However, uPA inhibitors also reduce tumor growth. We used a mutated form of plasminogen activator inhibitor type 1 to conform a correlation between the inactivation of uPA and tumor size; we have compared these results with the action of p-aminobenzamidine and amiloride, known inhibitors of uPA. Our results show that blocking uPA by uPA inhibitors reduces tumor size in experimental animals. Our molecular simulation of docking inhibitors to the urokinase reveals that all tested small molecule inhibitors bind in proximity of uPA's specificity pocket, a critical site for future search of novel anticancer uPA inhibitors.

Amiloride↗

The liverwort contains a lectin that is structurally and evolutionary related to the monocot mannose-binding lectins.

A mannose (Man)-binding lectin has been isolated and characterized from the thallus of the liverwort Marchantia polymorpha. N-terminal sequencing indicated that the M. polymorpha agglutinin (Marpola) shares sequence similarity with the superfamily of monocot Man-binding lectins. Searches in the databases yielded expressed sequence tags encoding Marpola. Sequence analysis, molecular modeling, and docking experiments revealed striking structural similarities between Marpola and the monocot Man-binding lectins. Activity and specificity studies further indicated that Marpola is a much stronger agglutinin than the Galanthus nivalis agglutinin and exhibits a preference for methylated Man and glucose, which is unprecedented within the family of monocot Man-binding lectins. The discovery of Marpola allows us, for the first time, to corroborate the evolutionary relationship between a lectin from a lower plant and a well-established lectin family from flowering plants. In addition, the identification of Marpola sheds a new light on the molecular evolution of the superfamily of monocot Man-binding lectins. Beside evolutionary considerations, the occurrence of a G. nivalis agglutinin homolog in a lower plant necessitates the rethinking of the physiological role of the whole family of monocot Man-binding lectins.

Amino Acid Sequence↗

Molecular determinants of the hanatoxin binding in voltage-gated K+-channel drk1.

The carboxyl terminus of S3 segment (S3(C)) in voltage-gated potassium channels was proposed to bear the binding site for gating modifier toxins like Hanatoxin and a helical secondary structural arrangement was suggested. Due to the lack of complete structure in high resolution for such a channel molecule, no further direct experimental data to elucidate the mechanism for their binding conformations could thus far be derived. In order to examine the putative three-dimensional structure of S3(C) and to illustrate the residues required for Hanatoxin binding, molecular simulation and docking were performed, based on the solution structure of Hanatoxin and the structural information from lysine-scanning results for S3(C) fragment. From our results, it is indicated that both hydrophobic and electrostatic interactions are utilized to stabilize the toxin binding. Detailed docking residues and appropriate orientation for binding regarding hydrophobic/-philic environments are also described. Compared with the functional data proposed by previous studies, the helical structural arrangement for the C-terminus of S3 segment in voltage-gated potassium channels can therefore be further emphasized.

Binding Sites↗

(S)-2-Amino-3-(3-hydroxy-7,8-dihydro-6H-cyclohepta[d]isoxazol-4-yl)propionic acid, a potent and selective agonist at the GluR5 subtype of ionotropic glutamate receptors. Synthesis, modeling, and molecular pharmacology.

We have previously described (RS)-2-amino-3-(3-hydroxy-7,8-dihydro-6H-cyclohepta[d]isoxazol-4-yl)propionic acid (4-AHCP) as a highly effective agonist at non-N-methyl-d-aspartate (non-NMDA) glutamate (Glu) receptors in vivo, which is more potent than (RS)-2-amino-3-(3-hydroxy-5-methylisoxazol-4-yl)propionic acid (AMPA) but inactive at NMDA receptors. However, 4-AHCP was found to be much weaker than AMPA as an inhibitor of [(3)H]AMPA binding and to have limited effect in a [(3)H]kainic acid binding assay using rat cortical membranes. To shed light on the mechanism(s) underlying this quite enigmatic pharmacological profile of 4-AHCP, we have now developed a synthesis of (S)-4-AHCP (6) and (R)-4-AHCP (7). At cloned metabotropic Glu receptors mGluR1alpha (group I), mGluR2 (group II), and mGluR4a (group III), neither 6 nor 7 showed significant agonist or antagonist effects. The stereoisomer 6, but not 7, activated cloned AMPA receptor subunits GluR1o, GluR3o, and GluR4o with EC(50) values in the range 4.5-15 microM and the coexpressed kainate-preferring subunits GluR6 + KA2 (EC(50) = 6.4 microM). Compound 6, but not 7, proved to be a very potent agonist (EC(50) = 0.13 microM) at the kainate-preferring GluR5 subunit, equipotent with (S)-2-amino-3-(5-tert-butyl-3-hydroxyisothiazol-4-yl)propionic acid [(S)-Thio-ATPA, 4] and almost 4 times more potent than (S)-2-amino-3-(5-tert-butyl-3-hydroxyisoxazol-4-yl)propionic acid [(S)-ATPA, 3]. Compound 6 thus represents a new structural class of GluR5 agonists. Molecular modeling and docking to a crystal structure of the extracellular binding domain of the AMPA subunit GluR2 has enabled identification of the probable active conformation and binding mode of 6. We are able to rationalize the observed selectivities by comparing the docking of 4 and 6 to subtype constructs, i.e., a crystal structure of the extracellular binding domain of GluR2 and a homology model of GluR5.

Alanine↗

Dibasic non-imidazole histamine H3 receptor antagonists with a rigid biphenyl scaffold.

A class of rigid, dibasic, non-imidazole H3 antagonists was developed, starting from a series of previously described flexible compounds. The original polymethylene chain between two tertiary amine groups was replaced by a rigid scaffold, composed by a phenyl ring or a biphenyl fragment. Modulation of the distance between the two amine groups, and of their alkyl substituents, was driven by superposition of molecular models and docking into a receptor model, resulting in the identification of 1,1'-[biphenyl-4,4'-diylbis(methylene)]bis-piperidine (5) as a subtype-selective H3 antagonist with high binding affinity (pKi=9.47) at human H3 histamine receptor.

Animals↗

Structural adaptability in the ligand-binding pocket of the ecdysone hormone receptor.

The ecdysteroid hormones coordinate the major stages of insect development, notably moulting and metamorphosis, by binding to the ecdysone receptor (EcR); a ligand-inducible nuclear transcription factor. To bind either ligand or DNA, EcR must form a heterodimer with ultraspiracle (USP), the homologue of retinoid-X receptor. Here we report the crystal structures of the ligand-binding domains of the moth Heliothis virescens EcR-USP heterodimer in complex with the ecdysteroid ponasterone A and with a non-steroidal, lepidopteran-specific agonist BYI06830 used in agrochemical pest control. The two structures of EcR-USP emphasize the universality of heterodimerization as a general mechanism common to both vertebrates and invertebrates. Comparison of the EcR structures in complex with steroidal and non-steroidal ligands reveals radically different and only partially overlapping ligand-binding pockets that could not be predicted by molecular modelling and docking studies. These findings offer new perspectives for the design of insect-specific, environmentally safe insecticides. The concept of a ligand-dependent binding pocket in EcR provides an insight into the moulding of nuclear receptors to their ligand, and has potential applications for human nuclear receptors.

Animals↗

Carbohydrate specificity of an insecticidal lectin isolated from the leaves of Glechoma hederacea (ground ivy) towards mammalian glycoconjugates.

Preliminary studies indicated that the potent insecticidal lectin, Gleheda, from the leaves of Glechoma hederacea (ground ivy) preferentially agglutinates human erythrocytes carrying the Tn (GalNAcalpha1-Ser/Thr) antigen. However, no details have been reported yet with respect to the fine specificity of the lectin. To corroborate the molecular basis of the insecticidal activity and physiological function of Gleheda, it is necessary to identify the recognition factors that are involved in the Gleheda-glycotope interaction. In the present study, the requirement of high-density multivalent carbohydrate structural units for Gleheda binding and a fine-affinity profile were evaluated using ELLSA (enzyme-linked lectinosorbent assay) with our extended glycan/ligand collections, a glycan array and molecular modelling. From the results, we concluded that a high-density of exposed multivalent Tn-containing glycoproteins (natural armadillo and asialo ovine salivary glycoproteins) were the most potent factors for Gleheda binding. They were, on a nanogram basis, 6.5x10(5), 1.5x10(4) and 3.1x10(3) times more active than univalent Gal (galactose), GalNAc (N-acetylgalactosamine) and Tn respectively. Among mono- and oligo-saccharides examined, simple clustered Tn (molecular mass <3000 Da) from ovine salivary glycoprotein was the best, being 37.5 and 1.7x10(3) times better than GalNAc and Gal respectively. GalNAc glycosides were significantly more active than Gal glycosides, indicating that the N-acetamido group at C-2 plays an important role in Gleheda binding. The results of glycan array support the conclusions drawn with respect to the specificity of Gleheda based on the ELLSA assays. These findings combined with the results of the molecular modelling and docking indicate the occurrence of a primary GalNAcalpha1-binding site in the Gleheda monomer. However, the extraordinary binding feature of Gleheda for glycoproteins demonstrates the importance of affinity enhancement by high-density multivalent glycotopes in the ligand-lectin interactions in biological processes.

Binding Sites↗

Identification of the cyclamate interaction site within the transmembrane domain of the human sweet taste receptor subunit T1R3.

The artificial sweetener cyclamate tastes sweet to humans, but not to mice. When expressed in vitro, the human sweet receptor (a heterodimer of two taste receptor subunits: hT1R2 + hT1R3) responds to cyclamate, but the mouse receptor (mT1R2 + mT1R3) does not. Using mixed-species pairings of human and mouse sweet receptor subunits, we determined that responsiveness to cyclamate requires the human form of T1R3. Using chimeras, we determined that it is the transmembrane domain of hT1R3 that is required for the sweet receptor to respond to cyclamate. Using directed mutagenesis, we identified several amino acid residues within the transmembrane domain of T1R3 that determine differential responsiveness to cyclamate of the human versus mouse sweet receptors. Alanine-scanning mutagenesis of residues predicted to line a transmembrane domain binding pocket in hT1R3 identified six residues specifically involved in responsiveness to cyclamate. Using molecular modeling, we docked cyclamate within the transmembrane domain of T1R3. Our model predicts substantial overlap in the hT1R3 binding pockets for the agonist cyclamate and the inverse agonist lactisole. The transmembrane domain of T1R3 is likely to play a critical role in the interconversion of the sweet receptor from the ground state to the active state.

Alanine↗

Cyclohexylcarbamic acid 3'- or 4'-substituted biphenyl-3-yl esters as fatty acid amide hydrolase inhibitors: synthesis, quantitative structure-activity relationships, and molecular modeling studies.

Fatty acid amide hydrolase (FAAH) is a promising target for modulating endocannabinoid and fatty acid ethanolamide signaling, which may have important therapeutic potential. We recently described a new class of O-arylcarbamate inhibitors of FAAH, including the cyclohexylcarbamic acid biphenyl-3-yl ester URB524 (half-maximal inhibitory concentration, IC(50) = 63 nM), which have significant anxiolytic-like properties in rats. In the present study, by introducing a selected group of substituents at the meta and para positions of the distal phenyl ring of URB524, we have characterized structure-activity profiles for this series of compounds and shown that introduction of small polar groups in the meta position greatly improves inhibitory potency. Most potent in the series was the m-carbamoyl derivative URB597 (4i, IC(50) = 4.6 nM). Furthermore, quantitative structure-activity relationship (QSAR) analysis of an extended set of meta-substituted derivatives revealed a negative correlation between potency and lipophilicity and suggested that small-sized substituents may undertake polar interactions with the binding pocket of the enzyme. Docking studies and molecular dynamics simulations, using the crystal structure of FAAH, indicated that the O-biphenyl scaffold of the carbamate inhibitors can be accommodated within a lipophilic region of the substrate-binding site, where their folded shape mimics the initial 10-12 carbon atoms of the arachidonyl moiety of anandamide (a natural FAAH substrate) and methyl arachidonyl fluorophosphonate (a nonselective FAAH inhibitor). Moreover, substituents at the meta position of the distal phenyl ring can form hydrogen bonds with atoms located on the polar section of a narrow channel pointing toward the membrane-associated side of the enzyme. The structure-activity characterization reported here should help optimize the pharmacodynamic and pharmacokinetic properties of this class of compounds.

Amidohydrolases↗

Targeting malaria with specific CDK inhibitors.

Cyclin-dependent protein kinases (CDKs) are attractive targets for drug discovery and efforts have led to the identification of novel CDK selective inhibitors in the development of treatments for cancers, neurological disorders, and infectious diseases. More recently, they have become the focus of rational drug design programs for the development of new antimalarial agents. CDKs are valid targets as they function as essential regulators of cell growth and differentiation. To date, several CDKs have been characterized from the genome of the malaria-causing protozoan Plasmodium falciparum. Our approach employs experimental and virtual screening methodologies to identify and refine chemical inhibitors of the parasite CDK Pfmrk, a sequence homologue of human CDK7. Chemotypes of Pfmrk inhibitors include the purines, quinolinones, oxindoles, and chalcones, which have sub-micromolar IC50 values against the parasite enzyme, but not the human CDKs. Additionally, we have developed and validated a pharmacophore, based on Pfmrk inhibitors, which contains two hydrogen bond acceptor functions and two hydrophobic sites, including one aromatic ring hydrophobic site. This pharmacophore has been exploited to identify additional compounds that demonstrate significant inhibitory activity against Pfmrk. A molecular model of Pfmrk designed using the crystal structure of human CDK7 highlights key amino acid substitutions in the ATP binding pocket. Molecular modeling and docking of the active site pocket with selective inhibitors has identified possible receptor-ligand interactions that may be responsible for inhibitor specificity. Overall, the unique biochemical characteristics associated with this protein, to include distinctive active site amino acid residues and variable inhibitor profiles, distinguishes the Pfmrk drug screen as a paradigm for CDK inhibitor analysis in the parasite.

Adenosine Triphosphate↗

Discovery of Helicobacter pylori shikimate kinase inhibitors: bioassay and molecular modeling.

Shikimate kinase (SK) is the fifth enzyme in the shikimate pathway and catalyzes the phosphate transfer from ATP to shikimate in generating shikimate 3-phosphate and ADP. SK has been developed as a promising target for the discovery of antibacterial agents. In this report, two small molecular inhibitors (compound 1, 3-methoxy-4-{[2-({2-methoxy-4-[(4-oxo-2-thioxo-1,3-thiazolidin-5-ylidene)methyl]phenoxy}methyl)benzyl]oxy}benzaldehyde; compound 2, 5-bromo-2-(5-{[1-(3,4-dichlorophenyl)-3,5-dioxo-4-pyrazolidinylidene]methyl}-2-furyl)benzoic acid) against Helicobacter pylori SK (HpSK) were successfully identified with IC(50) values of 5.5+/-1.2 and 6.4+/-0.4 microM, respectively. The inhibition kinetics shows that compound 1 is a noncompetitive inhibitor with respect to both shikimate and MgATP, and compound 2 is a competitive inhibitor toward shikimate and noncompetitive inhibitor with respect to MgATP. The surface plasmon resonance (SPR) technology based analysis reveals that the equilibrium dissociation constants (K(D)s) of compounds 1 and 2 with HpSK enzyme are 4.39 and 3.74 microM, respectively. The molecular modeling and docking of two inhibitors with HpSK reveals that the active site of HpSK is rather roomy and deep, forming an L-shape channel on the surface of the protein, and compound 1 prefers the corner area of L-shape channel, while compound 2 binds the short arm of the channel of SK in the binding interactions. It is expected that our current work might supply useful information for the development of novel SK inhibitors.

Algorithms↗