SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

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 937 records · Page 52Linked to original sources

PubMed2001

DNA repair mechanism by photolyase: electron transfer path from the photolyase catalytic cofactor FADH(-) to DNA thymine dimer.

Photolyase is an enzyme that catalyses photorepair of thymine dimers in UV damaged DNA by electron transfer reaction. The structure of the photolyase/DNA complex is unknown at present. Using crystal structure coordinates of the substrate-free enzyme from E. coli, we have recently built a computer molecular model of a thymine dimer docked to photolyase catalytic site and studied molecular dynamics of the system. In this paper, we present analysis of the electronic coupling and electron transfer pathway between the catalytic cofactor FADH(-) and the pyrimidine dimer by the method of interatomic tunneling currents. Electronic structure is treated in the extended Hückel approximation. The root mean square transfer matrix element is about 6 cm(-1), which is consistent with the experimentally determined rate of transfer. We find that electron transfer mechanism responsible for the repair utilizes an unusual folded conformation of FADH(-) in photolyases, in which the isoalloxazine ring of the flavin and the adenine are in close proximity, and the peculiar features of the docked orientation of the dimer. The tunneling currents show explicitly that despite of the close proximity between the donor and acceptor complexes, the electron transfer mechanism between the flavin and the thymine bases is not direct, but indirect, with the adenine acting as an intermediate. These calculations confirm the previously made conclusion based on an indirect evidence for such mechanism.

Animals↗

PubMed1995

Molecular orbital-generated QSARs in a homologous series of alkoxyresorufins and studies of their interactive docking with P450s.

1. Molecular and electronic structural parameters have been determined, by molecular orbital (MO) calculations, for a homologous series of 8 alkoxyresorufins (methoxy- to octoxy-). 2. Quantitative structure-activity relationships (QSARs) between these structural parameters and the rates of metabolism of the alkoxyresorufins in hepatic microsomes from the 3-methylcholanthrene (MC)-, and phenobarbital (PB)-pretreated mouse, and the beta-naphthoflavone (beta NF)-pretreated rat have been established. 3. The most significant single relationship is between beta NF-induction of cytochrome P4501 (CYP1A) and the total nucleophilic superdelocalizability (sigma SN) for the eight compounds in the series. 4. For double regressions, the electronic charge on the alkoxy oxygen, Q(O), or alpha-carbon Q(C), is important when combined with the hydrophobic substituent constant (pi). 5. These findings indicate that the rates of metabolism of these alkoxyresorufins are dependent upon their ability to cross cellular membranes, to fit the relevant CYP1A binding site, and on their ability to accept electrons from a donor nucleophilic species. 6. A different set of parameters correlated with CYP2B activity, namely, parameters of overall shape, which indicates that the way in which the alkoxyresorufins fit the CYP2B site, determines their differences in specificity. 7. Computer graphic interactive docking studies of the alkoxyresorufins with their affinity-specific cytochromes P450, namely, methoxy- with CYP1A2; ethoxy- with CYP1A1; pentoxy- with CYP2B1; and benzyloxy- with CYP3A, have also been undertaken to show the specific interactions of the alkoxyresorufins with the binding sites of the individual P450s.

Amino Acid Sequence↗

PubMed2006

Indolyl aryl sulphones as HIV-1 non-nucleoside reverse transcriptase inhibitors: synthesis, biological evaluation and binding mode studies of new derivatives at indole-2-carboxamide.

New non-nucleoside reverse transcriptase inhibitors (NNRTIs) that are active against the commonly occurring mutations of HIV are urgently needed for the treatment of AIDS. We synthesized new NNRTIs of the indolyl aryl sulphone (IAS) family, which are endowed with high antiviral potency against HIV-1 wt (wild-type), and the Y181C and K103N-Y181C drug resistant mutant strains. Several new compounds were highly active in lymphocytes infected with primary isolates carrying the K103N-V1081-M184V and L1001-V1081 mutations. The design of new IASs was based on three-dimensional quantitative structure-activity relationship (3D QSAR) studies and docking simulations. A cross-docking study was also undertaken to gain some insights in to the binding mode of the newly synthesized IASs in the wt and mutated isoforms of reverse transcriptase.

Cell Line↗

PubMed2004

CORCEMA refinement of the bound ligand conformation within the protein binding pocket in reversibly forming weak complexes using STD-NMR intensities.

We describe an intensity-restrained optimization procedure for refining approximate structures of ligands within the protein binding pockets using STD-NMR intensity data on reversibly forming weak complexes. In this approach, the global minimum for the bound-ligand conformation is obtained by a hybrid structure refinement method involving CORCEMA calculation of intensities and simulated annealing optimization of torsion angles of the bound ligand using STD-NMR intensities as experimental constraints and the NOE R-factor as the pseudo-energy function to be minimized. This method is illustrated using simulated STD data sets for typical carbohydrate and peptide ligands. Our procedure also allows for the optimization of side chain torsion angles of protein residues within the binding pocket. This procedure is useful in refining and improving initial models based on crystallography or computer docking or other algorithms to generate models for the bound ligand (e.g., a lead compound) within the protein binding pocket compatible with solution STD-NMR data. This method may facilitate structure-based drug design efforts.

Algorithms↗

PubMed2003

Membrane-docking loops of the cPLA2 C2 domain: detailed structural analysis of the protein-membrane interface via site-directed spin-labeling.

C2 domains are protein modules found in numerous eukaryotic signaling proteins, where their function is to target the protein to cell membranes in response to a Ca(2+) signal. Currently, the structure of the interface formed between the protein and the phospholipid bilayer is inaccessible to high-resolution structure determination, but EPR site-directed spin-labeling can provide a detailed medium-resolution view of this interface. To apply this approach to the C2 domain of cytosolic phospholipase A(2) (cPLA(2)), single cysteines were introduced at all 27 positions in the three Ca(2+)-binding loops and labeled with a methanethiosulfonate spin-label. Altogether, 24 of the 27 spin-labeled domains retained Ca(2+)-activated phospholipid binding. EPR spectra of these 24 labeled domains obtained in the presence and absence of Ca(2+) indicate that Ca(2+) binding triggers subtle changes in the dynamics of two localized regions within the Ca(2+)-binding loops: one face of the loop 1 helix and the junction between loops 1 and 2. However, no significant changes in loop structure were detected upon Ca(2+) binding, nor upon Ca(2+)-triggered docking to membranes. EPR depth parameters measured in the membrane-docked state allow determination of the penetration depth of each residue with respect to the membrane surface. Analysis of these depth parameters, using an improved, generalizable geometric approach, provides the most accurate picture of penetration depth and angular orientation currently available for a membrane-docked peripheral protein. Finally, the observation that Ca(2+) binding does not trigger large rearrangements of the membrane-docking loops favors the electrostatic switch model for Ca(2+) activation and disfavors, or places strong constraints on, the conformational switch model.

Amino Acid Substitution↗

PubMed2001

Selective inhibition of 6-phosphogluconate dehydrogenase from Trypanosoma brucei.

A number of triphenylmethane derivatives have been screened against 6-phosphogluconate dehydrogenase from Trypanosoma brucei and sheep liver. Some of these compounds show good inhibition of the enzymes and also selectivity towards the parasite enzyme. Modelling was undertaken to dock the compounds into the active sites of both enzymes. Using a combination of DOCK 3.5 and FLEXIDOCK a correlation was obtained between docking score and both activity for the enzymes and selectivity. Visualisation of the docked structures of the inhibitors in the active sites of the enzymes yielded a possible explanation of the selectivity for the parasite enzyme.

Amino Acid Sequence↗

PubMed2004

GAsDock: a new approach for rapid flexible docking based on an improved multi-population genetic algorithm.

Based on an improved multi-population genetic algorithm, a new fast flexible docking program, GAsDock, was developed. The docking accuracy, screening efficiency, and docking speed of GAsDock were evaluated by the docking results of thymidine kinase (TK) and HIV-1 reverse transcriptase (RT) enzyme with 10 available inhibitors of each protein and 990 randomly selected ligands. Nine of the ten known inhibitors of TK were accurately docked into the protein active site, the root-mean-square deviation (RMSD) values between the docking and X-ray crystal structures are less than 1.7A; binding poses (conformation and orientation) of 9 of the 10 known inhibitors of RT were reproduced by GAsDock with RMSD values less than 2.0A. The docking time is approximately in proportion to the number of rotatable bonds of ligands; GAsDock can finish a docking simulation within 60s for a ligand with no more than 20 rotatable bonds. Results indicate that GAsDock is an accurate and remarkably faster docking program in comparison with other docking programs, which is applausive in the application of virtual screening.

Algorithms↗

PubMed2002

Brownian dynamics simulations of the recognition of the scorpion toxin P05 with the small-conductance calcium-activated potassium channels.

The recognition of the scorpion toxin P05 and the small-conductance, calcium-activated potassium (SK) channels, rsk1, rsk2, and rsk3, has been studied by means of the Brownian dynamics (BD) method. All of the 25 available structures of P05 in the RCSB Protein Data Bank determined by NMR were considered during the simulation, which indicated that the conformation of P05 affects both the recognition and binding between the two proteins significantly. Comparing the top four high-frequency structures of P05 binding to the SK channels, we found that the rsk2 channel, with high frequencies and lowest electrostatic interaction energies (E (int)(ES)), is the most favorable for P05 binding, while rsk3 is intermediate, and rsk1 is the least favorable. Among the 25 structures of P05, the 13th structure docks into the binding site of the rsk2 channel with the highest probability and most favorable electrostatic interactions. From the P05-rsk2 channel binding model, we identified the residues critical for the recognition of these two proteins through triplet contact analyses. P05 locates around the extracellular mouth of the SK channels and contacts the SK channels using its alpha-helix rather than beta-sheets. The critical triplet contacts for recognition between P05 and the rsk2 channel are Arg6 (P05)-Asp364 (SK), Arg7 (P05)-Asn368 (SK), and Arg13 (P05)-Asp341 (SK). The structure of the P05-rsk2 complex with the most favorable electrostatic interaction energy was further refined by molecular mechanics, showing that six hydrogen bonding interactions exist between P05 and the rsk2 channel: one hydrogen bond is formed between Arg6 (P05) and Asp364(D) (rsk2); Arg7 (P05) forms three hydrogen bonds with Asp341(B) (rsk2)) and Asp364(C) (rsk2); two hydrogen bonds are formed by Arg13 (P05) with Asp341(A) (rsk2) and Asp364(B) (rsk2). The simulation results are in good agreement with the previous molecular biological experiments and can explain the binding phenomena between P05 and SK channels at the level of molecular structure. The consistency between the results of the BD simulations and the experimental data indicated that our 3D model of the P05-rsk2 channel complex is reasonable and can be employed in further biological studies, such as rational design of the novel therapeutic agents blocking the small-conductance, calcium-activated and apamin-sensitive potassium channels, and for mutagenesis studies in both toxins and SK channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp364 of the rsk2 channel is critical for its recognition and binding functionality towards P05. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating that this might be a new clue for further functional study of SK channels.

Amino Acid Sequence↗

PubMedDate not supplied

3D QSAR studies on protein tyrosine phosphatase 1B inhibitors: comparison of the quality and predictivity among 3D QSAR models obtained from different conformer-based alignments.

A set of 65 flexible peptidomimetic competitive inhibitors (52 in the training set and 13 in the test set) of protein tyrosine phosphatase 1B (PTP1B) has been used to compare the quality and predictive power of 3D quantitative structure-activity relationship (QSAR) comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models for the three most commonly used conformer-based alignments, namely, cocrystallized conformer-based alignment (CCBA), docked conformer-based alignment (DCBA), and global minima energy conformer-based alignment (GMCBA). These three conformers of 5-[(2S)-2-({(2S)-2-[(tert-butoxycarbonyl)amino]-3-phenylpropanoyl}amino)3-oxo-3-pentylamino)propyl]-2-(carboxymethoxy)benzoic acid (compound number 66) were obtained from the X-ray structure of its cocrystallized complex with PTP1B (PDB ID: 1JF7), its docking studies, and its global minima by simulated annealing. Among the 3D QSAR models developed using the above three alignments, the CCBA provided the optimal predictive CoMFA model for the training set with cross-validated r2 (q2)=0.708, non-cross-validated r2=0.902, standard error of estimate (s)=0.165, and F=202.553 and the optimal CoMSIA model with q2=0.440, r2=0.799, s=0.192, and F=117.782. These models also showed the best test set prediction for the 13 compounds with predictive r2 values of 0.706 and 0.683, respectively. Though the QSAR models derived using the other two alignments also produced statistically acceptable models in the order DCBA>GMCBA in terms of the values of q2, r2, and predictive r2, they were inferior to the corresponding models derived using CCBA. Thus, the order of preference for the alignment selection for 3D QSAR model development may be CCBA>DCBA>GMCBA, and the information obtained from the CoMFA and CoMSIA contour maps may be useful in designing specific PTP1B inhibitors.

Chemistry, Pharmaceutical↗

PubMed2001

Mutagenesis and computer modelling approach to study determinants for recognition of signal peptides by the mitochondrial processing peptidase.

Determinants for the recognition of a mitochondrial presequence by the mitochondrial processing peptidase (MPP) have been investigated using mutagenesis and bioinformatics approaches. All plant mitochondrial presequences with a cleavage site that was confirmed by experimental studies can be grouped into three classes. Two major classes contain an arginine residue at position -2 or -3, and the third class does not have any conserved arginines. Sequence logos revealed loosely conserved cleavage motifs for the first two classes but no significant amino acid conservation for the third class. Investigation of processing determinants for a class III precursor, Nicotiana plumbaginifolia F1beta precursor of ATP synthase (pF1beta), was performed using a series of pF1beta presequence mutants and mutant presequence peptides derived from the C-terminal portion of the presequence. Replacement of -2 Gln by Arg inhibited processing, whereas replacement of either the most proximally located -5 Arg or -15 Arg by Leu had only a low inhibitory effect. The C-terminal portion of the pF1beta presequence forms a helix-turn-helix structure. Mutations disturbing or prolonging the helical element upstream of the cleavage site inhibited processing significantly. Structural models of potato MPP and the C-terminal pF1beta presequence peptide were built by homology modelling and empirical conformational energy search methods, respectively. Molecular docking of the pF1beta presequence peptide to the MPP model suggested binding of the peptide to the negatively charged binding cleft formed by the alpha-MPP and beta-MPP subunits in close proximity to the H111XXE114H115X(116-190)E191 proteolytic active site on beta-MPP. Our results show for the first time that the amino acid at the -2 position, even if not an arginine, as well as structural properties of the C-terminal portion of the presequence are important determinants for the processing of a class III precursor by MPP.

Amino Acid Sequence↗

PubMed2006

Molecular modeling of purinergic receptor P2Y12 and interaction with its antagonists.

Purinergic receptors are a class of cell surface receptors for purines that prefer ATP or ADP over adenosine. The surface receptors for extracellular nucleotides are called P2 receptors. They are activated by both pyrimidine and purine nucleotides. ADP initiates platelet aggregation by 'simultaneous activation of two G protein-coupled receptors, P2Y1 and P2Y12. P2Y12 has been shown to be the target of the thienopyridine drugs, ticlopidine and clopidogrel. Here, the active sites of P2Y12 for ATP as well as ADP are predicted by bioinformatics and molecular modeling. First, the three-dimensional (3D) structure of P2Y12 was constructed by InsightII/Homology module using the corresponding bovine rhodopsin (PDB code: 1HZX) as the template. Then the primary structures were optimized by energy minimization that has been successfully accepted by the Protein Data Bank (PDB code: 1VZ1). Second, a simple scoring matrix was built up based on the analysis of 13 known ATP-binding proteins. And the most probable active sites of P2Y12 were predicted using the scoring matrix, which include three distant areas: "head area" (LGTGPLRTFV, 87-96), "middle area" (VGLITNGLAM, 38-47, and LGAKILSVVI, 139-148), and "bottom area" (RTRGVGKVPR, 222-231). Subsequently the structural model of P2Y12 was docked with ATP/ADP in comparison with P2Y1 (PDB code 1ddd). As a comparison, we docked its antagonists, such as ticlopidine and clopidogrel, to the most probable sites and calculated their intermolecular energy. Our results imply that P2Y12 has the potential to be inhibited by ADP/ATP analogs, and it suggests that P2Y12 acts as a target of new drugs that inhibit platelet aggregation.

Adenosine Diphosphate↗

PubMed2003

Design and synthesis of orally bioavailable inhibitors of inducible nitric oxide synthase. synthesis and biological evaluation of dihydropyridin-2(1H)-imines and 1,5,6,7-tetrahydro-2H-azepin-2-imines.

The process of discovery and biological evaluation of alpha,beta-unsaturated cyclic amidines, as selective inhibitors of inducible nitric oxide synthase (iNOS), is reported. Dihydropyridin-2(1H)-imines and 1,5,6,7-tetrahydro-2H-azepin-2-imines were synthesized and biologically evaluated both in vitro and in vivo using a nitric oxide synthase inhibition assay. Compounds 1, 5, 6, 8-12 and 16 exhibited potent inhibition of iNOS. Among these, compounds 6, 7, 10, 11 and 16 showed 5- to 19-fold isoform selectivity. Compounds 1, 6, 10, 11 and 16 also showed potent inhibitory activity in the NOx accumulation assay in mice. Compounds 1 and 6 showed excellent bioavailability (BA) in rats when administered orally. Full details are presented here, including the structure-activity relationship (SAR) studies, the chemistry of these compounds, and the pharmacokinetic data and the computer-aided docking study of 10 with hiNOS.

Animals↗

PubMed2003

A molecular modeling study of inhibitors of nuclear factor kappa-B (p50)--DNA binding.

Nuclear Factor-kappa B (NF-kappaB) is an inducible transcription factor of the Rel family, and is sequestered in the cytoplasm by the IkappaB family of proteins. NF-kappaB can exist in several dimeric forms, but the p50/p65 heterodimer is the predominant one. Activation of NF-kappaB by a range of stimuli including viral products, and oxidative stress, leads to phosphorylation and proteasome dependent degradation of IkappaB, leading to the release of free NF-kappaB. This free NF-kappaB then binds to its target sites (KB sites in the DNA) to initiate transcription. These kappaB sites are also present in the Long Terminal Repeat (LTR) of HIV-1, and hence NF-kappaB (p50 subunit) binding to LTR-DNA is critical in viral replication. Targeting direct p50-DNA binding, in this regard, is a novel approach to design anti-HIV gene expression inhibitors, which do not have the problem of resistance unlike in other anti-HIV strategies. The present study is a part of our search for leads for the specific inhibition of p50-DNA binding. We have been experimentally studying different types of these inhibitors, and in this work, we attempted to get a common definition of their structural mechanism onto p50-DNA binding. Using three different classes of inhibitors, we modelled their association with the DNA-Binding Region (DBR) of the p50 subunit of NF-kappaB. Docking studies were carried out using a genetic algorithm based program (GOLD). Further, to compare electrostatic complementarity in the association of the inhibitors with the DBR, Molecular Electrostatic Potentials (MEPs) were generated for the DBR and each inhibitor. The results of docking revealed a strong network of hydrogen bonding interactions for every active inhibitor, and the contrary for the less active ones. Further, the MEPs revealed that the DBR of p50 represents a surface of electropositive potential, and the active inhibitors represent a complementary electronegative surface. With the present modelling study we conclude that the principal properties to be possessed by the new leads against p50-DNA binding should be that of having the ability to make a strong network of hydrogen bonds with the DBR of p50, and preferably, having electronegative potentials in their peripheral surface.

Computer Simulation↗

PubMed2006

Docking of human interleukin-15 to its specific receptor alpha chain: correlation between molecular modeling and mutagenesis experimental data.

A structural model of the sushi domain of IL-15Ralpha was first obtained by homology modeling to study its interactions with IL-15 by means of molecular modeling, peptide scanning, and site-directed mutagenesis. From these experimental data, a putative interacting surface of IL-15Ralpha with a previously published IL-15 model was inferred: Leu25, Leu44, and Glu46 of IL-15 and Arg35 of IL-15Ralpha were found to be key interfacial residues and were subsequently used as filters for the construction of docking solutions. Human IL-15/IL-15Ralpha complexes were constructed in two stages, with a preliminary docking procedure, treating the two partners as rigid bodies and using these filters. In this first stage, two classes of docking solutions were characterized. From a topological point of view, each solution could be derived from the other by reverse orientation of one partner in relation to the other. In a second stage, several further energy refinements clearly favored one solution. Moreover, this unique docking solution was confirmed by molecular modeling of IL-15 mutants previously built and tested in our laboratory. Finally, this complex model, which is a useful tool to study the IL-15/IL-15Ralpha interface, was topologically compared to IL-2/IL-2Ralpha complexes (previous model in the literature and recent crystal structure).

Amino Acid Sequence↗

PubMed2003

Allosteric site on muscarinic acetylcholine receptors: identification of two amino acids in the muscarinic M2 receptor that account entirely for the M2/M5 subtype selectivities of some structurally diverse allosteric ligands in N-methylscopolamine-occupied receptors.

Two epitopes have been identified recently to be responsible for the high-affinity binding of alkane-bisammonium and caracurine V type allosteric ligands to N-methylscopolamine (NMS)-occupied M2 muscarinic acetylcholine receptors, relative to M5 receptors: the amino acid M2-Thr423 at the top of transmembrane region (TM) 7 and an epitope comprising the second extracellular loop (o2) of the M2 receptor including the flanking regions of TM4 and TM5. We aimed to find out whether a single amino acid could account for the contribution of this epitope to binding affinity. Allosteric interactions were investigated in wild-type and mutant receptors in which the orthosteric binding site was occupied by [3H]NMS (5 mM Na,K,Pi buffer, pH 7.4, 23 degrees C). Using M2/M5 chimeric and point-mutated receptors, the relevant epitope was narrowed down to M2-Tyr177. A double point-mutated M2 receptor in which both M2-Tyr177 and M2-Thr423 were replaced by the corresponding amino acids of M5 revealed that these two amino acids account entirely for the (approximately 100-fold) M2/M5 selectivity of the alkane-bisammonium and the caracurine V type allosteric ligands. At NMS-free M2 receptors, the caracurine V derivative also displayed approximately 100-fold M2/M5 selectivity, but the double point mutation reduced the M2 affinity by only approximately 10-fold; thus, additional epitopes may influence selectivity for the free receptors. A three-dimensional model of the M2 receptor was used to simulate allosteric agent docking to NMS-occupied receptors. M2-Tyr177 and M2-Thr423 seem to be located near the junction of the allosteric and the orthosteric areas of the M2 receptor ligand binding cavity.

Allosteric Regulation↗

PubMed2001

Theoretical models of catalytic domains of protein phosphatases 1 and 2A with Zn2+ and Mn2+ metal dications and putative bioligands in their catalytic centers.

The oligomeric metalloenzymes protein phosphatases dephosphorylate OH groups of Ser/Thr or Tyr residues of proteins whose actions depend on the phosphorus signal. The catalytic units of Ser/Thr protein phosphatases 1, 2A and 2B (PP1c, PP2Ac and PP2Bc, respectively), which exhibit about 45% sequence similarity, have their active centers practically identical. This feature strongly suggests that the unknown structure of PP2Ac could be successfully homology-modeled from the known structures of PP1c and/or PP2Bc. Initially, a theoretical model of PP1c was built, including a phosphate and a metal dication in its catalytic site. The latter was modeled, together with a structural hydroxyl anion, as a triangular pseudo-molecule (Zno or Mno), composed of two metal cations (double Zn2+ or Mn2+, respectively) and the OH- group. To the free PP1c two inhibitor sequences R29RRRPpTPAMLFR40 of DARPP-32 and R30RRRPpTPATLVLT42 of Inhibitor-1, and two putative substrate sequences LRRApSVA and QRRQRKpRRTI were subsequently docked. In the next step, a free PP2Ac model was built via homology re-modeling of the PP1c template and the same four sequences were docked to it. Thus, together, 20 starting model complexes were built, allowing for combination of the Zno and Mno pseudo-molecules, free enzymes and the peptide ligands docked in the catalytic sites of PP1c and PP2Ac. All models were subsequently subjected to 250-300 ps molecular dynamics using the AMBER 5.0 program. The equilibrated trajectories of the final 50 ps were taken for further analyses. The theoretical models of PP1c complexes, irrespective of the dication type, exhibited increased mobilities in the following residue ranges: 195-200, 273-278, 287-209 for the inhibitor sequences and 21-25, 194-200, 222-227, 261, 299-302 for the substrate sequences. Paradoxically, the analogous PP2Ac models appeared much more stable in similar simulations, since only their "prosegment" residues 6-10 and 14-18 exhibited an increased mobility in the inhibitor complexes while no areas of increased mobility were found in the substrate complexes. Another general observation was that the complexes with Mn dications were more stable than those with Zn dications for both PP1c and PP2Ac units.

Animals↗

PubMed2004

Predicting the three-dimensional structure of the human facilitative glucose transporter glut1 by a novel evolutionary homology strategy: insights on the molecular mechanism of substrate migration, and binding sites for glucose and inhibitory molecules.

The glucose transporters (GLUT/SLC2A) are members of the major facilitator superfamily. Here, we generated a three-dimensional model for Glut1 using a two-step strategy: 1), GlpT structure as an initial homology template and 2), evolutionary homology using glucose-6-phosphate translocase as a template. The resulting structure (PDB No. 1SUK) exhibits a water-filled passageway communicating the extracellular and intracellular domains, with a funnel-like exofacial vestibule (infundibulum), followed by a 15 A-long x 8 A-wide channel, and a horn-shaped endofacial vestibule. Most residues which, by mutagenesis, are crucial for transport delimit the channel, and putative sugar recognition motifs (QLS, QLG) border both ends of the channel. On the outside of the structure there are two positively charged cavities (one exofacial, one endofacial) delimited by ATP-binding Walker motifs, and an exofacial large side cavity of yet unknown function. Docking sites were found for the glucose substrate and its inhibitors: glucose, forskolin, and phloretin at the exofacial infundibulum; forskolin, and phloretin at an endofacial site next to the channel opening; and cytochalasin B at a positively charged endofacial pocket 3 A away from the channel. Thus, 1SUK accounts for practically all biochemical and mutagenesis evidence, and provides clues for the transport process.

Binding Sites↗

PubMed2005

Protein-Protein Docking Benchmark 2.0: an update.

We present a new version of the Protein-Protein Docking Benchmark, reconstructed from the bottom up to include more complexes, particularly focusing on more unbound-unbound test cases. SCOP (Structural Classification of Proteins) was used to assess redundancy between the complexes in this version. The new benchmark consists of 72 unbound-unbound cases, with 52 rigid-body cases, 13 medium-difficulty cases, and 7 high-difficulty cases with substantial conformational change. In addition, we retained 12 antibody-antigen test cases with the antibody structure in the bound form. The new benchmark provides a platform for evaluating the progress of docking methods on a wide variety of targets. The new version of the benchmark is available to the public at http://zlab.bu.edu/benchmark2.

Algorithms↗