SEARCH · Search PubMed
Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Recent improvements in flexible docking technology may lead to a bigger role for computational methods in lead discovery. Although fast and accurate computational prediction of binding affinities for an arbitrary molecule is still beyond the limits of current methods, the docking and screening procedures can select small sets of likely lead candidates from large libraries of either commercially or synthetically available compounds.
We have developed a microfluidic device for on-chip monitoring of cellular reactions. The device consists of two primary analytical functions: control of cell transport and immobilization, and dilution of an analyte solution to generate a concentration gradient. In this device, a dam structure in parallel to the fluid flow was constructed for docking and alignment of biological cells, which allows the fragile cells to move in the microfluidic channels and to be immobilized with controllable numbers in desired locations. The cells docked on the parallel dam structure are exposed to minimal stress caused by fluidic pressure. Additionally, a network of microfluidic channels was designed to generate a concentration gradient by controlled fluid distribution and diffusive mixing. An analyte solution could be diluted to different gradients as a function of distance along the dam. We used the ATP-dependent calcium uptake reaction of HL-60 cells as a model for on-chip measurement of the threshold ATP concentration that induces significant intracellular calcium signal. The results have demonstrated the feasibility of using the microchip for real-time monitoring of cellular processes upon treatment of a concentration gradient of a test solution. The integration of cell manipulation and solution manipulation on a microchip allows the measurement of concentration-dependent biological responses within a confined microscale feature.
The unique properties of fullerenes have raised the interest of using them for biomedical applications. Within this framework, the interactions of fullerenes with proteins have been an exciting research target, yet little is known about how native proteins can bind fullerenes, and what is the nature of these interactions. Moreover, though some proteins have been shown to interact with fullerenes, up to date, no crystal structure of such complexes was obtained. Here we report docking studies aimed at examining the interactions of fullerene in two forms (C60 nonsubstituted fullerene and carboxyfullerene) with four proteins that are known to bind fullerene derivatives: HIV protease, fullerene-specific antibody, human serum albumin, and bovine serum albumin. Our work provides docking models with detailed binding pockets information, which closely match available experimental data. We further compare the predicted binding sites using a novel multiple binding site alignment method. A high similarity between the physicochemical properties and surface geometry was found for fullerene's binding sites of HIV protease and the human and bovine serum albumins.
Random conjugation of therapeutic or diagnostic payloads to targeting proteins generates functionally heterogeneous products. Conjugation of payloads to an adapter that binds to a peptide tag engineered into a targeting protein provides an alternative strategy. To progress into clinical development, an adapter/docking tag system should include humanized components and be stable in circulation. We describe here an adapter/docking tag system based on mutated fragments of human RNase I that spontaneously bind to each other and form a conjugate with a disulfide bond between complimentary cysteine residues. This self-assembled "dock and lock" system utilizes the previously described fusion C-tag, a 1-15 aa fragment of human RNase I with the R4C amino acid substitution, and a newly engineered adapter protein (Ad-C), a 21-127-aa fragment of human RNase I with the V118C substitution. Two vastly different C-tagged recombinant proteins, human vascular endothelial growth factor (VEGF) and a 254-aa long N-terminal fragment of anthrax lethal factor (LFn), retain functional activities after spontaneous conjugation of Ad-C to N-terminal or C-terminal C-tag, respectively. Ad-C modified with pegylated phospolipid and inserted into the lipid membrane of drug-loaded liposomes (Doxil) retained the ability to conjugate C-tagged proteins, yielding targeted liposomes decorated with functionally active proteins. To further optimize the system, we engineered an adapter with an additional cysteine residue at position 88 for site-specific modification, conjugated it to C-tagged VEGF, and labeled with a near-infrared fluorescent dye Cy5.5, yielding a unique functionally active probe for in vivo molecular imaging. We expect that this self-assembled "dock and lock" system will provide new opportunities for using functionally active proteins for biomedical purposes.
We have recently reported that, in contrast to the glucocorticoid receptor, the thyroid hormone receptor does not bind to hsp90 when the receptor is translated in rabbit reticulocyte lysate [Dalman, F. C., Koenig, R. J., Perdew, G. H., Massa, E., & Pratt, W. B. (1990) J. Biol. Chem. 265, 3615-3618]. All of the steroid receptors that are known to bind hsp90 are recovered in the cytosolic fraction when hormone-free cells are ruptured in hypotonic buffer. In contrast, unliganded thyroid hormone receptors and retinoic acid receptors are tightly associated with nuclear components. In this paper, we translated the human estrogen receptor and the human retinoic acid receptor in reticulocyte lysate and then immunoadsorbed the [35S]methionine-labeled translation products with the 8D3 monoclonal antibody against hsp90. The estrogen receptor is bound to hsp90, as indicated by coimmunoadsorption, but the retinoic acid receptor is not. Translation and immunoadsorption of chimeric proteins containing the DNA binding domain of one receptor and the N-terminal and COOH-terminal segments of the other show that the DNA binding finger region of the estrogen receptor is neither necessary nor sufficient for hsp90 binding. These observations suggest that there are two classes within the steroid receptor family. In one class (e.g., glucocorticoid, mineralocorticoid, sex hormone, and dioxin receptors), the receptors bind to hsp90 and remain in some kind of inactive "docking" mode until hormone-triggered release of hsp90 occurs. In the retinoic acid/thyroid hormone class, the unligated receptors do not bind to hsp90, and the receptors appear to proceed directly to their high-affinity nuclear acceptor sites without entering the "docking" state.
Cytochrome c2 (cyt) is the mobile electron donor to the reaction center (RC) in photosynthetic bacteria. The electrostatic interactions involved in the dynamics of docking of cyt onto the RC were examined by double mutant studies of the rates of electron transfer between six modified Rhodobacter sphaeroides RCs in which negatively charged acid residues were replaced with Lys and five modified Rhodobacter capsulatus Cyt c2 molecules in which positively charged Lys residues were replaced with Glu. We measured the second-order rate constant, k2, for electron transfer from the reduced cyt to the oxidized primary donor on the RC, which reflects the energy of the transition state for the formation of the active electron transfer complex. Strong interactions were found between Lys C99 and Asp M184/Glu M95, and between Lys C54 and Asp L261/Asp L257. The interacting residues were found to be located close to each other in the recently determined crystal structure of the cyt-RC complex [Axelrod, H., et al. (2002) J. Mol. Biol. (in press)]. The interaction energies were approximately inversely proportional to the distances between charges. These results support earlier suggestions [Tetreault, M., et al. (2001) Biochemistry 40, 8452-8462] that the structure of the transition state in solution resembles the structure of the cyt-RC complex in the cocrystal and indicate that specific electrostatic interactions facilitate docking of the cyt onto the RC in a configuration optimized for both binding and electron transfer. The specific interaction between Asp M184 and Lys C99 may help to nucleate short-range hydrophobic contacts.
Steroid sulfatase (STS) is a new target for the endocrine therapy of breast cancer. To ascertain some of the requirements for inhibition of estrone sulfatase activity, a number of novel analogues of estrone 3-O-sulfate possessing sulfate surrogates were synthesized and evaluated as inhibitors of estrone sulfatase (STS) in comparison to a lead inhibitor, estrone-3-O-methylthiophosphonate (E1-3-MTP). Using a selective enzyme digestion, one of the diastereoisomers of this compound, (R(p))-E1-3-MTP, could be prepared and evaluated. From structure-activity studies, we show that chirality at the phosphorus atom, hydrophobicity, basicity, size, and charge all influence the ability of a compound to inhibit estrone sulfatase activity. Of these, hydrophobicity seems to be the most important since simple, active nonsteroidal inhibitors, based on 5,6,7,8-tetrahydronaphth-2-ol (THN), can be prepared, provided that they are lipophilic enough to partition into a nonpolar environment. Also, a negatively charged group is favorable for optimal binding, although it appears that the presence of a potentially cleavable group can compensate for lack of charge in certain cases. A homology model of STS has been constructed from the STS sequence, and molecular docking studies of inhibitors have been performed to broaden the understanding of enzyme/inhibitor interactions. This model clearly shows the positions of the key amino acid residues His136, His290, Lys134, and Lys368 in the putative catalytic region of the formylglycine at position 75, with residues Asp35, Asp36, Asp342, and Gln343 as ligands in the coordination sphere of the magnesium ion. Docking studies using the substrate and estrone-3-sulfate mimics that are active inhibitors indicate they are positioned in the area of proposed catalysis, confirming the predictive power of the model.
Angiogenin (ANG) is a potent inducer of angiogenesis and an RNase A homologue whose ribonucleolytic activity is essential for its biological action. Recently, we reported the identification of small non-nucleotide inhibitors of the enzymatic activity of ANG by high-throughput screening (HTS) [Kao, R. Y. T., et al. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 10066-10071]. Two of the inhibitors that were obtained, National Cancer Institute compound NSC-65828 [8-amino-5-(4'-hydroxybiphenyl-4-ylazo)naphthalene-2-sulfonate] and ChemBridge compound C-181431 [4,4'-dicarboxy-3,3'-bis(naphthylamido)diphenylmethanone], were judged to be suitable for further development, and one of these (NSC-65828) was shown to possess antitumor activity in mice. Here we have used computational docking as a guide for the identification of available NSC-65828 and C-181431 analogues that bind more tightly to ANG, and for the characterization of inhibitor binding modes. Numerous analogues were found to have greater avidity than the HTS compounds or any small nucleotide inhibitors; four were considered to be of interest as potential leads (K(i) = 5-25 microM). Two of these analogues bind more tightly to ANG than to RNase A, and are the first small molecules shown to exhibit this selectivity. The predicted binding orientations of the HTS compounds and the new lead inhibitors were evaluated by determining the effects of ANG active site mutations on inhibitory potency. The results with ANG variants R5A, H8A, N68A, and des(121-123) are highly consistent with the docking models. Affinity changes observed with Q12A and Q117G reveal aspects of active site function that are not apparent from the free ANG crystal structure or from the modeled complexes. These findings should prove to be useful in the design of more effective and specific ANG antagonists.
Serine/threonine/tyrosine (STY) protein kinase from peanut is developmentally regulated and is induced by abiotic stresses. In addition, STY protein kinase activity is regulated by tyrosine phosphorylation. Kinetic mechanism of plant dual specificity protein kinases is not studied so far. Recombinant STY protein kinase occurs as a monomer in solution as shown by gel filtration chromatography. The relative phosphorylation rate of kinase against increasing enzyme concentrations follows a first-order kinetics indicating an intramolecular phosphorylation mechanism. Moreover, the active recombinant STY protein kinase could not transphosphorylate a kinase-deficient mutant of STY protein kinase. Molecular docking studies revealed that the tyrosine kinase inhibitors bind the protein kinase at the same region as ATP. STY protein kinase activity was inhibited by the tyrosine kinase inhibitors, and the inhibitor potency series against the recombinant STY protein kinase was tyrphostin > genistein > staurosporine. The inhibition constant (K(i)), and the IC(50) value of STY protein kinase for tyrosine kinase inhibitors with ATP and histone are discussed. All the inhibitors competed with ATP. Genistein was an uncompetitive inhibitor with histone, whereas staurosporine and tyrphostin were linear mixed type noncompetitive inhibitors with histone. Molecular docking and kinetic analysis revealed that Y148F mutant of the "ATP-binding loop" and Y297F mutant of the "activation loop" showed a dramatic increase in K(i) values for genistein and tyrphostin with respect to wild-type STY protein kinase. Data presented here provide the direct evidence on the mechanism of inhibition of plant protein kinases by tyrosine kinase inhibitors. This study also suggests that tyrosine kinase inhibitors may be useful in unraveling the plant tyrosine phosphorylation signaling cascades.
Amyloid plaques composed of the peptide Abeta are an integral part of Alzheimer's disease (AD) pathogenesis. We have modeled the process of amyloid plaque growth by monitoring the deposition of soluble Abeta onto amyloid in AD brain tissue or synthetic amyloid fibrils and show that it is mediated by two distinct kinetic processes. In the first phase, "dock", Abeta addition to the amyloid template is fully reversible (dissociation t(1/2) approximately 10 min), while in the second phase, "lock", the deposited peptide becomes irreversibly associated (dissociation t(1/2) >> 1000 min) with the template in a time-dependent manner. The most recently deposited peptide dissociates first while Abeta previously deposited becomes irreversibly "locked" onto the template. Thus, the transition from monomer to neurotoxic amyloid is mediated by interaction with the template, a mechanism that has also been proposed for the prion diseases. Interestingly, two Abeta peptides bearing primary sequence alterations implicated in heritable Abeta amyloidoses displayed faster lock-phase kinetics than wild-type Abeta. Inhibiting the initial weak docking interaction between depositing Abeta and the template is a viable therapeutic target to prevent the critical conformational transition in the conversion of Abeta((solution)) to Abeta((amyloid)) and thus prevent stable amyloid accumulation. While thermodynamics suggest that inhibiting amyloid assembly would be difficult, the present study illustrates that the protein misfolding diseases are kinetically vulnerable to intervention.
In the current work, three-dimensional QSAR studies for one large set of quinazoline type epidermal growth factor receptor (EGF-R) inhibitors were conducted using two types of molecular field analysis techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). These compounds belonging to six different structural classes were randomly divided into a training set of 122 compounds and a test set of 13 compounds. The statistical results showed that the 3D-QSAR models derived from CoMFA were superior to those generated from CoMSIA. The most optimal CoMFA model after region focusing bears significant cross-validated r(2)(cv) of 0.60 and conventional r(2) of 0.92. The predictive power of the best CoMFA model was further validated by the accurate estimation to these compounds in the external test set, and the mean agreement of experimental and predicted log(IC(50)) values of the inhibitors is 0.6 log unit. Separate CoMFA models were conducted to evaluate the influence of different partial charges (Gasteiger-Marsili, Gasteiger-Hückel, MMFF94, ESP-AM1, and MPA-AM1) on the statistical quality of the models. The resulting CoMFA field map provides information on the geometry of the binding site cavity and the relative weights of various properties in different site pockets for each of the substrates considered. Moreover, in the current work, we applied MD simulations combined with MM/PBSA (Molecular mechanics/Possion-Boltzmann Surface Area) to determine the correct binding mode of the best inhibitor for which no ligand-protein crystal structure was present. To proceed, we define the following procedure: three hundred picosecond molecular dynamics simulations were first performed for the four binding modes suggested by DOCK 4.0 and manual docking, and then MM/PBSA was carried out for the collected snapshots. The most favorable binding mode identified by MM/PBSA has a binding free energy about 10 kcal/mol more favorable than the second best one. The most favorable binding mode identified by MM/PBSA can give satisfactory explanation of the SAR data of the studied molecules and is in good agreement with the contour maps of CoMFA. The most favorable binding mode suggests that with the quinazoline-based inhibitor, the N3 atom is hydrogen-bonded to a water molecule which, in turn, interacts with Thr 766, not Thr 830 as proposed by Wissner et al. (J. Med. Chem. 2000, 43, 3244). The predicted complex structure of quinazoline type inhibitor with EGF-R as well as the pharmacophore mapping from CoMFA can interpret the structure activities of the inhibitors well and afford us important information for structure-based drug design.
The 'model-free' screening engine ProPose implements a general method for performing simultaneous protein-ligand docking, ligand-ligand alignment, pharmacophore queries-and combinations thereof-in order to incorporate a priori information into screening protocols. In this manuscript we describe a case study on herpes simplex virus thymidine kinase, an important antiviral drug target, where we evaluate different approaches for handling a specific type of a priori information, i.e., multiple target structures. We demonstrate that a simultaneous alignment on two target structures--in conjunction with logic operations on interactions and docking constraints derived from protein structure--is an effective means of (i) improving the enrichment of chemical substructures that are compatible with the a priori known ligands, (ii) ensuring the steric fit into the target protein, and (iii) handling target flexibility. The combination of ligand- and receptor-based methods steers the virtual screening by ranking molecules according to the similarity of their interaction pattern with known ligands, thereby--to some extent--outweighing the deficiencies of simple scoring functions often used in initial virtual screening.
We have developed a novel iterative approach for calculation of partial charges in proteins within the framework of the 'molecular capacitance' model. The method operates by an effective 'inductive' electronegativity scale derived from a number of the conventional charge systems including CHARMM, AMBER, MMFF, OPLS, and PEOE among others. Our novel 'inductive' electronegativity equalization procedure allows rapid and conformation sensitive computation of adequate partial charges in proteins. Accuracy of the 'inductive' values was confirmed by their correlation with DFT-computed partial charges in common amino acids. A comparative docking study with an extended steroid data set not only illustrated the adequacy of 'inductive' protein charges but also demonstrated their superior performance compared to several conventional protein charging systems. Subsequent docking with 'inductive' charges resulted in identification of five potential leads as human Sex Hormone Binding Globulin (SHBG) ligands from a commercial library of natural compounds. When the selected substances were evaluated for their ability to bind SHBG in vitro, three of them displaced testosterone from the SHBG steroid-binding site, and with one compound this was achieved at micromolar concentrations.
Traditionally, algorithms for binding site characterization or identification focus on the problem of identifying atoms within a macromolecule that might be responsible for ligand binding. In this manuscript, we focus on the binding pocket problem from a different perspective as a challenge of calculating an artificial binding pocket boundary that is sufficient to isolate binding pocket volume. The approach involves the calculation of a macromolecule encapsulating surface (MES) that separates binding pocket volume from outside space. We show that the MES can be used to increase the efficiency of flexible docking as implemented in AutoDock 3.0. The most significant improvement in docking efficiency is seen when the entire protein is searched and results show additional support for the use of AutoDock, in and of itself, as a feasible tool for binding-site identification for cases in which a protein ligand is known.
Key pharmacophoric elements for the (aminoalkyl)indole (AAI) CB1 cannabinoid receptor agonists are the aminoalkyl moiety, the lipophilic aroyl group, and the heterocyclic indole ring. In the present study, the docking space allowed for (R)-[2,3-dihydro-5-methyl-3-[(4-morpholinyl)methyl]pyrrolo[1,2,3-de]-1,4-benzoxazin-6-yl](1-naphthalenyl)methanone (WIN55212-2; 1) within the CB1 receptor was extensively explored by a docking approach that combines Monte Carlo (MC) and molecular dynamics (MD) simulations. The goals were to understand the key binding interactions of AAIs within the CB1 receptor and to examine the role of the ligand in inducing a receptor conformational change. From the findings of extensive SAR studies on the cannabinoid compounds and correlation between AAI binding affinity data and calculated binding energies, we proposed two alternative binding conformations, aroyl-up1 and aroyl-up2. These denote the directionality of the ligand naphthyl ring within the receptor upward with respect to the extracellular side. A comprehensive structural analysis of 1 demonstrated that the aroyl ring moiety could be important as the steric trigger for inducing CB1 receptor conformational change. Thus, it appears that aromatic-aromatic interactions are important not only for the binding of 1 but also for inducing receptor conformational change. It is possible that differences in the nature of the ligand binding could contribute to ligand-specific conformational changes in the receptor.
Memapsin 2 (BACE) is an aspartyl protease known as beta-secretase that acts on the production of the beta-amyloid peptide in the human brain, a key event in the pathogenesis of Alzheimer's disease. Although it is expected that the net charge of the catalytic Asp diad would be -1 as in other kinds of aspartyl proteases, the exact protonation states of Asp32 and Asp228 have not been known without ambiguity. Two independent molecular dynamics (MD) simulations of BACE in complex with the potent inhibitor OM99-2 are carried out to determine the preferred protonation state of the Asp diad in the context that is consistent with the previous X-ray crystal structure. The results show that a strong hydrogen bond between the inhibitor hydroxyl group and Asp228 can be maintained only when Asp32 is neutral and Asp228 is ionized. The preference of this protonation state is further supported from the energetic and structural features found in the docking experiment of a novel potent inhibitor with the BACE active site. Thus, both MD and docking studies suggest that the role of hydrogen bond acceptor for the hydroxyl and piperazine groups of the inhibitors should be played by Asp228 instead of Asp32. This may be a key piece of information for the structure-based design/discovery of new inhibitor drugs.
Racemic CPCCOEt ((1aRS,7aRS)-2-hydroxyimino-1a, 2-dihydro-1H-7-oxacyclopropa[b]naphthalene-7a-carboxylic acid ethyl ester, (+/-)-1) derivatives have been shown to be subtype-selective metabotropic glutamate (mGlu) 1 receptor antagonists (Annoura et al. Bioorg. Med. Chem. Lett. 1996, 6, 763-766). The optical isomers of (+/-)-1 have been separated by chromatography on a chiral stationary phase. The absolute configuration at the C-1a and C-7a positions was determined using X-ray crystallography of an amide derivative with the methyl ester of L-phenylalanine (L-PheOMe) ((+)-6). In a phosphoinositol (PI) turnover assay at the cloned human mGlu1b receptor, (-)-1 and the new amide derivatives (-)-5 and (-)-6, all of which have (1aS,7aS)-stereochemistry on the chromane ring system, showed IC(50) values of 1.5, 0.43, and 0.93 microM, respectively. In contrast, (+)-1 and the new amide derivatives (+)-5 and (+)-6were found to be inactive up to a concentration of 30 microM indicating a selectivity for the (-)-enantiomers of at least 70-fold. In a previous study (Litschig et al. Mol. Pharmacol. 1999, 55, 453-461) we demonstrated using site-directed mutagenesis that the interaction site of (+/-)-1 is located in the transmembrane (TM) domain of hmGlu1b. To suggest a plausible binding mode of (-)-1, we have built a molecular mechanics model of the putative seven TM domain of hmGlu1 based on the alpha-carbon template of the TM helices of rhodopsin. A receptor docking hypothesis suggests that the OH of T815 (TMVII) comes in close contact with the oxime OH of (-)-1 and (-)-5, whereas no such close interactions could be demonstrated by docking of (+)-1.
Poly(ADP-ribose)polymerase-1 (PARP-1) is a nuclear enzyme that has recently emerged as an important player in the mechanisms leading to postischemic neuronal death, and PARP inhibitors have been proposed as potential neuroprotective agents. With the aim of clarifying the structural basis responsible for PARP inhibition, we carried out a computational study on 46 inhibitors available through the literature. Our computational approach is composed of three parts. In the first one, representative PARP inhibitors have been docked into the crystallographic structure of the catalytic domain of PARP by using the Autodock 2.4 program. The docking studies thus carried out have provided an alignment scheme that has been instrumental for superimposing all the remaining inhibitors. Upon the basis of this alignment scheme, a quantitative structure-activity relationship (QSAR) analysis has been carried out after electrostatic and steric interaction energies have been computed with the RECEPTOR program. The QSAR analysis yielded a predictive model able to explain much of the variance of the 46-compound data set. The inspection of the QSAR coefficients revealed that the major driving force for potent inhibition is given by the extension of the contact surface between enzyme and inhibitors while electrostatic energy and hydrogen bonding capability play a minor role. Finally, the projection of the QSAR coefficients back onto the X-ray structure of the catalytic domain of PARP provides insights into the role played by specific amino acid residues. This information will be useful to address the design of new selective and potent PARP inhibitors.