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 613 records · Page 34Linked to original sources

PubMed2002

Structural investigation of the HIV-1 envelope glycoprotein gp160 cleavage site.

The selective proteolytic activation of the HIV-1 envelope glycoprotein gp160 by furin and other precursor convertases (PCs) occurs at the carboxyl side of the sequence Arg508-Glu-Lys-Arg511 (site 1), in spite of the presence of another consensus sequence: Lys500-Ala-Lys-Arg503 (site 2). We report on the solution structural analysis of a 19-residue synthetic peptide, p498, which spans the two gp160-processing sites 1 and 2, and is properly digested by furin at site 1. A molecular model is obtained for p498, by means of molecular dynamics simulations, from NMR data collected in trifluoroethanol/water. The peptide N-terminal side presents a 9-residue helical segment, enclosing the processing site 2; the C-terminal segment can be described as a loop exposing the processing site 1. A hypothesis for the docking of p498 onto the catalytic domain of human furin, modeled by homology and fitting previous site-directed mutagenesis studies, is also presented. p498 site 1 is shown to have easy access to the furin catalytic site, unlike the nonphysiological site 2. Finally, on the basis of available data, we suggest a possible structural motif required for the gp160-PCs recognition.

Amino Acid Sequence↗

PubMed2006

Insight into the inhibition of human choline kinase: homology modeling and molecular dynamics simulations.

A homology model of human choline kinase (CK-alpha) based on the X-ray crystallographic structure of C. elegans choline kinase (CKA-2) is presented. Molecular dynamics simulations performed on CK-alpha confirm the quality of the model, and also support the putative ATP and choline binding sites. A good correlation between the MD results and reported CKA-2 mutagenesis assays has been found for the main residues involved in catalytic activity. Preliminary docking studies performed on the CK-alpha model indicate that inhibitors can bind to the binding sites of both substrates (ATP and choline). A possible reason for inhibition of choline kinase by Ca(2+) ion is also proposed.

Amino Acid Sequence↗

PubMed1999

Thermodynamics and kinetics of ligand-protein binding studied with the weighted histogram analysis method and simulated annealing.

The thermodynamics of ligand-protein molecular recognition is investigated by the energy landscape approach for two systems: methotrexate(MTX)--dihydrofolate reductase(DHFR) and biotin-streptavidin. The temperature-dependent binding free energy profile is determined using the weighted histogram analysis method. Two different force fields are employed in this study: a simplified model of ligand-protein interactions and the AMBER force field with a soft core smoothing component, used to soften the repulsive part of the potential. The results of multiple docking simulations are rationalized from the shape of the binding free energy profile that characterizes the thermodynamics of the binding process.

Biotin↗

PubMed1999

Role of the extracellular loops of G protein-coupled receptors in ligand recognition: a molecular modeling study of the human P2Y1 receptor.

The P2Y1 receptor is a G protein-coupled receptor (GPCR) and is stimulated by extracellular ADP and ATP. Site-directed mutagenesis of the three extracellular loops (ELs) of the human P2Y1 receptor indicates the existence of two essential disulfide bridges (Cys124 in EL1 and Cys202 in EL2; Cys42 in the N-terminal segment and Cys296 in EL3) and several specific ionic and H-bonding interactions (involving Glu209 and Arg287). Through molecular modeling and molecular dynamics simulations, an energetically sound conformational hypothesis for the receptor has been calculated that includes transmembrane (TM) domains (using the electron density map of rhodopsin as a template), extracellular loops, and a truncated N-terminal region. ATP may be docked in the receptor, both within the previously defined TM cleft and within two other regions of the receptor, termed meta-binding sites, defined by the extracellular loops. The first meta-binding site is located outside of the TM bundle, between EL2 and EL3, and the second higher energy site is positioned immediately underneath EL2. Binding at both the principal TM binding site and the lower energy meta-binding sites potentially affects the observed ligand potency. In meta-binding site I, the side chain of Glu209 (EL2) is within hydrogen-bonding distance (2.8 A) of the ribose O3', and Arg287 (EL3) coordinates both alpha- and beta-phosphates of the triphosphate chain, consistent with the insensitivity in potency of the 5'-monophosphate agonist, HT-AMP, to mutation of Arg287 to Lys. Moreover, the selective reduction in potency of 3'NH2-ATP in activating the E209R mutant receptor is consistent with the hypothesis of direct contact between EL2 and nucleotide ligands. Our findings support ATP binding to at least two distinct domains of the P2Y1 receptor, both outside and within the TM core. The two disulfide bridges present in the human P2Y1 receptor play a major role in the structure and stability of the receptor, to constrain the loops within the receptor, specifically stretching the EL2 over the opening of the TM cleft and thus defining the path of access to the binding site.

Amino Acid Sequence↗

PubMed2003

Computational docking of L-arginine and its structural analogues to C-terminal domain of Escherichia coli arginine repressor protein (ArgRc).

The arginine repressor (ArgR) of Escherichia coli binds to six L-arginine molecules that act as its co-repressor in order to bind to DNA. The binding of L-arginine molecules as well as its structural analogues is compared by means of computational docking. A grid-based energy evaluation method combined with a Monte Carlo simulated annealing process was used in the automated docking. For all ligands, the docking procedure proposed more than one binding site in the C-terminal domain of ArgR (ArgRc). Interaction patterns of ArgRc with L-arginine were also observed for L-canavanine and L-citrulline. L-lysine and L-homoarginine, on the other hand, were shown to bind poorly at the binding site. Figure A general overview of the sites found from docking the various ligands into ArgRc ( grey ribbons). Red coloured sticks: residues in binding site H that was selected for docking

Arginine↗

PubMed2001

Computer modeling of human angiogenin-dinucleotide substrate interaction.

Structures of substrate bound human angiogenin complexes have been obtained for the first time by computer modeling. The dinucleotides CpA and UpA have been docked onto human angiogenin using a systematic grid search procedure in torsion and Eulerian angle space. The docking was guided throughout by the similarity of angiogenin-substrate interactions with interactions of RNase A and its substrate. The models were subjected to 1 nanosecond of molecular dynamics to access their stability. Structures extracted from MD simulations were refined by simulated annealing. Stable hydrogen bonds that bridged protein and ligand residues during the MD simulations were taken as restraints for simulated annealing. Our analysis on the MD structures and annealed models explains the substrate specificity of human angiogenin and is in agreement with experimental results. This study also predicts the B2 binding site residues of angiogenin, for which no experimental information is available so far. In the case of one of the substrates, CpA, we have also identified the presence of a water molecule that invariantly bridges the B2 base with the protein. We have compared our results to the RNase A-substrate complex and highlight the similarities and differences.

Amino Acid Sequence↗

PubMed2001

Molecular dynamics simulations of the mononuclear zinc-beta-lactamase from Bacillus cereus complexed with benzylpenicillin and a quantum chemical study of the reaction mechanism.

Herein, we present results from MD simulations of the Michaelis complex formed between the B. cereus zinc-beta-lactamase enzyme and benzylpenicillin. The structural and dynamical effects induced by substrate-binding, the specific role of the conserved residues, and the near attack conformers of the Michaelis complex are discussed. Quantum chemical methods (HF/6-31G* and B3LYP/6-31G*) are also applied to study the hydrolysis reaction of N-methylazetidinone catalyzed by a monozinc system consisting of the side chains of the histidine residues (His86, His88, and His149) complexed with Zn-OH and the side chains of Asp90 and His210. From this model system, we built molecular-mechanics representations of the prereactive complex and transition state configurations docked into the active site. Linear-scaling semiempirical calculations coupled with a continuum solvent model were then performed on these static models. We propose that the experimental rate data for the B. cereus enzyme is compatible with a one-step mechanism for the hydrolysis of beta-lactam substrates in which His210 acts as a proton donor.

Bacillus cereus↗

PubMed2000

DARWIN: a program for docking flexible molecules.

A new program named "DARWIN" has been developed to perform docking calculations with proteins and other biological molecules. The program uses the Genetic Algorithm to optimize the molecule's conformation and orientation under the selective pressure of minimizing the potential energy of the complex. A unique feature of DARWIN is that it communicates with the molecular mechanics program CHARMM to make the energy calculations. A second important feature is its parallel interface, which allows simultaneous use of multiple stand-alone copies of CHARMM to rapidly evaluate large numbers of potential solutions. This permits an "accuracy first" approach to docking, which avoids many of the common assumptions and shortcuts often made to reduce computation time. The method was applied to three protein-carbohydrate complexes: the crystallographically determined structures of Concanavalin A and Fab Se155-4; and a model structure for Fab ME36.1. Conformations close to the crystal structures were obtained with this approach, but some "false positive" solutions were also selected. Many of these could be eliminated by introducing different methods for simulating solvent effects. An effective screening method for docking a database of compounds to a single target enzyme using DARWIN is also presented.

Algorithms↗

PubMedDate not supplied

Enrichment of ligands for the serotonin receptor using the Shape Signatures approach.

Shape Signatures, a new 3-dimensional molecular comparison method, has been adapted to rank ligands of the serotonin receptors. A set of 825 agonists and 400 antagonists together with approximately 10,000 randomly chosen compounds from the NCI database were used in this study. Both 1D and 2D Shape Signature databases were created, and enrichment studies were carried out. Results from these studies reveal that the 1D Shape Signature approach is highly efficient in separating agonists from a mixture of molecules which includes compounds randomly selected from the NCI database taken as inactives. It is also equally effective at separating agonists and antagonists from a pool of active ligands for the serotonin receptor. Parallel enrichment studies using 2D shape signatures showed high selectivity with more restricted coverage due to the high specificity of 2D signatures. The influence of conformational variation of the shape signature on enrichment was explored by docking a subset of ligands into the crystal structure of serotonin N-acetyltransferase. Enrichment studies on the resulting "docked" conformations produced only slightly improved results compared with the CORINA-generated conformations.

Computer Simulation↗

PubMed1997

Development of a structural model for the cytoplasmic domain of an integrin.

The cytoplasmic tails of integrin heterodimers play central roles in controlling the activation states of integrins and in transmitting intracellular signals. Despite their short length, no structure of any integrin cytoplasmic domain has been determined. Therefore, molecular models for the cytoplasmic domain of alpha(IIb)beta3, the major platelet integrin, were generated, including models for the individual cytoplasmic tails, the binary alphaIIb-calcium complex, and the ternary alphaIIb-beta3-calcium complex. Structural analysis of circular dichroism spectra were compiled with data obtained from short homologous sequences within crystallized proteins, and with secondary structural predictions to develop starting models for each subunit. These models were subjected to a series of energy minimization and molecular dynamic simulations to generate final models. AlphaIIb was predicted to be ordered at its N-terminus and its C-terminus could accommodate a cation in a multicoordinated complex. The structure of beta3 was dominated by a beta-turn at its NPXY motif (beta3 744-747). In docking of alphaIIb to different sites within beta3, the conformation of the beta3 juxta-transmembrane (beta3 716-721) was greatly altered. This region was confirmed to be a conformational 'hot-spot' by circular dichroism. The conformational flexibility of this juxta-transmembrane region, which is highly conserved amongst integrins, is ideally located to regulate signaling.

Amino Acid Sequence↗

PubMed2003

Molecular docking of substrates and inhibitors in the catalytic site of CYP6B1, an insect cytochrome p450 monooxygenase.

Furanocoumarins represent plant toxins that are used in the treatment of a variety of skin diseases and are metabolized by cytochrome p450 monooxygenases (p450s) existing in insects such as Papilio polyxenes (the black swallowtail). To elucidate the active site in the CYP6B1 protein that is the principal p450 existing in this species, we have constructed a homology model of it based on sequence and structure alignments with the bacterial CYP102 protein whose crystal structure has been defined and with the insect CYP6B4 protein that also metabolizes furanocoumarins. In the derived CYP6B1 model, Phe116 and His117 in SRS1, Phe371 in SRS5 and Phe484 in SRS6 contribute to the formation of a resonant network that stabilizes the p450's catalytic site and allows for interactions with its furanocoumarin substrates. The first two of these residues are absolutely conserved in all members of the insect CYP6B subfamily and the last two are variable in different members of the CYP6B subfamily. A combination of theoretical and experimental docking analyses of two substrates (xanthotoxin and bergapten) and two inhibitors (coumarin and pilocarpine) of this p450 provide significant information on the positioning of furanocoumarins within this catalytic pocket. Molecular replacement models based on the results of variations at two of these critical amino acids provide support for our furanocoumarin-docked model and begin to rationalize the altered substrate reactivities observed in experimental analyses.

5-Methoxypsoralen↗

PubMed2004

Computational analysis of tyrosine phosphatase inhibitor selectivity for the virulence factors YopH and SptP.

Bacterial pathogens such as Yersinia and Salmonella represent an important medical concern, causing human diseases ranging from gastrointestinal disease to the plague. The development of novel treatments of these bacterial infections has gained high priority recently due to the emergence of antibiotic resistance in these pathogens and the threat of the use of microbial agents as biological weapons. YopH of Yersinia and SptP of Salmonella are virulence factors that belong to the family of protein tyrosine phosphatases (PTPs). A great challenge remains in the design of selective PTPs inhibitors due to their highly conserved active site. In this paper, we present a comparative docking study to probe the selective inhibition of YopH and SptP with PTP1B in order to better understand their binding interactions with the bacterial tyrosine phosphates. Characterized binding sites in PTP1B were compared with YopH and SptP. Molecular dynamics simulations were used to incorporate ligand-induced conformational changes in the binding sites. These results, together with those binding modes and binding affinities distinguished in individual PTPs, provide insight into the structure-based design of inhibitors for YopH and SptP.

Bacterial Outer Membrane Proteins↗

PubMed1997

Comparative molecular modeling study of the three-dimensional structures of prostaglandin endoperoxide H2 synthase 1 and 2 (COX-1 and COX-2).

To understand the structural features that dictate the selectivity of diverse nonsteroidal antiinflammatory drugs for the two isoforms of the human prostaglandin H2 synthase (PGHS), the three-dimensional (3D) structure of human COX-2 was assessed by means of sequence homology modeling. The ovine COX-1 structure, solved by X-ray diffraction methods and sharing a 61% sequence identity with human COX-2, was used as template. Both structures were energy minimized using the AMBER 4.0 force field with a dielectric constant of 4r. (S)-Flurbiprofen, a nonselective COX inhibitor, and SC-558, a COX-2-selective ligand, were docked at the cyclooxygenase binding site in both isozymes, evidencing the role of different residues in the ligand-protein interaction. The 3D structures of the constructed four ligand-enzyme complexes were refined by energy minimization. Molecular dynamics simulations were also carried out, to understand more deeply the structural origins of the selectivity. Distances calculated during the dynamics process between the different ligands and the interacting residues of the two PGHS isozymes provided evidence of the flexible nature of the cyclooxygenase active site, permitting the identification of different conserved and nonconserved residues as responsible for ligand selectivity.

Amino Acid Sequence↗

PubMed2001

Investigations on inhibitors of human 17 alpha-hydroxylase-17,20-lyase and their interactions with the enzyme. Molecular modelling of 17 alpha-hydroxylase-17,20-lyase, Part II.

New methods in treatment of hormone-dependent diseases like prostate or breast cancer have become a major subject in medical and pharmaceutical research. Because of the direct correlation of cancer growth and hormone concentration, inhibition of hormone biosynthesis presents a promising strategy in cancer therapy. The key enzyme in androgen biosynthesis is the 17 alpha-hydroxylase-17,20-lyase a cytochrome P450 system, which specifically converts gestagens to androgens. Because the 3D-structure of the enzyme is still unknown most recently a ligand-based design was used to gain deeper insights into protein structure and function. In this paper we present molecular modelling studies on compounds acting as competitive inhibitors of the human 17 alpha-hydroxylase-17,20-lyase. The compounds developed by Hartmann et al. belong to two different structural classes and show a wide range of inhibitory potency. The physico-chemical properties of the molecules were investigated and compared by studying structural flexibility and by calculating molecular interactions fields. The superimposition of all inhibitors in a low energy conformation yielded in the common pharmacophore. In the second part of the paper individual inhibitors were docked into the active site of the enzyme model of CYP17 developed in our group. The dynamic behaviour and stability of the protein-inhibitor-complexes was studied. The protein ligand interactions observed in course of the molecular dynamics simulations correspond well with the experimental data.

Computer Simulation↗

PubMed1998

Computer simulations of carbon monoxide photodissociation in myoglobin: structural interpretation of the B states.

The early diffusion processes of a photodissociated ligand (carbon monoxide) in sperm whale myoglobin and its Phe29 mutant are studied computationally. An explicit solvent model is employed in which the protein is embedded in a box of at least 2300 water molecules. Electrostatic interactions are accounted for by using the particle mesh Ewald. Two hundred seventy molecular dynamics trajectories are computed for 10 ps. Different models of solvation and the ligand, and their influence on the diffusion are examined. The two B states of the CO are identified as "docking" sites in the heme pocket. The sites have a similar angle with respect to the heme normal, but differ in the orientation in the plane. The computational detection of the B states is stable under a reasonable variation of simulation conditions. However, in some trajectories only one of the states is observed. It is therefore necessary to use extensive simulation data to probe these states. Comparison to diffraction experiments and spectroscopy is performed. The shape of the experimental infrared spectra is computed. The overall linewidth is in an agreement with experiment. The contributions to the linewidth (van der Waals and electrostatic interactions) are discussed.

Animals↗

PubMed2006

Virtual screening for finding natural inhibitor against cathepsin-L for SARS therapy.

Recently Simmons et al. reported a new mechanism for SARS virus entry into target cells, where MDL28170 was identified as an efficient inhibitor of CTSL-meditated substrate cleavage with IC(50) of 2.5 nmol/l. Based on the molecule fingerprint searching method, 11 natural molecules were found in the Traditional Chinese Medicines Database (TCMD). Molecular simulation indicates that the MOL376 (a compound derived from a Chinese medicine herb with the therapeutic efficacy on the human body such as relieving cough, removing the phlegm, and relieving asthma) has not only the highest binding energy with the receptor but also the good match in geometric conformation. It was observed through docking studies that the van der Waals interactions made substantial contributions to the affinity, and that the receptor active pocket was too large for MDL21870 but more suitable for MOL736. Accordingly, MOL736 might possibly become a promising lead compound for CTSL inhibition for SARS therapy.

Algorithms↗

PubMed2005

Hierarchical docking of databases of multiple ligand conformations.

Ligand flexibility is an important problem in molecular docking and virtual screening. To address this challenge, we investigate a hierarchical pre-organization of multiple conformations of small molecules. Such organization of pre-calculated conformations removes the exploration of ligand conformational space from the docking calculation and allows for concise representation of what can be thousands of conformations. The hierarchy also recognizes and prunes incompatible conformations early in the calculation, eliminating redundant calculations of fit. We investigate the method by docking the MDL Drug Data Report (MDDR), an annotated database of 100,000 molecules, into apo and holo forms of seven unrelated targets. This annotated database allows us to track the ranking of tens to hundreds of annotated ligands in each of the docking systems. The binding sites and database are prepared in an automated fashion in an attempt to remove some human bias from the calculations. Many thousands of explicit and implicit ligand conformations may be docked in calculations not much longer than required for single conformer docking. As long as internal energies are not considered, recombination with the hierarchy is additive as the number of degrees of freedom is increased. Molecules with even millions of conformations can be docked in a few minutes on a single desktop computer.

Computer Simulation↗

PubMed2005

[In-silico prediction of pharmacokinetic properties].

In silico methods for predicting pharmacokinetic properties range from data-based approaches such as quantitative structure-activity relationships (QSARs), similarity searches, and 3-dimensional QSAR, to structure-based methods such as ligand-protein docking and pharmacophore modelling. Data-based modelling approaches are effective for many drug absorption, distribution, metabolism, and excretion (ADME) processes such as passive membrane permeation, where their molecular mechanism is barely delineated. Therefore QSAR approaches have been applied to simulate the relationships between ADME parameters and molecular structure and properties. In the present investigation, we describe the application of the genetic algorithm-combined partial least-squares (GA-PLS) method to QSAR modelling of various ADME properties. By selecting an appropriate set of molecular descriptors automatically using the genetic algorithm, many ADME properties could be well explained by simple molecular descriptors derived from the 2-dimensional chemical structure.

Algorithms↗