Sequence and structural analysis of two designed proteins with 88% identity adopting different folds

Protein folding is a natural phenomenon by which a sequence of amino acids folds into a unique functional three-dimensional structure. Although the sequence code that governs folding remains a mystery, one can identify key inter-residue contacts responsible for a given topology. In nature, there are many pairs of proteins of a given length that share little or no sequence identity. Similarly, there are many proteins that share a common topology but lack significant evidence of homology. In order to tackle this problem, protein engineering studies have been used to determine the minimal number of amino acid residues that codes for a particular fold. In recent years, the coupling of theoretical models and experiments in the study of protein folding has resulted in providing some fruitful clues. He et al. have designed two proteins with 88% sequence identity, which adopt different folds and functions. In this work, we have systematically analysed these two proteins by performing pentapeptide search, secondary structure predictions, variation in inter-residue interactions and residue-residue pair preferences, surrounding hydrophobicity computations, conformational switching and energy computations. We conclude that the local secondary structural preference of the two designed proteins at the Nand C-terminal ends to adopt either coil or strand conformation may be a crucial factor in adopting the different folds. Early on during the process of folding, both proteins may choose different energetically favourable pathways to attain the different folds.     

Origin of Increased Solvent Accessibility of Peptide Bonds in Mutual Synergetic Folding Proteins

Mutual Synergetic Folding (MSF) proteins belong to a recently discovered class of proteins. These proteins are disordered in their monomeric but ordered in their oligomeric forms. Their amino acid composition is more similar to globular proteins than to disordered ones. Our preceding work shed light on important structural aspects of the structural organization of these proteins, but the background of this behavior is still unknown. We suggest that solvent accessibility is an important factor, especially solvent accessibility of the peptide bonds can be accounted for this phenomenon. The side chains of the amino acids which form a peptide bond have a high local contribution to the shielding of the peptide bond from the solvent. During the oligomerization step, other non-local residues contribute to the shielding. We investigated these local and non-local effects of shielding based on Shannon information entropy calculations. We found that MSF and globular homodimeric proteins have different local contributions resulting from different amino acid pair frequencies. Their non-local distribution is also different because of distinctive inter-subunit contacts.     

Classification of G-protein coupled receptors at four levels

G-protein coupled receptors (GPCRs) are transmembrane proteins which via G-proteins initiate some of the important signaling pathways in a cell and are involved in various physiological processes. Thus, computational prediction and classification of GPCRs can supply significant information for the development of novel drugs in pharmaceutical industry. In this paper, a nearest neighbor method has been introduced to discriminate GPCRs from non-GPCRs and subsequently classify GPCRs at four levels on the basis of amino acid composition and dipeptide composition of proteins. Its performance is evaluated on a non-redundant dataset consisted of 1406 GPCRs for six families and 1406 globular proteins using the jackknife test. The present method based on amino acid composition achieved an overall accuracy of 96.4% and Matthew's correlation coefficient (MCC) of 0.930 for correctly picking out the GPCRs from globular proteins. The overall accuracy and MCC were further enhanced to 99.8% and 0.996 by dipeptide composition-based method. On the other hand, the present method has successfully classified 1406 GPCRs into six families with an overall accuracy of 89.6 and 98.8% using amino acid composition and dipeptide composition, respectively. For the subfamily prediction of 1181 GPCRs of rhodopsin-like family, the present method achieved an overall accuracy of 76.7 and 94.5% based on the amino acid composition and dipeptide composition, respectively. Finally, GPCRs belonging to the amine subfamily and olfactory subfamily of rhodopsin-like family were further analyzed at the type level. The overall accuracy of dipeptide composition-based method for the classification of amine type and olfactory type of GPCRs reached 94.5 and 86.9%, respectively, while the overall accuracy of amino acid composition-based method was very low for both subfamilies. In comparison with existing methods in the literature, the present method also displayed great competitiveness. These results demonstrate the effectiveness of our method on identifying and classifying GPCRs correctly. GPCRsIdentifier, a corresponding stand-alone executable program for GPCR identification and classification was also developed, which can be acquired freely on request from the authors for academic purposes.     

Detection and reduction of evolutionary noise in correlated mutation analysis

Direct or indirect inter-residue interactions in proteins are often reflected by mutations at one site that compensate for mutations at another site. Various bioinformatic methods have been developed for detecting such correlated mutations in order to obtain information about intra- and inter-protein interactions. Here, we show by carrying out a correlated mutation analysis for non-interacting proteins that the signal due to inter-residue interactions is of similar magnitude to the 'noise' that arises from other evolutionary processes related to common ancestry. A new method for detecting correlated mutations is presented that reduces this evolutionary noise by taking into account evolutionary distances in the protein family. It is shown that this method yields better signal-to-noise ratios and, therefore, can much better resolve, for example, correlated mutations that reflect true inter-residue interactions.     

Orientation-dependent coarse-grained potentials derived by statistical analysis of molecular structural databases

We present results obtained for anisotropic potentials for protein simulations extracted from the continually growing databases of protein structures. This work is based on the assumption that the detailed information on molecular conformations can be used to derive statistical (a.k.a. ‘knowledge-based’) potentials that can describe on a coarse-grained level the side chain-side chain interactions in peptides and proteins. The complexity of inter-residue interactions is reflected in a high degree of orientational anisotropy for the twenty amino acids. By including in this coarse-grained interaction model the possibility of quantifying the backbone-backbone and backbone-side chain interactions, important improvements are obtained in characterizing the native protein states. Results obtained from tests that involve the identification of native-like conformations from large sets of decoy structures are presented. The method for deriving orientation-dependent statistical potentials is also applied to obtain water-water interactions. Monte Carlo simulations using the new coarse-grained water model show that the locations of the minima and maxima of the oxygen-oxygen radial distribution function correspond well with experimental measurements.     

In Silico Modeling of the Influence of Environment on Amyloid Folding Using FOD-M Model

The role of the environment in amyloid formation based on the fuzzy oil drop model (FOD) is discussed here. This model assumes that the hydrophobicity distribution within a globular protein is consistent with a 3D Gaussian (3DG) distribution. Such a distribution is interpreted as the idealized effect of the presence of a polar solvent—water. A chain with a sequence of amino acids (which are bipolar molecules) determined by evolution recreates a micelle-like structure with varying accuracy. The membrane, which is a specific environment with opposite characteristics to the polar aquatic environment, directs the hydrophobic residues towards the surface. The modification of the FOD model to the FOD-M form takes into account the specificity of the cell membrane. It consists in “inverting” the 3DG distribution (complementing the Gaussian distribution), which expresses the exposure of hydrophobic residues on the surface. It turns out that the influence of the environment for any protein (soluble or membrane-anchored) is the result of a consensus factor expressing the participation of the polar environment and the “inverted” environment. The ratio between the proportion of the aqueous and the “reversed” environment turns out to be a characteristic property of a given protein, including amyloid protein in particular. The structure of amyloid proteins has been characterized in the context of prion, intrinsically disordered, and other non-complexing proteins to cover a wider spectrum of molecules with the given characteristics based on the FOD-M model.     

Physicochemical factors for discriminating between soluble and membrane proteins: hydrophobicity of helical segments and protein length

The average hydrophobicity of a polypeptide segment is consideredto be the most important factor in the formation of transmembranehelices, and the partitioning of the most hydrophobic (MH) segmentinto the alternative nonpolar environment, a membrane or hydrophobiccore of a globular protein may determine the type of proteinproduced. In order to elucidate the importance of the MH segmentin determining which of the two types of protein results froma given amino acid sequence, we statistically studied the characteristicsof MH helices, longer than 19 residues in length, in 97 membraneproteins whose three-dimensional structure or topology is known,as well as 397 soluble proteins selected from the Protein DataBank. The average hydrophobicity of MH helices in membrane proteinshad a characteristic relationship with the length of the protein.All MH helices in membrane proteins that were longer than 500residues had a hydrophobicity greater than 1.75 (Kyte and Doolittlescale), while the MH helices in membrane proteins smaller than100 residues could be as hydrophilic as 0.1. The possibilityof developing a method to discriminate membrane proteins fromsoluble ones, based on the effect of size on the type of proteinproduced, is discussed.     

Folding and Insertion of Transmembrane Helices at the ER

In eukaryotic cells, the endoplasmic reticulum (ER) is the entry point for newly synthesized proteins that are subsequently distributed to organelles of the endomembrane system. Some of these proteins are completely translocated into the lumen of the ER while others integrate stretches of amino acids into the greasy 30 Å wide interior of the ER membrane bilayer. It is generally accepted that to exist in this non-aqueous environment the majority of membrane integrated amino acids are primarily non-polar/hydrophobic and adopt an α-helical conformation. These stretches are typically around 20 amino acids long and are known as transmembrane (TM) helices. In this review, we will consider how transmembrane helices achieve membrane integration. We will address questions such as: Where do the stretches of amino acids fold into a helical conformation? What is/are the route/routes that these stretches take from synthesis at the ribosome to integration through the ER translocon? How do these stretches ‘know’ to integrate and in which orientation? How do marginally hydrophobic stretches of amino acids integrate and survive as transmembrane helices?     

Factors influencing the thermal stability of buried protein mutants

Hydrophobic interaction is believed to be the most important factor for the stability of proteins upon buried mutations. In this work, we have analyzed the influence of different interactions to the stability of buried protein mutants by means of 49 various physical-chemical, energetic and conformational properties of amino acid residues. We found that the mutant stability is attributed with several factors including hydrophobicity. In lysozyme T4, the properties reflecting hydrophobicity, flexibility, turn and coil tendency, and long-range interactions show a strong correlation with stability. Entropy plays an important role and the contribution of hydrophobicity is minimal in barnase. The stability of human lysozyme is attributed with both hydrophobicity and secondary structure. The stability of buried mutants in staphylococcal nuclease is influenced by hydrophobicity and physical properties. Our results indicate that the stability of buried protein mutants are influenced not only with hydrophobicity but also other factors, such as, secondary structure, shape, flexibility, entropy and inter-residue contacts play an important role to the stability. We obtained the highest single property correlation of 0.83 between amino acid properties and thermal stability of buried protein mutants. The properties showing high correlation coefficient with thermal stability agree very well with experimental observations. Further, multiple regression technique combining three properties leads to the correlation in the range of 0.83-0.92 in the considered proteins.     

Assessment and comparison of distributed model predictive control schemes: Application to a natural gas refrigeration plant

A number of decentralized and distributed control schemes based on model predictive control (MPC) have been introduced in the last years. They have been proposed as viable solutions to the computational, transmission and robustness issues arising in the centralized context in case of large-scale and/or distributed plants. Such MPC-based control schemes are very heterogeneous, based on different model structures and realizations, with different features and infrastructural/memory/computational requirements.In this paper, we test and compare, with a realistic case study, a robust non-cooperative scheme and a cooperative iterative one. The main scope is to analyze and unravel, in a fair comparison scenario, these methods from different viewpoints, spanning from the model realization issues to the communication and computational requirements, to the control performances. The benchmark case study consists of an existing natural gas refrigeration plant. Realistic simulations and validation tests are obtained through in the DynSim industrial process simulation environment.     

Protein Interactions with Nanoparticle Surfaces: Highlighting Solution NMR Techniques

In the last decade, nanoparticles (NPs) have become a key tool in medicine and biotechnology as drug delivery systems, biosensors and diagnostic devices. The composition and surface chemistry of NPs vary based on the materials used: typically organic polymers, inorganic materials, or lipids. Nanoparticle classes can be further divided into sub‐categories depending on the surface modification and functionalization. These surface properties matter when NPs are introduced into a physiological environment, as they will influence how nucleic acids, lipids, and proteins will interact with the NP surface. While small‐molecule interactions are easily probed using NMR spectroscopy, studying protein‐NP interactions using NMR introduces several challenges. For example, globular proteins may have a perturbed conformation when attached to a foreign surface, and the size of NP‐protein conjugates can lead to excessive line broadening. Many of these challenges have been addressed, and NMR spectroscopy is becoming a mature technique for in situ analysis of NP binding behavior. It is therefore not surprising that NMR has been applied to NP systems and has been used to study biomolecules on NP surfaces. Important considerations include corona composition, protein behavior, and ligand architecture. These features are difficult to resolve using classical surface and material characterization strategies, and NMR provides a complementary avenue of characterization. In this review, we examine how solution NMR can be combined with other analytical techniques to investigate protein behavior on NP surfaces.     

Assessing the Role of Lipids in the Molecular Mechanism of Membrane Proteins

Membrane proteins have evolved to work optimally within the complex environment of the biological membrane. Consequently, interactions with surrounding lipids are part of their molecular mechanism. Yet, the identification of lipid–protein interactions and the assessment of their molecular role is an experimental challenge. Recently, biophysical approaches have emerged that are compatible with the study of membrane proteins in an environment closer to the biological membrane. These novel approaches revealed specific mechanisms of regulation of membrane protein function. Lipids have been shown to play a role in oligomerization, conformational transitions or allosteric coupling. In this review, we summarize the recent biophysical approaches, or combination thereof, that allow to decipher the role of lipid–protein interactions in the mechanism of membrane proteins.     

Local polarity analysis: a sensitive method that discriminates between native proteins and incorrectly folded models

The evaluation of calculated protein structures is an importantstep in the protein design cycle. Known criteria for this assessmentof proteins are the polar and apolar, accessible and buriedsurface area, electrostatic interactions and other interactionsbetween the protein atoms (e.g. HO, S-S),atomic packing, analysisof amino acid environment and surface charge distribution. Weshow that a powerful test of accuracy of protein structure canbe derived by analysing the water contact of atoms and additionallytaking into account their polarity. On the basis of estimatedreference values of the polar fraction of typical globular proteinswith known structure (mean, SD and distribution), the evaluationof misfolded structures can be improved significantly. The referencevalues are derived by moving windows of different length (3–99amino acid residues) over the amino acid sequence. Model proteins,which are included in the Brookhaven protein structure databank,deliberately misfolded proteins, hypothetical proteins and predictedprotein structures are diagnosed as at least partially incorrectlyfolded. The local fault, mostly observed, is that polar groupsare buried too frequently in the interior of the protein. Thedatabase-derived quantities are useful in screening the designedproteins prior to experimentation and may also be useful inthe assessment of errors in the experimentally determined proteinstructures.     

Separation sequencing. use of information distance

A fast process synthesis procedure is developed for generating preliminary separation schemes. The main approach is based on similarity indices obtained from chemical, physical and operating properties between two chemical species in a separation scheme. For the whole set of species, a similarity matrix results. On this basis, algorithms for classifying chemical species and generating separation sequences are developed. A distance between flowsheets is used in order to compare different separation schemes with one another or to some reference cases. The weights of the properties and heuristics are adjusted to account for well-known optimal cases. Changes in the number of considered separation properties and in their hierarchy are also used to improve the procedure. The algorithm requires very little calculation and can even be implemented by hand in the simplest cases. Extensions of the algorithm include the use of actual values of the separation properties instead of boolean ones and the generation of non-sharp flowsheets.     

Contacting modes and behaviour classification of gas—solid and other two-phase suspensions

An attempt is made to extend schemes for classifying the behaviour of gas—solid contacting modes and other two-phase systems. The regime diagram approach of Reh (1971) is modified and extended to cover the operating regions of common reactors and contactors where a gas flows upwards through a bed of solids including fixed and moving packed beds, conventional fluidized beds, circulating beds, spouted beds, and pneumatically conveyed suspensions. New boundaries are proposed between groups A and B and between groups B and D of the Geldart (1972, 1973) powder classification scheme. These boundaries reflect more recent data and allow the classification scheme to be used for gases other than air and for temperatures and pressures other than atmospheric.     

Engineering Pyrrolysyl-tRNA Synthetase for the Incorporation of Non-Canonical Amino Acids with Smaller Side Chains

Site-specific incorporation of non-canonical amino acids (ncAAs) into proteins has emerged as a universal tool for systems bioengineering at the interface of chemistry, biology, and technology. The diversification of the repertoire of the genetic code has been achieved for amino acids with long and/or bulky side chains equipped with various bioorthogonal tags and useful spectral probes. Although ncAAs with relatively small side chains and similar properties are of great interest to biophysics, cell biology, and biomaterial science, they can rarely be incorporated into proteins. To address this gap, we report the engineering of PylRS variants capable of incorporating an entire library of aliphatic “small-tag” ncAAs. In particular, we performed mutational studies of a specific PylRS, designed to incorporate the shortest non-bulky ncAA (S-allyl-l-cysteine) possible to date and based on this knowledge incorporated aliphatic ncAA derivatives. In this way, we have not only increased the number of translationally active “small-tag” ncAAs, but also determined key residues responsible for maintaining orthogonality, while engineering the PylRS for these interesting substrates. Based on the known plasticity of PylRS toward different substrates, our approach further expands the reassignment capacities of this enzyme toward aliphatic amino acids with smaller side chains endowed with valuable functionalities.     

Engineering thermostability in serine protease inhibitors

Unlike most globular proteins, the native form of serine protease inhibitors (serpins) is strained. Previous studies of human alpha(1)-antitrypsin, a prototype plasma serpin, revealed that various unfavorable interactions, such as overpacking of side chains, buried polar groups and cavities, are the structural basis of the strain. The local strain could be relieved by various stabilizing single amino acid substitutions, which appeared to remove these unfavorable interactions. To improve the stability of other clinically important serpin members, here we examined whether the rules found in alpha(1)-antitrypsin studies are applicable to other serpins. Amino acid substitutions were introduced at various positions in human alpha(1)-antichymotrypsin and human antithrombin III that were equivalent to the sites of stabilizing substitutions of alpha(1)-antitrypsin. Two-thirds of the substitutions increased thermostability in all serpins tested. Mutational analysis and structural examination suggest that serpins are suboptimally folded using common structural strategies at many sites, even though some structural details can vary in individual members. The results suggest that schemes discovered with alpha(1)-antitrypsin, an easily manipulative serpin, are a useful basis for engineering conformational characteristics of other clinically important serpins.