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 739 records · Page 41Linked to original sources

PubMed2004

A computational study of the binding of propidium to the peripheral anionic site of human acetylcholinesterase.

Combined docking and molecular dynamics (MD) simulations were carried out in order to investigate the binding mode of propidium at the human acetylcholinesterase (HuAChE) peripheral site. Two different docking protocols followed by cluster analyses were performed, allowing the identification of five high-populated and low-energy configuration families. To dynamically explore the behavior of the ligand at the peripheral HuAChE binding site, six complexes (five low-energy and one high-energy) were submitted to 2.5 ns of MD simulations. The representative propidium/HuAChE binding modes were chosen on the basis of both the docking energy score and the dynamic stability of the complexes throughout the MD simulations. The most stable poses of propidium at HuAChE PAS were similar to those experimentally determined with the murine enzyme. We therefore suggest that the present modeling protocol might be used in the dynamic investigation of the interactions of a small-molecule inhibitor with a surface-like binding site of a protein. Finally, because of the biological relevance of the target studied here, the present results can be of interest for the rational design of molecules potentially useful in the treatment of the Alzheimer's disease.

Acetylcholinesterase↗

PubMed2003

Molecular recognition between 4a S/R-galanthamine diastereoisomers and alpha-cyclodextrin.

Molecular recognition between 4a S/R-galanthamine diastereoisomers (1: 4a S-galanthamine;2: 4a R-galanthamine) and alpha-cyclodextrin (alpha-CD) were studied by use of docking and molecular dynamics (MD) simulation approaches. The binding energy of constructed 2...alpha-CD complexes is approximately 17 kcal mol(-1) lower than that of 1...alpha-CD, implying a stronger binding ability of 2 with alpha-CD than that of1. The theoretical modeling result is consistent with our previous CZE result, which demonstrated that alpha-CD is an efficient chiral additive for separating 1 and 2. The modeling result also indicates that both hydrophobic interaction and H-bond force may work as major factors for molecular recognition between the galanthamine diastereoisomers and alpha-CD.

Computer Simulation↗

PubMed2003

RNA aptamer evaluation using comparative molecular interaction profile analysis system.

In this paper, we introduce a system named Comparative Molecular Interaction Profile Analysis system (CoMIPA) that utilizes virtual docking study for comprehensive analysis of molecular interactions. Research of RNA aptamer is very useful to creation of functional molecules, because RNA aptamer is excellent in binding selectivity and the stability of structure. By using 3-bases RNA fragment, the system can use computational molecular docking results to help for RNA aptamer screening. We believe that the system has the potential to be a useful for RNA aptamer design.

Binding Sites↗

PubMed2000

Interaction of pyrimethamine, cycloguanil, WR99210 and their analogues with Plasmodium falciparum dihydrofolate reductase: structural basis of antifolate resistance.

The nature of the interactions between Plasmodium falciparum dihydrofolate reductase (pfDHFR) and antimalarial antifolates, i.e., pyrimethamine (Pyr), cycloguanil (Cyc) and WR99210 including some of their analogues, was investigated by molecular modeling in conjunction with the determination of the inhibition constants (Ki). A three-dimensional structural model of pfDHFR was constructed using multiple sequence alignment and homology modeling procedures, followed by extensive molecular dynamics calculations. Mutations at amino acid residues 16 and 108 known to be associated with antifolate resistance were introduced into the structure, and the interactions of the inhibitors with the enzymes were assessed by docking and molecular dynamics for both wild-type and mutant DHFRs. The Ki values of a number of analogues tested support the validity of the model. A 'steric constraint' hypothesis is proposed to explain the structural basis of the antifolate resistance.

Amino Acid Sequence↗

PubMed2005

Regiospecificity of human cytochrome P450 1A1-mediated oxidations: the role of steric effects.

Cytochrome P450 1A1 oxidizes a diverse range of substrates, including the procarcinogenic xenobiotic benzo[a]pyrene (B[a]P) and endogenous fatty acid precursors of prostaglandins, such as arachidonic acid (AA) and eicosapentaenoic acid (EA). We have investigated the extent to which enzyme-substrate interactions govern regio- and stereoselectivity of oxidation of these compounds by using docking and molecular dynamics (MD) simulations to examine the likelihood of substrate oxidation at various sites. Due to structural differences between the substrates analyzed, B[a]P and its diols (planar, rigid), and the fatty acids AA and EA (long, flexible), different docking strategies were required. B[a]P, B[a]P-7,8-diols, (+) 7S,8S- and (-) 7R,8R-diols, were docked into the active site of a homology model of P450 1A1 using an automated routine, Affinity (Accelrys, San Diego, CA). AA and EA, on the other hand, required a series of restrained MD simulations to obtain a variety of productive binding modes. All complexes were evaluated by MD-based in silico site scoring to predict product profiles based on certain geometric criteria, such as angle and distance of a given substrate atom from the ferryl oxygen. For all substrates studied, the in vitro profiles were generally reflected by the in silico scores, which suggests that steric factors play a key role in determining regiospecificity in P450 1A1-mediated oxidations. We have also shown that molecular dynamics simulations may be very useful in determination of product profiles for structurally diverse substrates of P450 enzymes.

Binding Sites↗

PubMed2003

Nonnucleoside inhibitors of HIV-1 reverse transcriptase: from the biology of reverse transcription to molecular design.

Replication of human immunodeficiency virus 1 (HIV-1) uses a viral reverse transcriptase (RT) to convert its positive-strand RNA into double stranded DNA, which is then integrated into host genome. Reverse transcription is a complex event involving p66 and p51 RT subunits but also several viral proteins including Nef, Tat, Vif, IN, NCp7 and p55gag. Viral RNA itself forms a primer/template complex by association with a cellular tRNA(Lys3) which is already present in mature virions. A RT initiation complex (RTIC) is thus formed which may also involve cellular protein upon viral entry. X rays diffraction and NMR studies of free or inhibitor-bound RT have led to the recognition of RT 3D structure, and allowed a thorough understanding of the mode of action of classical competitive nucleoside RT inhibitors (NRTIs) and of the binding of allosteric, non NRTIs (NNRTIs) inhibitors. This also opened an access to computer-aided drug design and modeling. Current NNRTIs represent, in terms of chemical structures, a heterogeneous class of inhibitors currently undergoing extensive development. By contrast with NRTIs, they seem to block initiation steps of reverse transcription. Molecular dynamics, detailed analysis of their interaction with RT as well as the incidence, in the series, of cases of non classical biological behavior, as illustrated here for a new family of compounds, suggest mechanisms of action which are not understandable without considering the involvement of the RTIC as a whole. This opens the exciting perspective of developing new compounds based on this integrated knowledge. Key Words: Nonnucleoside reverse transcriptase inhibitors (NNRTIs); Reverse transcriptase initiation complex (RTIC); Human immunodeficiency virus (HIV); Non classical nonnucleoside reverse transcriptase inhibitors; Molecular modeling; Docking; QSAR; Natural endogenous reverse transcription (NERT).

Drug Delivery Systems↗

PubMed2006

Reverse-docking study of the TADDOL-catalyzed asymmetric hetero-Diels-Alder reaction.

The reverse-docking of a TADDOL catalyst to rigid transition-state (TS) representations of an asymmetric hetero-Diels-Alder reaction is described. The resulting docking poses represent a simplified geometric model of the TS for the catalyzed reaction. The conformational space of the catalyst in proximity to the catalyst-free TS models is sampled stochastically and the energetically favored poses are subjected to a clustering procedure to highlight structural attributes compatible with organocatalysis. Each pose is scored and ranked based on its molecular-mechanics docking energy. The reverse-docking procedure reveals a clear energetic trend in favor of the experimentally preferred product enantiomers. Analysis of the best poses suggests a geometric model that is consistent with principles of molecular recognition, catalysis, and experimental data.

Catalysis↗

PubMed2006

ADMET properties, database screening, molecular dynamics, density functional, and docking studies of novel potential anti-cancer compounds.

Database screening was performed in a large database (hundreds of thousands of molecules which we optimized at the AM1 level) yielding a set of potential bioactive ligands. One new ligand was selected among the top solutions and optimized at the B3LYP/6-31G* level, yielding also NBO (Natural Bond Order) charges. A flexible docking program was used to investigate the interactions between the receptor and the new ligand. The stability as well as the main protein-ligand contacts of our proposed novel ligand as well as the crystallographic RAR ligand was investigated by molecular dynamics. The ADMET (Absorption, Distribution, Metabolism, Excretion and Toxicity) properties as well as the parameters of the Rule of Five were investigated. The result of this work is compared with a crystallographic ligand of RAR. Our novel proposed anti-cancer ligand indicates hydrophobic interactions and strong polar interactions with the receptor.

Animals↗

PubMed2026

Computational prediction of a multi-epitope Human Metapneumovirus vaccine candidate through integrated reverse vaccinology and pan-genomic approaches.

Human metapneumovirus (HMPV) is a primary cause of global respiratory infections yet no approved vaccine currently exists. This study computationally predicts a multi-epitope vaccine candidate using a diverse dataset of 65 HMPV sequences spanning five continents. Following the screening of lead proteins for antigenicity and virulence, fifteen highly conserved MHC-I, MHC-II and B-cell epitopes were prioritized. These were integrated with a putative L7/L12 adjuvant using optimized AAY, GPGPG, and KK linkers to design three constructs (HMPV_V1-V3). Structural validation identified HMPV-V2 as the lead candidate that exhibits a Z-score of-5.24 and 87.7% of residues in favored Ramachandran regions indicating excellent stereochemical quality and structural stability. In silico docking indicated a strong predicted binding affinity between HMPV-V2 and the TLR4 receptor (energy: -969.2). Immune simulations predicted a robust adaptive response characterized by high IgG1 titers, memory B-cell maturation, and a Th1-dominant cytokine profile. Furthermore, molecular dynamics simulations suggested exceptional structural integrity for HMPV-V2, maintaining a low RMSD of 8.213 and RMSF of 0.737 throughout the simulation. Optimized in silico cloning into the pET28a (+) vector indicated a high potential for protein expression in E. coli systems. While these findings provide a theoretically grounded blueprint for vaccine development, this study is entirely computational and lacks experimental validation. Further in vitro and in vivo testing is required to confirm the actual safety and immunogenicity of the proposed candidate.

Metapneumovirus↗

PubMed1995

Docking of congo red to the surface of crystalline cellulose using molecular mechanics.

A method for the docking of ligands to crystalline surfaces using molecular mechanics was developed and applied to study the molecular features of the adsorption of Congo Red to cellulose crystals. The docking procedure was based on a grid search exploring the surface repeat unit of the cellulose crystals using several orientations of Congo Red at each grid point. Results suggest a preferential adsorption of Congo Red onto certain faces of the cellulose crystals and that the lower energy conformers have a similar positioning and orientation with respect to the cellulose chains at the surface repeat unit.

Adsorption↗

PubMed2005

Development of molecular probes for the identification of extra interaction sites in the mid-gorge and peripheral sites of butyrylcholinesterase (BuChE). Rational design of novel, selective, and highly potent BuChE inhibitors.

Tacrine heterobivalent ligands were designed as novel and reversible inhibitors of cholinesterases. On the basis of the investigation of the active site gorge topology of butyrylcholinesterase (BuChE) and acetylcholinesterase (AChE) and by using flexible docking procedures, molecular modeling studies formulated the hypothesis of extra interaction sites in the active gorge of hBuChE, namely, a mid-gorge interaction site and a peripheral interaction site. The design strategy led to novel BuChE inhibitors, balancing potency and selectivity. Among the compounds identified, the heterobivalent ligand 4m, containing an amide nitrogen and a sulfur atom at the 8-membered tether level, is one of the most potent and selective BuChE inhibitors described to date. The novel inhibitors, bearing postulated key features, validated the hypothesis of the presence of extra interaction sites within the hBuChE active site gorge.

Acetylcholinesterase↗

PubMed2005

Dynamic mechanism for the autophosphorylation of CheA histidine kinase: molecular dynamics simulations.

The two-component system (TCS) is an important signal transduction component for most bacteria. This signaling pathway is mediated by histidine kinases via autophosphorylation between P1 and P4 domains. Taking chemotaxis protein CheA as a model of TCS, the autophosphorylation mechanism of the TCS histidine kinases has been investigated in this study by using a computational approach integrated homology modeling, ligand-protein docking, protein-protein docking, and molecular dynamics (MD) simulations. Four nanosecond-scale MD simulations were performed on the free P4 domain, P4-ATP, P4-TNPATP, and P1-P4-ATP complexes, respectively. Upon its binding to the binding pocket of P4 with a folded conformation, ATP gradually extends to an open state with help from a water molecule. Meanwhile, ATP forms two hydrogen bonds with His413 and Lys494 at this state. Because of the lower energy of the folded conformations, ATP shrinks back to its folded conformations, leading to the rupture of the hydrogen bond between ATP and Lys494. Consequently, Lys494 moves away from the pocket entrance, resulting in an open of the ATP lid of P4. It is the open state of P4 that can bind tightly to P1, where the His45 of P1 occupies a favorable position for its autophosphorylation from ATP. This indicates that ATP is not only a phosphoryl group donor but also an activator for CheA phosphorylation. Accordingly, a mechanism of the autophosphorylation of CheA is proposed as that the ATP conformational switch triggers the opening of the ATP lid of P4, leading to P1 binding tightly, and subsequently autophosphorylation from ATP to P1.

Adenosine Triphosphate↗

PubMed2004

Modelling and simulation: a computational perspective in anticancer drug discovery.

The availability of high-quality molecular graphics tools in the public domain is changing the way macromolecular structure is perceived by researchers, educators and students alike. Computational methods have become increasingly important in a number of areas such as comparative or homology modelling, functional site location, characterization of ligand-binding sites in proteins, docking of small molecules into protein binding sites, protein-protein docking, and molecular dynamics simulations. The results obtained yield information that sometimes is beyond current experimental possibilities and can be used to guide and improve a vast array of experiments. On the basis of our improved level of understanding of molecular recognition and the widespread availability of target structures, it is reasonable to assume that computational methods will continue aiding not only in the design and interpretation of hypothesis-driven experiments in the field of cancer research but also in the rapid generation of new hypotheses.

Amino Acid Sequence↗

PubMed2005

Structure of Thermus thermophilus type 2 isopentenyl diphosphate isomerase inferred from crystallography and molecular dynamics.

Crystal structures of Thermus thermophilus and Bacillus subtilis type 2 IPP isomerases were combined to generate an almost complete model of the FMN-bound structure of the enzyme. In contrast to previous studies, positions of flexible loops were obtained and carefully analyzed by molecular dynamics. Docking simulations find a unique putative binding site for the IPP substrate.

Carbon-Carbon Double Bond Isomerases↗

PubMed2005

Searching for protein-protein interaction sites and docking by the methods of molecular dynamics, grid scoring, and the pairwise interaction potential of amino acid residues.

In CAPRI Rounds 1 and 2, we assumed that because there are many ionic charges that weaken electrostatic interaction forces in living cells, the hydrophobic interaction force might be important entropically. As a result of Rounds 1 and 2, the predictions for binding sites and geometric centers were acceptable, but those of the binding axes were poor, because only the largest benzene cluster was used for generating the initial docking structures. These were generated by fitting of benzene clusters formed on the surface of receptor and ligand. In CAPRI Rounds 3-5, the grid-scoring sum on the protein-protein interaction surface and the pairwise potential of the amino acid residues, which were indicated as coming easily into the protein-protein interaction regions, were used as the calculation methods, along with the smaller benzene clusters that participated in benzene cluster fitting. Good predicted models were obtained for Targets 11 and 12. When the modeled receptor proteins were superimposed on the experimental structures, the smallest ligand root-mean-square deviation (RMSD) values corresponding to the RMSD between the model and experimental structures were 6.2 A and 7.3 A, respectively.

Algorithms↗

PubMed1995

Protein docking combining symbolic descriptions of molecular surfaces and grid-based scoring functions.

With the growing number of known 3D protein structures, computing systems, that can predict where two protein molecules interact with each other is becoming of increasing interest. A system is presented, integrating preprocessing like the computation of molecular surfaces, segmentation, and searching for complementarity in the general framework of a pattern analyzing semantic network (ERNEST). The score of coarse symbolic computations is used by the problem independent control strategy of ERNEST to guide a more detailed analysis considering steric clash and judgements based on grid-based surface representations. Successful examples of the docking system are discussed that compare well with other approaches.

Binding Sites↗

PubMed2000

The third-dimensional structure of the complex between an Fv antibody fragment and an analogue of the main immunogenic region of the acetylcholine receptor: a combined two-dimensional NMR, homology, and molecular modeling approach.

Binding of autoantibodies to the acetylcholine receptor (AChR) plays a major role in the autoimmune disease Myasthenia gravis (MG). In this paper, we propose a structure model of a putative immunocomplex that gives rise to the reduction of functional AChR molecules during the course of MG. The model complex consists of the [G(70), Nle(76)] decapeptide analogue of the main immunogenic region (MIR), representing the major antigenic epitope of AChR, and the single chain Fv fragment of monoclonal antibody 198, a potent MG autoantibody. The structure of the complexed decapeptide antigen [G(70), Nle(76)]MIR was determined using two-dimensional nmr, whereas the antibody structure was derived by means of homology modeling. The final complex was constructed using calculational docking and molecular dynamics. We termed this approach "directed modeling," since the known peptide structure directs the prestructured antibody binding site to its final conformation. The independently derived structures of the peptide antigen and antibody binding site already showed a high degree of surface complementarity after the initial docking calculation, during which the peptide was conformationally restrained. The docking routine was a soft algorithm, applying a combination of Monte Carlo simulation and energy minimization. The observed shape complementarity in the docking process suggested that the structure assessments already led to anti-idiotypic conformations of peptide antigen and antibody fragment. Refinement of the complex by dynamic simulation yielded improved surface adaptation by small rearrangements within antibody and antigen. The complex presented herein was analyzed in terms of antibody-antigen interactions, properties of contacting surfaces, and segmental mobility. The structural requirements for AChR complexation by autoantibodies were explored and compared with experimental data from alanine scans of the MIR peptides. The analysis revealed that the N-terminal loop of the peptide structure, which is indispensable for antibody recognition, aligns three hydrophobic groups in a favorable arrangement leading to the burial of 40% of the peptide surface in the binding cleft upon complexation. These data should be valuable in the rational design of an Fv mutant with much improved affinity for the MIR and AChR to be used in therapeutic approaches in MG.

Amino Acid Sequence↗

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↗