Multiple Unfolding Intermediates Obtained by Molecular Dynamic Simulations under Stretching for Immunoglobulin-Binding Domain of Protein G

We have studied the mechanical properties of the immunoglobulin-binding domain of protein G at the atomic level under stretching at constant velocity using molecular dynamics simulations. We have found that the unfolding process can occur either in a single step or through intermediate states. Analysis of the trajectories from the molecular dynamic simulations showed that the mechanical unfolding of the immunoglobulin-binding domain of protein G is triggered by the separation of the terminal β-strands and the order in which the secondary-structure elements break is practically the same in two- and multi-state events and at the different extension velocities studied. It is seen from our analysis of 24 trajectories that the theoretical pathway of mechanical unfolding for the immunoglobulin-binding domain of protein G does not coincide with that proposed in denaturant studies in the absence of force

Development and testing of new force fields for molecular dynamics simulations

Recent progress in modeling of protein folding in Dr. Shaw laboratory has been achieved only after some improvements of potentials of covalent forces, taken from the standard AMBER force field; and still, the force field used is not quite satisfactory to reproduce folded structures of some larger proteins, having significant, about 5A, RMS deviation between the computed and experimentally determined 3D structures. The objective of this research is to develop and test new polarizable atomic force fields (FFs) for "in-vacuum" and "in-water" non-bonded interactions based on AMBER ff99SBILDN force fields, improved by inclusion of new terms. FFs parameter optimization will be done using our set of molecular crystals with crystallographic data from the Cambridge Structural Database and sublimation/solvation thermodynamics characteristics from various sources

To be folded, to be unfolded or to be aggregated with important functions: application of the directed coaggregation mechanism to combat bacterial communities

One of the reasons for the mortal danger to humans is the ability of pathogenic bacteria to form biofilms. The formation of biofilms is an evolutionarily conservative defense mechanism against adverse conditions. The use of this protection by pathogenic bacteria reduces the effectiveness of the main means of combating them - antibiotics, which complicates the production of new types of drugs. There are two types of antimicrobial agents that are not known antibiotics: nanoparticles and antimicrobial peptides. We demonstrated that peptides synthesized based on the amino acid sequence of proteins and capable of amyloid formation and coaggregation with the whole protein exhibit antimicrobial activity. The ability of peptides to coaggregate with target proteins can help combat biofilm-forming bacterial communities. We evaluated the antimicrobial effects of ten synthesized hybrid peptides, which were obtained based on the sequences of the S1 ribosomal protein of P. aeruginosa and S. aureus. It is important that some peptides demonstrated high antimicrobial activity comparable to the antibiotic gentamicin sulfate against pathogenic strains of MRSA, S. aureus, and P. aeruginosa. These peptides showed no toxicity to eukaryotic cells. Our study demonstrates the promise of hybrid peptides based on the amyloidogenic regions of the S1 ribosomal protein for the development of new antimicrobials against Gram-positive and Gram-negative bacteria resistant to traditional antibiotic.Belgrade : Institute of molecular genetics and genetic engineerin

Prots: A fragment based protein thermo‐stability potential

Designing proteins with enhanced thermo‐stability has been a main focus of protein engineering because of its theoretical and practical significance. Despite extensive studies in the past years, a general strategy for stabilizing proteins still remains elusive. Thus effective and robust computational algorithms for designing thermo‐stable proteins are in critical demand. Here we report PROTS, a sequential and structural four‐residue fragment based protein thermo‐stability potential. PROTS is derived from a nonredundant representative collection of thousands of thermophilic and mesophilic protein structures and a large set of point mutations with experimentally determined changes of melting temperatures. To the best of our knowledge, PROTS is the first protein stability predictor based on integrated analysis and mining of these two types of data. Besides conventional cross validation and blind testing, we introduce hypothetical reverse mutations as a means of testing the robustness of protein thermo‐stability predictors. In all tests, PROTS demonstrates the ability to reliably predict mutation induced thermo‐stability changes as well as classify thermophilic and mesophilic proteins. In addition, this white‐box predictor allows easy interpretation of the factors that influence mutation induced protein stability changes at the residue level. Proteins 2012; © 2011 Wiley Periodicals, Inc.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/89526/1/23163_ftp.pd

Proteome sequence features carry signatures of the environmental niche of prokaryotes

<p>Abstract</p> <p>Background</p> <p>Prokaryotic environmental adaptations occur at different levels within cells to ensure the preservation of genome integrity, proper protein folding and function as well as membrane fluidity. Although specific composition and structure of cellular components suitable for the variety of extreme conditions has already been postulated, a systematic study describing such adaptations has not yet been performed. We therefore explored whether the environmental niche of a prokaryote could be deduced from the sequence of its proteome. Finally, we aimed at finding the precise differences between proteome sequences of prokaryotes from different environments.</p> <p>Results</p> <p>We analyzed the proteomes of 192 prokaryotes from different habitats. We collected detailed information about the optimal growth conditions of each microorganism. Furthermore, we selected 42 physico-chemical properties of amino acids and computed their values for each proteome. Further, on the same set of features we applied two fundamentally different machine learning methods, Support Vector Machines and Random Forests, to successfully classify between bacteria and archaea, halophiles and non-halophiles, as well as mesophiles, thermophiles and mesothermophiles. Finally, we performed feature selection by using Random Forests.</p> <p>Conclusions</p> <p>To our knowledge, this is the first time that three different classification cases (domain of life, halophilicity and thermophilicity) of proteome adaptation are successfully performed with the same set of 42 features. The characteristic features of a specific adaptation constitute a signature that may help understanding the mechanisms of adaptation to extreme environments.</p

Translational Selection Is Ubiquitous in Prokaryotes

Codon usage bias in prokaryotic genomes is largely a consequence of background substitution patterns in DNA, but highly expressed genes may show a preference towards codons that enable more efficient and/or accurate translation. We introduce a novel approach based on supervised machine learning that detects effects of translational selection on genes, while controlling for local variation in nucleotide substitution patterns represented as sequence composition of intergenic DNA. A cornerstone of our method is a Random Forest classifier that outperformed previous distance measure-based approaches, such as the codon adaptation index, in the task of discerning the (highly expressed) ribosomal protein genes by their codon frequencies. Unlike previous reports, we show evidence that translational selection in prokaryotes is practically universal: in 460 of 461 examined microbial genomes, we find that a subset of genes shows a higher codon usage similarity to the ribosomal proteins than would be expected from the local sequence composition. These genes constitute a substantial part of the genome—between 5% and 33%, depending on genome size—while also exhibiting higher experimentally measured mRNA abundances and tending toward codons that match tRNA anticodons by canonical base pairing. Certain gene functional categories are generally enriched with, or depleted of codon-optimized genes, the trends of enrichment/depletion being conserved between Archaea and Bacteria. Prominent exceptions from these trends might indicate genes with alternative physiological roles; we speculate on specific examples related to detoxication of oxygen radicals and ammonia and to possible misannotations of asparaginyl–tRNA synthetases. Since the presence of codon optimizations on genes is a valid proxy for expression levels in fully sequenced genomes, we provide an example of an “adaptome” by highlighting gene functions with expression levels elevated specifically in thermophilic Bacteria and Archaea

Revisiting the Myths of Protein Interior: Studying Proteins with Mass-Fractal Hydrophobicity-Fractal and Polarizability-Fractal Dimensions

A robust marker to describe mass, hydrophobicity and polarizability distribution holds the key to deciphering structural and folding constraints within proteins. Since each of these distributions is inhomogeneous in nature, the construct should be sensitive in describing the patterns therein. We show, for the first time, that the hydrophobicity and polarizability distributions in protein interior follow fractal scaling. It is found that (barring ‘all-α’) all the major structural classes of proteins have an amount of unused hydrophobicity left in them. This amount of untapped hydrophobicity is observed to be greater in thermophilic proteins, than that in their (structurally aligned) mesophilic counterparts. ‘All-β’(thermophilic, mesophilic alike) proteins are found to have maximum amount of unused hydrophobicity, while ‘all-α’ proteins have been found to have minimum polarizability. A non-trivial dependency is observed between dielectric constant and hydrophobicity distributions within (α+β) and ‘all-α’ proteins, whereas absolutely no dependency is found between them in the ‘all-β’ class. This study proves that proteins are not as optimally packed as they are supposed to be. It is also proved that origin of α-helices are possibly not hydrophobic but electrostatic; whereas β-sheets are predominantly hydrophobic in nature. Significance of this study lies in protein engineering studies; because it quantifies the extent of packing that ensures protein functionality. It shows that myths regarding protein interior organization might obfuscate our knowledge of actual reality. However, if the later is studied with a robust marker of strong mathematical basis, unknown correlations can still be unearthed; which help us to understand the nature of hydrophobicity, causality behind protein folding, and the importance of anisotropic electrostatics in stabilizing a highly complex structure named ‘proteins’

How Quickly Do Proteins Fold and Unfold, and What Structural Parameters Correlate with These Values?

The correlations between the logarithm of the unfolding rate of 108 proteins and their structural parameters were calculated. We showed that there is a good correlation between the logarithm of folding rates (in native conditions) and unfolding rates (in denaturing conditions) (0.79) and protein stability and unfolding rate (0.79). Thus, the faster the protein folds, the faster it unfolds. Folding and unfolding rates are higher for the proteins with two-state kinetics, in comparison with the proteins with multi-state kinetics. At the same time, two-state bacterial proteins folds and unfolds two orders of magnitude faster than two-state eukaryotic proteins, and multi-state bacterial proteins folds and unfolds slower than multi-state eukaryotic proteins. Despite the fact that the folding rates of thermophilic and mesophilic proteins are close, the unfolding rates of thermophilic proteins is about two orders of magnitude lower than for mesophilic proteins. The correlation between unfolding rate and stability of thermophilic proteins is high (0.90). We also found that the unfolding rate correlates with such structural parameters as: size of the protein, radius of the cross-section, logarithm of absolute contact order, and radius of gyration. This information will be useful for engineering and designing new proteins with desired properties

Dataset of the molecular dynamics simulations of bilayers consisting of short amyloidogenic peptide VDSWNVLVAG from Bgl2p–glucantransferase of S. cerevisiae cell wall

The amyloidogenic peptide VDSWNVLVAG from Bgl2p–glucantransferase of Saccharomyces cerevisiae cell wall and its modifying analog VESWNVLVAG were taken for the construction of four types of bilayers which differ by orientation of the peptides in the layers and of the layers relative to each other. These bilayers were used as starting models for the molecular dynamics (MD) at three charge states (neutral, pH3, and pH5). The changes of the fraction of secondary structure during 1 ns simulations were received for 96 MD trajectories. The data article contains the necessary information for the construction of models of β-strands organization in the oligomer structure. These results were used in the associated research article “Structural model of amyloid fibrils for amyloidogenic peptide from Bgl2p–glucantransferase of S. cerevisiae cell wall and its modifying analog. New morphology of amyloid fibrils” (Selivanova et al., 2016) [1]