Chih-Chieh Chen

(PS)2: protein structure prediction server.

Protein structure prediction provides valuable insights into function, and comparative modeling is one of the most reliable methods to predict 3D structures directly from amino acid sequences. However, critical problems arise during the selection of the correct templates and the alignment of query sequences therewith. We have developed an automatic protein structure prediction server, (PS)2, which uses an effective consensus strategy both in template selection, which combines PSI-BLAST and IMPALA, and target–template alignment integrating PSI-BLAST, IMPALA and T-Coffee. (PS)2 was evaluated for 47 comparative modeling targets in CASP6 (Critical Assessment of Techniques for Protein Structure Prediction). For the benchmark dataset, the predictive performance of (PS)2, based on the mean GTD_TS score, was superior to 10 other automatic servers. Our method is based solely on the consensus sequence and thus is considerably faster than other methods that rely on the additional structural consensus of templates. Our results show that (PS)2, coupled with suitable consensus strategies and a new similarity score, can significantly improve structure prediction. Our approach should be useful in structure prediction and modeling. The (PS)2 is available through the website at .

(PS)2-v2: template-based protein structure prediction server

BackgroundTemplate selection and target-template alignment are critical steps for template-based modeling (TBM) methods. To identify the template for the twilight zone of 15~25% sequence similarity between targets and templates is still difficulty for template-based protein structure prediction. This study presents the (PS)2-v2 server, based on our original server with numerous enhancements and modifications, to improve reliability and applicability.ResultsTo detect homologous proteins with remote similarity, the (PS)2-v2 server utilizes the S2A2 matrix, which is a 60 × 60 substitution matrix using the secondary structure propensities of 20 amino acids, and the position-specific sequence profile (PSSM) generated by PSI-BLAST. In addition, our server uses multiple templates and multiple models to build and assess models. Our method was evaluated on the Lindahl benchmark for fold recognition and ProSup benchmark for sequence alignment. Evaluation results indicated that our method outperforms sequence-profile approaches, and had comparable performance to that of structure-based methods on these benchmarks. Finally, we tested our method using the 154 TBM targets of the CASP8 (Critical Assessment of Techniques for Protein Structure Prediction) dataset. Experimental results show that (PS)2-v2 is ranked 6th among 72 severs and is faster than the top-rank five serves, which utilize ab initio methods.ConclusionExperimental results demonstrate that (PS)2-v2 with the S2A2 matrix is useful for template selections and target-template alignments by blending the amino acid and structural propensities. The multiple-template and multiple-model strategies are able to significantly improve the accuracies for target-template alignments in the twilight zone. We believe that this server is useful in structure prediction and modeling, especially in detecting homologous templates with sequence similarity in the twilight zone.

Computational Analysis of KRAS Mutations: Implications for Different Effects on the KRAS p.G12D and p.G13D Mutations

Background The issue of whether patients diagnosed with metastatic colorectal cancer who harbor KRAS codon 13 mutations could benefit from the addition of anti-epidermal growth factor receptor therapy remains under debate. The aim of the current study was to perform computational analysis to investigate the structural implications of the underlying mutations caused by c.38G>A (p.G13D) on protein conformation. Methods Molecular dynamics (MD) simulations were performed to understand the plausible structural and dynamical implications caused by c.35G>A (p.G12D) and c.38G>A (p.G13D). The potential of mean force (PMF) simulations were carried out to determine the free energy profiles of the binding processes of GTP interacting with wild-type (WT) KRAS and its mutants (MT). Results Using MD simulations, we observed that the root mean square deviation (RMSD) increased as a function of time for the MT c.35G>A (p.G12D) and MT c.38G>A (p.G13D) when compared with the WT. We also observed that the GTP-binding pocket in the c.35G>A (p.G12D) mutant is more open than that of the WT and the c.38G>A (p.G13D) proteins. Intriguingly, the analysis of atomic fluctuations and free energy profiles revealed that the mutation of c.35G>A (p.G12D) may induce additional fluctuations in the sensitive sites (P-loop, switch I and II regions). Such fluctuations may promote instability in these protein regions and hamper GTP binding. Conclusions Taken together with the results obtained from MD and PMF simulations, the present findings implicate fluctuations at the sensitive sites (P-loop, switch I and II regions). Our findings revealed that KRAS mutations in codon 13 have similar behavior as KRAS WT. To gain a better insight into why patients with metastatic colorectal cancer (mCRC) and the KRAS c.38G>A (p.G13D) mutation appear to benefit from anti-EGFR therapy, the role of the KRAS c.38G>A (p.G13D) mutation in mCRC needs to be further investigated.

Characteristics and prevalence of KRAS, BRAF, and PIK3CA mutations in colorectal cancer by high-resolution melting analysis in Taiwanese population

BACKGROUND The identification of KRAS, BRAF, and PIK3CA mutations before the administration of anti-epidermal growth factor receptor therapy of colorectal cancer has become important. The aim of the present study was to investigate the occurrence of KRAS, BRAF, and PIK3CA mutations in the Taiwanese population with colorectal cancer. This study was undertaken to identify BRAF and PIK3CA mutations in patients with colorectal cancer by high-resolution melting (HRM) analysis. HRM analysis is a new gene scan tool that quickly performs the PCR and identifies sequence alterations without requiring post-PCR treatment. METHODS In the present study, DNAs were extracted from 182 cases of formalin-fixed, paraffin-embedded (FFPE) colorectal cancer samples for clinical KRAS mutational analysis by direct sequencing. All the samples were also tested for mutations within BRAF V600E and PIK3CA (exons 9 and 20) by HRM analysis. RESULTS The results were confirmed by direct sequencing. The frequency of BRAF and PIK3CA mutations is 1.1%, and 7.1%, respectively. Intriguingly, we found that nine patients (4.9%) with the KRAS mutation were coexistent with the PIK3CA mutation. Four patients (2.2%) without the KRAS mutation were existent with the PIK3CA mutation. Two patients (1.1%) without the KRAS mutation were existent with the BRAF mutation. CONCLUSIONS In the current study, we suppose that HRM analysis is rapid, feasible, and powerful diagnostic tool for the detection of BRAF and PIK3CA mutations in a clinical setting. Additionally, our results indicated the prevalence of KRAS, BRAF, and PIK3CA mutational status in the Taiwanese population.

Interplay between CRP, Atherogenic LDL, and LOX-1 and Its Potential Role in the Pathogenesis of Atherosclerosis

BACKGROUND Studies have shown that the classic acute-phase protein C-reactive protein (CRP) has proinflammatory effects on vascular cells and may play a causal role in the pathogenesis of coronary artery disease. A growing body of evidence has suggested that interplay between CRP, lectin-like oxidized LDL receptor-1 (LOX-1), and atherogenic LDL may underlie the mechanism of endothelial dysfunction that leads to atherosclerosis. CONTENT We review the biochemical evidence for an association of CRP, LOX-1, and either oxidized LDL (OxLDL) or electronegative L5 LDL with the pathogenesis of coronary artery disease. Artificially oxidized OxLDL has been studied extensively for its role in atherogenesis, as has electronegative L5 LDL, which is present at increased levels in patients with increased cardiovascular risks. OxLDL and L5 have been shown to stimulate human aortic endothelial cells to produce CRP, indicating that CRP is synthesized locally in the endothelium. The ligand-binding face (B-face) of CRP has been shown to bind the LOX-1 scavenger receptor and increase LOX-1 expression in endothelial cells, thereby promoting the uptake of OxLDL or L5 by LOX-1 into endothelial cells to induce endothelial dysfunction. SUMMARY CRP and LOX-1 may form a positive feedback loop with OxLDL or L5 in atherogenesis, whereby increased levels of atherogenic LDL in patients with cardiovascular risks induce endothelial cells to express CRP, which may in turn increase the expression of LOX-1 to promote the uptake of atherogenic LDL into endothelial cells. Further research is needed to confirm a causal role for CRP in atherogenesis.

Clinical relevance of KRAS mutations in codon 13: Where are we?

Recent advances in molecular diagnosis and the trend towards personalized medicine have made colorectal cancer one of the tumors where therapies have significantly improved patient survival after metastasis development. KRAS mutations in codon 12 and 13 are recognized biomarkers that are analyzed in clinic previously for anti-EGFR therapies administration. Since originally mutations in both codons were considered as a predictor of lack of response to cetuximab or panitumumab, the European Medicines Agency and the US Food and Drug Administration suggested that patients harboring any of those mutations should be excluded from the treatment. However, subsequent retrospective analysis has shown that mutations in codon 12 and codon 13 of KRAS gene could be different in their biological characteristics and as a result could confer variable effects in patients. In addition and increasing and sometimes contradictory number of solutions have been published demonstrating that patients with mutations in codon 13 could have worse outcome but could obtain a significant clinical benefit from anti-EGFR therapies. Here, we review and update the latest data on the biological role leading to a predictive outcome and benefit from anti-EGFR antibodies in patients with specific KRAS mutations in codon 13.

pKNOT v.2: the protein KNOT web server

Knotted proteins have recently received lots of attention due to their interesting topological novelty as well as its puzzling folding mechanisms. We previously published a pKNOT server, which provides a structural database of knotted proteins, analysis tools for detecting and analyzing knotted regions from structures as well as a Java-based 3D graphics viewer for visualizing knotted structures. However, there lacks a convenient platform performing similar tasks directly from ‘protein sequences’. In the current version of the web server, referred to as pKNOT v.2, we implement a homology modeling tool such that the server can now accept protein sequences in addition to 3D structures or Protein Data Bank (PDB) IDs and return knot analysis. In addition, we have updated the database of knotted proteins from the current PDB with a combination of automatic and manual procedure. We believe that the updated pKNOT server with its extended functionalities will provide better service to biologists interested in the research of knotted proteins. The pKNOT v.2 is available from http://pknot.life.nctu.edu.tw/.

Computational analysis of a novel mutation in ETFDH gene highlights its long-range effects on the FAD-binding motif

BackgroundMultiple acyl-coenzyme A dehydrogenase deficiency (MADD) is an autosomal recessive disease caused by the defects in the mitochondrial electron transfer system and the metabolism of fatty acids. Recently, mutations in electron transfer flavoprotein dehydrogenase (ETFDH) gene, encoding electron transfer flavoprotein:ubiquinone oxidoreductase (ETF:QO) have been reported to be the major causes of riboflavin-responsive MADD. To date, no studies have been performed to explore the functional impact of these mutations or their mechanism of disrupting enzyme activity.ResultsHigh resolution melting (HRM) analysis and sequencing of the entire ETFDH gene revealed a novel mutation (p.Phe128Ser) and the hotspot mutation (p.Ala84Thr) from a patient with MADD. According to the predicted 3D structure of ETF:QO, the two mutations are located within the flavin adenine dinucleotide (FAD) binding domain; however, the two residues do not have direct interactions with the FAD ligand. Using molecular dynamics (MD) simulations and normal mode analysis (NMA), we found that the p.Ala84Thr and p.Phe128Ser mutations are most likely to alter the protein structure near the FAD binding site as well as disrupt the stability of the FAD binding required for the activation of ETF:QO. Intriguingly, NMA revealed that several reported disease-causing mutations in the ETF:QO protein show highly correlated motions with the FAD-binding site.ConclusionsBased on the present findings, we conclude that the changes made to the amino acids in ETF:QO are likely to influence the FAD-binding stability.

(PS)2: protein structure prediction server version 3.0

Protein complexes are involved in many biological processes. Examining coupling between subunits of a complex would be useful to understand the molecular basis of protein function. Here, our updated (PS)2 web server predicts the three-dimensional structures of protein complexes based on comparative modeling; furthermore, this server examines the coupling between subunits of the predicted complex by combining structural and evolutionary considerations. The predicted complex structure could be indicated and visualized by Java-based 3D graphics viewers and the structural and evolutionary profiles are shown and compared chain-by-chain. For each subunit, considerations with or without the packing contribution of other subunits cause the differences in similarities between structural and evolutionary profiles, and these differences imply which form, complex or monomeric, is preferred in the biological condition for the subunit. We believe that the (PS)2 server would be a useful tool for biologists who are interested not only in the structures of protein complexes but also in the coupling between subunits of the complexes. The (PS)2 is freely available at http://ps2v3.life.nctu.edu.tw/.

PGAdb-builder: A web service tool for creating pan-genome allele database for molecular fine typing

With the advance of next generation sequencing techniques, whole genome sequencing (WGS) is expected to become the optimal method for molecular subtyping of bacterial isolates. To use WGS as a general subtyping method for disease outbreak investigation and surveillance, the layout of WGS-based typing must be comparable among laboratories. Whole genome multilocus sequence typing (wgMLST) is an approach that achieves this requirement. To apply wgMLST as a standard subtyping approach, a pan-genome allele database (PGAdb) for the population of a bacterial organism must first be established. We present a free web service tool, PGAdb-builder (http://wgmlstdb.imst.nsysu.edu.tw), for the construction of bacterial PGAdb. The effectiveness of PGAdb-builder was tested by constructing a pan-genome allele database for Salmonella enterica serovar Typhimurium, with the database being applied to create a wgMLST tree for a panel of epidemiologically well-characterized S. Typhimurium isolates. The performance of the wgMLST-based approach was as high as that of the SNP-based approach in Leekitcharoenphon’s study used for discerning among epidemiologically related and non-related isolates.