Search PubMed⌕ Search

Biomedical subjects

Chung-Jung Tsai

Publications and source records attributed to Chung-Jung Tsai.

At least 19 recordsLinked to original sources
PubMed2007

A hierarchical protein folding scheme based on the building block folding model.

The building block protein folding model states that the native protein structure is the product of a combinatorial assembly of relatively structurally independent contiguous parts of the protein that possess a hydrophobic core, i.e., building blocks (BBs). According to this model, our group proposed a three-stage scheme for a feasible time-wise semi ab-intio protein structure prediction. Given a protein sequence, at the first stage of the prediction scheme, we propose cutting the sequence into structurally assigned BBs. Next, we perform a combinatorial assembly and attempt to predict the relative three-dimensional arrangement of the BBs. In the third stage, we refine and rank the assemblies. The scheme has proven to be very promising in reducing the complexity of the protein folding problem and gaining insight into the protein folding process. In this chapter, we describe the different stages of the scheme and discuss a possible application of the model to protein design.

Algorithms↗
PubMed2006

Structure by design: from single proteins and their building blocks to nanostructures.

Nanotechnology realizes the advantages of naturally occurring biological macromolecules and their building-block nature for design. Frequently, assembly starts with the choice of a "good" molecule that is synthetically optimized towards the desired shape. By contrast, we propose starting with a pre-specified nanostructure shape, selecting candidate protein building blocks from a library and mapping them onto the shape and, finally, testing the stability of the construct. Such a shape-based, part-assembly strategy is conceptually similar to protein design through the combinatorial assembly of building blocks. If the conformational preferences of the building blocks are retained and their interactions are favorable, the nanostructure will be stable. The richness of the conformations, shapes and chemistries of the protein building blocks suggests a broad range of potential applications; at the same time, it also highlights their complexity. In this Opinion article, we focus on the first step: validating such a strategy against experimental data.

Drug Design↗
PubMed2006

Critical components of odors in evaluating the performance of food waste composting plants.

The current Taiwan government policy toward food waste management encourages composting for resource recovery. This study used olfactometry, gas chromatography-mass spectrometry (GC-MS) and gas detector tubes to evaluate the ambient air at three of the largest food waste composting plants in Taiwan. Ambient air inside the plants, at exhaust outlets and plant boundaries was examined to determine the comprehensive odor performance, critical components, and odor elimination efficiencies of various odor control engineering. Analytical results identified 29 compounds, including ammonia, amines, acetic acid, and multiple volatile organic compounds (VOCs) (hydrocarbons, ketones, esters, terpenes and S-compounds) in the odor from food waste composting plants. Concentrations of six components--ammonia, amines, dimethyl sulfide, acetic acid, ethyl benzene and p-Cymene--exceeded human olfactory thresholds. Ammonia, amines, dimethyl sulfide and acetic acid accounted for most odors compared to numerous VOCs. The results also show that the biotrickling filter was better at eliminating the concentrations of odor, NH(3), amines, S-compounds and VOCs than the chemical scrubber and biofilters. All levels measured by olfactometry at the boundaries of food waste composting plants (range, 74-115 Odor Concentration (OC)) exceeded Taiwan's EPA standard of 50 OC. This study indicated that the malodor problem continued to be a significant problem for food waste recovery.

Air Pollutants↗
PubMed2006

Structural stability and dynamics of an amyloid-forming peptide GNNQQNY from the yeast prion sup-35.

A seven amino acid yeast prion sup-35 fragment (GNNQQNY) forms amyloid fibrils. The availability of its detailed atomic oligomeric structure makes it a good model for studying the early stage of aggregation. Here we perform long all-atom explicit solvent molecular simulations of various sizes and arrangements of oligomer seeds of the wild-type and its mutants to study its stability and dynamics. Previous studies have suggested that the early stage rate-limiting step of oligomer formation occurs in high-order oligomers. Our simulations show that with the increase in the number of strands even from a dimer to a trimer, oligomer stability increases dramatically. This suggests that the minimal nucleus seed for GNNQQNY fibril formation could be small and is likely three or four peptides, in agreement with experiment, and that higher-order oligomers do not dissociate quickly since they have small diffusion coefficients and thus slow kinetics. Further, for the hydrophilic polar GNNQQNY, there are no hydrogen bonds and no hydrophobic interactions between adjacent beta-sheets. Simulations suggest that within the sheet, the driving forces to associate and stabilize are interstrand backbone-backbone and side chain-side chain hydrogen bonds, whereas between the sheets, shape-complementary by the dry polar steric zipper via the side chains of Asn-2, Gln-4, and Asn-6 holds the sheets together, as proposed in an earlier study. Since the polar side chains of Asn-2, Gln-4, and Asn-6 act as a hook to bind two neighboring sheets together, these geometric restraints reduce the conformational search for the correct side chain packing to a two-dimensional problem of intersheet side chain interactions. Mutant simulations show that substitution of Asn-2, Gln-4, or Asn-6 by Ala would disrupt this steric zipper, leading to unstable oligomers.

Amino Acid Sequence↗
PubMed2006

Designing a nanotube using naturally occurring protein building blocks.

Here our goal is to carry out nanotube design using naturally occurring protein building blocks. Inspection of the protein structural database reveals the richness of the conformations of proteins, their parts, and their chemistry. Given target functional protein nanotube geometry, our strategy involves scanning a library of candidate building blocks, combinatorially assembling them into the shape and testing its stability. Since self-assembly takes place on time scales not affordable for computations, here we propose a strategy for the very first step in protein nanotube design: we map the candidate building blocks onto a planar sheet and wrap the sheet around a cylinder with the target dimensions. We provide examples of three nanotubes, two peptide and one protein, in atomistic model detail for which there are experimental data. The nanotube models can be used to verify a nanostructure observed by low-resolution experiments, and to study the mechanism of tube formation.

HIV-1↗
PubMed2005

Energy landscape of amyloidogenic peptide oligomerization by parallel-tempering molecular dynamics simulation: significant role of Asn ladder.

Recent evidence suggests that amyloidogenic oligomers may be the toxic species in cell cultures. Thus, it is crucial to understand their structure and oligomerization mechanism in atomistic detail. By employing tens of fast central processing units and an advanced phase-space sampling algorithm, parallel-tempering molecular dynamics, we have explored the energy landscape of amyloidogenic peptide oligomerization in explicit water. A pentapeptide, DFNKF, derived from human calcitonin and its mutant, DFAKF, was simulated with a total simulation time of approximately 500 ns. The detailed oligomerization process of a DFNKF parallel beta-sheet formation at 300 K has been characterized. The assembly of a parallel beta-sheet from the amorphous state mainly occurs via a "bottleneck" channel where the interstrand Asn-Asn stacking is the major interaction. The interactions of Asn-Asn stacking include both backbone and side-chain hydrogen bonds. The Asn-Asn interactions work like "glue" by sticking the DFNKF strands together and assist the "on-pathway" oligomerization. The Asn-Asn stacking observed here is similar to the Asn ladder commonly found in globular beta-helix proteins. A control run shows that when Asn is mutated to Ala, the stability and population of the DFAKF parallel beta-sheet is decreased. Furthermore, our in vitro mutagenesis experiments show that the ability of DFAKF peptides to form amyloid fibrils is significantly reduced, in agreement with the simulations. Knowledge of the energy landscape of oligomerization may provide hints for rational drug design, preventing amyloid-associated diseases.

Algorithms↗
PubMed2004

Analysis of ordered and disordered protein complexes reveals structural features discriminating between stable and unstable monomers.

Most proteins exist in the cell as multi-component assemblies. However, which proteins need to be present simultaneously in order to perform a given function is frequently unknown. The first step toward this goal would be to predict proteins that can function only when in a complexed form. Here, we propose a scheme to distinguish whether the protein components are ordered (stable) or disordered when separated from their complexed partners. We analyze structural characteristics of several types of complexes, such as natively unstructured proteins, ribosomal proteins, two-state and three-state complexes, and crystal-packing dimers. Our analysis makes use of the fact that natively unstructured proteins, which undergo a disorder-to-order transition upon binding their partner, and stable monomeric proteins, which exist as dimers only in their crystal form, provide examples of two vastly different scenarios. We find that ordered monomers can be distinguished from disordered monomers on the basis of the per-residue surface and interface areas, which are significantly smaller for ordered proteins. With this scale, two-state dimers (where the monomers unfold upon dimer separation) and ribosomal proteins are shown to resemble disordered proteins. On the other hand, crystal-packing dimers, whose monomers are stable in solution, fall into the ordered protein category. While there should be a continuum in the distributions, nevertheless, the per-residue scale measures the confidence in the determination of whether a protein can exist as a stable monomer. Further analysis, focusing on the chemical and contact preferences at the interface, interior and exposed surface areas, reveals that disordered proteins lack a strong hydrophobic core and are composed of highly polar surface area. We discuss the implication of our results for de novo design of stable monomeric proteins and peptides.

Macromolecular Substances↗
PubMed2004

Insights into amyloid structural formation and assembly through computational approaches.

Amyloids are long, insoluble ordered fibers. Due to their insolubility, to date the determination of an amyloid structure with an atomic scale resolution has proven to be a difficult task. Under such circumstances, computational approaches are a preferred option, providing the means to build likely models, test their stabilities and figure out the chemistry of their prevailing interactions. Computational models can be validated by targeted experiments, such as introducing mutations and testing for amyloid formation. Computations further provide vehicles for the comprehension of the mechanisms of amyloid seed formation and oligomer toxicity. Nevertheless, computations face an immense hurdle, the outcome of the time scales involved in amyloid formation and the immense sizes of the systems. In an attempt to overcome these, we adopt a strategy that encompasses (1) bioinformatics studies of native proteins containing beta-sheet structures; (2) simulations of shorter peptides; and finally (3) construction of potential oligomeric models and tests of their stabilities. The results are correlated with experimental data where available. Here, we describe the computational methods in simple terms and present an overview of the results. The systems derive from amyloidogenic, disease-related proteins, including gelsolin, beta2-microglobulin, and peptides derived from the prion, Alzheimer's Abeta, IAPP and human calcitonin. Ultimately, obtaining molecular structures should facilitate efforts to therapy and drug design.

Amyloid↗
PubMed2004

A new, structurally nonredundant, diverse data set of protein-protein interfaces and its implications.

Here, we present a diverse, structurally nonredundant data set of two-chain protein-protein interfaces derived from the PDB. Using a sequence order-independent structural comparison algorithm and hierarchical clustering, 3799 interface clusters are obtained. These yield 103 clusters with at least five nonhomologous members. We divide the clusters into three types. In Type I clusters, the global structures of the chains from which the interfaces are derived are also similar. This cluster type is expected because, in general, related proteins associate in similar ways. In Type II, the interfaces are similar; however, remarkably, the overall structures and functions of the chains are different. The functional spectrum is broad, from enzymes/inhibitors to immunoglobulins and toxins. The fact that structurally different monomers associate in similar ways, suggests "good" binding architectures. This observation extends a paradigm in protein science: It has been well known that proteins with similar structures may have different functions. Here, we show that it extends to interfaces. In Type III clusters, only one side of the interface is similar across the cluster. This structurally nonredundant data set provides rich data for studies of protein-protein interactions and recognition, cellular networks and drug design. In particular, it may be useful in addressing the difficult question of what are the favorable ways for proteins to interact. (The data set is available at http://protein3d.ncifcrf.gov/~keskino/ and http://home.ku.edu.tr/~okeskin/INTERFACE/INTERFACES.html.)

Algorithms↗
PubMed2004

In silico protein design by combinatorial assembly of protein building blocks.

Utilizing concepts of protein building blocks, we propose a de novo computational algorithm that is similar to combinatorial shuffling experiments. Our goal is to engineer new naturally occurring folds with low homology to existing proteins. A selected protein is first partitioned into its building blocks based on their compactness, degree of isolation from the rest of the structure, and hydrophobicity. Next, the protein building blocks are substituted by fragments taken from other proteins with overall low sequence identity, but with a similar hydrophobic/hydrophilic pattern and a high structural similarity. These criteria ensure that the designed protein has a similar fold, low sequence identity, and a good hydrophobic core compared with its native counterpart. Here, we have selected two proteins for engineering, protein G B1 domain and ubiquitin. The two engineered proteins share approximately 20% and approximately 25% amino acid sequence identities with their native counterparts, respectively. The stabilities of the engineered proteins are tested by explicit water molecular dynamics simulations. The algorithm implements a strategy of designing a protein using relatively stable fragments, with a high population time. Here, we have selected the fragments by searching for local minima along the polypeptide chain using the protein building block model. Such an approach provides a new method for engineering new proteins with similar folds and low homology.

Amino Acid Sequence↗
PubMed2004

The stability and dynamics of the human calcitonin amyloid peptide DFNKF.

The stability and dynamics of the human calcitonin-derived peptide DFNKF (hCT(15-19)) are studied using molecular dynamics (MD) simulations. Experimentally, this peptide is highly amyloidogenic and forms fibrils similar to the full length calcitonin. Previous comparative MD studies have found that the parallel beta-stranded sheet is a stable organization of the DFNKF protofibril. Here, we probe the stability and dynamics of the small parallel DFNKF oligomers. The results show that even small DFNKF oligomers, such as trimers and tetramers, are stable for a sufficient time in the MD simulations, indicating that the crucial nucleus seed size for amyloid formation can be quite small. The simulations also show that the stability of DFNKF oligomers increases with their sizes. The small but stable seed may reflect the experimental rapid formation of the DFNKF fibrils. Further, a noncooperative process of parallel beta-sheet formation from the out-of-register trimer is observed in the simulations. In general, the residues of DFNKF peptides near the N-/C-termini are more flexible, whereas the interior residues are more stable. Simulations of mutants and capped peptides show that both interstrand hydrophobic and electrostatic interactions play important roles in stabilizing the DFNKF parallel oligomers. This study provides insights into amyloid formation.

Amyloid↗
PubMed2003

Hierarchical protein folding pathways: a computational study of protein fragments.

We have previously presented a building block folding model. The model postulates that protein folding is a hierarchical top-down process. The basic unit from which a fold is constructed, referred to as a hydrophobic folding unit, is the outcome of combinatorial assembly of a set of "building blocks." Results obtained by the computational cutting procedure yield fragments that are in agreement with those obtained experimentally by limited proteolysis. Here we show that as expected, proteins from the same family give very similar building blocks. However, different proteins can also give building blocks that are similar in structure. In such cases the building blocks differ in sequence, stability, contacts with other building blocks, and in their 3D locations in the protein structure. This result, which we have repeatedly observed in many cases, leads us to conclude that while a building block is influenced by its environment, nevertheless, it can be viewed as a stand-alone unit. For small-sized building blocks existing in multiple conformations, interactions with sister building blocks in the protein will increase the population time of the native conformer. With this conclusion in hand, it is possible to develop an algorithm that predicts the building block assignment of a protein sequence whose structure is unknown. Toward this goal, we have created sequentially nonredundant databases of building block sequences. A protein sequence can be aligned against these, in order to be matched to a set of potential building blocks.

Algorithms↗
PubMed2003

Extended disordered proteins: targeting function with less scaffold.

It has been estimated that a large fraction of cellular proteins are natively disordered. Current opinion largely holds that natively disordered proteins are more 'adaptive', leading to advantages in regulation and in binding diverse ligands. Here, we argue for another, simple, physically based reason. Disordered proteins often have large intermolecular interfaces, the size of which is dictated by protein function. For proteins to be stable as monomers with extensive interfaces, protein size would need to be 2-3 times larger. This would either increase cellular crowding or enlarge the size of the cell by 15-30%, owing to the increase in the sequence length. Smaller sizes of cells, proteins, DNA and RNA conserve energy. Thus, disordered proteins provide a simple yet elegant solution to having large intermolecular interfaces, but with smaller protein, genome and cell sizes.

Models, Molecular↗
PubMed2003

Reducing the computational complexity of protein folding via fragment folding and assembly.

Understanding, and ultimately predicting, how a 1-D protein chain reaches its native 3-D fold has been one of the most challenging problems during the last few decades. Data increasingly indicate that protein folding is a hierarchical process. Hence, the question arises as to whether we can use the hierarchical concept to reduce the practically intractable computational times. For such a scheme to work, the first step is to cut the protein sequence into fragments that form local minima on the polypeptide chain. The conformations of such fragments in solution are likely to be similar to those when the fragments are embedded in the native fold, although alternate conformations may be favored during the mutual stabilization in the combinatorial assembly process. Two elements are needed for such cutting: (1) a library of (clustered) fragments derived from known protein structures and (2) an assignment algorithm that selects optimal combinations to "cover" the protein sequence. The next two steps in hierarchical folding schemes, not addressed here, are the combinatorial assembly of the fragments and finally, optimization of the obtained conformations. Here, we address the first step in a hierarchical protein-folding scheme. The input is a target protein sequence and a library of fragments created by clustering building blocks that were generated by cutting all protein structures. The output is a set of cutout fragments. We briefly outline a graph theoretic algorithm that automatically assigns building blocks to the target sequence, and we describe a sample of the results we have obtained.

Algorithms↗
PubMed2002

Maximal stabilities of reversible two-state proteins.

The hydrophobic effect is the major force driving protein folding. Around room temperature, small organic solutes and hydrophobic amino acids have low solubilities in water and the hydrophobic effect is the strongest. These facts suggest that globular proteins should be maximally stable around room temperature. While this fundamental paradigm has been expected, it has not actually been shown to hold. Toward this goal, we have collected and analyzed experimental thermodynamic data for 31 proteins that show reversible two-state folding <--> unfolding transitions at or near neutral pH. Twenty-six of these are unique, and 20 of the 26 are maximally stable around room temperature irrespective of their structural properties, the melting temperature, or the living temperatures of their source organisms. Their average temperature of maximal stability is 293 +/- 8 K (20 +/- 8 degrees C). These proteins differ in size, fold, and number of domains, hydrophobic folding units, and oligomeric states. They derive from the cold-loving psychrophiles, from mesophiles, and from thermophiles. Analysis of the single-domain proteins present in this set shows that the variations in their thermodynamic parameters are correlated in a way which may explain the adaptation of the proteins to the living temperatures of the organisms from which they derive. The average energetic contribution of the individual amino acids toward protein stability decreases with an increase in protein size, suggesting that there may be an upper limit for protein maximal thermodynamic stability. For the remaining proteins, deviation of the maximal stability temperatures from room temperature may be due to greater uncertainties in their heat capacity change (DeltaC(p)) values, a weaker hydrophobic effect, and/or a stronger electrostatic contribution.

Animals↗
PubMed2002

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway.

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.

Computer Simulation↗
PubMed2002

The hydrophobic effect: a new insight from cold denaturation and a two-state water structure.

Herein we provide a new insight into the hydrophobic effect in protein folding. Our proposition explains the molecular basis of cold denaturation, and of intermediate states in heat and their absence in cold denaturation. The exposure of non-polar surface reduces the entropy and enthalpy of the system, at low and at high temperatures. At low temperatures the favorable reduction in enthalpy overcomes the unfavorable reduction in entropy, leading to cold denaturation. At high temperatures, folding/unfolding is a two-step process: in the first, the entropy gain leads to hydrophobic collapse, in the second, the reduction in enthalpy due to protein-protein interactions leads to the native state. The different entropy and enthalpy contributions to the Gibbs energy change at each step at high, and at low, temperatures can be conveniently explained by a two-state model of the water structure. The model provides a clear view of the dominant factors in protein folding and stability. Consequently, it appears to provide a microscopic view of the hydrophobic effect and is consistently linked to macroscopic thermodynamic parameters.

Cold Temperature↗