Dataset Information

Homology modeling, molecular docking, and molecular dynamics simulations elucidated ?-fetoprotein binding modes.

ABSTRACT: BACKGROUND: An important mechanism of endocrine activity is chemicals entering target cells via transport proteins and then interacting with hormone receptors such as the estrogen receptor (ER). ?-Fetoprotein (AFP) is a major transport protein in rodent serum that can bind and sequester estrogens, thus preventing entry to the target cell and where they could otherwise induce ER-mediated endocrine activity. Recently, we reported rat AFP binding affinities for a large set of structurally diverse chemicals, including 53 binders and 72 non-binders. However, the lack of three-dimensional (3D) structures of rat AFP hinders further understanding of the structural dependence for binding. Therefore, a 3D structure of rat AFP was built using homology modeling in order to elucidate rat AFP-ligand binding modes through docking analyses and molecular dynamics (MD) simulations. METHODS: Homology modeling was first applied to build a 3D structure of rat AFP. Molecular docking and Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) scoring were then used to examine potential rat AFP ligand binding modes. MD simulations and free energy calculations were performed to refine models of binding modes. RESULTS: A rat AFP tertiary structure was first obtained using homology modeling and MD simulations. The rat AFP-ligand binding modes of 13 structurally diverse, representative binders were calculated using molecular docking, (MM-GBSA) ranking and MD simulations. The key residues for rat AFP-ligand binding were postulated through analyzing the binding modes. CONCLUSION: The optimized 3D rat AFP structure and associated ligand binding modes shed light on rat AFP-ligand binding interactions that, in turn, provide a means to estimate binding affinity of unknown chemicals. Our results will assist in the evaluation of the endocrine disruption potential of chemicals.

SUBMITTER: Shen J

PROVIDER: S-EPMC3851483 | biostudies-literature | 2013

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Homology modeling, molecular docking, and molecular dynamics simulations elucidated α-fetoprotein binding modes.

Shen Jie J Zhang Wenqian W Fang Hong H Perkins Roger R Tong Weida W Hong Huixiao H

BMC bioinformatics 20131009

<h4>Background</h4>An important mechanism of endocrine activity is chemicals entering target cells via transport proteins and then interacting with hormone receptors such as the estrogen receptor (ER). α-Fetoprotein (AFP) is a major transport protein in rodent serum that can bind and sequester estrogens, thus preventing entry to the target cell and where they could otherwise induce ER-mediated endocrine activity. Recently, we reported rat AFP binding affinities for a large set of structurally di ...[more]

PMID: 24266910

Similar Datasets

Project description:We have modeled the structure of human lymphotactin (hLpnt), by homology modeling and molecular dynamics simulations. This chemokine is unique in having a single disulfide bond and a long C-terminal tail. Because other structural classes of chemokines have two pairs of Cys residues, compared to one in Lpnt, and because it has been shown that both disulfide bonds are required for stability and function, the question arises how the Lpnt maintains its structural integrity. The initial structure of hLpnt was constructed by homology modeling. The first 63 residues in the monomer of hLpnt were modeled using the structure of the human CC chemokine, RANTES, whose sequence appeared most similar. The structure of the long C-terminal tail, missing in RANTES, was taken from the human muscle fatty-acid binding protein. In a Protein Data Bank search, this protein was found to contain a sequence that was most homologous to the long tail. Consequently, the modeled hLpnt C-terminal tail consisted of both alpha-helical and beta-motifs. The complete model of the hLpnt monomer consisted of two alpha-helices located above the five-stranded beta-sheet. Molecular dynamics simulations of the solvated initial model have indicated that the stability of the predicted fold is related to the geometry of Pro78. The five-stranded beta-sheet appeared to be preserved only when Pro78 was modeled in the cis conformation. Simulations were also performed both for the C-terminal truncated forms of the hLpnt that contained one or two (CC chemokine-like) disulfide bonds, and for the chicken Lpnt (cLpnt). Our MD simulations indicated that the turn region (T30-G34) in hLpnt is important for the interactions with the receptor, and that the long C-terminal region stabilizes both the turn (T30-G34) and the five-stranded beta-sheet. The major conclusion from our theoretical studies is that the lack of one disulfide bond and the extension of the C-terminus in hLptn are mutually complementary. It is very likely that removal of two Cys residues sufficiently destabilizes the structure of a chemokine molecule, particularly the core beta-sheet, to abolish its biological function. However, this situation is rectified by the long C-terminal segment. The role of this long region is most likely to stabilize the first beta-turn region and alpha-helix H1, explaining how this chemokine can function with a single disulfide bond.

Project description:The transient receptor potential vanilloid type 1 (TRPV1) is a heat-activated cation channel protein, which contributes to inflammation, acute and persistent pain. Antagonists of human TRPV1 (hTRPV1) represent a novel therapeutic approach for the treatment of pain. Developing various antagonists of hTRPV1, however, has been hindered by the unavailability of a 3D structure of hTRPV1. Recently, the 3D structures of rat TRPV1 (rTRPV1) in the presence and absence of ligand have been reported as determined by cryo-EM. rTRPV1 shares 85.7% sequence identity with hTRPV1. In the present work, we constructed and reported the 3D homology tetramer model of hTRPV1 based on the cryo-EM structures of rTRPV1. Molecular dynamics (MD) simulations, energy minimizations, and prescreen were applied to select and validate the best model of hTRPV1. The predicted binding pocket of hTRPV1 consists of two adjacent monomers subunits, which were congruent with the experimental rTRPV1 data and the cyro-EM structures of rTRPV1. The detailed interactions between hTRPV1 and its antagonists or agonists were characterized by molecular docking, which helped us to identify the important residues. Conformational changes of hTRPV1 upon antagonist/agonist binding were also explored by MD simulation. The different movements of compounds led to the different conformational changes of monomers in hTRPV1, indicating that TRPV1 works in a concerted way, resembling some other channel proteins such as aquaporins. We observed that the selective filter was open when hTRPV1 bound with an agonist during MD simulation. For the lower gate of hTRPV1, we observed large similarities between hTRPV1 bound with antagonist and with agonist. A five-point pharmacophore model based on several antagonists was established, and the structural model was used to screen in silico for new antagonists for hTRPV1. By using the 3D TRPV1 structural model above, the pilot in silico screening has begun to yield promising hits with activity as hTRPV1 antagonists, several of which showed substantial potency.

Dataset Information

Homology modeling, molecular docking, and molecular dynamics simulations elucidated ?-fetoprotein binding modes.

Publications

Homology modeling, molecular docking, and molecular dynamics simulations elucidated α-fetoprotein binding modes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets