Optimizing nucleotide sequence ensembles for combinatorial protein libraries using a genetic algorithm

Protein libraries are essential to the field of protein engineering. Increasingly, probabilistic protein design is being used to synthesize combinatorial protein libraries, which allow the protein engineer to explore a vast space of amino acid sequences, while at the same time placing restrictions on the amino acid distributions. To this end, if site-specific amino acid probabilities are input as the target, then the codon nucleotide distributions that match this target distribution can be used to generate a partially randomized gene library. However, it turns out to be a highly nontrivial computational task to find the codon nucleotide distributions that exactly matches a given target distribution of amino acids. We first showed that for any given target distribution an exact solution may not exist at all. Formulated as a constrained optimization problem, we then developed a genetic algorithm-based approach to find codon nucleotide distributions that match as closely as possible to the target amino acid distribution. As compared with the previous gradient descent method on various objective functions, the new method consistently gave more optimized distributions as measured by the relative entropy between the calculated and the target distributions. To simulate the actual lab solutions, new objective functions were designed to allow for two separate sets of codons in seeking a better match to the target amino acid distribution

Modeling study on the validity of a possibly simplified representation of proteins

The folding characteristics of sequences reduced with a possibly simplified representation of five types of residues are shown to be similar to their original ones with the natural set of residues (20 types or 20 letters). The reduced sequences have a good foldability and fold to the same native structure of their optimized original ones. A large ground state gap for the native structure shows the thermodynamic stability of the reduced sequences. The general validity of such a five-letter reduction is further studied via the correlation between the reduced sequences and the original ones. As a comparison, a reduction with two letters is found not to reproduce the native structure of the original sequences due to its homopolymeric features.Comment: 6 pages with 4 figure

Protein design in a lattice model of hydrophobic and polar amino acids

A general strategy is described for finding which amino acid sequences have native states in a desired conformation (inverse design). The approach is used to design sequences of 48 hydrophobic and polar aminoacids on three-dimensional lattice structures. Previous studies employing a sequence-space Monte-Carlo technique resulted in the successful design of one sequence in ten attempts. The present work also entails the exploration of conformations that compete significantly with the target structure for being its ground state. The design procedure is successful in all the ten cases.Comment: RevTeX, 12 pages, 1 figur

Protein sequence and structure: Is one more fundamental than the other?

We argue that protein native state structures reside in a novel "phase" of matter which confers on proteins their many amazing characteristics. This phase arises from the common features of all globular proteins and is characterized by a sequence-independent free energy landscape with relatively few low energy minima with funnel-like character. The choice of a sequence that fits well into one of these predetermined structures facilitates rapid and cooperative folding. Our model calculations show that this novel phase facilitates the formation of an efficient route for sequence design starting from random peptides.Comment: 7 pages, 4 figures, to appear in J. Stat. Phy

Comparative modeling of DNA and RNA polymerases from Moniliophthora perniciosa mitochondrial plasmid

<p>Abstract</p> <p>Background</p> <p>The filamentous fungus <it>Moniliophthora perniciosa </it>(Stahel) Aime & Phillips-Mora is a hemibiotrophic Basidiomycota that causes witches' broom disease of cocoa (<it>Theobroma cacao </it>L.). This disease has resulted in a severe decrease in Brazilian cocoa production, which changed the position of Brazil in the market from the second largest cocoa exporter to a cocoa importer. Fungal mitochondrial plasmids are usually invertrons encoding DNA and RNA polymerases. Plasmid insertions into host mitochondrial genomes are probably associated with modifications in host generation time, which can be involved in fungal aging. This association suggests activity of polymerases, and these can be used as new targets for drugs against mitochondrial activity of fungi, more specifically against witches' broom disease. Sequencing and modeling: DNA and RNA polymerases of <it>M. perniciosa </it>mitochondrial plasmid were completely sequenced and their models were carried out by Comparative Homology approach. The sequences of DNA and RNA polymerase showed 25% of identity to 1XHX and 1ARO (pdb code) using BLASTp, which were used as templates. The models were constructed using Swiss PDB-Viewer and refined with a set of Molecular Mechanics (MM) and Molecular Dynamics (MD) in water carried out with AMBER 8.0, both working under the ff99 force fields, respectively. Ramachandran plots were generated by Procheck 3.0 and exhibited models with 97% and 98% for DNA and RNA polymerases, respectively. MD simulations in water showed models with thermodynamic stability after 2000 ps and 300 K of simulation.</p> <p>Conclusion</p> <p>This work contributes to the development of new alternatives for controlling the fungal agent of witches' broom disease.</p

The crystal structure of human Rogdi provides insight into the causes of Kohlschutter-Tonz Syndrome

Kohlschutter-T&amp;ouml;nz syndrome (KTS) is a rare autosomal-recessive disorder of childhood onset characterized by global developmental delay, spasticity, epilepsy, and amelogenesis imperfecta. Rogdi, an essential protein, is highly conserved across metazoans, and mutations in Rogdi are linked to KTS. However, how certain mutations in Rogdi abolish its physiological functions and cause KTS is not known. In this study, we determined the crystal structure of human Rogdi protein at atomic resolution. Rogdi forms a novel elongated curved structure comprising the ?? domain, a leucine-zipper-like four-helix bundle, and a characteristic ??-sheet domain. Within the ?? domain, the N-terminal H1 helix (residues 19-45) pairs with the C-terminal H6 helix (residues 252-287) in an antiparallel manner, indicating that the integrity of the four-helix bundle requires both N- and C-terminal residues. The crystal structure, in conjunction with biochemical data, indicates that the ?? domain might undergo a conformational change and provide a structural platform for protein-protein interactions. Disruption of the four-helix bundle by mutation results in significant destabilization of the structure. This study provides structural insights into how certain mutations in Rogdi affect its structure and cause KTS, which has important implications for the development of pharmaceutical agents against this debilitating neurological disease

C-Terminal Extension of the Yeast Mitochondrial DNA Polymerase Determines the Balance between Synthesis and Degradation

Saccharomyces cerevisiae mitochondrial DNA polymerase (Mip1) contains a C-terminal extension (CTE) of 279 amino acid residues. The CTE is required for mitochondrial DNA maintenance in yeast but is absent in higher eukaryotes. Here we use recombinant Mip1 C-terminal deletion mutants to investigate functional importance of the CTE. We show that partial removal of the CTE in Mip1Δ216 results in strong preference for exonucleolytic degradation rather than DNA polymerization. This disbalance in exonuclease and polymerase activities is prominent at suboptimal dNTP concentrations and in the absence of correctly pairing nucleotide. Mip1Δ216 also displays reduced ability to synthesize DNA through double-stranded regions. Full removal of the CTE in Mip1Δ279 results in complete loss of Mip1 polymerase activity, however the mutant retains its exonuclease activity. These results allow us to propose that CTE functions as a part of Mip1 polymerase domain that stabilizes the substrate primer end at the polymerase active site, and is therefore required for efficient mitochondrial DNA replication in vivo

Inversion of the balance between hydrophobic and hydrogen bonding interactions in protein folding and aggregation.

Identifying the forces that drive proteins to misfold and aggregate, rather than to fold into their functional states, is fundamental to our understanding of living systems and to our ability to combat protein deposition disorders such as Alzheimer's disease and the spongiform encephalopathies. We report here the finding that the balance between hydrophobic and hydrogen bonding interactions is different for proteins in the processes of folding to their native states and misfolding to the alternative amyloid structures. We find that the minima of the protein free energy landscape for folding and misfolding tend to be respectively dominated by hydrophobic and by hydrogen bonding interactions. These results characterise the nature of the interactions that determine the competition between folding and misfolding of proteins by revealing that the stability of native proteins is primarily determined by hydrophobic interactions between side-chains, while the stability of amyloid fibrils depends more on backbone intermolecular hydrogen bonding interactions