Biomedical subjects

Hsuan-Liang Liu

Publications and source records attributed to Hsuan-Liang Liu.

18 recordsLinked to original sources

PubMed2005

Improving disulfide connectivity prediction with sequential distance between oxidized cysteines.

SUMMARY: Predicting disulfide connectivity precisely helps towards the solution of protein structure prediction. In this study, a descriptor derived from the sequential distance between oxidized cysteines (denoted as DOC) is proposed. An approach using support vector machine (SVM) method based on weighted graph matching was further developed to predict the disulfide connectivity pattern in proteins. When DOC was applied, prediction accuracy of 63% for our SVM models could be achieved, which is significantly higher than those obtained from previous approaches. The results show that using the non-local descriptor DOC coupled with local sequence profiles significantly improves the prediction accuracy. These improvements demonstrate that DOC, with a proper scaling scheme, is an effective feature for the prediction of disulfide connectivity. The method developed in this work is available at the web server PreCys (prediction of cys-cys linkages of proteins).

Chymotrypsinogen↗

PubMed2005

Recent developments in structural proteomics for protein structure determination.

The major challenges in structural proteomics include identifying all the proteins on the genome-wide scale, determining their structure-function relationships, and outlining the precise three-dimensional structures of the proteins. Protein structures are typically determined by experimental approaches such as X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. However, the knowledge of three-dimensional space by these techniques is still limited. Thus, computational methods such as comparative and de novo approaches and molecular dynamic simulations are intensively used as alternative tools to predict the three-dimensional structures and dynamic behavior of proteins. This review summarizes recent developments in structural proteomics for protein structure determination; including instrumental methods such as X-ray crystallography and NMR spectroscopy, and computational methods such as comparative and de novo structure prediction and molecular dynamics simulations.

Aspergillus↗

PubMed2005

Homology models of the tetramerization domain of six eukaryotic voltage-gated potassium channels Kv1.1-Kv1.6.

The homology models of the tetramerization (T1) domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were constructed based on the crystal structure of the Shaker T1 domain. The results of amino acid sequence alignment indicate that the T1 domains of these K+ channels are highly conserved, with the similarities varying from 77% between Shaker and Kv1.6 to 93% between Kv1.2 and Kv1.3. The homology models reveal that the T1 domains of these Kv channels exhibit similar folds as those of Shaker K+ channel. These models also show that each T1 monomer consists of three distinct layers, with N-terminal layer 1 and C-terminal layer 3 facing the cytoplasm and the membrane, respectively. Layer 2 exhibits the highest structural conservation because it is located around the central hydrophobic core. For each Kv channel, four identical subunits assemble into the homotetramer architecture around a four-fold axis through the hydrogen bonds and salt bridges formed by 15 highly conserved polar residues. The narrowest opening of the pore is formed by the four conserved residues corresponding to R115 of the Shaker T1 domain. The homology models of these Kv T1 domains provide particularly attractive targets for further structure-based studies.

Amino Acid Sequence↗

PubMed2004

Cysteine separations profiles on protein sequences infer disulfide connectivity.

MOTIVATION: Disulfide bonds play an important role in protein folding. A precise prediction of disulfide connectivity can strongly reduce the conformational search space and increase the accuracy in protein structure prediction. Conventional disulfide connectivity predictions use sequence information, and prediction accuracy is limited. Here, by using an alternative scheme with global information for disulfide connectivity prediction, higher performance is obtained with respect to other approaches. RESULT: Cysteine separation profiles have been used to predict the disulfide connectivity of proteins. The separations among oxidized cysteine residues on a protein sequence have been encoded into vectors named cysteine separation profiles (CSPs). Through comparisons of their CSPs, the disulfide connectivity of a test protein is inferred from a non-redundant template set. For non-redundant proteins in SwissProt 39 (SP39) sharing less than 30% sequence identity, the prediction accuracy of a fourfold cross-validation is 49%. The prediction accuracy of disulfide connectivity for proteins in SwissProt 43 (SP43) is even higher (53%). The relationship between the similarity of CSPs and the prediction accuracy is also discussed. The method proposed in this work is relatively simple and can generate higher accuracies compared to conventional methods. It may be also combined with other algorithms for further improvements in protein structure prediction. AVAILABILITY: The program and datasets are available from the authors upon request. CONTACT: cykao@csie.ntu.edu.tw.

Algorithms↗

PubMed2004

A set of homology models of pore loop domain of six eukaryotic voltage-gated potassium channels Kv1.1-Kv1.6.

Homology models of the pore loop domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were generated based on the crystallographic structure of KcsA. The results of amino acid sequence alignment indicate that these Kv channels are composed of two structurally and functionally independent domains: the N-terminal 'voltage sensor' domain and the C-terminal 'pore loop' domain. The homology models reveal that the pore loop domains of these Kv channels exhibit similar folds to those of KcsA. The structural features and specific packing of aromatic residues around the selectivity filter of these Kv channels are nearly identical to those of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. The distribution of polar and nonpolar side chains on the surfaces of the KcsA and Kv channels reveals that they exhibit a segregation of side chains common to most integral membrane proteins. As the hydrogen bond between Glu71 and Asp80 in KcsA plays an important role in stabilizing the channel, the substituted Val residue in the Kv family corresponding to Glu71 of KcsA stabilizes the channel by making hydrophobic contact with Tyr residue from the signature sequence of the selectivity filter. The homology models of these Kv channels provide particularly attractive subjects for further structure-based studies.

Amino Acid Sequence↗

PubMed2004

Molecular docking of the scorpion toxin Tc1 to the structural model of the voltage-gated potassium channel Kv1.1 from human Homo sapiens.

In this study, structural model of the pore loop region of the voltage-gated potassium channel Kv1.1 from human Homo sapiens was constructed based on the crystallographic structure of KcsA by structural homology. The pore loop region of Kv1.1 exhibits similar folds as that of KcsA. The structural feature of the selectivity filter of Kv1.1 is nearly identical to that of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. Molecular docking experiments of the scorpion toxin Tc1 from Tityus cambridgei to the outer vestibule of KcsA as well as Kv1.1 were subsequently performed with various initial Tc1 orientations. Tc1 was found to form the most stable complexes with these two K+ channels when the side chain of Lys14 occupies the pore of the selectivity filter through electrostatic interaction. Tc1 binds preferentially towards Kv1.1 than KcsA due to stronger hydrophobic and electrostatic interactions formed between the toxin and the selectivity filter and outer vestibule of Kv1.1. Furthermore, surface complementarity of the outer vestibules of the channels to the Tc1 spatial conformations also plays an important role in stabilizing both the Tc1/KcsA and Tc1/Kv1.1 complexes.

Amino Acid Sequence↗

PubMed2004

The stabilizing effects of O-glycosylation on the secondary structural integrity of the designed alpha-loop-alpha motif by molecular dynamics simulations.

In this study, various 400 ps molecular dynamics simulations were conducted to determine the stabilizing effect of O-glycosylation on the secondary structural integrity of the design alpha-loop-alpha motif, which has the optimal loop length of 7 Gly residues (denoted as N-A16G7A16-C). In general, O-glycosylation stabilizes the structural integrity of the model peptide regardless of the length and position of glycosylation sites because it decreases the opportunity for water molecules to compete for the intramolecular hydrogen bonds. The designed peptide exhibits the highest helicity when residues 11 and 31 are replaced with Ser residues followed by O-linked with 3 galactose residues, representing the "face-to-face" glycosylation near the loop. In this case, the loop exhibits an extended conformation and several new hydrogen bonds are observed between the main chain of the loop and the galactose residues, resulting in decreasing the fluctuation and increasing the stability of the entire peptide. When the glycosylation are made close to the loop, the secondary structural integrity of the alpha-loop-alpha motif increases with the number of galactose residues. In addition, "face-to-face" glycosylation increases the structural integrity of this motif to a greater extent than "back-to-back" glycosylation. However, when the glycosylation are created away from the loop and near the N- and C-termini, no general rule is found for the stabilizing effect.

Amino Acid Motifs↗

PubMed2003

The effects of metabolites from the indigenous Acidithiobacillus thiooxidans and temperature on the bioleaching of cadmium from soil.

The effect of metabolites from the indigenous Acidithiobacillus thiooxidans and temperature on the bioleaching of cadmium from soil was investigated in the present study. Bioleaching was found to be more effective than chemical leaching of cadmium. The metabolite, mainly sulfuric acid, which was shown to be growth-associated in the exponential phase, plays a major role in bioleaching. The maximum amount of cadmium leached was obtained after 8 days of precultivation when cells were directly involved in the leaching process. It indicates that cells in the exponential growth phase exhibit higher activity toward bioleaching. In contrast, the maximum amount of cadmium leached and the maximum initial rate for bioleaching were reached after 16 days of precultivation when only metabolites were involved in the bioleaching process. It implies that higher sulfuric acid concentration results in higher leaching efficiency. In addition, higher temperature leads to higher leaching efficiency. The optimal operation condition for bioleaching was determined to be a two-stage process: The first stage involves the precultivation of the indigenous A. thiooxidans at 30 degrees C for 8 days followed by 20 minutes of centrifugation to discard cells. The second stage involves the bioleaching with the subsequent supernatant at 50 degrees C.

Acidithiobacillus thiooxidans↗

PubMed2003

Molecular dynamics simulations of the unfolding mechanism of the catalytic domain from Aspergillus awamori var. X100 glucoamylase.

In this study, 200 ps molecular dynamics simulations were conducted to investigate the unfolding mechanism of the catalytic domain of glucoamylase from Aspergillus awamori var. X100. The unfolding of this domain was suggested to follow a putative hierarchical manner, in which the heavily O-glycosylated belt region from residues T440 to A471 acted as the initiation site, followed by the alpha-helix secondary structure destruction, and then the collapse of the catalytic center pocket. The O-glycosylated belt region surrounded the surface of the catalytic domain in its native state at low temperature, whereas it was extended and is more suitable to be classified as part of the subsequent linker domain at high temperatures due to its high flexibility. The inner set helices of the (alpha/alpha)(6)-barrel seemed to exhibit higher helical content than the outer set ones at all temperatures examined. The distances between the C(alpha) of the three Cys residue pairs fluctuated rapidly at higher temperatures, indicating that these disulfide bonds have little effect on the structural stabilization. The melting temperature, at which the residual total helicity of the catalytic domain is 50%, is much lower than the critical temperature, at which the catalytic center pocket has lost its structural integrity.

Aspergillus↗

PubMed2003

Molecular dynamics simulations of the unfolding of the starch binding domain from Aspergillus niger glucoamylase.

The 600 ps molecular dynamics simulations to investigate the unfolding of the starch binding domain from Aspergillus niger glucoamylase were conducted in vacuum as well as in an external field with the dielectric constant of 80 with temperature jump technique. Electrostatic interactions play an important role in determining the stability of the beta-strands in this domain. The starch binding site 1 is less stable than site 2 since it is more exposed to the surface. The disulfide bond between C509 and C604 is unstable since these two residues are located near the flexible linker domain and in the mobile loop region between beta-strands 6 and 7, respectively. The melting temperature, at which the total residual beta-strand content is 50% that of the solution structure, is about 544K for the simulations with dielectric constant of 80, leading to the estimated unfolding timescale of 0.48 ms in vitro. In addition, the unfolding of the starch binding domain is proposed to initiate from the interior region by the lost of the integrity of the secondary structure.

Animals↗

PubMed2003

Molecular dynamics simulations to determine the optimal loop length in the helix-loop-helix motif.

In this study, 100 ps molecular dynamics simulations were conducted to determine the optimal length of the Gly loop in the helix-loop-helix motif. The results revealed that the length of the loop significantly affected the stability of the two alpha-helices. In addition, the optimal loop length to maintain the highest helicity of the target peptide was found to be the one corresponding to 7 Gly residues. Longer loop length resulted in structural destabilization due to the entropic effect, while shorter loop length was not suitable to achieve intrachain packing of the two alpha-helices.

Amino Acid Sequence↗

PubMed2003

Molecular simulations to determine the chelating mechanisms of various metal ions to the His-tag motif: a preliminary study.

In the present study, molecular simulations were performed to investigate the chelating mechanisms of various metal ions to the His-tag motifs with various His residues. The chelation mostly involved the i and i+2 His residues for Ni(2+), Zn(2+), Cu(2+), and Co(2+), while the cooperation of 3 His residues was necessary when Fe(3+) was involved in chelation with His-tags having more than 4 His residues. Metal ion was best fitted into the pocket formed by the imidazole nitrogens while it was about equally located among these nitrogen atoms. His-tag6 was found to have little effect on the structural integrity while the target protein contains more than 68 amino acid residues. Ni(2+) interacted with the imidazole nitrogen of His3 in the beginning of chelation, and then entered into the pocket formed by His3 and His5 at 4 ns during the 10 ns molecular dynamics simulations. The fast chelating process resulted in successful application of IMAC techniques in efficient protein purification.

Amino Acid Motifs↗

PubMed2003

Protein engineering to improve the thermostability of glucoamylase from Aspergillus awamori based on molecular dynamics simulations.

Twelve mutations were constructed to improve the thermostability of glucoamylase from Aspergillus awamori based on the results of molecular dynamics simulations. The thermal unfolding of the catalytic domain followed a putative hierarchical behavior. In addition, the unfolding of the 13 alpha-helices obeyed the random ordered mechanism, in which the alpha-helices 8, 1 and 11 unfolded more rapidly than the others. The catalytic center was well protected by the (alpha/alpha)(6)-barrel at simulation temperatures up to 600 K, whereas the catalytic base, E400, migrated from its original interior pocket to the surface of the catalytic domain by surmounting the hydrophobic barrier provided by alpha-helices 12 and 13 at 800 K. The disulfide bonds engineered to 'lock' the alpha-helix 11 on the surface of the catalytic domain dramatically increased the thermostability. Substituting G396 and G407 with Ala residues slightly increased the thermostability, whereas their specific activity and catalytic efficiency were reduced. This indicates that the introduced residues with higher hydrophobicity were favorable in the loop between alpha-helices 12 and 13, whereas they partially destroyed the hydrogen bond and salt linkage network in the catalytic center. Alpha-helices 12 and 13 can be stabilized by introducing residues with higher hydrophobicity, except for the H391M mutation.

Amino Acid Substitution↗

PubMedDate not supplied

The effect of metal ions on the binding of ethanol to human alcohol dehydrogenase beta2beta2.

Molecular docking simulations were performed in this study to investigate the importance of both structural and catalytic zinc ions in the human alcohol dehydrogenase beta(2)beta(2) on substrate binding. The structural zinc ion is not only important in maintaining the structural integrity of the enzyme, but also plays an important role in determining substrate binding. The replacement of the catalytic zinc ion or both catalytic and structural zinc ions with Cu(2+) results in better substrate binding affinity than with the wild-type enzyme. The width of the bottleneck formed by L116 and V294 in the substrate binding pocket plays an important role for substrate entrance. In addition, unfavorable contacts between the substrate and T48 and F93 prevent the substrate from moving too close to the metal ion. The optimal binding position occurs between 1.9 and 2.4 A from the catalytic metal ion.

Alcohol Dehydrogenase↗

PubMedDate not supplied

Production, IMAC purification, and molecular modeling of N-carbamoyl-D-amino acid amidohydrolase C-terminally fused with a six-his peptide.

A six-His peptide was genetically engineered to the C-terminus of Agrobacterium radiobacter N-carbamoyl-D-amino acid amidohydrolase monomer to facilitate the protein purification with immobilized metal affinity chromatography (IMAC). The fusion enzyme, named as DCaseH, was overexpressed in Escherichia coli and one-step IMAC-purified. The production study showed that DCaseH was optimally produced at 15 degrees C for 25 h by the induction of 0.05 mM IPTG. Both Co(2+)-chelated TANOL gels and Ni(2+)-chelated nitriloacetic acid agarose gels efficiently purified DCaseH, with the former yielding purer enzyme than the latter. Highly pure DCaseH was obtained in the former purification with the addition of 5 mM imidazole in the washing buffer, and the specific enzyme activity was increased more than 11-fold. Denaturing IMAC purification successfully purified DCaseH from inclusion bodies that were mostly composed of the overexpressed DCaseH, while the attempt to refold the purified enzyme by either dialysis or solid-state refolding was not achieved. The purified native enzyme was optimally active at pH 6.5 and 50 degrees C, and the presence of 10% glycerol increased the activity. The molecular modeling of dimeric DCaseH indicated that the six-His tags were freely exposed to the protein surface, resulting in the selective and effective IMAC purification of DCaseH.

Amidohydrolases↗

PubMedDate not supplied

Predicted unfolding order of the 13 alpha-helices in the catalytic domain of glucoamylase from Aspergillus awamori var. X100 by molecular dynamics simulations.

The unfolding mechanism of the 13 alpha-helices in the catalytic domain of Aspergillus awamori var. X100 glucoamylase was investigated by 200 ps molecular dynamics simulations in explicit water with temperature jump technique. Rather than a simultaneous event, the unfolding of these 13 alpha-helices followed a random ordered mechanism as alpha8-->alpha1-->alpha11-->alpha7-->alpha10-->alpha3-->alpha12-->alpha13-->alpha4-->alpha5-->alpha9-->alpha6-->alpha2. No significant relationships were found between the unfolding order and the length and the hydrophobicity of the helix. alpha-Helix 8 located in the inner region of the catalytic domain was predicted to be the first helix to unfold, indicating that the destruction of the secondary structure motif was initiated from the inner region of the catalytic domain. The dynamic behavior of these alpha-helices induced by increased kinetic energy during the unfolding process is considered to be similar to the expansion and compression of a series of springs under the influence of mechanical stress.

Amino Acid Sequence↗

PubMedDate not supplied

Molecular dynamics simulations to determine the effect of supercritical carbon dioxide on the structural integrity of hen egg white lysozyme.

In this study, various molecular dynamics simulations were conducted to investigate the effect of supercritical carbon dioxide on the structural integrity of hen egg white lysozyme. The analyses of backbone root-mean-square deviation, radius of gyration, and secondary structure stability all show that supercritical CO(2) exhibits the ability to increase the stability of this protein, probably as a result of the solvent with less polarity, where hydrophobic interactions stabilizing the native structure are weakened and simultaneously the local hydrogen bonds are strengthened, resulting in stabilization of the secondary structures. The hydrophobic cores in the alpha- and beta-domains also play an important role in preventing this protein from thermal unfolding. As supercritical CO(2) has been attractive for biomedical applications because of the advantages of mild critical condition, nonflammability, nontoxity, and the purity of the resulting products, the structural stabilizing effect found in this study strongly suggests that it is possible to increase the thermostability of hen egg white lysozyme by pretreatment with supercritical CO(2), leading to better industrial applications of this protein.

Carbon Dioxide↗

PubMedDate not supplied

Homology model and molecular dynamics simulation of carp ovum cystatin.

In this study, a homology model of carp ovum cystatin was constructed based on the crystal structure of chicken egg white cystatin. The results of amino acid sequence alignment indicate that these two proteins exhibit 36.11% of sequence identity. The resultant homology model reveals that carp ovum cystatin shares similar folds as chicken egg white cystatin, particularly in the conserved regions of Q48-V49-G52 and P98-W99 and the locations of two disulfide bonds, C67-C76 and C90-C110. However, the results of 1 ns molecular dynamics simulations show that carp ovum cystatin exhibits less structural integrity than chicken egg white cystatin in explicit water at 300 K. The relatively hydrophilic Met62 of carp ovum cystatin, corresponding to the hydrophobic Leu68 of human cystatin C and Ile66 of chicken egg white cystatin, may destabilize the hydrophobic core and form a dimeric structure more easily through domain swapping. A total of 16 positively charged residues are equally distributed on the surface of carp ovum cystatin, resulting in agglutination with the negatively charged spermatozoa via electrostatic interaction. Thus, carp ovum cystatin is considered to be important in preventing carp eggs from polyspermy.

Agglutination↗