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 919 records · Page 51Linked to original sources

PubMed2004

Predicting the cell differentiation activity of 1 alpha,25-dihydroxyvitamin D3 side chain analogues from docking simulations.

We present a receptor-based protocol for the prediction of the cell differentiation activities of a series of side chain analogues of 1 alpha,25-dihydroxyvitamin D(3), a compound that exhibits a very large variety of biological functions. Our protocol is able to reproduce the activity of the compounds studied here. It also sheds light on the relative importance of binding site residues in the biological activity and on the mechanism behind it.

Binding Sites↗

PubMed2001

X-ray structure of HPr kinase: a bacterial protein kinase with a P-loop nucleotide-binding domain.

HPr kinase/phosphatase (HprK/P) is a key regulatory enzyme controlling carbon metabolism in Gram- positive bacteria. It catalyses the ATP-dependent phosphorylation of Ser46 in HPr, a protein of the phosphotransferase system, and also its dephosphorylation. HprK/P is unrelated to eukaryotic protein kinases, but contains the Walker motif A characteristic of nucleotide-binding proteins. We report here the X-ray structure of an active fragment of Lactobacillus casei HprK/P at 2.8 A resolution, solved by the multiwavelength anomalous dispersion method on a seleniated protein (PDB code 1jb1). The protein is a hexamer, with each subunit containing an ATP-binding domain similar to nucleoside/nucleotide kinases, and a putative HPr-binding domain unrelated to the substrate-binding domains of other kinases. The Walker motif A forms a typical P-loop which binds inorganic phosphate in the crystal. We modelled ATP binding by comparison with adenylate kinase, and designed a tentative model of the complex with HPr based on a docking simulation. The results confirm that HprK/P represents a new family of protein kinases, first identified in bacteria, but which may also have members in eukaryotes.

Adenosine Diphosphate↗

PubMed2006

Molecular determination by electron microscopy of the actin filament end structure.

In eukaryotic cells, actin filaments play various crucial roles by altering their spatial and temporal distributions in the cell. The distribution of actin filaments is regulated by the binding of end-binding proteins, including capping protein (CapZ in muscle), the Arp2/3 complex, gelsolin, formin and tropomodulin, to the end of the actin filament. In order to determine the nature of these regulations, structural elucidations of actin filament-end-binding protein complexes are crucially important. Here, we have developed new procedures on the basis of single-particle analysis to determine the structure of the end of actin filaments from electron micrographs. In these procedures, the polarity of the actin filament image, as well as the azimuth orientation and the axial position of each actin protomer within a short stretch near the filament end, were determined accurately. This improved both the stability and accuracy of the structural determination dramatically. We tested our procedures by reconstructing structures from simulated filament images, which were obtained from 24 model structures for the actin-CapZ complex. These model structures were generated by random docking of the atomic structure of CapZ to the barbed end of an atomic model of the actin filament. Of the 24 model structures, 23 were recovered correctly by the present procedures. We found that our analysis was robust against local aberrations of the helical twist near the end of the actin filament. Finally, the procedures were applied successfully to determine the structure of the actin-CapZ complex from real cryo-electron micrographs of the complex. This is the first method for elucidating the detailed 3D structures at the end of the actin filament.

Actin Cytoskeleton↗

PubMed1996

Computer-aided molecular modeling of cathepsin E, a possible endothelin-converting enzyme.

A three-dimensional model of human cathepsin E, a possible endothelin-converting enzyme, is constructed using computer-aided molecular modeling techniques. The structure of porcine pepsin, another aspartic protease, was used as a template. The final structure, after all gaps and deletions were made, was optimized using the AMBER-4 package. A dipeptide (Trp-Val) representing the substrate was docked in the putative active site and the whole structure was optimized after several runs of minimization and dynamics calculations. The result of this modeling study showed that the structure of cathepsin E is similar to that of porcine pepsin and has three disulfide bonds that are conserved in both enzymes. There are two Asp-Thr-Gly sequences at the active site of enzyme. The active site cavity is large enough to accommodate its substrate.

Amino Acid Sequence↗

PubMed2000

Domain motions of EF-G bound to the 70S ribosome: insights from a hand-shaking between multi-resolution structures.

Molecular modeling and information processing techniques were combined to refine the structure of translocase (EF-G) in the ribosome-bound form against data from cryoelectron microscopy (cryo-EM). We devised a novel multi-scale refinement method based on vector quantization and force-field methods that gives excellent agreement between the flexibly docked structure of GDP. EF-G and the cryo-EM density map at 17 A resolution. The refinement reveals a dramatic "induced fit" conformational change on the 70S ribosome, mainly involving EF-G's domains III, IV, and V. The rearrangement of EF-G's structurally preserved regions, mediated and guided by flexible linkers, defines the site of interaction with the GTPase-associated center of the ribosome.

Computer Simulation↗

PubMed2004

Modeling correlated main-chain motions in proteins for flexible molecular recognition.

We describe a new method for modeling protein and ligand main-chain flexibility, and show its ability to model flexible molecular recognition. The goal is to sample the full conformational space, including large-scale motions that typically cannot be reached in molecular dynamics simulations due to the computational intensity, as well as conformations that have not been observed yet by crystallography or NMR. A secondary goal is to assess the degree of flexibility consistent with protein-ligand recognition. Flexibility analysis of the target protein is performed using the graph-theoretic algorithm FIRST, which also identifies coupled networks of covalent and noncovalent bonds within the protein. The available conformations of the flexible regions are then explored with ROCK by random-walk sampling of the rotatable bonds. ROCK explores correlated motions by only sampling dihedral angles that preserve the coupled bond networks in the protein and generates conformers with good stereochemistry, without using a computationally expensive potential function. A representative set of the conformational ensemble generated this way can be used as targets for docking with SLIDE, which handles the flexibility of protein and ligand side-chains. The realism of this protein main-chain conformational sampling is assessed by comparison with time-resolved NMR studies of cyclophilin A motions. ROCK is also effective for modeling the flexibility of large cyclic and polycyclic ligands, as demonstrated for cyclosporin and zearalenol. The use of this combined approach to perform docking with main-chain flexibility is illustrated for the cyclophilin A-cyclosporin complex and the estrogen receptor in complex with zearalenol, while addressing the question of how much flexibility is allowed without hindering molecular recognition.

Biochemical Phenomena↗

PubMedDate not supplied

Comparative molecular modeling analysis of-5-amidinoindole and benzamidine binding to thrombin and trypsin: specific H-bond formation contributes to high 5-amidinoindole potency and selectivity for thrombin and factor Xa.

The coagulation cascade enzymes thrombin and factor Xa are known to have specificity pockets very similar to those of trypsin and plasmin. However, comparative molecular modeling analysis of the crystal structures of benzamidine-thrombin and benzamidine-trypsin, in conjunction with a docking analysis of 5-amidinoindole and related inhibitors in both enzymes reveals subtle differences between the specificity sites of the two types of enzymes. Specifically, thrombin and factor Xa, which have an alanine residue at position 190, exhibit increased activities for the rigid and more bulky bicyclic derivatives of benzamidine (e.g. amidinobenzofuran, amidinothiophene and amidinoindole), because of additional hydrophobic and H-bond interactions between the inhibitors and the specificity sites, whereas enzymes with a serine residue at position 190, like trypsin and plasmin, exhibit little difference in activity among the same set of compounds because of the orientational restriction imposed on the inhibitors by Ser190, which forms an additional H-bond with the amidino group of the inhibitors. Enzymes of both groups show similar responses to the flexible aminobenzamidine since the smaller size and the rotatable anilino group of the inhibitor would allow the inhibitor to achieve favorable electrostatic interactions with both groups of enzymes. 5-amidinoindole is the most dramatic example of the rigid bicyclic type inhibitor. Based on our docking analysis, we propose that a selective H-bond with the hydroxyl group of the catalytic Ser195 and the subtle differences in steric fit imposed by Ala/Ser at position 190 explain the high potency and selectivity of 5-amidinoindole for thrombin and factor Xa.

Amino Acid Substitution↗

PubMed2006

Selective pharmacophore design for alpha1-adrenoceptor subtypes.

Alpha1-adrenoceptors are G-protein coupled receptors found in a variety of vascular tissues and responsible for vasoconstriction. Selectivity for each of the three subtypes is an important consideration in drug design in order to minimise the possibility of side effects. Using Catalyst we developed ligand-based pharmacophores from alpha(1a,b,d)-selective antagonists available in the literature using three separate training sets. Four-feature pharmacophores were developed for the alpha(1a) and alpha(1b) subtype-selective antagonists and a five-feature pharmacophore was developed for the alpha(1d) subtype-selective antagonists. The alpha(1a) pharmacophore represents both class I and II compounds with good predictivity for other compounds outside the training set as well. The alpha(1b) pharmacophore best predicts the activity of prazosin analogues as these make up the majority of alpha(1b)-selective antagonists. Unexpectedly, no positive ionisable feature was incorporated in the alpha(1b) pharmacophore. The alpha(1d) pharmacophore was based primarily on one structural class of compounds, but has good predictivity for a heterogeneous test set. Preliminary docking studies using AutoDock and optimised alpha1-adrenoceptor homology models, conducted with the antagonists prazosin (32) and 66, showed good agreement with the findings from the pharmacophores.

Adrenergic alpha-1 Receptor Antagonists↗

PubMed2004

Structural mining: self-consistent design on flexible protein-peptide docking and transferable binding affinity potential.

A flexible protein-peptide docking method has been designed to consider not only ligand flexibility but also the flexibility of the protein. The method is based on a Monte Carlo annealing process. Simulations with a distance root-mean-square (dRMS) virtual energy function revealed that the flexibility of protein side chains was as important as ligand flexibility for successful protein-peptide docking. On the basis of mean field theory, a transferable potential was designed to evaluate distance-dependent protein-ligand interactions and atomic solvation energies. The potential parameters were developed using a self-consistent process based on only 10 known complex structures. The effectiveness of each intermediate potential was judged on the basis of a Z score, approximating the gap between the energy of the native complex and the average energy of a decoy set. The Z score was determined using experimentally determined native structures and decoys generated by docking with the intermediate potentials. Using 6600 generated decoys and the Z score optimization criterion proposed in this work, the developed potential yielded an acceptable correlation of R(2) = 0.77, with binding free energies determined for known MHC I complexes (Class I Major Histocompatibility protein HLA-A(*)0201) which were not present in the training set. Test docking on 25 complexes further revealed a significant correlation between energy and dRMS, important for identifying native-like conformations. The near-native structures always belonged to one of the conformational classes with lower predicted binding energy. The lowest energy docked conformations are generally associated with near-native conformations, less than 3.0 Angstrom dRMS (and in many cases less than 1.0 Angstrom) from the experimentally determined structures.

Algorithms↗

PubMed2000

Bivalent sequential binding model of a Bacillus thuringiensis toxin to gypsy moth aminopeptidase N receptor.

Specificity for target insects of Bacillus thuringiensis insecticidal Cry toxins is largely determined by toxin affinity for insect midgut receptors. The mode of binding for one such toxin-receptor complex was investigated by extensive toxin mutagenesis, followed by real-time receptor binding analysis using an optical biosensor (BIAcore). Wild-type Cry1Ac, a three-domain, lepidopteran-specific toxin, bound purified gypsy moth (Lymantria dispar) aminopeptidase N (APN) biphasically. Site 1 displayed fast association and dissociation kinetics, while site 2 possessed slower kinetics, yet tighter affinity. We empirically determined that two Cry1Ac surface regions are involved in in vivo toxicity and APN binding. Mutations within domain III affected binding rates to APN site 1, whereas mutations in domain II affected binding rates to APN site 2. Furthermore, domain III contact is completely inhibited in the presence of N-acetylgalactosamine, indicating loss of domain III binding eliminates all APN binding. Based upon these observations, the following model is proposed. A cavity in lectin-like domain III initiates docking through recognition of an N-acetylgalactosamine moiety on L. dispar APN. Following primary docking, a higher affinity domain II binding mechanism occurs, which is critical for insecticidal activity.

Animals↗

PubMed1999

Homology model of juvenile hormone esterase from the crop pest, Heliothis virescens.

Juvenile Hormone Esterase (JHE) plays an essential role in the development of insects since it is partially responsible for clearing juvenile hormone (JH), one of the hormones that is responsible for insect metamorphosis. JHE is a 60 kDa enzyme that selectively hydrolyzes the alpha/beta unsaturated ester of JH. Because of its pivotal role in insect development, we have targeted JHE for use as a biopesticide. In this study, we have constructed a homology-based molecular model of JHE from the agricultural crop pest, Heliothis virescens. JHE is a member of the alpha/beta hydrolase fold family of enzymes and was built according to two structures in the same family: acetylcholinesterase from Torpedo californica and lipase from Geotrichum candidum. Analysis of the active site region reveals extensive conservation between JHE and its templates. A surprise was the presence of a conserved Ser near the catalytic triad. Docking of JH III into the active site has provided insight into protein-substrate interactions that are corroborated by experimental observation. The model is being used as a predictive basis to design biopesticides. In this regard, we have identified a site on the protein surface that is suggestive of a recognition site for the putative JHE receptor.

Acetylcholinesterase↗

PubMed1999

Homology modeling of gelatinase catalytic domains and docking simulations of novel sulfonamide inhibitors.

Three-dimensional models for the catalytic domain of gelatinases (MMP-9 and -2) have been constructed based on the X-ray crystal structure of MMP-3. Conformations of the loop segment which forms the bottom half of the S1' subsite but shows conformational diversity among the crystal structures of other MMPs have been explored by simulated annealing of each gelatinase model complexed with two highly potent "probe" inhibitors. Representative catalytic domain models have been selected for each gelatinase from the set of generated conformations based on shape complementarity of the loop to the probe inhibitors. The single model selected for MMP-9 was utilized to explain the structure-activity relationship of our novel sulfonamide inhibitors. Molecular dynamics (MD) simulations of the complex models revealed important features of the binding mechanism of our inhibitors: (i) the ligand carboxylate group coordinating to the catalytic zinc ion and hydrogen bonding to the Glu219 side chain, (ii) one of the sulfonyl oxygens forming hydrogen bonds with the main chain NHs (Leu181 and Ala182), (iii) the sulfonyl substituent making extensive hydrophobic contact with the S1' subsite. The gauche conformation exclusively adopted by the sulfonamide C-N-S-C torsion plays an important role in achieving the third binding feature by properly directing the substituent into the S1' subsite. Improvement of the inhibitory activity according to straight elongation of the sulfonyl substituent was attributed to an increase of the hydrophobic contact between the substituent and the S1' subsite. Structural modifications which alter the straight shape of the substituent lead to deterioration of the activity. On the other hand, the two candidate models selected for MMP-2 differ in the bottom shape of the S1' subsite: one with a channel-like subsite and the other with a pocket-like subsite resembling that of the MMP-9 model. The bottom shape was experimentally probed by chemical synthesis of inhibitors having elongated sulfonyl substituents whose terminal alkyl groups were shown by MD simulations to protrude from the S1' subsite bottom into the solvent. Gelatinase assays of these inhibitors showed that elongation of the substituent significantly reduces activity against MMP-9 while retaining activity against MMP-2, consequently increasing the selectivity between MMP-2 and -9. The results confirm that MMP-9 has a pocket-like S1' subsite with a floorboard and MMP-2 has a channel-like S1' subsite.

Amino Acid Sequence↗

PubMed2005

The molecular basis of resilience to the effect of the Lys103Asn mutation in non-nucleoside HIV-1 reverse transcriptase inhibitors studied by targeted molecular dynamics simulations.

A series of targeted molecular dynamics simulations have been carried out in an attempt to assess the effect that the common Lys103Asn mutation in HIV-1 reverse transcriptase (RT) has on the binding of three representative non-nucleoside RT inhibitors (NNRTI), nevirapine, efavirenz, and etravirine. We have shown previously that, in the absence of an incoming inhibitor, creation of the NNRTI binding pocket is hampered due to the existence of a hydrogen bond between the side chains of Asn103 and Tyr188 for which no equivalent exists in the wild-type enzyme. As an extension of this work, we now apply the same methodology to drive the enzyme's conformation from the unbound state to the drug-bound state in the presence of the NNRTI. The location of each drug outside the binding pocket was determined by an automated docking program, and steering into the binding pocket followed a route that is likely to represent the actual entrance pathway. The additional hurdle to inhibitor entry imposed by the extra Asn103-Tyr188 hydrogen bond is seen to affect each NNRTI differently, with the ability to disrupt this interaction increasing in the order etravirine >> efavirenz > or = nevirapine, in good accord with the experimental findings. This coherent picture strongly suggests that attempts to overcome resistance through structure-based drug design may be considerably more successful if dynamic structural aspects of the type studied here are considered, particularly in cases where binding energy-based structure-activity relationship methods are unable to provide the required information.

Alkynes↗

PubMed1995

Modeled structure of the 75-kDa neurotrophin receptor.

Motifs in ligand-binding domains of the neurotrophin (NTR) and lymphotoxin (TNFR-I) receptors define a family of receptors that mediates programmed cell death. We have explored relationships of architecture and function in this family through a molecular model of NTR, also called p75NGFR or LANR. Modeling by homology took advantage of four modular subdomains in the crystal structure of TNFR-I that also occur in NTR. Hypothetical complexes between the model and a ligand structure (for nerve growth factor, NGF) were then examined using docking software. NTR appears to bind in the dimer interface of NGF, making two sets of contacts. NTR subdomains III and IV provide the ligand-contact surfaces, in contrast to TNFR, in which subdomains II and III contact TNF-beta. NTR subdomain II appears to have been evolutionarily modified, potentially contributing to an interface between receptor subunits. These and other specific predictions of the model will require experimental confirmation.

Amino Acid Sequence↗

PubMed2003

In silico mutations and molecular dynamics studies on a winged bean chymotrypsin inhibitor protein.

Winged bean chymotrypsin inhibitor (WCI) has an intruding residue Asn14 that plays a crucial role in stabilizing the reactive site loop conformation. This residue is found to be conserved in the Kunitz (STI) family of serine protease inhibitors. To understand the contribution of this scaffolding residue on the stability of the reactive site loop, it was mutated in silico to Gly, Ala, Ser, Thr, Leu and Val and molecular dynamics (MD) simulations were carried out on the mutants. The results of MD simulations reveal the conformational variability and range of motions possible for the reactive site loop of different mutants. The N-terminus side of the scissile bond, which is close to a beta-barrel, is conformationally less variable, while the C-terminus side, which is relatively far from any such secondary structural element, is more variable and needs stability through hydrogen-bonding interactions. The simulated structures of WCI and the mutants were docked in the peptide-binding groove of the cognate enzyme chymotrypsin and the ability to form standard hydrogen-bonding interactions at P3, P1 and P2' residues were compared. The results of the MD simulations coupled with docking studies indicate that hydrophobic residues like Leu and Val at the 14th position are disruptive for the integrity of the reactive site loop, whereas a residue like Thr, which can stabilize the C-terminus side of the scissile bond, can be predicted at this position. However, the size and charge of the Asn residue made it most suitable for the best maintenance of the integrity of the reactive site loop, explaining its conserved nature in the family.

Binding Sites↗

PubMed1993

Molecular modelling of the specific interactions involved in the amylose complexation by fatty acids.

Comprehensive modelling of a fatty acid molecule inside a VH amylose helix is described. In a first step, the docking of an acetic acid molecule near the helix entry was performed. The low energy solutions were propagated by an iterative procedure involving the sequential addition of single CH2 groups up to a C12 fatty acid followed by energy minimizations. The main result is the superposition of the aliphatic and the helix axes. For the low-energy complexes, the mean plane of the aliphatic carbons has three potential orientations. In each, the aliphatic hydrogens point towards the less crowded regions near the glycosidic oxygens of the amylose. The close packing is due to the related symmetries of both the helix and aliphatic chain. In a second step, the relative roles of the aliphatic part and the polar group were studied separately. For the aliphatic chain, a map based on the two major internal parameters (translation and rotation) along the helix axis shows that the isolated docking solutions are related by a combination of a 60 degrees (360 degrees/6) rotation and a translation of p/6 (p = 0.804 nm corresponds to the pitch of Vhydrate amylose). The H5 glucopyranose atoms participate in close contacts and are responsible for steric conflicts in structures intermediate to the stable docking solutions. The four possible low-energy arrangements of the carboxylic group were added to the calculated amylose/aliphatic structures. Two stable conformations of the total fatty acid molecule were determined. For both stable solutions, the polar group is located near the entrance of the helix cavity.(ABSTRACT TRUNCATED AT 250 WORDS)

Amylose↗

PubMed2003

Catalysis, specificity, and ACP docking site of Streptomyces coelicolor malonyl-CoA:ACP transacylase.

Malonyl-CoA:ACP transacylase (MAT), the fabD gene product of Streptomyces coelicolor A3(2), participates in both fatty acid and polyketide synthesis pathways, transferring malonyl groups that are used as extender units in chain growth from malonyl-CoA to pathway-specific acyl carrier proteins (ACPs). Here, the 2.0 A structure reveals an invariant arginine bound to an acetate that mimics the malonyl carboxylate and helps define the extender unit binding site. Catalysis may only occur when the oxyanion hole is formed through substrate binding, preventing hydrolysis of the acyl-enzyme intermediate. Macromolecular docking simulations with actinorhodin ACP suggest that the majority of the ACP docking surface is formed by a helical flap. These results should help to engineer polyketide synthases (PKSs) that produce novel polyketides.

Amino Acid Sequence↗

PubMedDate not supplied

3D-QSAR and docking studies of selective GSK-3beta inhibitors. Comparison with a thieno[2,3-b]pyrrolizinone derivative, a new potential lead for GSK-3beta ligands.

The three-dimensional structures of 3-anilino-4-arylmaleimides, selective GSK-3beta inhibitors, were correlated to their biological affinities by 3D-QSAR studies (CoMFA method). The cocrystallographic data of GSK-3beta vs 3-anilino-4-arylmaleimide allowed us to compare 3D-QSAR results to experimental intermolecular interactions. The results of the CoMFA analysis did not really correspond to the interactions recorded in the active site, but they characterized fundamental features (areas of the active site) of the interactions ligand-receptor. These studies were the starting point to analyze a new GSK-3beta ligand, a thieno[2,3-b]pyrrolizinone derivative. This comparison based on docking and simulation approaches allowed us to confirm one preferential orientation of this ligand inside the active site, explaining the relationship with the reference 3-anilino-4-arylmaleimide derivatives and its biological affinity.

Enzyme Inhibitors↗