PubMed Health⌕ Search

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 361 records · Page 20Linked to original sources
PubMed2001

Monosaccharide-linked inhibitors of O(6)-methylguanine-DNA methyltransferase (MGMT): synthesis, molecular modeling, and structure-activity relationships.

A series of potential inhibitors of the human DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) were synthesized, characterized in detail by NMR, and tested for their ability to deplete MGMT activity in vitro. The new compounds, omega-[O(6)-R-guan-9-yl]-(CH(2))(n)-beta-d-glucosides with R = benzyl or 4-bromothenyl and omega = n = 2, 4,. 12, were compared with the established inhibitors O(6)-benzylguanine (O(6)-BG), 8-aza-O(6)-benzylguanine (8-aza-BG), and O(6)-(4-bromothenyl)guanine (4-BTG), which exhibit in an in vitro assay IC(50) values of 0.62, 0.038, and 0.009 microM, respectively. Potential advantages of the glucosides are improved water solubility and selective uptake in tumor cells. The 4-BTG glucosides with n = 2, 4, 6 show moderate inhibition with an IC(50) of ca. 0.5 microM, while glucosides derived from BG and 8-aza-BG showed significantly poorer inhibition compared to the parent compounds. The 4-BTG glucosides with n = 8, 10, 12 were effective inhibitors with IC(50) values of ca. 0.03 microM. To understand this behavior, extensive molecular modeling studies were performed using the published crystal structure of MGMT (PDB entry: ). The inhibitor molecules were docked into the BG binding pocket, and molecular dynamics simulations with explicit water molecules were carried out. Stabilization energies for the interactions of specific regions of the inhibitor and individual amino acid residues were calculated. The alkyl spacer is located in a cleft along helix 6 of MGMT. With increasing spacer length there is increasing interaction with several amino acid residues which play an important role in the proposed nucleotide flipping mechanism required for DNA repair.

Cell-Free System↗
PubMed2006

Understanding the bacterial polysaccharide antigenicity of Streptococcus agalactiae versus Streptococcus pneumoniae.

Bacterial surface capsular polysaccharides (CPS) that are similar in carbohydrate sequence may differ markedly in immunogenicity and antigenicity. The structural origin of these phenomena is poorly understood. Such a case is presented by the Gram-positive bacteria Streptococcus agalactiae (Group B Streptococcus; GBS) type III (GBSIII) and Streptococcus pneumoniae (Pn) type 14 (Pn14), which share closely related CPS sequences. Nevertheless, antibodies (Abs) against GBSIII rarely cross-react with the CPS from Pn14. To establish the origin for the variation in CPS antigenicity, models for the immune complexes of CPS fragments from GBSIII and Pn14, with the variable fragment (Fv) of a GBS-specific mAb (mAb 1B1), are presented. The complexes are generated through a combination of comparative Ab modeling and automated ligand docking, followed by explicitly solvated 10-ns molecular dynamics simulations. The relationship between carbohydrate sequence and antigenicity is further quantified through the computation of interaction energies using the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method, augmented by conformational entropy estimates. Despite the electrostatic differences between Pn14 and GBSIII CPS, analysis indicates that entropic penalties are primarily responsible for the loss of affinity of the highly flexible Pn14 CPS for mAb 1B1. The similarity of the solution conformation of the relatively rigid GBSIII CPS with that in the immune complex characterizes the previously undescribed 3D structure of the conformational epitope. The analysis provides a comprehensive interpretation for a large body of biochemical and immunological data related to Ab recognition of bacterial polysaccharides and should be applicable to other Ab-carbohydrate interactions.

Amino Acid Sequence↗
PubMed2002

The X-ray crystallographic study of long chain left handed Z-DNA and polyamine complex and structural chemistry study of effection by polyamine for the Z-DNA.

We cleared the stabilization mechanism of Z-DNA by using the molecular dynamics calculation and the grid map simulated annealing docking method. Then, we clarified Z-DNA was stabilized by the metal counter ions whose electrical charge and the number of electrons are too many and by polyamine to go into the minor groove. We succeeded in stabilizing Z-DNA of the long chain by using these methods, and succeeded in the crystallization and the data collection up to 1.7 A resolution.

Crystallography, X-Ray↗
PubMed1996

Exploring the energy landscapes of molecular recognition by a genetic algorithm: analysis of the requirements for robust docking of HIV-1 protease and FKBP-12 complexes.

Energy landscapes of molecular recognition are explored by performing "semi-rigid" docking of FK-506 and rapamycin with the Fukisawa binding protein (FKBP-12), and flexible docking simulations of the Ro-31-8959 and AG-1284 inhibitors with HIV-1 protease by a genetic algorithm. The requirements of a molecular recognition model to meet thermodynamic and kinetic criteria of ligand-protein docking simultaneously are investigated using a family of simple molecular recognition energy functions. The critical factor that determines the success rate in predicting the structure of ligand-protein complexes is found to be the roughness of the binding energy landscape, in accordance with a minimal frustration principle. The results suggest that further progress in structure prediction of ligand-protein complexes can be achieved by designing molecular recognition energy functions that generate binding landscapes with reduced frustration.

Algorithms↗
PubMed2005

Molecular mechanics PBSA ligand binding energy and interaction of Efavirenz derivatives with HIV-1 reverse transcriptase.

In order to evaluate the properties of several HIV-1 reverse transcripase(RT) inhibitors, Efavirenz (SUSTIVA) and a set of its derivatives (benzoxazinones) have been placed into the nonnucleoside analogue binding site of the enzyme by molecular docking. The resulting geometries were used for a molecular dynamics simulation and binding energy calculations. The enzyme-inhibitor binding energies were estimated from experimental inhibitory activities (IC90). The correlation of the predicted and experimental binding energies were satisfactory acceptable as indicated by r2 = 0.865. Based on MD simulations, the obtained results indicate that the tight association of the ligand to the HIV-1 RT binding pocket was based on hydrogen bonding between Efavirenz's N1 and the oxygen of the backbone of Lys 101, with an estimated average distance of 1.88 A. Moreover, electrostatic interaction was mainly contributed by two amino acid residues in the binding site; Lys 101 and His 235. MD simulations open the possibility to study the reaction of the flexible enzyme to those substances as well as the overall affinity.

Adipates↗
PubMed2006

Protein-protein interactions: modeling the hepatitis C virus ion channel p7.

The p7 protein is a small ion-channel-forming membrane polypeptide encoded by the hepatitis C virus which consists of two transmembrane alpha-helices, TM1 and TM2, and can be blocked by long-alkyl-chain iminosugar derivatives. The length of TM1 and TM2 was estimated by employing different secondary structure prediction algorithms and is proposed to span from Ala-10 to Leu-32 for TM1 and from Trp-36 to Pro-58 for TM2. A configurational search protocol based on simulated annealing combined with short restrained molecular dynamics simulations is used in addition to protein-protein docking to investigate the packing of TM1/TM2. Full p7 oligomeric bundles were generated, and in the most plausible models serines and threonines are facing the hydrophilic pore. In these models, His-17 would be a pore-facing residue, suggesting that p7 may be sensitive to pH in respect to its function.

Amantadine↗
PubMed2002

Kappa-opioid receptor model in a phospholipid bilayer: molecular dynamics simulation.

A three-dimensional molecular model of the transmembrane domain of the kappa-opioid receptor in a phospholipid bilayer is presented. The endogenous ligand, dynorphin A (1), and synthetic ligands, benzomorphan-based compounds (2a, 2b) (Figure 1), are docked into the model. We report the results of a 500 ps molecular dynamics simulation of these protein-ligand complexes in a simplified bilayer of 97 molecules of the lipid dipalmitoylphosphatidylcholine and 26 water molecules per lipid. The simulations explore the stability and conformational dynamics of the model in a phospholipid bilayer; we also investigate the interactions of the protein with its ligands. Molecular simulation of the receptor-ligand complexes, endogenous and synthetic, has confirmed the existence of different binding domains for peptide and non-peptide ligands. Similarities are found in the dynamics and binding mode of all conformations of the synthetic ligands studied. The protonated hydrogen of the benzomorphan is always involved in an H-bond with Asp138, and other potentially stabilizing receptor-ligand interactions found involve the hydroxyl substituent on the benzomorphan, which may form an H-bond with Tyr139 or Gly190 according to the different molecules. The ester group of 2a may therefore form an H-bond with Ile316, while the carbonyl group of 2b forms an H-bond with Gln115 and Tyr312. The remaining part of the ligand is located in the extracellular portion of the pocket. It is surrounded by hydrophobic residues in the transmembrane region (TM), and it interacts with different sets of residues. The results obtained are in general agreement with site-directed mutagenesis data that have highlighted the importance of all TM regions for synthetic-ligand affinity with the kappa-opioid receptor.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2007

Photoselected electron transfer pathways in DNA photolyase.

Cyclobutane dimer photolyases are proteins that bind to UV-damaged DNA containing cyclobutane pyrimidine dimer lesions. They repair these lesions by photo-induced electron transfer. The electron donor cofactor of a photolyase is a two-electron-reduced flavin adenine dinucleotide (FADH(-)). When FADH(-) is photo-excited, it transfers an electron from an excited pi --> pi* singlet state to the pyrimidine dimer lesion of DNA. We compute the lowest excited singlet states of FADH(-) using ab initio (time-dependent density functional theory and time-dependent Hartree-Fock), and semiempirical (INDO/S configuration interaction) methods. The calculations show that the two lowest pi --> pi* singlet states of FADH(-) are localized on the side of the flavin ring that is proximal to the dimer lesion of DNA. For the lowest-energy donor excited state of FADH(-), we compute the conformationally averaged electronic coupling to acceptor states of the thymine dimer. The coupling calculations are performed at the INDO/S level, on donor-acceptor cofactor conformations obtained from molecular dynamics simulations of the solvated protein with a thymine dimer docked in its active site. These calculations demonstrate that the localization of the (1)FADH(-)* donor state on the flavin ring enhances the electronic coupling between the flavin and the dimer by permitting shorter electron-transfer pathways to the dimer that have single through-space jumps. Therefore, in photolyase, the photo-excitation itself enhances the electron transfer rate by moving the electron towards the dimer.

Computer Simulation↗
PubMed2002

Further development and validation of empirical scoring functions for structure-based binding affinity prediction.

New empirical scoring functions have been developed to estimate the binding affinity of a given protein-ligand complex with known three-dimensional structure. These scoring functions include terms accounting for van der Waals interaction, hydrogen bonding, deformation penalty, and hydrophobic effect. A special feature is that three different algorithms have been implemented to calculate the hydrophobic effect term, which results in three parallel scoring functions. All three scoring functions are calibrated through multivariate regression analysis of a set of 200 protein-ligand complexes and they reproduce the binding free energies of the entire training set with standard deviations of 2.2 kcal/mol, 2.1 kcal/mol, and 2.0 kcal/mol, respectively. These three scoring functions are further combined into a consensus scoring function, X-CSCORE. When tested on an independent set of 30 protein-ligand complexes. X-CSCORE is able to predict their binding free energies with a standard deviation of 2.2 kcal/mol. The potential application of X-CSCORE to molecular docking is also investigated. Our results show that this consensus scoring function improves the docking accuracy considerably when compared to the conventional force field computation used for molecular docking.

Binding Sites↗
PubMed2006

Towards the design of highly selective recognition sites into molecular imprinting polymers: a computational approach.

A computational approach to simulate the formation of possible imprinted polymers in acetonitrile solution for theophylline (THO) is proposed, using combined molecular dynamics (MD), molecular mechanics (MM), docking and site mapping computational techniques. Methacrylic acid (MAA) and methylmethacrylate (MMA) monomers are used to simulate possible homo and copolymer structures. The model is able predict binding affinity and selectivity when considering THO analogues, such as caffeine, theobromine, xanthine and 3-methylxanthine. Comparison with available experimental data is proposed.

Biosensing Techniques↗
PubMed2006

Toxin insights into nicotinic acetylcholine receptors.

Venomous species have evolved cocktails of bioactive peptides to facilitate prey capture. Given their often exquisite potency and target selectivity, venom peptides provide unique biochemical tools for probing the function of membrane proteins at the molecular level. In the field of the nicotinic acetylcholine receptors (nAChRs), the subtype specific snake alpha-neurotoxins and cone snail alpha-conotoxins have been widely used to probe receptor structure and function in native tissues and recombinant systems. However, only recently has it been possible to generate an accurate molecular view of these nAChR-toxin interactions. Crystal structures of AChBP, a homologue of the nAChR ligand binding domain, have now been solved in complex with alpha-cobratoxin, alpha-conotoxin PnIA and alpha-conotoxin ImI. The orientation of all three toxins in the ACh binding site confirms many of the predictions obtained from mutagenesis and docking simulations on homology models of mammalian nAChR. The precise understanding of the molecular determinants of these complexes is expected to contribute to the development of more selective nAChR modulators. In this commentary, we review the structural data on nAChR-toxin interactions and discuss their implications for the design of novel ligands acting at the nAChR.

Animals↗
PubMed2002

Boar salivary lipocalin. Three-dimensional X-ray structure and androsterol/androstenone docking simulations.

The X-ray structure of variant A of authentic boar salivary lipocalin (SAL), a pheromone-binding protein specifically expressed in the submaxillary glands of the boar, has been solved and refined at 2.1 A resolution. The structure displays a classical lipocalin fold with a nine-stranded sandwiched beta barrel and an alpha helix. A putative glycosylation site, at position 53, has been found to carry a GlcNAc sugar residue. In contrast with what was expected on the basis of mass spectroscopy reports, the internal cavity was found to be devoid of bound pheromonal compound (androstenone or androstenol). Instead, a small electron density volume could be satisfied by a glycerol molecule, a component of the cryoprotecting liquor. The internal cavity was revealed to be very small for steroid compound accommodation. Therefore, docking and molecular dynamics experiments were performed with both pheromonal compounds. These simulations clearly demonstrate a volume increase of the cavity upon steroid binding and the adaptation of the amino-acid side chains to the steroid molecules. This explains the higher affinity of SAL for both steroid molecules compared to other smaller molecules, although no specific interaction is established with either compound.

Amino Acid Sequence↗
PubMed1999

A proposed molecular model for the interaction of calcineurin with the cyclosporin A-cyclophilin A complex.

Cyclosporin A (CsA) and FK506 are potent natural product immunosuppressants that induce their biological effects by forming an initial complex with cytosolic proteins termed immunophilins. These drug immunophilin complexes then bind to and inhibit the serine/threonine protein phosphatase calcineurin (CN). Two classes of immunophilin have been identified with cyclophilins (CyP's) being proteins specifically binding CsA and FKBPs specifically binding FK506. Solution and crystal structures of various CsA-CyP and FK506-FKBP complexes have been determined and show no apparent structural similarity between the two classes of drug protein complexes. These findings raise the question as to how, given their structural differences, these two complexes can both inhibit CN. While the crystal structure of the FK506-FKBP12-CN complex has been reported, no structure for a CsA-CyP CN complex has been determined. Here are reported studies that use various modelling strategies to construct a model for the interaction of the cyclosporin A- cyclophilin A complex with calcineurin. The first stage of constructing this model consisted of using conformational comparison of CsA and FK506, GRID and GROUP analysis and restrained molecular dynamics to dock CsA into the FK506 binding site of the FK506-FKBP12-CN structure. An initial model for the CsA-CyPA-CN complex was then constructed by superimposing the structure of the CsA-CyPA complex onto the docked CsA molecule. This model was then optimised with molecular dynamics simulations run on sterically clashing regions. The validity of the model for the CsA-CyPA-CN complex was then examined with respect to the effect of chemical modifications to CsA and amino acid substitutions within CyPA on the ability of the drug-immunophilin complex to inhibit calcineurin.

Amino Acid Substitution↗
PubMedDate not supplied

WIN55212-2 docking to the CB1 cannabinoid receptor and multiple pathways for conformational induction.

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.

Benzoxazines↗
PubMed2006

Computational study of antagonist/alpha1A adrenoceptor complexes--observations of conformational variations on the formation of ligand/receptor complexes.

As selective antagonist inhibition may relieve the symptoms of benign prostatic hyperplasia, we have examined the interactions of antagonists including quinazoline and imidazolidinium/guanidinium compounds complexed with a homology model of the alpha(1A) adrenoceptor. Our approach involves docking of ligands of various structural classes followed by molecular dynamics simulations of antagonist/receptor complexes, which demonstrates that different structural classes of antagonist induce different receptor conformations upon binding with particular variations noted in the conformation of TM-V. Subsequently, we examined the interactions and the conformational flexibility of alpha(1) and alpha(1A) adrenoceptor antagonists, with the ligand-induced receptor conformers. This study indicated that a receptor conformation induced by one structural class of antagonist is not suitable for direct screening of another class. Our analysis indicates that computational high-throughput screening is likely to give inaccurate data on binding and selectivity and such studies need to consider conformational changes in the receptor.

Adrenergic alpha-1 Receptor Antagonists↗
PubMed2004

Molecular recognition in purinergic receptors. 1. A comprehensive computational study of the h-P2Y1-receptor.

P2Y receptors (P2Y-Rs) are attractive pharmaceutical targets due to their involvement in the modulation of many tissues and organs. The lack of experimental structural data on P2Y-Rs impedes structure-based drug design. The need to elucidate the receptor's molecular recognition, together with the limitations of previous receptor models, triggered the construction of a new molecular model for the h-P2Y1-R. Therefore, a h-P2Y1-R model was constructed by homology modeling using the 2.6 A crystal structure of bovine rhodopsin as a template and subsequently refined by constrained molecular dynamics (MD) simulations in a fully hydrated lipid bilayer environment. ATP was docked into the receptor binding site, followed by binding site refinement using Monte Carlo and MD simulations. Analysis of the h-P2Y1-R-ATP complex suggests that the triphosphate moiety is tightly bound by a multitude of interactions possibly including a Mg2+ ion, the ribose ring is not involved in specific interactions, and the adenine ring is bound via N1, N7, and N6. The molecular recognition of the h-P2Y1-R was further probed by ATP derivatives modified on the adenine ring, and correlated with EC50 values for these derivatives. Analysis of receptor:ligand complexes and quantum mechanical studies on model compounds support the role of both steric and electronic effects in improving H-bonding (via N1 and N6) and pi-stacking interactions. The computed h-P2Y1-R model was validated with respect to our previous biochemical results. We believe that this new model of the h-P2Y1-R provides the means for understanding phenomena such as the ligand's potency and receptor subtype selectivity.

Adenosine Triphosphate↗
PubMed1998

Molecular anchors with large stability gaps ensure linear binding free energy relationships for hydrophobic substituents.

Ligand-protein docking simulations are employed to analyze the binding energy landscape of the pipecolinyl fragment that serves as a recognition core of the FK506 ligand in binding with the FKBP12 protein. This fragment acts as a molecular anchor that specifically binds within the protein active site in a unique binding mode, in harmony with the structure of the FK506-FKBP12 complex. Molecular anchors are characterized by a large stability gap, defined to be the free energy of a ligand bound in the native binding mode relative to the free energy of alternative binding modes. For ligands that share a common anchor fragment, a linear binding free energy relationship may be expected for hydrophobic substituents provided they do not abrogate the anchor binding mode. Changes in solvent-accessible surface area for these peripheral groups are used to rationalize the relative binding affinities of a series of FKBP12-ligand complexes which share the pipecolinyl anchor fragment. A series of benzene derivatives that bind to a mutant form of T4 lysozyme is also analyzed, and implications for structure-based drug design are discussed.

Allosteric Site↗