Yuno Lee

RNA-Seq approach for genetic improvement of meat quality in pig and evolutionary insight into the substrate specificity of animal carbonyl reductases.

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 C21-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 C21-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.

Synthesis, characterization and catalytic properties of TS-1 monoliths

Titanium silicalite-1 monolith (TS-1M) was prepared using polyurethane foam as the template. The characterization of the TS-1M was carried out using XRD, SEM, UV–vis and FT-IR spectroscopies, and water vapor adsorption. Catalytic activity was measured for 1-hexene and 2,5-dihydrofuran epoxidation using H2O2 as an oxidant. These studies revealed that TS-1 monoliths have essentially identical chemical properties to the TS-1 powders, but the crystal morphology was different and diffusion plays an important role in the catalysis of the liquid-phase epoxidation reactions.

AtObgC, a plant ortholog of bacterial Obg, is a chloroplast-targeting GTPase essential for early embryogenesis.

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.

Probing possible egress channels for multiple ligands in human CYP3A4: A molecular modeling study

Human cytochrome P450 (CYP) 3A4 extensively contributes to metabolize 50% of the marketed drugs. Recently, a CYP3A4 structure with two molecules of ketoconazole (2KT) was identified. However, channels for egresses of these inhibitors are unexplored. Thus, we applied molecular dynamics simulations followed by channel analyses. Two simulations of empty and 2KT-bound CYP3A4 results revealed the multiple ligand-induced conformational changes in channel forming regions, which appear to be important for the regulation of channels. In addition, we observed that the channel-3 entrance is closed due to the large structural deviation of the key residues from Phe-cluster. F215 and F220 are known as entrance blockers of channel-2 in metyrapone-bound CYP3A4. Currently, F220 blocks the channel-3 along with F213 and F241. Therefore, it suggested that channel-1 and 2 could potentially serve as egress routes for 2KT. It is also supported by the results from MOLAxis analyses, in which the frequency of channel occurrence and bottleneck radius during simulation favor channel-1 and 2. Several bottleneck residues of these channels may have critical roles in 2KT egresses, especially S119. Our modeling study for multiple ligand-channeling of CYP3A4 could be very helpful to gain new insights into channel selectivity of CYP3A4.

Dynamic Structure-Based Pharmacophore Model Development: A New and Effective Addition in the Histone Deacetylase 8 (HDAC8) Inhibitor Discovery

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.

Binding Mode Analyses and Pharmacophore Model Development for Stilbene Derivatives as a Novel and Competitive Class of α-Glucosidase Inhibitors

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.

Molecular Modeling Study on Tunnel Behavior in Different Histone Deacetylase Isoforms

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.

DISCOVERY AND EVALUATION OF POTENTIAL SONIC HEDGEHOG SIGNALING PATHWAY INHIBITORS USING PHARMACOPHORE MODELING AND MOLECULAR DYNAMICS SIMULATIONS

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.

A Novel Competitive Class of α‐Glucosidase Inhibitors: (E)‐1‐Phenyl‐3‐(4‐Styrylphenyl)Urea Derivatives

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 (IC50=8.4 μM, Ki=3.2 μM), manifested a simple slow‐binding inhibition profile for α‐glucosidase with the kinetic parameters k3=0.005256 μM−1 min−1, k4=0.003024 min−1, and

Novel chemical scaffolds of the tumor marker AKR1B10 inhibitors discovered by 3D QSAR pharmacophore modeling

{K{{{\rm app}\hfill \atop {\rm i}\hfill}}}