Search

Modulation of γ-aminobutyric acid (GABA) levels has been required in various disorders. GABA itself cannot be directly introduced into central nervous system (CNS) because of the blood brain barrier; inhibition of GABA aminotransferase (GABA-AT), which degrades GABA in CNS, has been the target for the modulation of GABA levels in CNS. Given that root extract of valerian (Valeriana officinalis) has been used for millennia as anti-anxiolytic and sedative, in silico approach was carried out to investigate valerian compounds exhibiting GABA-AT inhibiting activity. The 3D structure of human GABA-AT was created from pig crystal structure via homology modeling. Inhibition of GABA-AT by 18 valerian compounds was analyzed using molecular docking and molecular dynamics simulations and compared with known GABA-AT inhibitors such as vigabatrin and valproic acid. Isovaleric acid and didrovaltrate exhibited GABA-AT inhibiting activity in computational analysis, albeit less potent compared with vigabatrin. However, multiple compounds with low activity may have additive effects when the total extract of valeriana root was used in traditional usage. In addition, isovaleric acid shares similar backbone structure to GABA, suggesting that isovaleric acid might be a valuable starting structure for the development of more efficient GABA-AT inhibitors for disorders related with low level of GABA in the CNS.

Peroxiredoxins (Prx) have received considerable attention during recent years. This study demonstrates that two typical Pseudomonas-derived 2-Cys Prx proteins, PpPrx and PaPrx can alternatively function as a peroxidase and chaperone. The amino acid sequences of these two Prx proteins exhibit 93% homology, but PpPrx possesses an additional cysteine residue, Cys112, instead of the alanine found in PaPrx. PpPrx predominates with a high molecular weight (HMW) complex and chaperone activity, whereas PaPrx has mainly low molecular weight (LMW) structures and peroxidase activity. Mass spectrometry and structural analyses showed the involvement of Cys112 in the formation of an inter-disulfide bond, the instability of LMW structures, the formation of HMW complexes, and increased hydrophobicity leading to functional switching of Prx proteins between peroxidase and chaperone.

In contrast to a number of studies on the humanization of non-human antibodies, the reshaping of a non-human antibody into a chicken antibody has never been attempted. Therefore, nothing is known about the animal species-dependent compatibility of the framework regions (FRs) that sustain the appropriate conformation of the complementarity-determining regions (CDRs). In this study, we attempted the reshaping of the variable domains of the mouse catalytic anti-nucleic acid antibody 3D8 (m3D8) into the FRs of a chicken antibody (“chickenization”) by CDR grafting, which is a common method for the humanization of antibodies. CDRs of the acceptor chicken antibody that showed a high homology to the FRs of m3D8 were replaced with those of m3D8, resulting in the chickenized antibody (ck3D8). ck3D8 retained the biochemical properties (DNA binding, DNA hydrolysis, and cellular internalizing activities) and three-dimensional structure of m3D8 and showed reduced immunogenicity in chickens. Our study demonstrates that CDR grafting can be applied to the chickenization of a mouse antibody, probably due to the interspecies compatibility of the FRs.

The modular auxiliary activity (AA) family of proteins is believed to cause amorphogenesis in addition to oxidative cleavage of crystalline cellulose although the supporting evidence is limited. HcAA10-2 is a modular AA10 family protein (58 kDa) composed of a AA10 module and a family two carbohydrate binding module (CBM2), joined by a long stretch of 222 amino acids of unknown function. The protein was expressed in Escherichia coli and purified to homogeneity. Scanning electron microscopy and X-ray diffraction analysis of Avicel treated with HcAA10-2 provided evidence for the disruption of the cellulose microfibrils ("amorphogenesis") and reduction of the crystallinity index, resulting in a twofold increase of cellulase adsorption on the polysaccharide surface. HcAA10-2 exhibited weak endoglucanase-like activity toward soluble cellulose and cello-oligosaccharides with an optimum at pH 6.5 and 45 °C. HcAA10-2 catalyzed oxidative cleavage of crystalline cellulose released native and oxidized cello-oligosaccharides in the presence of copper and an electron donor such as ascorbic acid. Multiple sequence alignment indicated that His1, His109, and Phe197 in the AA10 module formed the conserved copper-binding site. The reducing sugar released from Avicel by the endoglucanase Cel5 and Celluclast accompanying HcAA10-2 was increased by four- and sixfold, respectively. Moreover, HcAA10-2 and Celluclast acted synergistically on pretreated wheat straw biomass resulting in a threefold increase in reducing sugar than Celluclast alone. Taken together, these results suggest that HcAA10-2 is a novel multifunctional modular AA10 protein possessing amorphogenesis, weak endoglucanase, and oxidative cleavage activities useful for efficient degradation of crystalline cellulose.

As previously reported, the activity of the large-conductance calcium (Ca(2+))-activated potassium (K(+)) (BK(Ca)) channel is strongly potentiated from the extracellular side of the cell membrane by certain benzofuroindole derivatives. Here, the mechanism of action of one of the most potent activators, 4-chloro-7-(trifluoromethyl)-10H-benzofuro[3,2-b]indole-1-carboxylic acid (CTBIC), is characterized. This compound, Compound 22 in the previous report (Chembiochem 6:1745-1748, 2005), potentiated the activity of the channel by shifting its conductance-voltage relationship toward the more negative direction. Cotreatment with CTBIC reduced the affinity of charybdotoxin, a peptide pore-blocker, whereas that of tetraethylammonium, a small pore-blocking quaternary ammonium, was not significantly altered. Guided by these results, scanning mutagenesis of the outer vestibule of the BK(Ca) channel was launched to uncover the molecular determinants that affect CTBIC binding. Alanine substitution of several amino acid residues in the turret region and the S6 helix of the channel decreased potentiation by CTBIC. Homology modeling and molecular dynamics simulation showed that some of these residues formed a CTBIC binding pocket between two adjacent α-subunits in the outer vestibule of the channel. Thus, it can be envisioned that benzofuroindole derivatives stabilize the open conformation of the channel by binding to the residues clustered across the extracellular part of the subunit interface. The present results indicate that the interface between different α-subunits of the BK(Ca) channel may play a critical role in the modulation of channel activity. Therefore, this interface represents a potential therapeutic target site for the regulation of K(+) channels.

Histone deacetylase 8 (HDAC8) is an enzyme involved in deacetylating the amino groups of terminal lysine residues, thereby repressing the transcription of various genes including tumor suppressor gene. The over expression of HDAC8 was observed in many cancers and thus inhibition of this enzyme has emerged as an efficient cancer therapeutic strategy. In an effort to facilitate the future discovery of HDAC8 inhibitors, we developed two pharmacophore models containing six and five pharmacophoric features, respectively, using the representative structures from two molecular dynamic (MD) simulations performed in Gromacs 4.0.5 package. Various analyses of trajectories obtained from MD simulations have displayed the changes upon inhibitor binding. Thus utilization of the dynamically-responded protein structures in pharmacophore development has the added advantage of considering the conformational flexibility of protein. The MD trajectories were clustered based on single-linkage method and representative structures were taken to be used in the pharmacophore model development. Active site complimenting structure-based pharmacophore models were developed using Discovery Studio 2.5 program and validated using a dataset of known HDAC8 inhibitors. Virtual screening of chemical database coupled with drug-like filter has identified drug-like hit compounds that match the pharmacophore models. Molecular docking of these hits reduced the false positives and identified two potential compounds to be used in future HDAC8 inhibitor design.

The acetylcholinesterase (AChE) is important to terminate acetylcholine-mediated neurotransmission at cholinergic synapses. The pivotal role of AChE in apoptosome formation through the interactions with cytochrome c (Cyt c) was demonstrated in recent study. In order to investigate the proper binding conformation between the human AChE (hAChE) and human Cyt c (hCyt c), macro-molecular docking simulation was performed using DOT 2.0 program. The hCyt c was bound to peripheral anionic site (PAS) on hAChE and binding mode of the docked conformation was very similar to the reported crystal structure of the AChE and fasciculin-II (Fas-II) complex. Two 10 ns molecular dynamics (MD) simulations were carried out to refine the binding mode of docked structure and to observe the differences of the binding conformations between the absent (Apo) and presence (Holo) of heme group. The key hydrogen bonding residues between hAChE and hCyt c proteins were found in Apo and Holo systems, as well as each Tyr341 and Trp286 residue of hAChE was participated in cation-pi (π) interactions with Lys79 of hCyt c in Apo and Holo systems, respectively. From the present study, although the final structures of the Apo and Holo systems have similar binding pattern, several differences were investigated in flexibilities, interface interactions, and interface accessible surface areas. Based on these results, we were able to predict the reasonable binding conformation which is indispensable for apoptosome formation.

Obg is a ribosome-associated GTPase essential for bacterial viability and is conserved in most organisms, from bacteria to eukaryotes. Obg is also expressed in plants, which predicts an important role for this molecule in plant viability; however, the functions of the plant Obg homologs have not been reported. Here, we first identified Arabidopsis AtObgC as a plant chloroplast-targeting Obg and elucidated its molecular biological and physiological properties. AtObgC encodes a plant-specific Obg GTPase that contains an N-terminal region for chloroplast targeting and has intrinsic GTP hydrolysis activity. A targeting assay using a few AtObgC N-terminal truncation mutants revealed that AtObgC localizes to chloroplasts and its transit peptide consists of more than 50 amino acid residues. Interestingly, GFP-fused full-length AtObgC exhibited a punctate staining pattern in chloroplasts of Arabidopsis protoplasts, which suggests a dimerization or multimerization of AtObgC. Moreover, its Obg fold was indispensable for the generation of the punctate staining pattern, and thus, was supposed to be important for such oligomerization of AtObgC by mediating the protein-protein interaction. In addition, the T-DNA insertion AtObgC null mutant exhibited an embryonic lethal phenotype that disturbed the early stage of embryogenesis. Altogether, our results provide a significant implication that AtObgC as a chloroplast targeting GTPase plays an important role at the early embryogenesis by exerting its function in chloroplast protein synthesis.

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.

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.

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.

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.

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.