BACKGROUND
the entry of virus into target cell represents one of the most attractive targets in the search for new drugs to treat hiv infection.
the entry of hiv- <dig> into target cells requires the cooperation of the viral envelope glycoproteins gp <dig> and gp <dig>  and of two host-cell proteins, the primary receptor cd <dig> and a chemokine co-receptor  <cit> .
several agents have been developed to target these key regulatory proteins that are essential for hiv replication.
several of them are in clinical trials and one of them has been approved by the fda for clinical use.
therefore, drugs targeting hiv- <dig> entry are an exiting prospect in terms of prevention of aids.
recently another cell-surface protein was found to be involved in hiv- <dig> entry, the oxidoreductase protein disulfide isomerase  which catalyzes thiol-disulfide interchange reactions  <cit> .
it is present mostly in the endoplasmic reticulum and act as oxidase to forms disulfide bonds in nascent proteins and assists in protein folding  <cit> .
it also occurs at the surface of mammalian cells, where it acts as a reductase to cleave disulfide bonds of proteins attached to the cell  <cit> .
its redox function is based on the presence of two cysteine residues in its active sites cys-gly-his-cys .
when the cysteine of cxxc bears two cysteinyl thiols, it breaks neighboring disulfide bonds.
in the event of hiv- <dig> entry, the viral glycoprotein gp <dig> attaches the virus to the cell by binding to its receptor cd <dig> which also contains disulfide bonds.
after cd <dig> binding, various gp <dig> domains interact with the enzyme pdi and the chemokine co-receptors forms a pdi-cd4-gp120-chemokine complex.
pdi can reach the complex and reduce disulfide bonds in gp <dig>  which causes key conformational changes in gp <dig> and activate gp <dig> for the fusogenic potential of the viral envelope  <cit> .
it has been shown that inhibition of hiv- <dig> entry can be brought about by introducing membrane impermeant sulfhydryl agents that can block the redox function of pdi  <cit> .
these agents will stop the generation of two free thiols in a gp <dig> and an oxidized form of cxxc motif in pdi.
it was reported that the membrane-impermeant thiol reagent dithionitrobenzoic acid  causes 100% inhibition of soluble pdi activity at  <dig>  mm concentration  <cit> .
the exact mode of binding interaction is yet to be elucidated and this would give more insights into development of new effective drugs that target pdi activity.
therefore, this necessitates a rational study on the mode of binding of the inhibitors to pdi.
this can be achieved by molecular docking studies to determine whether two molecules interact and to find the orientation that maximizes this interaction as well as minimizing the total energy of the interaction complex.
predicting the mode of protein interaction with other molecules promises deduction of protein function and the enhancement of drug discovery.
a tangible example can be seen with hiv- <dig> protease  <cit> .
the current study attempts to find the mode of binding of dtnb and its related compounds on pdi.
the accelrys discovery studio and autodock  <dig>   <cit>  were used to study the interaction.
therefore, by inhibiting pdi, one can not only prevent the viral entry, but also circumvent the problem of viral resistance because pdi is a host protein  <cit> .
methods
the set of ligand molecules studied in this work include dithionitrobenzoic acid   and its structurally similar bioactive compounds such as nsc <dig>  sulfonyl-2-nitrobenzoic acid) , thionitrobenzoic acid , 2-nitro-5-thiocyanobenzoic acid , 2-nitro-5-sulfo-sulfonyl-benzoic acid  and nsc <dig>  disulfonyl-5-nitrobenzoic acid) .
these ligand molecules were retrieved from ncbi-pubchem compound database  <cit> .
the structure of these compounds is shown in figure  <dig>  the three dimensional structure of second thioredoxin-like domain of human protein disulfide-isomerase  was obtained from protein data bank   <cit> .
this structure was determined using nmr experimental method  <cit> .
the energy of the ligand molecules and target pdi were minimized in steepest descent followed by conjugate gradient method using accelrys discovery studio , the most comprehensive suite of modeling and simulation solutions for drug discovery available.
each of the minimization methods were carried out with charmm force field.
in order to carry out the docking simulation, we used the autodock  <dig>  suite as molecular-docking tool  <cit> .
it is suitable software for performing automated docking of ligands to their macromolecular receptors.
typically, the ligands are substrates or drug candidates and the macromolecule is a protein of known three dimensional structures.
in this docking simulation, we used semi-flexible docking protocols.
in which the target protein pdi was kept as rigid.
the ligands being docked were kept flexible, in order to explore an arbitrary number of torsional degrees of freedom in addition to the six spatial degrees of freedom spanned by the translational and rotational parameters.
the graphical user interface program "autodock tools" was used to prepare, run, and analyze the docking simulations.
kollman united atom charges, solvation parameters and polar hydrogens were added into the receptor pdb file for the preparation of protein in docking simulation.
this pdi enzyme structure does not have any water molecules and/or ligands to remove from its pdb file and make a free receptor.
since ligands are not peptides, gasteiger charge was assigned and then non-polar hydrogens were merged.
the rigid roots of each ligand were defined automatically instead of picking manually.
the amide bonds were made non-rotatable.
all rotable dihedrals in the ligands were assigned with the help of auxiliary program auto-tors and were allowed to rotate freely.
autodock requires pre-calculated grid maps, one for each atom type, present in the ligand being docked and its stores the potential energy arising from the interaction with macromolecule.
this grid must surround the region of interest in the macromolecule.
in the present study, the binding site was selected based on the amino acid residues, which are involved in redox function and is required for hiv entry into the host cell.
therefore, the grid was centered in the catalytic active region of redox-active cys-gly-his-cys motif and includes all amino acid residues that surround active site.
the grid box size was set at  <dig>   <dig>  and  <dig> a° , though it was changed depending on the ligand size.
autogrid  <dig>  program, supplied with autodock  <dig>  was used to produce grid maps.
the spacing between grid points was  <dig>  angstroms.
the lamarckian genetic algorithm  was chosen to search for the best conformers.
during the docking process, a maximum of  <dig> conformers was considered for each compound.
the population size was set to  <dig> and the individuals were initialized randomly.
maximum number of energy evaluation was set to  <dig>  maximum number of generations  <dig>  maximum number of top individual that automatically survived set to  <dig>  mutation rate of  <dig> , crossover rate of  <dig> , step sizes were  <dig>  Å for translations,  <dig> ° for quaternions and  <dig> ° for torsions.
cluster tolerance  <dig> a°, external grid energy  <dig> , max initial energy  <dig> , max number of retries  <dig> and  <dig> lga runs were performed.
all the autodock docking runs were performed in intel pentium pd- <dig> cpu @  <dig>  ghz of hcl infosystem origin, with  <dig> gb ddr ram.
autodock  <dig>  was compiled and run under windows xp operating system.
RESULTS
existing nmr structure of protein disulfide isomerase was downloaded from pdb .
an examination of this nmr structure for the domain of pdi showed the cysteines in their cxxc motifs  are in the form of free thiols.
at the cell surface, this isomerase predominantly acts as reductase leads to the cleavage of disulfide bonds and production of free thiols in the interacting protein .
the criteria which we have followed in the selection of ligands includes only targeting cell surface pdi by using membrane impermeant sulfhydryl inhibitors.
investigational membrane impermeant pdi inhibitors include thiol blockers such as dtnb, nsc <dig>  thionitrobenzoic acid, 2-nitro-5-thiocyanobenzoic acid, 2-nitro-5-sulfo-sulfonyl-benzoic acid and nsc <dig> were retrieved from pubchem database.
all the ligands  and pdi structures were confirmed to have minimized energy based on the final potential energy obtained from the results as shown in the table  <dig>  the potential energy value of all ligands was below - <dig> kcal/mol.
these ligand molecules have the potential to change the original states of it is into a system of binding state in a target pdi molecule when the energy is released.
the potential energy of pdi was - <dig>   <dig> kcal/mol.
molecular docking simulations were conducted with the autodock  <dig>  software suite.
<dig> docking runs were performed for each of the six different ligands into redox-active site  of pdi protein.
for each docking simulation, the activity inferred by redox inhibitory docked models, in which the docked conformers of ligands and pdi would permit redox active interaction  and matched with the wet lab experimental observation results of inhibition  <cit> .
to verify the accuracy of the autodock  <dig>  results, we also considered some top clusters of conformations/orientations in addition to the best scored one.
here, the docking accuracy was evaluated in terms of the root mean square deviation  between the docked position and the experimentally determined position for the ligand.
in this molecular docking study, prediction was considered successful if the rmsd value is less then  <dig>  Å for the best-scored conformation  <cit> .
the redox inhibitory conformations were energetically and statistically validated.
docking of dtnb into pdi
as shown in table  <dig>  the docking results are ranked according to the ascent of the binding energies for each of the ligands investigated.
docking simulation of dtnb into the redox active site of the pdi produced eight clusters of conformers using rmsd-tolerance of  <dig>  Å out of  <dig> docking runs .
the conformation of the # <dig> ranked cluster was favored in that structure and repeated two times out of  <dig> runs.
inspection of this conformation showed that it could explain the non-redox inhibitory mode of pdi , since docked region was out of redox active site.
the conformation of the # <dig> ranked cluster of dtnb in the pdi was also highly energetically and statistically favored.
in which, it was observed that two redox inhibitory docking modes were found in same cluster rank, which would explain the inhibitor effect of the enzyme.
in these two conformations, there were, three, two hydrogen bonds formed, respectively .
the thiol group of cys <dig> and backbone-oxygen atom of phe <dig> in pdi enzyme interacts with the two different atom of sulfur  in dtnb to form multiple hydrogen bonds in both the two conformations with a bond length of approximately  <dig>  Å, respectively.
in addition, side chain nitrogen atom of arg <dig> was also involved in the hydrogen bond interaction with o <dig> atom of dtnb in the lowest binding energy complex of # <dig> ranked clusters.
the dtnb was also made additional van der waals forces with the residues ala <dig>  trp <dig>  his <dig>  cys <dig>  and pro <dig> in the scaling factor around  <dig>  Å to further stables the interaction .
the lowest binding energy of the docked complex was  <dig>  kcal/mol.
the lowest cluster root mean square deviation  between the docked conformation and the conformation in the crystal structure was  <dig>  Å. the conformations of other clusters of dtnb in pdi were all non-redox inhibitory model.
in these instances, dtnb was bound in similar fashion to that seen with cluster rank # <dig>
docking of nsc <dig> into pdi
docking simulation of nsc <dig>  sulfonyl-2-nitrobenzoic acid) into pdi resulted in five clusters of conformers .
the docked models of all the clusters appeared non-redox inhibitory mode, since docked region of the compound in out side of the redox active site.
in detailed inspection of sub ranked cluster conformations, it was observed that the second sub-rank of # <dig> cluster represented redox inhibitory mode.
in this docked complex, side chain nh group of trp <dig> and oh <dig> group of nsc <dig> acts as hydrogen bond donors to make multiple hydrogen bonds with sulfur atom of cys <dig> and o <dig> atom of nsc <dig>  respectively.
the bond distance between donors and acceptors atoms of hydrogen bond was in approximately  <dig>  Å .
in addition to hydrogen bonds, residues ala  <dig>  his <dig>  thr <dig>  ser <dig>  and phe <dig> were involved in van der waals interactions further stabilized the ligand .
among these residues ala <dig> and his <dig> were in similar binding mode as compared to # <dig> ranked cluster of dtnb docking.
the conformations of the # <dig> ranked cluster of nsc <dig> in pdi were generated five more sub-ranked complexes with non-redox inhibitory modes.
these conformations had significantly high frequency of non-redox docking occurrences than those of other ranked clusters.
docking of thionitrobenzoic acid into pdi
the ten runs of docking simulations of thionitrobenzoic acid to pdi resulted six clusters of conformers.
the conformations of the ranked and sub-ranked # <dig> cluster were highly favored statistically and energetically and repeated four times  with the high frequency of occurrences.
analysis of the positions and orientations of the ligand in the docking models in relationship to the location of the redox-avtice site indicating that the first and second sub-ranked conformation of the # <dig> cluster represented redox inhibitory mode.
in these two conformations, there were same atoms involved in hydrogen bond formation.
in which, thiol group of cys <dig> participated in hydrogen bond formation with the o <dig> atom of the ligand .
as a result, conformations were having similar binding energies with similar bond length of  <dig>  Å. in addition to cys <dig> residue, the thionitrobenzate was held in position by van der waals forces with the residues trp <dig> and phe <dig> .
these same residues were participated in the redox inhibitory mode of dtnb in pdi enzyme.
the cluster conformers of the other ranked clusters were in non-redox inhibitory modes.
docking of 2-nitro-5-thiocyanobenzoic acid into pdi
the results for docking 2-nitro-5-thiocyanobenzoic acid into pdi enzyme are given in table  <dig>  one redox inhibitory mode was found  with the 2-nitro-5-thiocyanobenzoic acid, which could explain the redox inhibitory activity of the enzyme.
this conformation was energetically favored.
the amino acid residues of pdi involved in the interaction with 2-nitro-5-thiocyanobenzoic acid, which were mostly similar to those residues involved in the redox-inhibitory mode of dtnb in pdi .
in the conformation of # <dig> ranked cluster, there were four hydrogen bonds formed between the 2-nitro-5-thiocyanobenzoic acid with residues cys <dig>  cys <dig>  phe <dig> of pdi target with binding energy of  <dig>  kcal/mol.
in is interesting to note that the redox-active sites of cys <dig> and cys <dig> were both involved in hydrogen bond formation.
in addition, s-h and s atoms of cys <dig> itself acts as hydrogen bond donor and acceptor to interact with o <dig> and o-h <dig> atoms of ligand in the formation of hydrogen bonds, respectively.
this cluster docked conformation was also stabled by van der waals with ala <dig>  trp <dig> and thr <dig>
docking of 2-nitro-5-sulfosulfanyl-benzoic acid into pdi
the number of docked clusters observed with 2-nitro-5-sulfo-sulfonyl-benzoic acid for pdi was less frequency of occurrences, since produced seven different clusters with single conformer and one cluster  with three conformers , indicating a decreased specificity of this ligand with pdi.
there were two clusters  of conformers appeared as redox-active mode out of eight clusters generated.
the lowest binding energy of #  <dig> ranked clusters was observed with  <dig>  kcal/mol and it was significantly lower than the other ligands docked.
there was a single hydrogen bond formed between electronegative atoms of sulfur in cys <dig> with o-h <dig> group of 2-nitro-5-sulfosulfanyl-benzoic acid .
the residues ala <dig>  trp <dig>  gly <dig>  his <dig> and phe <dig> were also taken part in van der waals interactions with quite same as redox-inhibitory mode of dtnb in pdi in the scaling factor around  <dig>  Å.
docking of nsc <dig> into pdi
the docked models of clusters of nsc <dig>  disulfanyl-5-nitrobenzoic acid) in pdi all appeared non-redox inhibitory mode, since the ligand was placed far away from the redox-active binding site.
in the second sub-ranked conformation of the # <dig> ranked cluster, nsc <dig> was bound to pdi in a fashion similar to the redox inhibitory modes of dtnb, 2-nitro-5-sulfosulfanyl-benzoic acid and nsc <dig> in pdi such that it was the quite same amino acid residues involved in the molecular interaction .
in the hydrogen bond formation, sulfur atom of cys <dig> acts as hydrogen bond acceptors to made interaction with o-h <dig> group of nsc <dig> with bond length of  <dig>  Å as approximately similar to 2-nitro-5-sulfo-sulfonyl-benzoic acid docking in concern to both numbers of hydrogen bond and bond length.
the van der waals interactions between the pdi and nsc <dig> were observed in the residues of ala <dig>  trp <dig>  his <dig>  ser <dig>  thr <dig> and phe <dig> in a maximum similar to those observed in the case of 2-nitro-5-sulfo-sulfonyl-benzoic acid in pdi.
as shown in table  <dig>  it was observed that all the six different ligand molecules participated in hydrogen bond formation with sulphur atom of the cys <dig> and/or cys <dig> in pdi.
in addition to this, sulfur atom in dithionitrobenzoic acid itself act as hydrogen bond acceptor.
sulfur atoms are larger and have a more diffuse electron cloud than oxygen and nitrogen.
the o-h...s- hydrogen bond strength in cysteine thiolate anion was estimated experimentally as  <dig>  ±  <dig>  kcal/mol  <cit> .
this makes binding energy of the complex in positive score.
lipinski's rule of five was calculated for all the six ligand molecules that satisfy the 'rule-of-5' and it was found that all the ligand molecules satisfied the rule for potent inhibitors .
CONCLUSIONS
in this work, a molecular docking simulation study was undertaken to investigate the binding mechanism of dtnb and its structurally similar membrane impermeable inhibitors such as dtnb, nsc <dig>  thionitrobenzoic acid, 2-nitro-5-thiocyanobenzoic acid, 2-nitro-5-sulfo-sulfonyl-benzoic acid and nsc <dig> to the protein disulfide isomerase enzyme to enable the finding of potential anti-hiv drugs .
we determined the energetically favored docking sites for the inhibitors.
our docking results explain that the numbers of clusters with each ligand vary in the range from five to eight for pdi enzyme , indicating that the binding specificity of each ligand is varying in pdi.
in this study, the hydrogen bond makes important contributions to the interactions between ligand and enzyme.
from the frequency of residue's occurrence in the formation of hydrogen bonding, it is evident that cys <dig> plays an important role due to the presence of s-h group and it acts as both hydrogen bond donor and acceptor.
and trp <dig>  cys <dig>  phe <dig>  ag <dig> also take part in the hydrogen bonding with comparatively high frequency.
this information could be supportive for new drug design for aids, in which these potential inhibitors should interact strongly with above-mentioned residues.
cys <dig> is of great concern in anti-hiv drug design, which is believed to be one of redox-active site in pdi  <cit> .
it is significant to have this cys <dig> s-h group locked such that its reduction function is depressed.
therefore the above-mentioned inhibitors will block the cleavage of disulfide bond in gp <dig>  leading to inactive gp <dig> and resistance to the entry of hiv.
it was observed that there were maximum of four hydrogen bonds formed between pdi and the inhibitors used .
besides hydrogen bonding, van der waals interactions were also taking part in the stabilizations of inhibitors binding with high frequency of residues such as ala <dig>  trp <dig>  cys <dig>  his <dig>  thr <dig> and phe <dig> in pdi.
the redox-inhibitory mode of all six inhibitors with pdi was consistent with the laboratory experimental result of ryser et al  <cit> .
competing interests
the authors declare that they have no competing interests.
authors' contributions
ug, mj and ds designed the methods and experimental setup.
ug and mj carried out the implementation of the various methods.
mj and ds wrote the manuscript.
all authors have read and approved the final manuscript.

 the entry of virus into target cell represents one of the most attractive targets in the search for new drugs to treat hiv infection. the entry of hiv- into target cells requires the cooperation of the viral envelope glycoproteins gp and gp and of two host-cell proteins, the primary receptor cd and a chemokine co-receptor. several agents have been developed to target these key regulatory proteins that are essential for hiv replication. several of them are in clinical trials and one of them has been approved by the fda for clinical use. therefore, drugs targeting hiv- entry are an exiting prospect in terms of prevention of aids. recently another cell-surface protein was found to be involved in hiv- entry, the oxidoreductase protein disulfide isomerase which catalyzes thiol-disulfide interchange reactions. it is present mostly in the endoplasmic reticulum and act as oxidase to forms disulfide bonds in nascent proteins and assists in protein folding. it also occurs at the surface of mammalian cells, where it acts as a reductase to cleave disulfide bonds of proteins attached to the cell. its redox function is based on the presence of two cysteine residues in its active sites cys-gly-his-cys. when the cysteine of cxxc bears two cysteinyl thiols, it breaks neighboring disulfide bonds. in the event of hiv- entry, the viral glycoprotein gp attaches the virus to the cell by binding to its receptor cd which also contains disulfide bonds. after cd binding, various gp domains interact with the enzyme pdi and the chemokine co-receptors forms a pdi-cd4-gp120-chemokine complex. pdi can reach the complex and reduce disulfide bonds in gp which causes key conformational changes in gp and activate gp for the fusogenic potential of the viral envelope. it has been shown that inhibition of hiv- entry can be brought about by introducing membrane impermeant sulfhydryl agents that can block the redox function of pdi. these agents will stop the generation of two free thiols in a gp and an oxidized form of cxxc motif in pdi. it was reported that the membrane-impermeant thiol reagent dithionitrobenzoic acid causes 100% inhibition of soluble pdi activity at mm concentration. the exact mode of binding interaction is yet to be elucidated and this would give more insights into development of new effective drugs that target pdi activity. therefore, this necessitates a rational study on the mode of binding of the inhibitors to pdi. this can be achieved by molecular docking studies to determine whether two molecules interact and to find the orientation that maximizes this interaction as well as minimizing the total energy of the interaction complex. predicting the mode of protein interaction with other molecules promises deduction of protein function and the enhancement of drug discovery. a tangible example can be seen with hiv- protease. the current study attempts to find the mode of binding of dtnb and its related compounds on pdi. the accelrys discovery studio and autodock were used to study the interaction. therefore, by inhibiting pdi, one can not only prevent the viral entry, but also circumvent the problem of viral resistance because pdi is a host protein. methods the set of ligand molecules studied in this work include dithionitrobenzoic acid and its structurally similar bioactive compounds such as nsc sulfonyl-2-nitrobenzoic acid), thionitrobenzoic acid, 2-nitro-5-thiocyanobenzoic acid, 2-nitro-5-sulfo-sulfonyl-benzoic acid and nsc disulfonyl-5-nitrobenzoic acid). these ligand molecules were retrieved from ncbi-pubchem compound database. the structure of these compounds is shown in figure the three dimensional structure of second thioredoxin-like domain of human protein disulfide-isomerase was obtained from protein data bank. this structure was determined using nmr experimental method. the energy of the ligand molecules and target pdi were minimized in steepest descent followed by conjugate gradient method using accelrys discovery studio, the most comprehensive suite of modeling and simulation solutions for drug discovery available. each of the minimization methods were carried out with charmm force field. in order to carry out the docking simulation, we used the autodock suite as molecular-docking tool. it is suitable software for performing automated docking of ligands to their macromolecular receptors. typically, the ligands are substrates or drug candidates and the macromolecule is a protein of known three dimensional structures. in this docking simulation, we used semi-flexible docking protocols. in which the target protein pdi was kept as rigid. the ligands being docked were kept flexible, in order to explore an arbitrary number of torsional degrees of freedom in addition to the six spatial degrees of freedom spanned by the translational and rotational parameters. the graphical user interface program "autodock tools" was used to prepare, run, and analyze the docking simulations. kollman united atom charges, solvation parameters and polar hydrogens were added into the receptor pdb file for the preparation of protein in docking simulation. this pdi enzyme structure does not have any water molecules and/or ligands to remove from its pdb file and make a free receptor. since ligands are not peptides, gasteiger charge was assigned and then non-polar hydrogens were merged. the rigid roots of each ligand were defined automatically instead of picking manually. the amide bonds were made non-rotatable. all rotable dihedrals in the ligands were assigned with the help of auxiliary program auto-tors and were allowed to rotate freely. autodock requires pre-calculated grid maps, one for each atom type, present in the ligand being docked and its stores the potential energy arising from the interaction with macromolecule. this grid must surround the region of interest in the macromolecule. in the present study, the binding site was selected based on the amino acid residues, which are involved in redox function and is required for hiv entry into the host cell. therefore, the grid was centered in the catalytic active region of redox-active cys-gly-his-cys motif and includes all amino acid residues that surround active site. the grid box size was set at and a°, though it was changed depending on the ligand size. autogrid program, supplied with autodock was used to produce grid maps. the spacing between grid points was angstroms. the lamarckian genetic algorithm was chosen to search for the best conformers. during the docking process, a maximum of conformers was considered for each compound. the population size was set to and the individuals were initialized randomly. maximum number of energy evaluation was set to maximum number of generations maximum number of top individual that automatically survived set to mutation rate of, crossover rate of, step sizes were Å for translations, ° for quaternions and ° for torsions. cluster tolerance a°, external grid energy, max initial energy, max number of retries and lga runs were performed. all the autodock docking runs were performed in intel pentium pd- cpu @ ghz of hcl infosystem origin, with gb ddr ram. autodock was compiled and run under windows xp operating system.  existing nmr structure of protein disulfide isomerase was downloaded from pdb. an examination of this nmr structure for the domain of pdi showed the cysteines in their cxxc motifs are in the form of free thiols. at the cell surface, this isomerase predominantly acts as reductase leads to the cleavage of disulfide bonds and production of free thiols in the interacting protein. the criteria which we have followed in the selection of ligands includes only targeting cell surface pdi by using membrane impermeant sulfhydryl inhibitors. investigational membrane impermeant pdi inhibitors include thiol blockers such as dtnb, nsc thionitrobenzoic acid, 2-nitro-5-thiocyanobenzoic acid, 2-nitro-5-sulfo-sulfonyl-benzoic acid and nsc were retrieved from pubchem database. all the ligands and pdi structures were confirmed to have minimized energy based on the final potential energy obtained from the results as shown in the table the potential energy value of all ligands was below - kcal/mol. these ligand molecules have the potential to change the original states of it is into a system of binding state in a target pdi molecule when the energy is released. the potential energy of pdi was - kcal/mol. molecular docking simulations were conducted with the autodock software suite. docking runs were performed for each of the six different ligands into redox-active site of pdi protein. for each docking simulation, the activity inferred by redox inhibitory docked models, in which the docked conformers of ligands and pdi would permit redox active interaction and matched with the wet lab experimental observation results of inhibition. to verify the accuracy of the autodock results, we also considered some top clusters of conformations/orientations in addition to the best scored one. here, the docking accuracy was evaluated in terms of the root mean square deviation between the docked position and the experimentally determined position for the ligand. in this molecular docking study, prediction was considered successful if the rmsd value is less then Å for the best-scored conformation. the redox inhibitory conformations were energetically and statistically validated. docking of dtnb into pdi as shown in table the docking results are ranked according to the ascent of the binding energies for each of the ligands investigated. docking simulation of dtnb into the redox active site of the pdi produced eight clusters of conformers using rmsd-tolerance of Å out of docking runs. the conformation of the # ranked cluster was favored in that structure and repeated two times out of runs. inspection of this conformation showed that it could explain the non-redox inhibitory mode of pdi, since docked region was out of redox active site. the conformation of the # ranked cluster of dtnb in the pdi was also highly energetically and statistically favored. in which, it was observed that two redox inhibitory docking modes were found in same cluster rank, which would explain the inhibitor effect of the enzyme. in these two conformations, there were, three, two hydrogen bonds formed, respectively. the thiol group of cys and backbone-oxygen atom of phe in pdi enzyme interacts with the two different atom of sulfur in dtnb to form multiple hydrogen bonds in both the two conformations with a bond length of approximately Å, respectively. in addition, side chain nitrogen atom of arg was also involved in the hydrogen bond interaction with o atom of dtnb in the lowest binding energy complex of # ranked clusters. the dtnb was also made additional van der waals forces with the residues ala trp his cys and pro in the scaling factor around Å to further stables the interaction. the lowest binding energy of the docked complex was kcal/mol. the lowest cluster root mean square deviation between the docked conformation and the conformation in the crystal structure was Å. the conformations of other clusters of dtnb in pdi were all non-redox inhibitory model. in these instances, dtnb was bound in similar fashion to that seen with cluster rank # docking of nsc into pdi docking simulation of nsc sulfonyl-2-nitrobenzoic acid) into pdi resulted in five clusters of conformers. the docked models of all the clusters appeared non-redox inhibitory mode, since docked region of the compound in out side of the redox active site. in detailed inspection of sub ranked cluster conformations, it was observed that the second sub-rank of # cluster represented redox inhibitory mode. in this docked complex, side chain nh group of trp and oh group of nsc acts as hydrogen bond donors to make multiple hydrogen bonds with sulfur atom of cys and o atom of nsc respectively. the bond distance between donors and acceptors atoms of hydrogen bond was in approximately Å. in addition to hydrogen bonds, residues ala his thr ser and phe were involved in van der waals interactions further stabilized the ligand. among these residues ala and his were in similar binding mode as compared to # ranked cluster of dtnb docking. the conformations of the # ranked cluster of nsc in pdi were generated five more sub-ranked complexes with non-redox inhibitory modes. these conformations had significantly high frequency of non-redox docking occurrences than those of other ranked clusters. docking of thionitrobenzoic acid into pdi the ten runs of docking simulations of thionitrobenzoic acid to pdi resulted six clusters of conformers. the conformations of the ranked and sub-ranked # cluster were highly favored statistically and energetically and repeated four times with the high frequency of occurrences. analysis of the positions and orientations of the ligand in the docking models in relationship to the location of the redox-avtice site indicating that the first and second sub-ranked conformation of the # cluster represented redox inhibitory mode. in these two conformations, there were same atoms involved in hydrogen bond formation. in which, thiol group of cys participated in hydrogen bond formation with the o atom of the ligand. as a result, conformations were having similar binding energies with similar bond length of Å. in addition to cys residue, the thionitrobenzate was held in position by van der waals forces with the residues trp and phe. these same residues were participated in the redox inhibitory mode of dtnb in pdi enzyme. the cluster conformers of the other ranked clusters were in non-redox inhibitory modes. docking of 2-nitro-5-thiocyanobenzoic acid into pdi the results for docking 2-nitro-5-thiocyanobenzoic acid into pdi enzyme are given in table one redox inhibitory mode was found with the 2-nitro-5-thiocyanobenzoic acid, which could explain the redox inhibitory activity of the enzyme. this conformation was energetically favored. the amino acid residues of pdi involved in the interaction with 2-nitro-5-thiocyanobenzoic acid, which were mostly similar to those residues involved in the redox-inhibitory mode of dtnb in pdi. in the conformation of # ranked cluster, there were four hydrogen bonds formed between the 2-nitro-5-thiocyanobenzoic acid with residues cys cys phe of pdi target with binding energy of kcal/mol. in is interesting to note that the redox-active sites of cys and cys were both involved in hydrogen bond formation. in addition, s-h and s atoms of cys itself acts as hydrogen bond donor and acceptor to interact with o and o-h atoms of ligand in the formation of hydrogen bonds, respectively. this cluster docked conformation was also stabled by van der waals with ala trp and thr docking of 2-nitro-5-sulfosulfanyl-benzoic acid into pdi the number of docked clusters observed with 2-nitro-5-sulfo-sulfonyl-benzoic acid for pdi was less frequency of occurrences, since produced seven different clusters with single conformer and one cluster with three conformers, indicating a decreased specificity of this ligand with pdi. there were two clusters of conformers appeared as redox-active mode out of eight clusters generated. the lowest binding energy of # ranked clusters was observed with kcal/mol and it was significantly lower than the other ligands docked. there was a single hydrogen bond formed between electronegative atoms of sulfur in cys with o-h group of 2-nitro-5-sulfosulfanyl-benzoic acid. the residues ala trp gly his and phe were also taken part in van der waals interactions with quite same as redox-inhibitory mode of dtnb in pdi in the scaling factor around Å. docking of nsc into pdi the docked models of clusters of nsc disulfanyl-5-nitrobenzoic acid) in pdi all appeared non-redox inhibitory mode, since the ligand was placed far away from the redox-active binding site. in the second sub-ranked conformation of the # ranked cluster, nsc was bound to pdi in a fashion similar to the redox inhibitory modes of dtnb, 2-nitro-5-sulfosulfanyl-benzoic acid and nsc in pdi such that it was the quite same amino acid residues involved in the molecular interaction. in the hydrogen bond formation, sulfur atom of cys acts as hydrogen bond acceptors to made interaction with o-h group of nsc with bond length of Å as approximately similar to 2-nitro-5-sulfo-sulfonyl-benzoic acid docking in concern to both numbers of hydrogen bond and bond length. the van der waals interactions between the pdi and nsc were observed in the residues of ala trp his ser thr and phe in a maximum similar to those observed in the case of 2-nitro-5-sulfo-sulfonyl-benzoic acid in pdi. as shown in table it was observed that all the six different ligand molecules participated in hydrogen bond formation with sulphur atom of the cys and/or cys in pdi. in addition to this, sulfur atom in dithionitrobenzoic acid itself act as hydrogen bond acceptor. sulfur atoms are larger and have a more diffuse electron cloud than oxygen and nitrogen. the o-h...s- hydrogen bond strength in cysteine thiolate anion was estimated experimentally as ± kcal/mol. this makes binding energy of the complex in positive score. lipinski's rule of five was calculated for all the six ligand molecules that satisfy the'rule-of-5' and it was found that all the ligand molecules satisfied the rule for potent inhibitors.  in this work, a molecular docking simulation study was undertaken to investigate the binding mechanism of dtnb and its structurally similar membrane impermeable inhibitors such as dtnb, nsc thionitrobenzoic acid, 2-nitro-5-thiocyanobenzoic acid, 2-nitro-5-sulfo-sulfonyl-benzoic acid and nsc to the protein disulfide isomerase enzyme to enable the finding of potential anti-hiv drugs. we determined the energetically favored docking sites for the inhibitors. our docking results explain that the numbers of clusters with each ligand vary in the range from five to eight for pdi enzyme, indicating that the binding specificity of each ligand is varying in pdi. in this study, the hydrogen bond makes important contributions to the interactions between ligand and enzyme. from the frequency of residue's occurrence in the formation of hydrogen bonding, it is evident that cys plays an important role due to the presence of s-h group and it acts as both hydrogen bond donor and acceptor. and trp cys phe ag also take part in the hydrogen bonding with comparatively high frequency. this information could be supportive for new drug design for aids, in which these potential inhibitors should interact strongly with above-mentioned residues. cys is of great concern in anti-hiv drug design, which is believed to be one of redox-active site in pdi. it is significant to have this cys s-h group locked such that its reduction function is depressed. therefore the above-mentioned inhibitors will block the cleavage of disulfide bond in gp leading to inactive gp and resistance to the entry of hiv. it was observed that there were maximum of four hydrogen bonds formed between pdi and the inhibitors used. besides hydrogen bonding, van der waals interactions were also taking part in the stabilizations of inhibitors binding with high frequency of residues such as ala trp cys his thr and phe in pdi. the redox-inhibitory mode of all six inhibitors with pdi was consistent with the laboratory experimental result of ryser et al. competing interests the authors declare that they have no competing interests. authors' contributions ug, mj and ds designed the methods and experimental setup. ug and mj carried out the implementation of the various methods. mj and ds wrote the manuscript. all authors have read and approved the final manuscript. 
we also identified that cys of the enzyme plays an important role in hydrogen bonding with inhibitors. we performed molecular docking simulation with six different inhibitors, which includes dtnb, nsc thionitrobenzoic acid, 2-nitro-5-thiocyanobenzoic acid, 2-nitro-5-sulfo-sulfonyl-benzoic acid and nsc into the redox-active site of the pdi enzyme and the activity was inferred by redox inhibitory models. this residue can be considered to being an active site for anti-hiv drug design. there have been studies with agents that inhibit pdi activity, but the exact mode of binding remains to be elucidated; this might provide insights to develop new drugs to target pdi. therefore, by inhibiting pdi, one can, not only prevent the viral entry but also circumvent the problem of viral resistance entry of hiv- into human lymphoid requires activities of viral envelope glycoproteins gp and gp and two host-cell proteins, the primary receptor cd and a chemokine co-receptor. this study attempts to perceive the mode of binding of dithionitrobenzoic acid, and its structurally related compounds on pdi enzyme. in addition, a third cell-surface protein called protein disulfide isomerase is found to play a major role in hiv- entry. in this scenario, inhibition of hiv- entry can be brought about by introducing agents that can block thiol-disulfide interchange reaction of cell surface pdi. all ligands showed favorable interactions and most of them seemed to bind to hydrophobic amino acids ala trp cys cys his thr and phe the redox inhibitory conformations were energetically and statistically favored and supported the evidence from wet laboratory experiments reported in the literature.
 entry of hiv- into human lymphoid requires activities of viral envelope glycoproteins gp and gp and two host-cell proteins, the primary receptor cd and a chemokine co-receptor. in addition, a third cell-surface protein called protein disulfide isomerase is found to play a major role in hiv- entry. pdi is capable of mediating thio-disulfide interchange reactions and could enable the reduction of gp disulfide bonds, which triggers the major conformational changes in gp and gp required for virus entry. in this scenario, inhibition of hiv- entry can be brought about by introducing agents that can block thiol-disulfide interchange reaction of cell surface pdi. there have been studies with agents that inhibit pdi activity, but the exact mode of binding remains to be elucidated; this might provide insights to develop new drugs to target pdi. this study attempts to perceive the mode of binding of dithionitrobenzoic acid, and its structurally related compounds on pdi enzyme.  we performed molecular docking simulation with six different inhibitors, which includes dtnb, nsc thionitrobenzoic acid, 2-nitro-5-thiocyanobenzoic acid, 2-nitro-5-sulfo-sulfonyl-benzoic acid and nsc into the redox-active site of the pdi enzyme and the activity was inferred by redox inhibitory models. all ligands showed favorable interactions and most of them seemed to bind to hydrophobic amino acids ala trp cys cys his thr and phe the redox inhibitory conformations were energetically and statistically favored and supported the evidence from wet laboratory experiments reported in the literature.  we demonstrated that in silico docking experiment can be effectively carried out to recognize the redox inhibitory models of pdi with inhibitor molecules. interestingly we found that number of docked clusters with each ligand varies in the range of five to eight and conveys that the binding specificity of each inhibitor varies for pdi. we also identified that cys of the enzyme plays an important role in hydrogen bonding with inhibitors. this residue can be considered to being an active site for anti-hiv drug design. therefore, by inhibiting pdi, one can, not only prevent the viral entry but also circumvent the problem of viral resistance 20– october asia pacific bioinformatics network seventh international conference on bioinformatics taipei, taiwan 