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 343 records · Page 19Linked to original sources

PubMed2006

Molecular docking studies on 4-thiazolidinones as HIV-1 RT inhibitors.

Flexible docking simulations were performed on two series of 4-thiazolidinones as HIV-1 reverse transcriptase (HIV-1 RT) inhibitors. This was done by analyzing the interaction of these compounds with the allosteric site of the HIV-1 reverse transcriptase enzyme. The binding scores for these compounds were also congruent with their anti-HIV activity. A good correlation between the predicted binding free energies and the experimentally observed inhibitory activities (EC(50)) suggest that the identified binding conformations of these inhibitors are reliable. The results of docking studies provide an insight into the pharmacophoric structural requirements for the HIV-1 RT inhibitory activity of this class of molecules.

Allosteric Site↗

PubMed2004

Computational prediction of structure, substrate binding mode, mechanism, and rate for a malaria protease with a novel type of active site.

The histo-aspartic protease (HAP) from the malaria parasite P. falciparum is one of several new promising targets for drug intervention. The enzyme possesses a novel type of active site, but its 3D structure and mechanism of action are still unknown. Here we use a combination of homology modeling, automated docking searches, and molecular dynamics/reaction free energy profile simulations to predict the enzyme structure, conformation of bound substrate, catalytic mechanism, and rate of the peptide cleavage reaction. We find that the computational tools are sufficiently reliable both for identifying substrate binding modes and for distinguishing between different possible reaction mechanisms. It is found that the favored pathway only involves direct participation by the catalytic aspartate, with the neighboring histidine providing critical stabilization (by a factor of approximately 10000) along the reaction. The calculated catalytic rate constant of about 0.1 s(-1) for a hexapeptide substrate derived from the alpha chain of human hemoglobin is in excellent agreement with experimental kinetic data for a similar peptide fragment.

Animals↗

PubMed2005

Experimental validation of the docking orientation of Cdc25 with its Cdk2-CycA protein substrate.

Cdc25 phosphatases are key activators of the eukaryotic cell cycle and compelling anticancer targets because their overexpression has been associated with numerous cancers. However, drug discovery targeting these phosphatases has been hampered by the lack of structural information about how Cdc25s interact with their native protein substrates, the cyclin-dependent kinases. Herein, we predict a docked orientation for Cdc25B with its Cdk2-pTpY-CycA protein substrate by a rigid-body docking method and refine the docked models with full-scale molecular dynamics simulations and minimization. We validate the stable ensemble structure experimentally by a variety of in vitro and in vivo techniques. Specifically, we compare our model with a crystal structure of the substrate-trapping mutant of Cdc25B. We identify and validate in vivo a novel hot-spot residue on Cdc25B (Arg492) that plays a central role in protein substrate recognition. We identify a hot-spot residue on the substrate Cdk2 (Asp206) and confirm its interaction with hot-spot residues on Cdc25 using hot-spot swapping and double mutant cycles to derive interaction energies. Our experimentally validated model is consistent with previous studies of Cdk2 and its interaction partners and initiates the opportunity for drug discovery of inhibitors that target the remote binding sites of this protein-protein interaction.

Binding Sites↗

PubMed2002

A model of the human M2 muscarinic acetylcholine receptor.

The M2 muscarinic acetylcholine receptor belongs to the family of rhodopsin like G-Protein Coupled Receptors. This subtype of muscarinic receptors is of special interest because it bears, aside from an orthosteric binding site, also an allosteric binding site. Based on the X-ray structure of bovine rhodopsin a complete homology model of the human M2 receptor was developed. For the orthosteric binding site point mutations and binding studies with different agonists and antagonists are available. This knowledge was utilized for an initial verification of the M2 model. Allosteric modulation of activity is mediated by structurally different ligands such as gallamine, caracurine V salts or W84 (a hexamethonium-derivative). Caracurine V derivatives with different affinities to M2 were docked using GRID-fields. Subsequent molecular dynamics simulations yielded different binding energies based on diverse electrostatic and lipophilic interactions. The calculated affinities are in good agreement to experimentally determined affinities.

Allosteric Site↗

PubMed2006

Modeling the structure of the StART domains of MLN64 and StAR proteins in complex with cholesterol.

Steroidogenic acute regulatory protein-related lipid transfer (StART) domains are ubiquitously involved in intracellular lipid transport and metabolism and other cell-signaling events. In this work, we use a flexible docking algorithm, comparative modeling, and molecular dynamics (MD) simulations to generate plausible three-dimensional atomic models of the StART domains of human metastatic lymph node 64 (MLN64) and steroidogenic acute regulatory protein (StAR) proteins in complex with cholesterol. Our results show that cholesterol can adopt a similar conformation in the binding cavity in both cases and that the main contribution to the protein-ligand interaction energy derives from hydrophobic contacts. However, hydrogen-bonding and water-mediated interactions appear to be important in the fine-tuning of the binding affinity and the position of the ligand. To gain insights into the mechanism of binding, we carried out steered MD simulations in which cholesterol was gradually extracted from within the StAR model. These simulations indicate that a transient opening of loop Omega1 may be sufficient for uptake and release, and they also reveal a pathway of intermediate states involving residues known to be crucial for StAR activity. Based on these observations, we suggest specific mutagenesis targets for binding studies of cholesterol and its derivatives that could improve our understanding of the structural determinants for ligand binding by sterol carrier proteins.

Algorithms↗

PubMed2004

SDOCKER: a method utilizing existing X-ray structures to improve docking accuracy.

This paper introduces a new strategy for structure-based drug design that combines high-quality docking with data from existing ligand-protein cocrystal X-ray structures. The main goal of SDOCKER, a new algorithm that implements this strategy, is docking accuracy improvement. In this new paradigm, simulated annealing molecular dynamics is used for conformational sampling and optimization and an additional similarity force is applied on the basis of the positions of ligands from X-ray data that focus the sampling on relevant regions of the active site. Because the structural information from both the ligand and protein active site is included, this approach is more effective in finding the optimal conformation for a ligand-protein complex than the classical docking or similarity overlays. Interestingly, it was found that a 3D similarity-only approach gives comparable docking accuracy to the regular force field approach used in classical docking, given the final structures are minimized in the presence of the protein. The combination of both, as implemented in SDOCKER, is shown here to be more accurate. A significant improvement in docking accuracy has been observed for three different test systems. Specifically an improvement of 10%, 17.5%, and 10% is seen for 37 HIV-1 protease, 32 thrombin, and 23 CDK2 ligands, respectively, compared to docking using the force field alone. In addition, SDOCKER's accuracy performance dependence on the similarity template is discussed. The strategy of utilizing existing ligand X-ray information should prove effective in light of the multitude of structures available from structural genomics approaches.

Algorithms↗

PubMed1998

Modeling of angiogenin - 3-NMP complex.

Angiogenin belongs to the Ribonuclease superfamily and has a weak enzymatic activity that is crucial for its biological function of stimulating blood vessel growth. Structural studies on ligand bound Angiogenin will go a long way in understanding the mechanism of the protein as well as help in designing drugs against it. In this study we present the first available structure of nucleotide ligand bound Angiogenin obtained by computer modeling. The importance of this study in itself notwithstanding, is a precursor to modeling a full dinucleotide substrate onto Angiogenin. Bovine Angiogenin, the structure of which has been solved at a high resolution, was earlier subjected to Molecular Dynamics simulations for a nanosecond. The MD structures offer better starting points for docking as they offer lesser obstruction than the crystal structure to ligand binding. The MD structure with the least serious short contacts was modeled to obtain a steric free Angiogenin - 3' mononucleotide complex structure. The structures were energetically minimized and subjected to a brief spell of Molecular Dynamics. The results of the simulation show that all the ligand-Angiogenin interactions and hydrogen bonds are retained, redeeming the structure and docking procedure. Further, following ligand - protein interactions in the case of the ligands 3'-CMP and 3'-UMP we were able to speculate on how Angiogenin, a predominantly prymidine specific ribonuclease prefers Cytosine to Uracil in the first base position.

Animals↗

PubMed2003

Detailed analysis of grid-based molecular docking: A case study of CDOCKER-A CHARMm-based MD docking algorithm.

The influence of various factors on the accuracy of protein-ligand docking is examined. The factors investigated include the role of a grid representation of protein-ligand interactions, the initial ligand conformation and orientation, the sampling rate of the energy hyper-surface, and the final minimization. A representative docking method is used to study these factors, namely, CDOCKER, a molecular dynamics (MD) simulated-annealing-based algorithm. A major emphasis in these studies is to compare the relative performance and accuracy of various grid-based approximations to explicit all-atom force field calculations. In these docking studies, the protein is kept rigid while the ligands are treated as fully flexible and a final minimization step is used to refine the docked poses. A docking success rate of 74% is observed when an explicit all-atom representation of the protein (full force field) is used, while a lower accuracy of 66-76% is observed for grid-based methods. All docking experiments considered a 41-member protein-ligand validation set. A significant improvement in accuracy (76 vs. 66%) for the grid-based docking is achieved if the explicit all-atom force field is used in a final minimization step to refine the docking poses. Statistical analysis shows that even lower-accuracy grid-based energy representations can be effectively used when followed with full force field minimization. The results of these grid-based protocols are statistically indistinguishable from the detailed atomic dockings and provide up to a sixfold reduction in computation time. For the test case examined here, improving the docking accuracy did not necessarily enhance the ability to estimate binding affinities using the docked structures.

Algorithms↗

PubMed2004

A hybrid method of molecular dynamics and harmonic dynamics for docking of flexible ligand to flexible receptor.

We have developed a new docking method to consider receptor flexibility, a hybrid method of molecular dynamics and harmonic dynamics. The global motions of the whole receptor were approximately introduced into those of the receptor in the docking simulation as harmonic dynamics. On the other hand, the local flexibility of the side chains was also considered by conventional molecular dynamics. We confirmed that this new method can reproduce the fluctuations of the whole receptor by making a comparison of the directions and amplitudes of the global fluctuations. Then this method was applied to the docking of HIV-1 protease and its ligand. As a result, we observed a docking process where the ligand enters into the binding pocket well, which implies that this method is effective enough to reproduce a molecular complex formation.

Algorithms↗

PubMed2004

Simulation of the interaction between ScyTx and small conductance calcium-activated potassium channel by docking and MM-PBSA.

Computational methods are employed to simulate interaction of scorpion toxin ScyTx in complex with the small conductance calcium-activated potassium channel rsk2. All of available 25 structures of ScyTx in the Protein Data Bank determined by NMR were considered for improving performance of rigid protein docking of ZDOCK. Four main binding modes were found among a large number of predicted complexes by using clustering analysis, screening with expert knowledge, energy minimization, and molecular dynamics simulations. The quality and validity of the resulting complexes were further evaluated by molecular dynamics simulations with the generalized Born solvation model and by calculation of relative binding free energies with the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) in the AMBER 7 suit of programs. The complex formed by the 22nd structure of the ScyTx and rsk2 channel was identified as the most favorable complex by using a combination of computational methods, which contain further introduction of flexibility without restraining residue side chain. From the resulted spatial structure of the ScyTx and rsk2 channel, ScyTx associates the mouth of the rsk2 channel with alpha-helix rather than beta-sheet. Structural analysis first revealed that Arg(13) played a novel and vital role of blocking the pore of the rsk2 channel, whose role is remarkably different from that of highly homologous scorpion toxin P05. Between the interfaces in the ScyTx-rsk2 complex, strong electrostatic interaction and hydrogen bonds exist between Arg(13) of ScyTx and Gly-Tyr-Gly-Asp sequential residues located in the four symmetrical chains of the pore region. Simultaneously, five hydrogen bonds between Arg(6) of ScyTx and Asp(341)(C), Val(366)(C), and Pro(367)(C), and electrostatic interaction between Arg(6) of ScyTx and Asp(364)(B) and Asp(341)(C) are also found by structural analysis. In addition, His(31) located at the C-terminal of ScyTx is surrounded by Val(342)(A), Asp(364)(A), Met(365)(A), Pro(367)(B), and Asn(366)(B) within a contact distance of 4.0 A. These simulation results are in good agreement with experimental data and can effectively explain the binding phenomena between ScyTx and the potassium channel at the level of molecular spatial structure. The consistency between results of molecular modeling and experimental data strongly suggests that our spatial structure model of the ScyTx-rsk2 complex is reasonable. Therefore, molecular docking combined with molecular dynamics simulations followed by molecular mechanics Poisson-Boltzmann surface area analysis is an attractive approach for modeling scorpion toxin-potassium channel complexes a priori for further biological studies.

Animals↗

PubMed2004

Molecular modeling study of the editing active site of Escherichia coli leucyl-tRNA synthetase: two amino acid binding sites in the editing domain.

Aminoacyl-tRNA synthetases (aaRSs) strictly discriminate their cognate amino acids. Some aaRSs accomplish this via proofreading and editing mechanisms. Mursinna and coworkers recently reported that substituting a highly conserved threonine (T252) with an alanine within the editing domain of Escherichia coli leucyl-tRNA synthetase (LeuRS) caused LeuRS to cleave its cognate aminoacylated leucine from tRNA(Leu) (Mursinna et al., Biochemistry 2001;40:5376-5381). To achieve atomic level insight into the role of T252 in LeuRS and the editing reaction of aaRSs, a series of molecular modeling studies including homology modeling and automated docking simulations were carried out. A 3D structure of E. coli LeuRS was constructed via homology modeling using the X-ray structure of Thermus thermophilus LeuRS as a template because the E. coli LeuRS structure is not available from X-ray or NMR studies. However, both the X-ray T. thermophilus and homology-modeled E. coli structures were used in our studies. Amino acid binding sites in the proposed editing domain, which is also called the connective polypeptide 1 (CP1) domain, were investigated by automated docking studies. The root mean square deviation (RMSD) for backbone atoms between the X-ray and homology-modeled structures was 1.18 A overall and 0.60 A for the editing (CP1) domain. Automated docking studies of a leucine ligand into the editing domain were performed for both structures: homology structure of E. coli LeuRS and X-ray structure of T. thermophilus LeuRS for comparison. The results of the docking studies suggested that there are two possible amino acid binding sites in the CP1 domain for both proteins. The first site lies near a threonine-rich region that includes the highly conserved T252 residue, which is important for amino acid discrimination. The second site is located in a flexible loop region surrounded by residues E292, A293, M295, A296, and M298. The important T252 residue is at the bottom of the first binding pocket.

Adenosine Monophosphate↗

PubMed2005

The 3D structure of the binding pocket of the human oxytocin receptor for benzoxazine antagonists, determined by molecular docking, scoring functions and 3D-QSAR methods.

Molecular docking and 3D-QSAR studies were performed to determine the binding mode for a series of benzoxazine oxytocin antagonists taken from the literature. Structural hypotheses were generated by docking the most active molecule to the rigid receptor by means of AutoDock 3.05. The cluster analysis yielded seven possible binding conformations. These structures were refined by using constrained simulated annealing, and the further ligands were aligned in the refined receptor by molecular docking. A good correlation was found between the estimated deltaG(bind) and the pKi values for complex F. The Connolly-surface analysis, CoMFA and CoMSIA models q2(CoMFA) = 0.653, q2(CoMSA) = 0.630 and r2(pred,CoMFA) = 0.852 , r2(pred,CoMSIA) = 0.815) confirmed the scoring function results. The structural features of the receptor-ligand complex and the CoMFA and CoMSIA fields are in closely connected. These results suggest that receptor-ligand complex F is the most likely binding hypothesis for the studied benzoxazine analogs.

Benzoxazines↗

PubMed2003

Homology model of the CB1 cannabinoid receptor: sites critical for nonclassical cannabinoid agonist interaction.

Association of cannabimimetic compounds such as cannabinoids, aminoalkylindoles (AAIs), and arachidonylethanolamide (anandamide) with the brain cannabinoid (CB(1)) receptor activates G-proteins and relays signals to regulate neuronal functions. A CB(1) receptor homology model was constructed using the published x-ray crystal structure of bovine rhodopsin (Palczewski et al., Science, 2000, Vol. 289, pp. 739-745) in the conformation most likely to represent the "high-affinity" state for agonist binding to G-protein coupled receptors (GPCRs). A molecular docking approach that combined Monte Carlo and molecular dynamics simulations was used to identify the putative binding conformations of nonclassical cannabinoid agonists, including AC-bicyclic CP47497 and CP55940, and ACD-tricyclic CP55244. Placement of these ligands was based upon the assumption of a critical hydrogen bond between the A-ring OH and the side chain N of Lys192 in transmembrane helix 3. We evaluated two alternative binding conformations, C3-in and C3-out, denoting the directionality of the ligand C3 side chain within the receptor with respect to the inside or the outside of the cell. Assuming both the C3-in or C3-out conformation, the calculated ligand-receptor binding energy (DeltaE(bind)) was correlated with the experimentally observed binding affinity (K(i)) for a series of nonclassical cannabinoid agonists. The C3-in conformation was marginally better than the alternative C3-out conformation in predicting the rank order of the tested nonclassical cannabinoid analogs. Adopting the C3-in conformation due to the greater number of receptor interactions with known pharmacophoric elements of the ligand, key residues were identified comprising the presumed hydrophobic pocket that interacts with the C3 side chain of cannabinoid agonists. Key hydrogen bonds would form between both K3.28(192) and E(258) and the A-ring OH, and between Q(261) and the C-ring C-12 hydroxypropyl. In summary, the present study represents one of the first attempts to construct a homology model of the CB(1) cannabinoid receptor based upon the published bovine rhodopsin x-ray crystal structure and to elucidate the putative ligand binding site for nonclassical cannabinoid agonists. We postulated sites of the CB(1) receptor critical for the ligand interaction, including the hydrophobic pocket interacting with the key pharmacophoric moiety, the C3 side chain. More work is needed to delineate between two alternative (and possibly other) binding conformations of the nonclassical cannabinoid ligands within the CB(1) receptor. The present study provides a consistent framework for further investigation of the CB(1) receptor-ligand interaction and for the study of CB(1) receptor activation.

Amino Acid Sequence↗

PubMed2004

Design, synthesis, and computational affinity prediction of ester soft drugs as inhibitors of dihydrofolate reductase from Pneumocystis carinii.

A series of dihydrofolate reductase (DHFR) inhibitors, where the methylenamino-bridge of non-classical inhibitors was replaced with an ester function, have been prepared as potential soft drugs intended for inhalation against Pneumocystis carinii pneumonia (PCP). Several of the new ester-based inhibitors that should serve as good substrates for the ubiquitous esterases and possibly constitute safer alternatives to metabolically stable DHFR inhibitors administered orally, were found to be potent inhibitors of P. carinii DHFR (pcDHFR). Although the objectives of the present program is to achieve a favorable toxicity profile by applying the soft drug concept, a high preference for inhibition of the fungal DHFR versus the mammalian DHFR is still desirable to suppress host toxicity at the site of administration. Compounds with a slight preference for the fungal enzyme were identified. The selection of the target compounds for synthesis was partly guided by an automated docking and scoring procedure as well as molecular dynamics simulations. The modest selectivity of the synthesized inhibitors was reasonably well predicted, although a correct ranking of the relative affinities was not successful in all cases.

Drug Design↗

PubMed2005

R-state hemoglobin bound to heterotropic effectors: models of the DPG, IHP and RSR13 binding sites.

We performed a docking study followed by a 500-ps molecular dynamics simulation of R-state human adult hemoglobin (HbA) complexed to different heterotropic effectors [2,3-diphosphoglycerate (DPG), inositol hexaphosphate (IHP), and 2-[4-[(3,5-dichlorophenylcarbamoyl)-]methyl]-phenoxy]-2-methylpropionic acid (RSR13)) to propose a molecular basis for recently reported interactions of effectors with oxygenated hemoglobin. The simulations were carried out with counterions and explicit solvation. As reported for T-state HbA, the effector binding sites are also located in the central cavity of the R-state and differ depending on effector anionic character. DPG and IHP bind between the alpha-subunits and the RSR13 site spans the alpha1-, alpha2- and beta2-subunits. The generated models provide the first report of the molecular details of R-state HbA bound to heterotropic effectors.

Aniline Compounds↗

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↗