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 685 records · Page 38Linked to original sources

PubMed2001

Theoretical evidence of a salt bridge disruption as the initiating process for the alpha1d-adrenergic receptor activation: a molecular dynamics and docking study.

This study reports the building of the three-dimensional structure of the rat alpha1d-adrenergic receptor through a topology approach based on the structure of the rhodopsin receptor from cryoelectron microscopy. The validity and reliability of the receptor model were assessed through exhaustive molecular dynamics and docking studies. Some interesting ligand-receptor interactions were identified along with significant differences between the binding mode of agonists and antagonists. The importance of the disruption of a salt bridge as a possible initial event leading to receptor activation is discussed on the basis of data from mutagenesis and molecular dynamics studies.

Amino Acid Sequence↗

PubMed2005

Hypocrea jecorina (Trichoderma reesei) Cel7A as a molecular machine: A docking study.

Hypocrea jecorina (formerly Trichoderma reesei) Cel7A has a catalytic domain (CD) and a cellulose-binding domain (CBD) separated by a highly glycosylated linker. Very little is known of how the 2 domains interact to degrade crystalline cellulose. Based on the interaction energies and forces on cello-oligosaccharides computationally docked to the CD and CBD, we propose a molecular machine model, where the CBD wedges itself under a free chain end on the crystalline cellulose surface and feeds it to the CD active site tunnel. Enzyme-substrate interactions produce the forces required to pull cellulose chains from the surface and also to help the enzyme move on the cellulose chain for processive hydrolysis. The energy to generate these forces is ultimately derived from the chemical energy of glycosidic bond breakage.

Binding Sites↗

PubMed1997

The cannabinoid receptor: computer-aided molecular modeling and docking of ligand.

A three-dimensional model of human cannabinoid receptor is constructed using computer-aided molecular modeling techniques. The helices of bacteriorhodopsin were used as the initial template to construct the transmembrane helices. The extracellular and intracellular loops were added using the SYBYL molecular modeling package. The extracellular N terminus was modeled on the basis of its similarity to rat oncomodulin. Similarly, the C terminus was constructed on the basis of similarity to bovine prothrombin fragment 1. The final structure was refined by several runs of minimization and dynamics calculation using the CHARMm package. delta 9-Tetra hydrocannabinol was docked into the internal cavity using the AUTODOCK program. Our study snows that there may be a calcium-binding site in the extracellular N terminus of this receptor. The ligand binds mainly to a hydrophobic site, which consists of residues Met-240, Trp-241 (TMH-4), Trp-356, Leu-359, Leu-360 (TMH-6), and Ala-283 (TMH-5). Its phenolic hydroxyl group forms a hydrogen bond with the carboxy group of Ala-198 (TMH-3). The results of modeling agree well with experimental QSAR studies.

Amino Acid Sequence↗

PubMed2002

Inhibition of trypsin by cowpea thionin: characterization, molecular modeling, and docking.

Higher plants produce several families of proteins with toxic properties, which act as defense compounds against pests and pathogens. The thionin family represents one family and comprises low molecular mass cysteine-rich proteins, usually basic and distributed in different plant tissues. Here, we report the purification and characterization of a new thionin from cowpea (Vigna unguiculata) with proteinase inhibitory activity. Cowpea thionin inhibits trypsin, but not chymotrypsin, binding with a stoichiometry of 1:1 as shown with the use of mass spectrometry. Previous annotations of thionins as proteinase inhibitors were based on their erroneous identification as homologues of Bowman-Birk family inhibitors. Molecular modeling experiments were used to propose a mode of docking of cowpea thionin with trypsin. Consideration of the dynamic properties of the cowpea thionin was essential to arrive at a model with favorable interface characteristics comparable with structures of trypsin-inhibitor complexes determined by X-ray crystallography. In the final model, Lys11 occupies the S1 specificity pocket of trypsin as part of a canonical style interaction.

Amino Acid Sequence↗

PubMed1992

A robust and efficient automated docking algorithm for molecular recognition.

A completely automated method is described for determining the most likely mode of binding of two (macro)molecules from the knowledge of their three-dimensional structures alone. The method is based on well-known graph theoretical techniques and has been used successfully to determine and rationalize the binding of a number of known macromolecular complexes. In this article we present results for a special case of the general molecular recognition problem--given the information concerning the particular atoms involved in the binding for one of the molecules, the algorithm can correctly identify the corresponding (contacting) atoms of the other molecule. The approach used can be easily extended to the general molecular recognition problem and requires the extraction of maximal common subgraphs. In these studies the docking of the macromolecules was achieved without the aid of computer graphics or other visual aids. The algorithm has been used to determine the correct mode of binding of a protein antigen to an antibody in approximately 100 min on a DEC micro VAX 3600.

Algorithms↗

PubMed2002

Regioselectivity of CYP2B6: homology modeling, molecular dynamics simulation, docking.

Human cytochrome P450 (CYP) 2B6 activates the anticancer prodrug cyclophosphamide (CPA) by 4-hydroxylation. In contrast, the same enzyme catalyzes N-deethylation of a structural isomer, the prodrug ifosfamide (IFA), thus causing severe adverse drug effects. To model the molecular interactions leading to a switch in regioselectivity, the structure of CYP2B6 was modeled based on the structure of rabbit CYP2C5. We modeled the missing 22-residue loop in CYP2C5 between helices F and G (the F-G loop), which is not resolved in the X-ray structure, by molecular dynamics (MD) simulations using a simulated annealing protocol. The modeled conformation of the loop was validated by unconstrained MD simulations of the complete enzymes (CYP2C5 and CYP2B6) in water for 70 and 120 ps, respectively. The simulations were stable and led to a backbone r.m.s. deviation of 1.7 A between the two CYPs. The shape of the substrate binding site of CYP2B6 was further analyzed. It consists of three well-defined hydrophobic binding pockets adjacent to the catalytic heme. Size, shape and hydrophobicity of these pockets were compared to the shapes of the two structurally isomeric substrates. In their preferred orientation in the binding site, both substrates fill all three binding pockets without repulsive interactions. The distance to the heme iron is short enough for 4-hydroxylation and N-deethylation to occur for CPA and IFA, respectively. However, if the substrates are docked in the non-preferred orientation (such that 4-hydroxylation and N-deethylation would occur for IFA and CPA, respectively), one pocket is left empty, and clashes were observed between the substrates and the enzyme.

Amino Acid Sequence↗

PubMed2003

Plausible interaction of an alpha-fetoprotein cyclopeptide with the G-protein-coupled receptor model GPR30: docking study by molecular dynamics simulated annealing.

In this manuscript, the procedure of molecular dynamics simulated annealing is applied to locate a probable receptor and binding site of a cyclicpeptide that inhibits estrogen-stimulated proliferation of breast cancer. The hydrophilic cyclopeptide EMTOVNOGQ (O = 4-hydroxyproline), derived from alpha-fetoprotein, is an inhibitor of estrogen-stimulated proliferation of human breast cancer. This peptide has been shown to act through a mechanism different from that of estrogen; however, its receptor is unknown. We report computer experiments that suggest that this peptide may execute its actions by interacting with GPR30, a G-protein-coupled receptor. The subject of this work is the simulation, by molecular dynamics simulated annealing, of the interaction of cyclopeptide EMTOVNOGQ with receptor GPR30 protein. A conformational analysis of the cyclopeptide was undertaken and the final structure was docked on several sites of the GPR30 3D model. Our results show that the cyclopeptide interacts on the pocket located between TM6 and TM7 transmembrane helices of the G-protein, triggering a slight conformational change in the secondary structure of the receptor in the complex. Based on differences in accessible surface areas between GPR30 and its ligand, the residues in the interaction zone were identified. The cyclopeptide is stabilized in the active site by forming a network of hydrogen bonds between Glu, Thr, (1)Pro(OH) and GLn residues of the ligand and Arg-259, Cys-271, Asn-316, Asn-320 and Tyr-324 of the G-protein. Moreover, the study of the electrostatic surface potential on the GPR30 receptor shows that the active site is more positively charged than the other sites. Our modeling indicates a plausible interaction of the cyclopeptide with the seven transmembrane GPR30 protein. This may have profound implications for the treatment of breast cancer.

Amino Acid Sequence↗

PubMed2005

Toward selective histone deacetylase inhibitor design: homology modeling, docking studies, and molecular dynamics simulations of human class I histone deacetylases.

Histone deacetylases (HDACs) play an important role in gene transcription. Inhibitors of HDACs induce cell differentiation and suppress cell proliferation in tumor cells. Although many HDAC inhibitors have been designed and synthesized, selective inhibition for class I HDAC isoforms is a goal that has yet to be achieved. To understand the difference between class I HDAC isoforms that could be exploited for the design of isoform-specific HDAC inhibitors, we have built three-dimensional models of four class I histone deacetylases, HDAC1, HDAC2, HDAC3, and HDAC8. Comparison of the homology model of HDAC8 with the recently published X-ray structure shows excellent agreement and validates the approach. A series of HDAC inhibitors were docked to the homology models to understand the similarities and differences between the binding modes. Molecular dynamic simulations of these HDAC-inhibitor complexes indicate that the interaction between the protein surface and inhibitor is playing an important role; also some active site residues show some flexibility, which is usually not included in routine docking protocols. The implications of these results for the design of isoform-selective HDAC inhibitors are discussed.

Amino Acid Sequence↗

PubMed2006

Enantioselective substrate binding in a monooxygenase protein model by molecular dynamics and docking.

The two-component flavoenzyme styrene monooxygenase (SMO) is an efficient alternative to several chemical epoxidation catalysts on a preparative scale. A first homology model of the catalytic domain (StyA) of SMO was constructed (Protein Data Bank ID 2HD8) based on the structure of para-hydroxybenzoate hydroxylase. The StyA protein structure was optimized by restrained molecular dynamics to reproduce specific pre-S binding orientations of styrene. Effects of all 10 point mutations examined were explained by the distance of the site to the styrene and FAD binding sites. Thirteen of 20 ligands could be accommodated in a catalytically active binding orientation, and predicted affinities correlated well with experimental turnover and inhibition. The binding cavity is almost completely hydrophobic except for a hydrogen-bonded network formed by three water molecules, the backbone of residues 300-302, and the flavin ribityl, similar to P293, and three crystal waters in para-hydroxybenzoate hydroxylase suggest that P302, T47, and the waters in StyA are a vital component of the catalytic mechanism. The current optimized and validated StyA model provides a good starting point for elucidation of the structural basis of StyA ligand binding and catalysis. Novel insights in the binding of ligands to SMO/StyA, provided by the current protein model, will aid the rational design of mutants with specific, altered enantioselective properties.

Binding Sites↗

PubMed2005

Molecular modeling and docking simulations of scorpion toxins and related analogs on human SKCa2 and SKCa3 channels.

The small-conductance Ca2+-activated K+ (SKCa) channels modulate cytosolic Ca2+ concentration in excitable and non-excitable tissues by regulating the membrane potential and are responsible of slow action potential after hyperpolarization that inhibits cell firing. Among these, human SKCa2 and SKCa3 channels differ in the pore region by only two residues: Ala331 and Asn367 (human small-conductance calcium-activated potassium channel, hSKCa2) instead of Val485 and His521 (hSKCa3). To design highly selective blockers of hSKCa channels, a number of known hSKCa2 and/or hSKCa3-active peptides (i.e. scorpion toxins and analogs thereof) were analyzed for their interactions and selectivities toward these channels. Molecular models of hSKCa2 and hSKCa3 channels (S5-H5-S6 portion) were generated, and scorpion toxins/peptides of unsolved three-dimensional (3D) structures were modeled. Models of toxin-channel complexes were generated by the bimolecular complex generation with global evaluation, and ranking (BiGGER) docking software and selected by using a screening method of the docking solutions. A high degree of correlation was found to exist between docking energies and experimental Kd values of peptides that blocked hSKCa2 and/or hSKCa3 channels, suggesting it could be appropriate to predict Kd values of other bioactive peptides. The best scoring complexes were also used to identify key residues of both interacting partners, indicating that such an approach should help the design of more active and/or selective peptide blockers of targeted ion channels.

Amino Acid Motifs↗

PubMedDate not supplied

A combination of molecular dynamics and docking calculations to explore the binding mode of ADS-J1, a polyanionic compound endowed with anti-HIV-1 activity.

The HIV-1 entry process is an important target for the design of new pharmaceuticals for the multidrug therapy of AIDS. A lot of polyanionic compounds, such as polysulfonated and polysulfated, are reported in the literature for their ability to block early stages of HIV-1 replication. Several studies have been performed to elucidate the mechanism of the anti-HIV-1 activity of sulfated polysaccharides and polyanions in general, including binding to cell surface CD4 and interfering with the gp120-coreceptor interaction. Here, we show molecular modeling investigations on ADS-J1, a polyanionic compound with anti-HIV activity that is able to interfere with gp120-coreceptor interactions. Agreeing with experimental data, computer simulations suggested that the V3 loop of gp120 was the preferential binding site for ADS-J1 onto HIV-1. Moreover, mutations induced by the inhibitor significantly changed the stereoelectronic properties of the gp120 surface, justifying a marked drop in the affinity of ADS-J1 toward an ADS-J1-resistant HIV-1 strain.

Anti-HIV Agents↗

PubMed2006

LysM domains of Medicago truncatula NFP protein involved in Nod factor perception. Glycosylation state, molecular modeling and docking of chitooligosaccharides and Nod factors.

The establishment of the symbiosis between legume plants and rhizobial bacteria depends on the production of rhizobial lipo-chitooligosaccharidic signals (the Nod factors) that are specifically recognized by roots of the host plant. In Medicago truncatula, specific recognition of Sinorhizobium meliloti and its Nod factors requires the NFP (Nod factor perception) gene, which encodes a putative serine/threonine receptor-like kinase (RLK). The extracellular region of this protein contains three tandem lysin motifs (LysMs), a short peptide domain that is implicated in peptidoglycan or chitin binding in various bacterial or eukaryotic proteins, respectively. We report here the homology modeling of the three LysM domains of M. truncatula NFP based on the structure of a LysM domain of the Escherichia coli membrane-bound lytic murein transglycosidase D (MltD). Expression of NFP in a homologous system (M. truncatula roots) revealed that the protein is highly N-glycosylated, probably with both high-mannose and complex glycans. Surface analysis and docking calculations performed on the models of the three domains were used to predict the most favored binding modes for chitooligosaccharides and Nod factors. A convergent model can be proposed where the sulfated, O-acetylated lipo-chitooligosaccharidic Nod factor of S. meliloti binds in similar orientation to the three LysM domains of M. truncatula NFP. N-glycosylation is not expected to interfere with Nod factor binding in this orientation.

Amino Acid Motifs↗

PubMed2003

Construction of molecular assemblies via docking: modeling of tetramers with D2 symmetry.

Comparative modeling methods are commonly used to construct models of homologous proteins or oligomers. However, comparative modeling may be inapplicable when the number of subunits in a modeled oligomer is different than in the modeling template. Thus, a dimer cannot be a template for a tetramer because a new monomer-monomer interface must be predicted. We present in this study a new prediction approach, which combines protein-protein docking with either of two tetramer-forming algorithms designed to predict the structures of tetramers with D2 symmetry. Both algorithms impose symmetry constraints. However, one of them requires identification of two of the C2 dimers within the tetramer in the docking step, whereas the other, less demanding algorithm, requires identification of only one such dimer. Starting from the structure of one subunit, the procedures successfully reconstructed 16 known D2 tetramers, which crystallize with either a monomer, a dimer or a tetramer in the asymmetric unit. In some cases we obtained clusters of native-like tetramers that differ in the relative rotation of the two identical dimers within the tetramer. The predicted structural pliability for concanavalin-A, phosphofructokinase, and fructose-1,6-bisphosphatase agrees semiquantitatively with the observed differences between the several experimental structures of these tetramers. Hence, our procedure identifies a structural soft-mode that allows regulation via relative rigid-body movements of the dimers within these tetramers. The algorithm also predicted three nearly correct tetramers from model structures of single subunits, which were constructed by comparative modeling from subunits of homologous tetrameric, dimeric, or hexameric systems.

Algorithms↗

PubMed1994

Hydrophobic docking: a proposed enhancement to molecular recognition techniques.

In the classical procedures for predicting the structure of protein complexes two molecules are brought in contact at multiple relative positions, the extent of complementarity (geometric and/or energy) at the surface of contact is assessed at each position, and the best fits are retrieved. In view of the higher occurrence of hydrophobic groups at contact sites, their contribution results in more intermolecular atom-atom contacts per unit area for correct matches than for false positive fits. The hydrophobic groups are also potentially less flexible at the surface. Thus, from a practical point of view, a partial representation of the molecules based on hydrophobic groups should improve the quality of the results in finding molecular recognition sites, as compared to full representation. We tested this proposal by applying the idea to an existing geometric fit procedure and compared the results obtained with full vs. hydrophobic representations of molecules in known molecular complexes. The hydrophobic docking yielded distinctly higher signal-to-noise ratio so that the correct match is discriminated better from false positive fits. It appears that nonhydrophobic groups contribute more to false matches. The results are discussed in terms of their relevance to molecular recognition techniques as compared to energy calculations.

Algorithms↗

PubMed2005

Molecular flexibility in ab initio drug docking to DNA: binding-site and binding-mode transitions in all-atom Monte Carlo simulations.

The dynamics of biological processes depend on the structure and flexibility of the interacting molecules. In particular, the conformational diversity of DNA allows for large deformations upon binding. Drug-DNA interactions are of high pharmaceutical interest since the mode of action of anticancer, antiviral, antibacterial and other drugs is directly associated with their binding to DNA. A reliable prediction of drug-DNA binding at the atomic level by molecular docking methods provides the basis for the design of new drug compounds. Here, we propose a novel Monte Carlo (MC) algorithm for drug-DNA docking that accounts for the molecular flexibility of both constituents and samples the docking geometry without any prior binding-site selection. The binding of the antimalarial drug methylene blue at the DNA minor groove with a preference of binding to AT-rich over GC-rich base sequences is obtained in MC simulations in accordance with experimental data. In addition, the transition between two drug-DNA-binding modes, intercalation and minor-groove binding, has been achieved in dependence on the DNA base sequence. The reliable ab initio prediction of drug-DNA binding achieved by our new MC docking algorithm is an important step towards a realistic description of the structure and dynamics of molecular recognition in biological systems.

Algorithms↗

PubMed2025

Genomics-informed drug-repurposing strategy identifies two therapeutic targets for preventing liver disease associated with metabolic dysfunction.

Identification of drug-repurposing targets with genetic and biological support is an economically and temporally efficient strategy for improving the treatment of diseases. We employed a cross-disciplinary approach to identify potential therapeutics for the prevention of metabolic-dysfunction-associated steatotic liver disease (MASLD) in at-risk individuals by using humans as a model organism. We identified 212 putative candidate genes associated with MASLD by using data from a large multi-ancestry genetic association study, of which 158 (74.5%) were previously unreported. From this set, we identified 57 genes that encode for druggable protein targets and for which the effects of increasing genetically predicted gene expression on MASLD risk align with the function of that drug on the protein target. We then used We then evaluated these potential targets for evidence of efficacy by using Mendelian randomization, pathway analysis, and protein structural modeling. Through these approaches, we present compelling evidence to suggest that the activation of FADS1 by icosapent ethyl, as well as S1PR2 by fingolimod, could be a promising therapeutic strategy for MASLD prevention.

Humans↗

PubMed1996

Molecular similarity based on DOCK-generated fingerprints.

An alternative method for defining molecular similarity is presented. By using the docking program DOCK and a reference panel of protein binding sites, fingerprints for a set of molecules have been generated, based on calculated interaction energies. These binding patterns allowed us to calculate matrices of similarity coefficients which subsequently were used for nearest-neighbor searches within the database. Our results indicate that the method is suitable for finding significant similarities of compounds of the same biological activity. Although the overall performance of a traditional 2D similarity method is better in the test systems investigated, our 3D approach can be regarded as complementary since it is able to detect similarities independent of the covalent structure of the compounds. Thus it should be a useful 3D database-searching tool for rational lead discovery.

Amino Acid Sequence↗

PubMed1993

Molecular recognition analyzed by docking simulations: the aspartate receptor and isocitrate dehydrogenase from Escherichia coli.

Protein docking protocols are used for the prediction of both small molecule binding to DNA and protein macromolecules and of complexes between macromolecules. These protocols are becoming increasingly automated and powerful tools for computer-aided drug design. We review the basic methodologies and strategies used for analyzing molecular recognition by computer docking algorithms and discuss recent experiments in which (i) substrate and substrate analogues are docked to the active site of isocitrate dehydrogenase and (ii) maltose binding protein is docked to the extracellular domain of the receptor, which signals maltose chemotaxis.

Amino Acid Sequence↗