SEARCH · PubMed Health

Results for “Molecular docking and molecular dynamics 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 91 records · Page 5Linked to original sources

PubMed2002

Novel computer program for fast exact calculation of accessible and molecular surface areas and average surface curvature.

New computer programs, SurfRace and FastSurf, perform fast calculations of the solvent accessible and molecular (solvent excluded) surface areas of macromolecules. Program SurfRace also calculates the areas of cavities inaccessible from the outside. We introduce the definition of average curvature of molecular surface and calculate average molecular surface curvatures for each atom in a structure. All surface area and curvature calculations are analytic and therefore yield exact values of these quantities. High calculation speed of this software is achieved primarily by avoiding computationally expensive mathematical procedures wherever possible and by efficient handling of surface data structures. The programs are written initially in the language C for PCs running Windows 2000/98/NT, but their code is portable to other platforms with only minor changes in input-output procedures. The algorithm is robust and does not ignore either multiplicity or degeneracy of atomic overlaps. Fast, memory-efficient and robust execution make this software attractive for applications both in computationally expensive energy minimization algorithms, such as docking or molecular dynamics simulations, and in stand-alone surface area and curvature calculations.

Algorithms↗

PubMed2004

Novel approaches for modeling of the A1 adenosine receptor and its agonist binding site.

The present work describes the building of a human A(1) adenosine receptor (hA(1)AR) model, based on the X-ray crystal structure of bovine rhodopsin, and its use as a basis for the investigation of some important structural characteristics of the receptor. One of the issues investigated was the protonation position of two histidine residues known to influence ligand binding, with protonation of His251 (6.52) in epsilon position and His278 (7.43) in delta position showing the best agreement with experimental evidence. The model was also used to study the position and structural role of water molecules present in the helical bundle. Finally, the binding site location and the ligand docking were investigated using an objective strategy. A suitable site for the binding of the ribose moiety of adenosine was first postulated and further confirmed by means of a novel chemometric strategy based on GRIND descriptors. Using this position as an anchor point, the binding of adenosine was studied by docking and molecular dynamics simulations obtaining two putative binding positions in good agreement with experimental data.

Adenosine↗

PubMed2026

Cross-tissue Mendelian randomization prioritizes RAB27B as a brain-derived candidate protein for postpartum depression.

OBJECTIVE: Postpartum depression (PPD) is one of the most common and debilitating complications of childbirth, yet the candidate proteins linking genetic risk to disease remain poorly defined. Building on recent genome-wide association studies (GWAS), we sought to integrate cross-tissue proteogenomic data to identify candidate proteins for PPD and explore therapeutic opportunities. METHODS: We conducted two-sample Mendelian randomization (MR) using genome-wide significant cis-protein QTLs from brain (n = 608 proteins), cerebrospinal fluid (CSF; n = 214), and plasma (n = 612). PPD summary statistics were obtained from FinnGen R8 (13,657 cases, 236,178 controls) and replicated in an independent GWAS. Phenome-wide association (PheWAS) was used to assess pleiotropy. Potential therapeutic targets were evaluated through DSigDB drug repurposing, molecular docking, and molecular dynamics simulations. RESULTS: Among all proteins tested, RAB27B was the only brain-derived protein surpassing Bonferroni correction (OR = 1.60; 95% CI: 1.30-1.96; P = 6.6 × 10⁻⁶), whereas no significant proteins were identified in CSF or plasma. This association was replicated in an independent GWAS (OR = 1.27; 95% CI: 1.02-1.58; P = 0.037). PheWAS identified no pleiotropic associations at genome-wide significance. In silico drug repurposing identified pregnenolone as a candidate ligand with computationally predicted stable binding to RAB27B, providing a starting point for future experimental validation. CONCLUSION: This study provides the first cross-tissue proteogenomic evidence that RAB27B is a brain-derived, reproducible candidate protein genetically associated with PPD. By extending GWAS signals to functional protein-level mechanisms and therapeutic inference, our findings nominate RAB27B and pregnenolone as promising directions for postpartum psychiatric research.

Humans↗

PubMed2003

Structural and functional characterization of the human CCR5 receptor in complex with HIV gp120 envelope glycoprotein and CD4 receptor by molecular modeling studies.

The entry of human immunodeficiency virus (HIV) into cells depends on a sequential interaction of the gp120 envelope glycoprotein with the cellular receptors CD4 and members of the chemokine receptor family. The CC chemokine receptor CCR5 is such a receptor for several chemokines and a major coreceptor for the entry of R5 HIV type-1 (HIV-1) into cells. Although many studies focus on the interaction of CCR5 with HIV-1, the corresponding interaction sites in CCR5 and gp120 have not been matched. Here we used an approach combining protein structure modeling, docking and molecular dynamics simulation to build a series of structural models of the CCR5 in complexes with gp120 and CD4. Interactions such as hydrogen bonds, salt bridges and van der Waals contacts between CCR5 and gp120 were investigated. Three snapshots of CCR5-gp120-CD4 models revealed that the initial interactions of CCR5 with gp120 are involved in the negatively charged N-terminus (Nt) region of CCR5 and positively charged bridging sheet region of gp120. Further interactions occurred between extracellular loop2 (ECL2) of CCR5 and the base of V3 loop regions of gp120. These interactions may induce the conformational changes in gp120 and lead to the final entry of HIV into the cell. These results not only strongly support the two-step gp120-CCR5 binding mechanism, but also rationalize extensive biological data about the role of CCR5 in HIV-1 gp120 binding and entry, and may guide efforts to design novel inhibitors.

CD4 Antigens↗

PubMed2005

Computational studies and drug design for HIV-1 reverse transcriptase inhibitors of 3',4'-di-O-(S)-camphanoyl-(+)-cis-khellactone (DCK) analogs.

Molecular docking and molecular dynamics simulation were applied to study the binding mode of 3',4'-di-O-(S)-camphanoyl-(+)-cis-khellactone (DCK) analogs anti-HIV inhibitors with HIV-1 RT. The results suggest that there is a strong hydrogen bond between DCK O16 and NH of Lys101, and that DCK analogues might act similarly as other types of HIV-1 RT inhibitors. The investigation about drug resistance for DCK shows no remarkable influence on the most frequently observed mutation K103N of HIV-1 RT. Based on the proposed mechanism, some new structures were designed and predicted by a SVM model. All compounds exhibited potent inhibitory activities against HIV replication in H9 lymphocytes with EC50 values lower than 1.95 microM. The rationality of the method was validated by experimental results.

Camphor↗

PubMed2005

Modelling the interaction of catecholamines with the alpha 1A adrenoceptor towards a ligand-induced receptor structure.

Adrenoceptors are members of the important G protein coupled receptor family for which the detailed mechanism of activation remains unclear. In this study, we have combined docking and molecular dynamics simulations to model the ligand induced effect on an homology derived human alpha1A adrenoceptor. Analysis of agonist/alpha1A adrenoceptor complex interactions focused on the role of the charged amine group, the aromatic ring, the N-methyl group of adrenaline, the beta hydroxyl group and the catechol meta and para hydroxyl groups of the catecholamines. The most critical interactions for the binding of the agonists are consistent with many earlier reports and our study suggests new residues possibly involved in the agonist-binding site, namely Thr-174 and Cys-176. We further observe a number of structural changes that occur upon agonist binding including a movement of TM-V away from TM-III and a change in the interactions of Asp-123 of the conserved DRY motif. This may cause Arg-124 to move out of the TM helical bundle and change the orientation of residues in IC-II and IC-III, allowing for increased affinity of coupling to the G-protein.

Adrenergic alpha-1 Receptor Agonists↗

PubMed2005

Identification of a drug target motif: an anti-tumor drug NK109 interacts with a PNxxxxP.

The synthetic compound NK109 shows anti-tumor effects against a number of human cancer cell lines. The mechanism of action is thought to involve the inhibition of DNA topoisomerase II. However, NK109 also exhibits potent anti-tumor activities against doxorubicin-, cisplatin- and etoposide-resistant human cell lines. This paper describes target validation of NK109 using biotinylated NK109 and a T7 phage library screening procedure. Phage particles displaying an affinity for NK109 were isolated and the DNA sequence determined. The amino acid sequences of selected peptides, and the results of mutation experiments by alanine scanning, confirmed that the binding target motif of NK109 is PNxxxxP. In silico analysis of the interaction between NK109 and the peptide, by docking simulation and molecular dynamics, supported this conclusion. The PNxxxxP motif exists in the C2 domain of protein kinase Calpha. NK109 was confirmed to bind the C2 domain from surface plasmon resonance analysis. Furthermore, NK109 moderately inhibited protein kinase C activity in vitro. Our results show that the anti-tumor activity of NK109 stems from interactions with multiple protein targets.

Amino Acid Motifs↗

PubMed2006

Stereoselectivity of 8-OH-DPAT toward the serotonin 5-HT1A receptor: biochemical and molecular modeling study.

The great majority of pharmacological investigations of 5-HT1A receptors' reactivity has been performed using racemic 8-OH-DPAT, therefore the biochemical as well as behavioral profiles of both 8-OH-DPAT enantiomers are not circumstantiated. In the biochemical study capability of racemic 8-OH-DPAT (0.05, 0.1 mg/kg s.c.) and its counterparts R-8-OH-DPAT (0.05, 0.1 mg/kg s.c.) and S-8-OH-DPAT (0.05, 0.1 mg/kg s.c.) to influence 5-HT synthesis rate in rats' prefrontal cortex, hypothalamus, hippocampus and brainstem was evaluated by HPLC/ED technique. Biochemical results are supported by the exhaustive computational study of possible differences between R- and S-enantiomer toward the 5-HT1A receptor. A reliable 3D model of the rat 5-HT1A receptor was constructed from the amino acid sequence using the crystal structure of bovine rhodopsin as a structural template. The structure of the receptor model was validated through docking studies and molecular dynamics simulations that gave results consistent with experimental data. Docking studies and the dynamics of ligand-receptor complexes emphasized different profiles of both enantiomers at the molecular level. The results of both biochemical and computational studies confirmed that R-enantiomer in contrast to S-8-OH-DPAT acts as full and potent agonist, whilst racemic form may display similar pharmacological profile to R-8-OH-DPAT.

8-Hydroxy-2-(di-n-propylamino)tetralin↗

PubMed2005

Macrocyclic inhibitors of the malarial aspartic proteases plasmepsin I, II, and IV.

The first macrocyclic inhibitor of the Plasmodium falciparum aspartic proteases plasmepsin I, II, and IV with considerable selectivity over the human aspartic protease cathepsin D has been identified. A series of macrocyclic compounds were designed and synthesized. Cyclizations were accomplished using ring-closing metathesis with the second generation Grubbs catalyst. These compounds contain either a 13-membered or a 16-membered macrocycle and incorporate a 1,2-dihydroxyethylene as transition state mimicking unit. The binding mode of this new class of compounds was predicted with automated docking and molecular dynamics simulations, with an estimation of the binding affinities through the linear interaction energy (LIE) method.

Animals↗

PubMed2005

Structural basis for the GSK-3beta binding affinity and selectivity against CDK-2 of 1-(4-aminofurazan-3yl)-5-dialkylaminomethyl-1H-[1,2,3] triazole-4-carboxylic acid derivatives.

A novel structural class of glycogen synthase kinase-3beta inhibitors is modeled using quantum mechanics, automated docking, and molecular dynamics simulations. The proposed binding modes identify important hydrogen bonds and salt-bridges with the ATP-binding pocket of the kinase. The modeled complexes justify the observed structure-activity relationships and provide a structural basis for the high selectivity of these inhibitors against cyclin dependent kinase-2.

Adenosine Triphosphate↗

PubMed2005

Structural insight into the interactions of SoxV, SoxW and SoxS in the process of transport of reductants during sulfur oxidation by the novel global sulfur oxidation reaction cycle.

Microbial redox reactions involving inorganic sulfur compounds, mainly the sulfur anions, are one of the vital reactions responsible for the environmental sulfur balance. These reactions are mediated by phylogenetically diverse prokaryotes, some of which also take part in the extraction of metal ions from their sulfur containing ores. These sulfur oxidizers oxidize inorganic sulfur compounds like sulfide, thiosulfate etc. to produce reductants that are used for carbon dioxide fixation or in respiratory electron transfer chains. The sulfur oxidizing gene cluster (sox) of alpha-Proteobacteria comprises of at least 15 genes, forming two transcriptional units, viz., soxSR and soxVWXYZABCDEFGH. SoxV is known to be a CcdA homolog involved in the transport of reductants from cytoplasm to periplasm. SoxW and SoxS are periplasmic thioredoxins, which (SoxW) interact with SoxV and thereby help in the redox reactions. We have employed homology modeling to construct the three-dimensional structures of the SoxV, SoxW and SoxS proteins from Rhodovulum sulfidophilum. With the help of docking and molecular dynamics simulations we have identified the amino acid residues of these proteins involved in the interaction. The probable biochemical mechanism of the transport of reductants through the interactions of these proteins has also been investigated. Our study provides a rational basis to interpret the molecular mechanism of the biochemistry of sulfur anion oxidation reactions by these ecologically important organisms.

Alphaproteobacteria↗

PubMed2003

Design, synthesis and biological evaluation of ambenonium derivatives as AChE inhibitors.

Ambenonium (1), an old AChE inhibitor, is endowed with an outstanding affinity and a peculiar mechanism of action that, taken together, make it a very promising pharmacological tool for the treatment of Alzheimer's disease (AD). Unfortunately, the bisquaternary structure of 1 prevents its passage through the blood brain barrier. In a search of centrally active ambenonium derivatives, we planned to synthesize tertiary amines of 1, such as 2 and 3. In addition, to add new insights into the binding mechanism of the inhibitor, we designed constrained analogues of ambenonium by incorporating the diamine functions into cyclic moieties (4-12). The biological evaluation of the new compounds has been assessed in vitro against human AChE and BChE. All tertiary amine derivatives resulted more than 1000-fold less potent than 1 and, unlike prototype, did not show any selectivity between the two enzymes. This result, because of recent findings concerning the role of BChE in AD, makes our compounds, endowed with a well-balanced profile of AChE/BChE inhibition, valuable candidates for further development. To better clarify the interactions that account for the high affinity of 1, docking simulations and molecular dynamics studies on the AChE-1 complex were also carried out.

Acetylcholinesterase↗

PubMed2002

Nitrophenyl derivatives as aldose reductase inhibitors.

Nitrophenyl derivatives were recently discovered as a new class of ALR2 inhibitors by means of docking and database screening of the National Cancer Institute database of organic molecules. The nitro group was predicted to bind to the Tyr48 and His110 active site residues of the enzyme, the site where acidic ALR2 inhibitors such as carboxylic acids bind in their anionic form. Given the novelty of these compounds, we decided to expand their structure-activity relationships by synthesizing and testing a series of derivatives and the corresponding compounds having a carboxylic group instead of the nitro moiety; the results obtained were rationalized by means of docking and molecular dynamics simulations. On the whole there is an agreement between inhibitory data and the results of molecular modeling experiments, supporting the hypothesized binding mode of these compounds.

Aldehyde Reductase↗

PubMed2004

Localization of the antimitotic peptide and depsipeptide binding site on beta-tubulin.

Several naturally occurring peptides and depsipeptides which include the cryptophycins, dolastatin 10, hemiasterlin, and phomopsin A have been found to be potent antimitotic agents, causing cell death at picomolar or low nanomolar concentrations. These compounds inhibit microtubule growth, modulate the dynamics of microtubules, and induce the self-association of tubulin dimers into single-walled rings and spirals. These peptides exhibit mutual competitive inhibition in binding to beta-tubulin, while noncompetitively inhibiting the binding of vinblastine and vincristine to beta-tubulin. Despite the abundance of biochemical information, the details of their molecular interactions with tubulin are not known. In this study, using a combination of molecular dynamics simulations and molecular docking studies, a common binding site for cryptophycin 1, cryptophycin 52, dolastatin 10, hemiasterlin, and phomopsin A on beta-tubulin has been identified. Application of these same methods to alpha-tubulin indicated no interaction between alpha-tubulin and any of the peptides. On the basis of the docking results, a model for the mechanism of microtubule disruption and formation of aberrant nonmicrotubule structures is proposed. Both the active site and mechanism of microtubule depolymerization predictions are in good agreement with experimental findings.

Amino Acid Sequence↗

PubMedDate not supplied

Modeling the structural basis of human CCR5 chemokine receptor function: from homology model building and molecular dynamics validation to agonist and antagonist docking.

This article describes the construction and validation of a three-dimensional model of the human CCR5 receptor using a homology-based approach starting from the X-ray structure of the bovine rhodopsin receptor. The reliability of the model is assessed through molecular dynamics and docking simulations using both natural agonists and a synthetic antagonist. Some important structural and functional features of the receptor cavity and the extracellular loops are identified, in agreement with data available from site-directed mutagenesis. The results of this study help to explain the structural basis for the recognition, activation, and inhibition processes of CCR5 and may provide fresh insights for the design of HIV-1 entry blockers.

CCR5 Receptor Antagonists↗

PubMed2002

Molecular docking studies of a group of hydroxamate inhibitors with gelatinase-A by molecular dynamics.

We have performed docking and molecular dynamics simulations of hydroxamates complexed with human gelatinase-A (MMP-2) to gain insight into the structural and energetic preferences of these inhibitors. The study was conducted on a selected set of eleven compounds with variation in structure and activity. Molecular dynamics simulations were performed at 300 K for 100 ps with equilibration for 50 ps. The structural analyses of the trajectories indicate that the coordinate bond interactions, the hydrogen bond interactions, the van der Waals interactions as well as the hydrophobic interactions between ligand and receptor are responsible simultaneously for the preference of inhibition and potency. The ligand hydroxamate group is coordinated to the catalytic zinc ion and form stable hydrogen bonds with the carbonyl oxygen of Gly 162. The P1' group makes extensive van der Waals and hydrophobic contacts with the nonpolar side chains of several residues in the S1' subsite, including Leu 197, Val 198, Leu 218 and Tyr 223. Moreover, four to eight hydrogen bonds between hydroxamates and MMP-2 are formed to stabilize the inhibitors in the active site. Compared with the P2' and P3' groups, the P1' groups of inhibitors are oriented regularly, which is produced by the restrain of the S1' subsite. From the relationship between the length of the nonpolar P1' group and the biological activity, we confirm that MMP-2 has a pocket-like S1' subsite, not a channel-like S1' subsite proposed by Kiyama (Kiyama, R. et al., J. Med. Chem. 42 (1999), 1723). The energetic analyses show that the experimental binding free energies can be well correlated with the interactions between the inhibitors and their environments, which could be used as a simple score function to evaluate the binding affinities for other similar hydroxamates. The validity of the force field parameters and the MD simulations can be fully testified by the satisfactory agreements between the experimental structure-activity relationship and the information from the structural and energetic analyses. The information generated from the predicted complexes should be useful for further work in the area of structure-based design of new compounds.

Binding Sites↗

PubMed2000

Novel carbohydrate binding site recognizing blood group A and B determinants in a hybrid of cholera toxin and Escherichia coli heat-labile enterotoxin B-subunits.

The B-subunits of cholera toxin (CTB) and Escherichia coli heat-labile enterotoxin (LTB) are structurally and functionally related. However, the carbohydrate binding specificities of the two proteins differ. While both CTB and LTB bind to the GM1 ganglioside, LTB also binds to N-acetyllactosamine-terminated glycoconjugates. The structural basis of the differences in carbohydrate recognition has been investigated by a systematic exchange of amino acids between LTB and CTB. Thereby, a CTB/LTB hybrid with a gain-of-function mutation resulting in recognition of blood group A and B determinants was obtained. Glycosphingolipid binding assays showed a specific binding of this hybrid B-subunit, but not CTB or LTB, to slowly migrating non-acid glycosphingolipids of human and animal small intestinal epithelium. A binding-active glycosphingolipid isolated from cat intestinal epithelium was characterized by mass spectrometry and proton NMR as GalNAcalpha3(Fucalpha2)Galbeta4(Fucalpha3)Glc NAcbeta3Galbeta4Glc NAcbeta3Galbeta4Glcbeta1Cer. Comparison with reference glycosphingolipids showed that the minimum binding epitope recognized by the CTB/LTB hybrid was Galalpha3(Fucalpha2)Galbeta4(Fucalpha3)GlcNAc beta. The blood group A and B determinants bind to a novel carbohydrate binding site located at the top of the B-subunit interfaces, distinct from the GM1 binding site, as found by docking and molecular dynamics simulations.

ABO Blood-Group System↗

PubMed2003

Characterization of endostatin binding to heparin and heparan sulfate by surface plasmon resonance and molecular modeling: role of divalent cations.

Endostatin (20 kDa) is a C-terminal proteolytic fragment of collagen XVIII that is localized in vascular basement membrane zones in various organs. It binds zinc, heparin/heparan sulfate, laminin, and sulfatides and inhibits angiogenesis and tumor growth. Here we determined the kinetics and affinity of the interaction of endostatin with heparin/heparan sulfate and investigated the effects of divalent cations on these interactions and on the biological activities of endostatin. The binding of human recombinant endostatin to heparin and heparan sulfate was studied by surface plasmon resonance using BIAcore technology and further characterized by docking and molecular dynamics simulations. Kinetic data, evaluated using a 1:1 interaction model, showed that heparan sulfate bound to and dissociated from endostatin faster than heparin and that endostatin bound to heparin and heparan sulfate with a moderate affinity (K(D) approximately 2 microm). Molecular modeling of the complex between endostatin and heparin oligosaccharides predicted that, compared with mutagenesis studies, two further arginine residues, Arg(47) and Arg(66), participated in the binding. The binding of endostatin to heparin and heparan sulfate required the presence of divalent cations. The addition of ZnCl(2) to endostatin enhanced its binding to heparan sulfate by approximately 40% as well as its antiproliferative effect on endothelial cells stimulated by fibroblast growth factor-2, suggesting that this activity is mediated by the binding of endostatin to heparan sulfate. In contrast, no increase in the antiangiogenic and anti-proliferative activities of endostatin promoted by vascular endothelial growth factor was observed upon the addition of zinc.

Cations, Divalent↗