Search

Introduction: GTPases known as translation factor play a vital role as ribosomal subunit assembly chaperone. The bacterial Obg proteins (-associated -binding protein) belong to the subfamily of P-loop GTPase proteins and now it is considered as one of the new target for antibacterial drug. The majority of bacterial Obgs have been commonly found to be associated with ribosome, implying that these proteins may play a fundamental role in ribosome assembly or maturation. In addition, one of the experimental evidences suggested that Bacillus subtilis Obg (BsObg) protein binds to the L13 ribosomal protein (BsL13) which is known to be one of the early assembly proteins of the 50S ribosomal subunit in Escherichia coli. In order to investigate binding mode between the BsObg and the BsL13, protein-protein docking simulation was carried out after generating 3D structure of the BsL13 structure using homology modeling method. Materials and Methods: Homology model structure of BsL13 was generated using the EcL13 crystal structure as a template. Protein-protein docking of BsObg protein with ribosomal protein BsL13 was performed by DOT, a macro-molecular docking software, in order to predict a reasonable binding mode. The solvated energy minimization calculation of the docked conformation was carried out to refine the structure. Results and Discussion: The possible binding conformation of BsL13 along with activated Obg fold in BsObg was predicted by computational docking study. The final structure is obtained from the solvated energy minimization. From the analysis, three important H-bond interactions between the Obg fold and the L13 were detected: Obg:Tyr27-L13:Glu32, Obg:Asn76-L13:Glu139, and Obg:Ala136-L13:Glu142. The interaction between the BsObg and BsL13 structures were also analyzed by electrostatic potential calculations to examine the interface surfaces. From the results, the key residues for hydrogen bonding and hydrophobic interaction between the two proteins were predicted. Conclusion and Prospects: In this study, we have focused on the binding mode of the BsObg protein with the ribosomal BsL13 protein. The interaction between the activated Obg and target protein was investigated with protein-protein docking calculations. The binding pattern can be further used as a base for structure-based drug design to find a novel antibacterial drug.

Adenosine kinase (AK) is a ubiquitous intracellular enzyme, which catalyzes the phosphorylation of adenosine (ADO) to adenosine monophosphate (AMP). AK inhibitors have therapeutic potential as analgesic and antiinflammatory agents. A chemical feature based pharmacophore model has been generated from known AK inhibitors (26 training set compounds) by HypoGen module implemented in CATALYST software. The top ranked hypothesis (Hypo1) contained four features of two hydrogen-bond acceptors (HBA) and two hydrophobic aromatics (Z). Hypo1 was validated by 124 test set molecules with a correlation coefficient of 0.905 between experimental and estimated activity. It was also validated by CatScramble method. Thus, the Hypo1 was exploited for searching new lead compounds over 238,819 chemical compounds in NCI database and then the selected compounds were screened based on restriction estimated activity and Lipinski's rules to evaluate their drug-like properties. Finally we could obtain 72 new lead candidates and the two best compound structures from them were posted.

The activation process of the redox-regulated chaperone heat shock protein 33 (Hsp33) is constituted by the oxidation-induced unfolding of the C-terminal zinc-binding domain and concomitant oligomerization of the N-terminal core domain. Herein, the semi-empirical solution structure of Escherichia coli Hsp33 in the reduced, inactive form was generated through conformational space annealing calculations, utilizing minimalistic NMR data and multiple homology restraints. The various conformations of oxidized Hsp33 and some mutant forms were also investigated in solution. Interestingly, a specific region concentrated around the interdomain linker stretch and its interacting counterparts, the N-terminal β-strand 1 and α-helix 1, hardly showed up as signals in the NMR measurements. The NMR spectra of an Hsp33 derivative with a six-residue deletion in the disordered N-terminus implied a plausible conformational exchange associated with the identified region, and the corresponding exchange rate appeared slower than that of the wild type. Subsequent mutations that destroyed the structure of the β1 or α1 elements resulted in the formation of a reduced but active monomer, without the unfolding of the zinc-binding domain. Collectively, structural insights into the inactive and active conformations, including wild-type and mutant proteins, suggest that the dynamic interactions of the N-terminal segments with their contacting counterpart, the interdomain linker stretch, in the reduced, inactive state are the structural determinants regulating the activation process of the post-translationally regulated chaperone, Hsp33.

AtTDX, a thioredoxin-like plant-specific protein present in Arabidospis is a thermo-stable and multi-functional enzyme. This enzyme is known to act as a thioredoxin and as a molecular chaperone depending upon its oligomeric status. The present study examines the effects of γ-irradiation on the structural and functional changes of AtTDX. Holdase chaperone activity of AtTDX was increased and reached a maximum at 10 kGy of γ-irradiation and declined subsequently in a dose-dependent manner, together with no effect on foldase chaperone activity. However, thioredoxin activity decreased gradually with increasing irradiation. Electrophoresis and size exclusion chromatography analysis showed that AtTDX had a tendency to form high molecular weight (HMW) complexes after γ-irradiation and γ-ray-induced HMW complexes were tightly associated with a holdase chaperone activity. The hydrophobicity of AtTDX increased with an increase in irradiation dose till 20 kGy and thereafter decreased further. Analysis of the secondary structures of AtTDX using far UV-circular dichroism spectra revealed that the irradiation remarkably increased the exposure of β-sheets and random coils with a dramatic decrease in α-helices and turn elements in a dose-dependent manner. The data of the present study suggest that γ-irradiation may be a useful tool for increasing holdase chaperone activity without adversely affecting foldase chaperone activity of thioredoxin-like proteins.

While structural and dynamical properties of water on a labile metal ion such as Na+ and Ca2+ are straightforwardly evaluated using molecular dynamics (MD) simulations, direct evaluation of water exchange kinetics on a relatively more inert metal ion, Mg2+, is computationally demanding. Here, we propose a universal relationship which allows us to determine the water exchange time on a metal ion with varying charges and hydration radii. The proposed relationship highlights the ultra-sensitivity of water lifetime to the ion size, as exemplified between two divalent ions, Ca2+ (∼100 ps) and Mg2+ (∼1.5 μs). Importantly, the current MD force fields underestimate the hydration radius of Mg2+ by ∼2.5 %, which results in ∼ 400 fold overestimation of water exchange time. Accurate determination of ion size is critical for correctly modeling the water dynamics around multivalent ions.

Repeat proteins have recently attracted much attention as alternative scaffolds to immunoglobulin antibodies due to their unique structural and biophysical features. In particular, repeat proteins show high stability against temperature and chaotic agents. Despite many studies, structural features for the stability of repeat proteins remain poorly understood. Here we present an interesting result from in silico analyses pursuing the factors which affect the stability of repeat proteins. Previously developed repebody structure based on variable lymphocytes receptors (VLRs) which consists of leucine-rich repeat (LRR) modules was used as initial structure for the present study. We constructed extra six repebody structures with varying numbers of repeat modules and those structures were used for molecular dynamics simulations. For the structures, the intramolecular interactions including backbone H-bonds, van der Waals energy, and hydrophobicity were investigated and then the radius of gyration, solvent-accessible surface area, ratio of secondary structure, and hydration free energy were also calculated to find out the relationship between the number of LRR modules and stability of the protein. Our results show that the intramolecular interactions lead to more compact structure and smaller surface area of the repebodies, which are critical for the stability of repeat proteins. The other features were also well compatible with the experimental results. Based on our observations, the repebody-5 was proposed as the best structure from the all repebodies in structure optimization process. The present study successfully demonstrated that our computer-based molecular modeling approach can significantly contribute to the experiment-based protein engineering challenge.

Porcine testicular carbonyl reductase, PTCR which is one of the short chain dehydrogenases/reductases (SDR) superfamily catalyzes the NADPH-dependent reduction of carbonyl compounds including steroids and prostaglandins. Previously we reported C-terminal tail of PTCR was deleted due to a nonsynonymous single nucleotide variation (nsSNV). Here we identified from kinetic studies that the enzymatic properties for 5 alpha-dihydrotestosterone (5 alpha-DHT) were different between wild-type and C-terminal-deleted PTCRs. Compared to wild-type PTCR, C-terminal-deleted PTCR has much higher reduction rate. To investigate structural difference between wild-type and C-terminal-deleted PTCRs upon 5a-DHT binding, we performed molecular dynamics simulations for two complexes. Using trajectories, molecular interactions including hydrogen bonding patterns, distance between 5a-DHT and catalytic Tyr193, and interaction energies are analyzed and compared. During the MD simulation time, the dynamic behavior of C-terminal tail in wild-type PTCR is also examined using essential dynamics analysis. The results of our simulations reveal that the binding conformation of 5 alpha-DHT in C-terminal-deleted PTCR is more favorable for reduction reaction in PTCR, which shows strong agreement with kinetic data. These structural findings provide valuable information to understand substrate specificity of PTCR and further kinetic properties of enzymes belonging to the SDR superfamily.

Multiple isoforms of Arabidopsis thaliana h-type thioredoxins (AtTrx-hs) have distinct structural and functional specificities. AtTrx-h3 acts as both a disulfide reductase and as a molecular chaperone. We prepared five representative AtTrx-hs and compared their protein structures and disulfide reductase and molecular chaperone activities. AtTrx-h2 with an N-terminal extension exhibited distinct functional properties with respect to other AtTrx-hs. AtTrx-h2 formed low-molecular-mass structures and exhibited only disulfide reductase activity, whereas the other AtTrx-h isoforms formed high-molecular-mass complexes and displayed both disulfide reductase and molecular chaperone activities. The domains that determine the unique structural and functional properties of each AtTrx-hs protein were determined by constructing a domain-swap between the N- and C-terminal regions of AtTrx-h2 and AtTrx-h3 (designated AtTrx-h-2N3C and AtTrx-h-3N2C respectively), an N-terminal deletion mutant of AtTrx-h2 [AtTrx-h2-N(∆19)] and site-directed mutagenesis of AtTrx-h3. AtTrx-h2-N(∆19) and AtTrx-h-3N2C exhibited similar properties to those of AtTrx-h2, but AtTrx-h-2N3C behaved more like AtTrx-h3, suggesting that the structural and functional specificities of AtTrx-hs are determined by their C-terminal regions. Hydrophobicity profiling and molecular modelling revealed that Ala100 and Ala106 in AtTrx-h3 play critical roles in its structural and functional regulation. When these two residues in AtTrx-h3 were replaced with lysine, AtTrx-h3 functioned like AtTrx-h2. The chaperone function of AtTrx-hs conferred enhanced heat-shock-resistance on a thermosensitive trx1/2-null yeast mutant.

The interaction of Obg (Spo0B-associated GTP-binding protein) GTPase and SpoT, which is a bifunctional ppGpp (guanosine 3',5'-bispyrophosphate) hydrolase/synthetase, is vital for the modulation of intracellular ppGpp levels during bacterial responses to environmental cues. It has been recently reported that the ppGpp level is also inducible by various stresses in the chloroplasts of plant cells. However, the function of the Obg-SpoT interaction in plants remains elusive. The results from the present and previous studies suggest that AtRSH1 is a putative bacterial SpoT homolog in Arabidopsis and that its transcription levels are responsive to wounding and salt stresses. In this study, we used a yeast two-hybrid analysis to map the regions required for the AtObgC-AtRSH1 interaction. Moreover, protein-protein docking simulations revealed reasonable geometric and electrostatic complementarity in the binding surfaces of the two proteins. The data support our experimental results, which suggest that the conserved domains in AtObgC and the N terminus of AtRSH1 containing the TGS domain contribute to their interaction. In addition, quantitative real-time PCR (qRT-PCR) analyses showed that the expression of AtObgC and AtRSH1 exhibit a similar inhibition pattern under wounding and salt-stress conditions, but the inhibition pattern was not greatly influenced by the presence or absence of light. Based on in vivo analyses, we further confirmed that the AtRSH1 and AtObgC proteins similarly localize in chloroplasts. Based on these results, we propose that the AtObgC-AtRSH1 interaction plays a vital role in ppGpp-mediated stress responses in chloroplasts.

Human chymase catalyzes the hydrolysis of peptide bonds. Three chymase inhibitors with very similar chemical structures but highly different inhibitory profiles towards the hydrolase function of chymase were selected with the aim of elucidating the origin of disparities in their biological activities. As a substrate (angiotensin-I) bound crystal structure is not available, molecular docking was performed to dock the substrate into the active site. Molecular dynamics simulations of chymase complexes with inhibitors and substrate were performed to calculate the binding orientation of inhibitors and substrate as well as to characterize conformational changes in the active site. The results elucidate details of the 3D chymase structure as well as the importance of K40 in hydrolase function. Binding mode analysis showed that substitution of a heavier Cl atom at the phenyl ring of most active inhibitor produced a great deal of variation in its orientation causing the phosphinate group to interact strongly with residue K40. Dynamics simulations revealed the conformational variation in region of V36-F41 upon substrate and inhibitor binding induced a shift in the location of K40 thus changing its interactions with them. Chymase complexes with the most active compound and substrate were used for development of a hybrid pharmacophore model which was applied in databases screening. Finally, hits which bound well at the active site, exhibited key interactions and favorable electronic properties were identified as possible inhibitors for chymase. This study not only elucidates inhibitory mechanism of chymase inhibitors but also provides key structural insights which will aid in the rational design of novel potent inhibitors of the enzyme. In general, the strategy applied in the current study could be a promising computational approach and may be generally applicable to drug design for other enzymes.

Activation of the peroxisome proliferator-activated receptor γ (PPARγ) is important for the treatment of type 2 diabetes and obesity through the regulation of glucose metabolism and fatty acid accumulation. Hence, the discovery of novel PPARγ agonists is necessary to overcome these diseases. In this study, a newly developed approach, multi-conformation dynamic pharmacophore modeling (MCDPM), was used for screening candidate compounds that can properly bind PPARγ. Highly populated structures obtained from molecular dynamics (MD) simulations were selected by clustering analysis. Based on these structures, pharmacophore models were generated from the ligand-binding pocket and then validated to check the rationality. Consequently, two hits were retrieved as final candidates by utilizing virtual screening and molecular docking simulations. These compounds can be used in the design of novel PPARγ agonists.

Hypertension is characterized with stress on the heart and blood vessels thus increasing the risk of heart attack and renal diseases. The Renin angiotensin system (RAS) plays a major role in blood pressure control. Renin is the enzyme that controls the RAS at the rate-limiting step. Our aim is to develop new drug-like leads which can inhibit renin and thereby emerge as therapeutics for hypertension. To achieve this, molecular dynamics (MD) simulation and receptor-based pharmacophore modeling were implemented, and three rennin-inhibitor complex structures were selected based on IC50 value and scaffolds of inhibitors. Three pharmacophore models were generated considering conformations induced by inhibitor. The compounds mapped to these models were selected and subjected to drug-like screening. The identified hits were docked into the active site of renin. Finally, hit1 satisfying the binding mode and interaction energy was selected as possible lead candidate to be used in novel renin inhibitors., {"references":["Murray CJL, Lopez AD, \"Alternative projections of mortality and disability by cause 1990-2020: Global Burden of Disease Study,\" The Lancet, vol.349, no.9064, pp.1498-1504, 1997.","Dzau\tV,\t\"The\tcardiovascular\tcontinuum\tand renin-angiotensin-aldosterone system blockade,\" Journal of hypertension,vol. 23, S9, 2005.","Kim S, Iwao H, \"Molecular and cellular mechanisms of angiotensin II-mediated cardiovascular and renal diseases,\" Pharmacological reviews,vol. 52, no. 1, pp. 11, 2000.","Cooper ME, \"The role of the renin-angiotensin-aldosterone system in diabetes and its vascular complications,\" American journal of hypertension,vol. 17, pp. S16-S20, 2004.","Norris K, Vaughn C, \"The role of reninangiotensin aldosterone system inhibition in chronic kidney disease,\" Expert review of cardiovascular therapy,vol. 1, no. 1, pp. 51-63, 2003.","Wood JM, Stanton JL, Hofbauer KG, \"Inhibitors of renin as potential therapeutic agents,\" Journal of enzyme inhibition,vol. 1, no. 3, pp. 169, 1987.","Sielecki AR, Hayakawa K, Fujinaga M, Murphy M, Fraser M, Muir AK, Carilli CT, Lewicki JA, Baxter JD, James M, \"Structure of recombinant human renin, a target for cardiovascular-active drugs, at 2.5 A resolution,\" Science,vol. 243, no. 4896, pp. 1346, 1989.","Rahuel J, Rasetti V, Maibaum J, Rueger H, Goschke R, Cohen N, Stutz S, Cumin F, Fuhrer W,\"Structure-based drug design: the discovery of novel nonpeptide orally active inhibitors of human renin,\"Chemistry & Biology,vol. 7, no. 7, pp. 493-504, 2000.","James MNG, Sielecki AR,\"Stereochemical analysis of peptide bond hydrolysis catalyzed by the aspartic proteinase penicillopepsin,\"Biochemistry,vol. 24, no. 14, pp. 3701-3713, 1985. \n[10]\tStanton A, \"Therapeutic potential of renin inhibitors in the management of cardiovascular disorders,\" American Journal of Cardiovascular Drugs,vol. 3, no. 6, pp. 389-394, 2003.\n[11]\tFisher NDL, Hollenberg NK, \"Is there a future for renin inhibitors?,\" Expert Opinion on Investigational Drugs,vol. 10, no. 3, pp. 417-426, 2001.\n[12]\tWood JM, Maibaum J, Rahuel J, Grutter MG, Cohen NC, Rasetti V, Ruger H, Goschke R, Stutz S, Fuhrer W, \"Structure-based design of aliskiren, a novel orally effective renin inhibitor,\" Biochemical and biophysical research communications,vol. 308, no. 4, pp. 698-705, 2003. \n[13]\tUdenigwe CC, Li H, Aluko RE, \"Quantitative structure-activity relationship modeling of renin-inhibiting dipeptides,\" Amino Acids, pp. 1-8, 2011.\n[14]\tClaude Cohen N, \"Structure-Based Drug Design and the Discovery of Aliskiren (Tektuma): Perseverance and Creativity to Overcome a RD Pipeline Challenge,\" Chemical Biology & Drug Design,vol. 70, no. 6, pp. 557-565, 2007.\n[15]\tBezenc on 0, Bur D, Weller T, Richard-Bildstein S, Remen L, Sifferlen T, Corminboeuf 0, Grisostomi C, Boss C, Prade L, \"Design and Preparation of Potent, Nonpeptidic, Bioavailable Renin Inhibitors,\" Journal of Medicinal Chemistry,vol. 52, no. 12, pp. 3689-3702, 2009. \n[16]\tDeng J, Lee KW, Sanchez T, Cui M, Neamati N, Briggs JM, \"Dynamic receptor-based pharmacophore model development and its application in designing novel HIV-1 integrase inhibitors,\" Journal of Medicinal Chemistry,vol. 48, no. 5, pp. 1496-1505, 2005.\n[17]\tTice CM, Xu Z, Yuan J, Simpson RD, Cacatian ST, Flaherty PT, Zhao W, Guo J, Ishchenko A, Singh SB, \"Design and optimization of renin inhibitors: Orally bioavailable alkyl amines,\" Bioorganic & medicinal chemistry letters,vol. 19, no. 13, pp. 3541-3545, 2009.\n[18]\tVan Der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJC, \"GROMACS: fast, flexible, and free,\" Journal of computational chemistry,vol. 26, no. 16, pp. 1701-1718, 2005.\n[19]\tSchuttelkopf AW, Van Aalten DMF, \"PRODRG: a tool for high-throughput crystallography of protein-ligand complexes,\" ActaCrystallographica Section D: Biological Crystallography,vol. 60, no. 8, pp. 1355-1363, 2004.\n[20]\tEssmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG, \"A smooth particle mesh Ewald method,\" Journal of Chemical Physics,vol. 103, no. 19, pp. 8577-8593, 1995.\n[21]\tHess B, Bekker H, Berendsen HJC, Fraaije JGEM, \"L1NCS: a linear constraint solver for molecular simulations,\" Journal of computational chemistry,vol. 18, no. 12, pp. 1463-1472, 1997.\n[22]\tMiyamoto S, Kollman PA, \"SETTLE: an analytical version of the SHAKE and RATTLE algorithm for rigid water models,\" Journal of computational chemistry,vol. 13, no. 8, pp. 952-962, 1992.\n[23]\tVerdonk ML, Cole JC, Hartshorn MJ, Murray CW, Taylor RD, \"Improved protein-ligand docking using GOLD,\" Proteins: Structure, Function, and Bioinformatics,vol. 52, no. 4, pp. 609-623, 2003."]}

In our study, a structure-based virtual screening study was conducted to identify potent ITK inhibitors, as ITK is considered to play an important role in the treatment of inflammatory diseases. We developed a structure-based pharmacophore model using the crystal structure (PDB ID: 3MJ2) of ITK complexed with BMS-50944. The most predictive model, SB-Hypo1, consisted of six features: three hydrogen-bond acceptors (HBA), one hydrogen-bond donor (HBD), one ring aromatic (RA), and one hydrophobic (HY). The statistical significance of SB-Hypo1 was validated using wide range of test set molecules and a decoy set. The resulting well-validated model could then be confidently used as a 3D query to screen for drug-like molecules in a database, in order to retrieve new chemical scaffolds that may be potent ITK inhibitors. The hits retrieved from this search were filtered based on the maximum fit value, drug-likeness, and ADMET properties, and the hits that were retained were used in a molecular docking study to find the binding mode and molecular interactions with crucial residues at the active site of the protein. These hits were then fed into a molecular dynamics simulation to study the flexibility of the activation loop of ITK upon ligand binding. This combination of methodologies is a valuable tool for identifying structurally diverse molecules with desired biological activities, and for designing new classes of selective ITK inhibitors.

Changes in meat quality traits are strongly associated with alterations in postmortem metabolism which depend on genetic variations, especially nonsynonymous single nucleotide variations (nsSNVs) having critical effects on protein structure and function. To selectively identify metabolism-related nsSNVs, next-generation transcriptome sequencing (RNA-Seq) was carried out using RNAs from porcine liver, which contains a diverse range of metabolic enzymes. The multiplex SNV genotyping analysis showed that various metabolism-related genes had different nsSNV alleles. Moreover, many nsSNVs were significantly associated with multiple meat quality traits. Particularly, ch7:g.22112616A>G SNV was identified to create a single amino acid change (Thr/Ala) at the 145th residue of H1.3-like protein, very close to the putative 147th threonine phosphorylation site, suggesting that the nsSNV may affect multiple meat quality traits by affecting the epigenetic regulation of postmortem metabolism-related gene expression. Besides, one nonsynonymous variation, probably generated by gene duplication, led to a stop signal in porcine testicular carbonyl reductase (PTCR), resulting in a C-terminal (E281-A288) deletion. Molecular docking and energy minimization calculations indicated that the binding affinity of wild-type PTCR to 5α-DHT, a C(21)-steroid, was superior to that of C-terminal-deleted PTCR or human carbonyl reductase, which was very consistent with experimental data, reported previously. Furthermore, P284 was identified as an important residue mediating the specific interaction between PTCR and 5α-DHT, and phylogenetic analysis showed that P284 is an evolutionarily conserved residue among animal carbonyl reductases, which suggests that the C-terminal tails of these reductases may have evolved under evolutionary pressure to increase the substrate specificity for C(21)-steroids and facilitate metabolic adaptation. Altogether, our RNA-Seq revealed that selective nsSNVs were associated with meat quality traits that could be useful for successful marker-assisted selection in pigs and also represents a useful resource to enhance understanding of protein folding, substrate specificity, and the evolution of enzymes such as carbonyl reductase.

2-Cys peroxiredoxins (Prxs) play important roles in the protection of chloroplast proteins from oxidative damage. Arabidopsis NADPH-dependent thioredoxin reductase isotype C (AtNTRC) was identified as efficient electron donor for chloroplastic 2-Cys Prx-A. There are three isotypes (A, B, and C) of thioredoxin reductase (TrxR) in Arabidopsis. AtNTRA contains only TrxR domain, but AtNTRC consists of N-terminal TrxR and C-terminal thioredoxin (Trx) domains. AtNTRC has various oligomer structures, and Trx domain is important for chaperone activity. Our previous experimental study has reported that the hybrid protein (AtNTRA-(Trx-D)), which was a fusion of AtNTRA and Trx domain from AtNTRC, has formed variety of structures and shown strong chaperone activity. But, electron transfer mechanism was not detected at all. To find out the reason of this problem with structural basis, we performed two different molecular dynamics (MD) simulations on AtNTRC and AtNTRA-(Trx-D) proteins with same cofactors such as NADPH and flavin adenine dinucleotide (FAD) for 50 ns. Structural difference has found from superimposition of two structures that were taken relatively close to average structure. The main reason that AtNTRA-(Trx-D) cannot transfer the electron from TrxR domain to Trx domain is due to the difference of key catalytic residues in active site. The long distance between TrxR C153 and disulfide bond of Trx C387-C390 has been observed in AtNTRA-(Trx-D) because of following reasons: i) unstable and unfavorable interaction of the linker region, ii) shifted Trx domain, and iii) different or weak interface interaction of Trx domains. This study is one of the good examples for understanding the relationship between structure formation and reaction activity in hybrid protein. In addition, this study would be helpful for further study on the mechanism of electron transfer reaction in NADPH-dependent thioredoxin reductase proteins.

Sonic hedgehog (Shh) plays an important role in the activation of Shh signaling pathway that regulates preservation and rebirth of adult tissues. An abnormal activation of this pathway has been identified in hyperplasia and various tumorigenesis. Hence the inhibition of this pathway using a Shh inhibitor might be an efficient way to treat a wide range of malignancies. This study was done in order to develop a lead chemical candidate that has an inhibitory function in the Shh signaling pathway. We have generated common feature pharmacophore models using three-dimensional (3D) structural information of robotnikinin, an inhibitor of the Shh signaling pathway, and its analogs. These models have been validated with fit values of robotnikinin and its analogs, and the best model was used as a 3D structural query to screen chemical databases. The hit compounds resulted from the screening docked into a proposed binding site of the Shh named pseudo-active site. Molecular dynamics (MD) simulations were performed to investigate detailed binding modes and molecular interactions between the hit compounds and functional residues of the pseudo-active site. The results of the MD simulation analyses revealed that the hit compounds can bind the pseudo-active site with high affinity than robotnikinin. As a result of this study, a candidate inhibitor (GK 03795) was selected as a potential lead to be employed in future Shh inhibitor design.

Competitive glycosidase inhibitors are generally sugar mimics that are costly and tedious to obtain because they require challenging and elongated chemical synthesis, which must be stereo- and regiocontrolled. Here, we show that readily accessible achiral (E)-1-phenyl-3-(4-strylphenyl)ureas are potent competitive α-glucosidase inhibitors. A systematic synthesis study shows that the 1-phenyl moiety on the urea is critical for ensuring competitive inhibition, and substituents on both terminal phenyl groups contribute to inhibition potency. The most potent inhibitor, compound 12 (IC(50)=8.4 μM, K(i)=3.2 μM), manifested a simple slow-binding inhibition profile for α-glucosidase with the kinetic parameters k(3)=0.005256 μM(-1) min(-1), k(4)=0.003024 min(-1), and K(i)(app) =0.5753 μM.

Bacterial Hfq is a highly conserved thermostable protein of about 10 kDa. The Hfq protein was discovered in 1968 as an E. coli host factor that was essential for replication of the bacteriophage Q. It is now clear that Hfq has many important physiological roles. In E. coli, Hfq mutants show a multiple stress response related phenotypes. Hfq is now known to regulate the translation of two major stress transcription factors RpoS and RpoE in Enterobacteria and mediates its plieotrophic effects through several mechanisms. It interacts with regulatory sRNA and facilitates their antisense interaction with their targets. It also acts independently to modulate mRNA decay and in addition acts as a repressor of mRNA translation. Recent paper from Arluison et al. [9] provided the first evidence indicating that Hfq is an ATP-binding protein. They determined a plausible ATP-binding site in Hfq and tested Hfq’s ATP-binding affinity and stoichiometry. Experimental data suggest that the ATP-binding by the Hfq–RNA complex results in its significant destabilization of the protein and the result also proves important role of Tyr25 that flanks the cleft and stabilizes the adenine portion of ATP, possibly via aromatic stacking. In our study, the ATP molecule was docked into the predicted binding cleft using GOLD docking software. The binding nature of ATP and its effect on Hfq–RNA complex was studied using molecular dynamics simulations. Importance of Tyr25 residue was monitored and revealed using mutational study on the modeled systems. Our data and the corresponding results point to one of Hfq functional structural consequences due to ATP binding and Tyr25Ala mutation.

The bacterial Obg proteins (Spo0B-associated GTP-binding protein) belong to the subfamily of P-loop GTPase proteins that contain two equally and highly conserved domains, a C-terminal GTP binding domain and an N-terminal glycine-rich domain which is referred as the "Obg fold" and now it is considered as one of the new targets for antibacterial drug. When the Obg protein is associated with GTP, it becomes activated, because conformation of Obg fold changes due to the structural changes of GTPase switch elements in GTP binding site. In order to investigate the effects and structural changes in GTP bound to Obg and GTPase switch elements for activation, four different molecular dynamics (MD) simulations were performed with/without the three different nucleotides (GTP, GDP, and GDP + Pi) using the Bacillus subtilis Obg (BsObg) structure. The protein structures generated from the four different systems were compared using their representative structures. The pattern of C(alpha)-C(alpha) distance plot and angle between the two Obg fold domains of simulated apo form and each system (GTP, GDP, and GDP+Pi) were significantly different in the GTP-bound system from the others. The switch 2 element was significantly changed in GTP-bound system. Also root-mean-square fluctuation (RMSF) analysis revealed that the flexibility of the switch 2 element region was much higher than the others. This was caused by the characteristic binding mode of the nucleotides. When GTP was bound to Obg, its gamma-phosphate oxygen was found to interact with the key residue (D212) of the switch 2 element, on the contrary there was no such interaction found in other systems. Based on the results, we were able to predict the possible binding conformation of the activated form of Obg with L13, which is essential for the assembly with ribosome.

Hfq is an abundant RNA-binding bacterial protein that was first identified in E. coli as a required host factor for phage Qbeta RNA replication. The pleiotrophic phenotype resulting from the deletion of Hfq predicates the importance of this protein. Two RNA-binding sites have been characterized: the proximal site which binds sRNA and mRNA and the distal site which binds poly(A) tails. Previous studies mainly focused on the key residues in the proximal site of the protein. A recent mutation study in E. coli Hfq showed that a distal residue Val43 is important for the protein function. Interestingly, when we analyzed the sequence and structure of Staphylococcus aureus Hfq using the CONSEQ server, the results elicited that more functional residues were located far from the nucleotide-binding portion (NBP). From the analysis seven individual residues Asp9, Leu12, Glu13, Lys16, Gln31, Gly34 and Asp40 were selected to investigate the conformational changes in Hfq-RNA complex due to point mutation effect of those residues using molecular dynamics simulations. Results showed a significant effect on Asn28 which is an already known highly conserved functionally important residue. Mutants D9A, E13A and K16A depicted effects on base stacking along with increase in RNA pore diameter, which is required for the threading of RNA through the pore for the post-translational modification. Further, the result of protein stability analysis by the CUPSAT server showed destabilizing effect in the most mutants. From this study we characterized a series of important residues located far from the NBP and provide some clues that those residues may affect sRNA binding in Hfq.

Stilbene urea derivatives as a novel and competitive class of non-glycosidic α-glucosidase inhibitors are effective for the treatment of type II diabetes and obesity. The main purposes of our molecular modeling study are to explore the most suitable binding poses of stilbene derivatives with analyzing the binding affinity differences and finally to develop a pharmacophore model which would represents critical features responsible for α-glucosidase inhibitory activity. Three-dimensional structure of S. cerevisiae α-glucosidase was built by homology modeling method and the structure was used for the molecular docking study to find out the initial binding mode of compound 12, which is the most highly active one. The initial structure was subjected to molecular dynamics (MD) simulations for protein structure adjustment at compound 12-bound state. Based on the adjusted conformation, the more reasonable binding modes of the stilbene urea derivatives were obtained from molecular docking and MD simulations. The binding mode of the derivatives was validated by correlation analysis between experimental Ki value and interaction energy. Our results revealed that the binding modes of the potent inhibitors were engaged with important hydrogen bond, hydrophobic, and π-interactions. With the validated compound 12-bound structure obtained from combining approach of docking and MD simulation, a proper four featured pharmacophore model was generated. It was also validated by comparison of fit values with the Ki values. Thus, these results will be helpful for understanding the relationship between binding mode and bioactivity and for designing better inhibitors from stilbene derivatives.

Histone deacetylases (HDACs) have emerged as effective therapeutic targets in the treatment of various diseases including cancers as these enzymes directly involved in the epigenetic regulation of genes. However the development of isoform-selective HDAC inhibitors has been a challenge till date since all HDAC enzymes possess conserved tunnel-like active site. In this study, using molecular dynamics simulation we have analyzed the behavior of tunnels present in HDAC8, 10, and 11 enzymes of class I, II, and IV, respectively. We have identified the equivalent tunnel forming amino acids in these three isoforms and found that they are very much conserved with subtle differences to be utilized in selective inhibitor development. One amino acid, methionine of HDAC8, among six tunnel forming residues is different in isoforms of other classes (glutamic acid (E) in HDAC10 and leucine (L) in HDAC 11) based on which mutations were introduced in HDAC11, the less studied HDAC isoform, to observe the effects of this change. The HDAC8-like (L268M) mutation in the tunnel forming residues has almost maintained the deep and narrow tunnel as present in HDAC8 whereas HDAC10-like (L268E) mutation has changed the tunnel wider and shallow as observed in HDAC10. These results explained the importance of the single change in the tunnel formation in different isoforms. The observations from this study can be utilized in the development of isoform-selective HDAC inhibitors.

TWIK (tandem-pore domain weak inward rectifying K(+))-related spinal cord K(+) channel, TRESK, a member of the tandem-pore domain K(+) channel family, is the most recently cloned K(2P) channel. TRESK is highly expressed in dorsal root ganglion neuron, a pain sensing neuron, which is a target for analgesics. In this study, a reliable 3D structure for transmembrane (TM) region of mouse TRESK (mTRESK) was constructed, and then the reasonable blocker binding mode of the protein was investigated. The 3D structure of the mTRESK built by homology modeling method was validated with recommend value of stereochemical quality. Based on the validated structure, K(+) channel blocker-bound conformation was obtained by molecular docking and 5ns MD simulation with DPPC lipid bilayer. Our docking study provides the plausible binding mode of known blockers with key interacting residues, especially, F156 and F364. Finally, these modeling results were verified by experimental study with mutation from phenylalanine to alanine (F156A, F364A and F156A/F364A) at the TM2 and TM4. This is the first modeling study for TRESK that can provide structural information of the protein including ligand binding information. These results can be useful in structure based drug design for finding new blockers of the TRESK as potential therapeutic target of pain treatment.