Hydrophobic folding units derived from dissimilar monomer structures and their interactions

Understanding the Folding Mediated Assembly of the Bacteriophage MS2 Coat Protein Dimers

The capsids of RNA viruses such as MS2 are great models for studying protein self-assembly because they are made almost entirely of multiple copies of a single coat protein (CP). Although CP is the minimal repeating unit of the capsid, previous studies have shown that CP exists as a homodimer (CP2) even in an acid-disassembled system, indicating that CP2 is an obligate dimer. Here, we investigate the molecular basis of this obligate dimerization using coarse-grained structure-based models and molecular dynamics simulations. We find that, unlike monomeric proteins of similar size, CP populates a single partially folded ensemble whose "foldedness" is sensitive to denaturing conditions. In contrast, CP2 folds similarly to single-domain proteins populating only the folded and the unfolded ensembles, separated by a prominent folding free energy barrier. Several intramonomer contacts form early, but the CP2 folding barrier is crossed only when the intermonomer contacts are made. A dissection of the structure of CP2 through mutant folding simulations shows that the folding barrier arises both from the topology of CP and the interface contacts of CP2. Together, our results show that CP2 is an obligate dimer because of kinetic stability, that is, dimerization induces a folding barrier and that makes it difficult for proteins in the dimer minimum to partially unfold and access the monomeric state without completely unfolding. We discuss the advantages of this obligate dimerization in the context of dimer design and virus stability.

Intragenic antimicrobial peptides (IAPs) from human proteins with potent antimicrobial and anti-inflammatory activity

Following the treads of our previous works on the unveiling of bioactive peptides encrypted in plant proteins from diverse species, the present manuscript reports the occurrence of four proof-of-concept intragenic antimicrobial peptides in human proteins, named Hs IAPs. These IAPs were prospected using the software Kamal, synthesized by solid phase chemistry, and had their interactions with model phospholipid vesicles investigated by differential scanning calorimetry and circular dichroism. Their antimicrobial activity against bacteria, yeasts and filamentous fungi was determined, along with their cytotoxicity towards erythrocytes. Our data demonstrates that Hs IAPs are capable to bind model membranes while attaining α-helical structure, and to inhibit the growth of microorganisms at concentrations as low as 1μM. Hs02, a novel sixteen residue long internal peptide (KWAVRIIRKFIKGFIS-NH₂) derived from the unconventional myosin 1h protein, was further investigated in its capacity to inhibit lipopolysaccharide-induced release of TNF-α in murine macrophages. Hs02 presented potent anti-inflammatory activity, inhibiting the release of TNF-α in LPS-primed cells at the lowest assayed concentration, 0.1 μM. A three-dimensional solution structure of Hs02 bound to DPC micelles was determined by Nuclear Magnetic Resonance. Our work exemplifies how the human genome can be mined for molecules with biotechnological potential in human health and demonstrates that IAPs are actual alternatives to antimicrobial peptides as pharmaceutical agents or in their many other putative applications.

Folding with a protein's native shortcut network

A complex network approach to protein folding is proposed, wherein a protein's contact map is reconceptualized as a network of shortcut edges, and folding is steered by a structural characteristic of this network. Shortcut networks are generated by a known message passing algorithm operating on protein residue networks. It is found that the shortcut networks of native structures (SCN0s) are relevant graph objects with which to study protein folding at a formal level. The logarithm form of their contact order (SCN0_lnCO) correlates significantly with folding rate of two-state and nontwo-state proteins. The clustering coefficient of SCN0s (C_SCN0 ) correlates significantly with folding rate, transition-state placement and stability of two-state folders. Reasonable folding pathways for several model proteins are produced when C_SCN0 is used to combine protein segments incrementally to form the native structure. The folding bias captured by C_SCN0 is detectable in non-native structures, as evidenced by Molecular Dynamics simulation generated configurations for the fast folding Villin-headpiece peptide. These results support the use of shortcut networks to investigate the role protein geometry plays in the folding of both small and large globular proteins, and have implications for the design of multibody interaction schemes in folding models. One facet of this geometry is the set of native shortcut triangles, whose attributes are found to be well-suited to identify dehydrated intraprotein areas in tight turns, or at the interface of different secondary structure elements.

Binding Free Energy and the structural changes determination in hGH protein with different concentrations of copper ions (A molecular dynamics simulation study)

In this study, we estimated the optimum concentration of copper ions that are effective in the stability and the structural changes of human growth hormone (hGH) protein in the combination of different concentrations of these ions at the molecular level using molecular dynamics simulation by Gromacs 4.6.5 software. Moreover, to estimate the binding affinity of copper ions to hGH protein, binding free energies is calculated by the molecular mechanics Poisson–Boltzmann Surface Area (MM-PBSA). The analysis of molecular dynamics (MD) trajectories as dictionary of the secondary structure of protein (DSSP), solvent accessible surface area (SASA) and binding free energy calculations show that hGH protein structure is more stabilized by increasing a limited concentration of copper ions. These findings align with our previous experimental studies.

Sparse Markov chain-based semi-supervised multi-instance multi-label method for protein function prediction

Automated assignment of protein function has received considerable attention in recent years for genome-wide study. With the rapid accumulation of genome sequencing data produced by high-throughput experimental techniques, the process of manually predicting functional properties of proteins has become increasingly cumbersome. Such large genomics data sets can only be annotated computationally. However, automated assignment of functions to unknown protein is challenging due to its inherent difficulty and complexity. Previous studies have revealed that solving problems involving complicated objects with multiple semantic meanings using the multi-instance multi-label (MIML) framework is effective. For the protein function prediction problems, each protein object in nature may associate with distinct structural units (instances) and multiple functional properties (class labels) where each unit is described by an instance and each functional property is considered as a class label. Thus, it is convenient and natural to tackle the protein function prediction problem by using the MIML framework. In this paper, we propose a sparse Markov chain-based semi-supervised MIML method, called Sparse-Markov. A sparse transductive probability graph is constructed to encode the affinity information of the data based on ensemble of Hausdorff distance metrics. Our goal is to exploit the affinity between protein objects in the sparse transductive probability graph to seek a sparse steady state probability of the Markov chain model to do protein function prediction, such that two proteins are given similar functional labels if they are close to each other in terms of an ensemble Hausdorff distance in the graph. Experimental results on seven real-world organism data sets covering three biological domains show that our proposed Sparse-Markov method is able to achieve better performance than four state-of-the-art MIML learning algorithms.

Multi-instance multilabel learning with weak-label for predicting protein function in electricigens

Nature often brings several domains together to form multidomain and multifunctional proteins with a vast number of possibilities. In our previous study, we disclosed that the protein function prediction problem is naturally and inherently Multi-Instance Multilabel (MIML) learning tasks. Automated protein function prediction is typically implemented under the assumption that the functions of labeled proteins are complete; that is, there are no missing labels. In contrast, in practice just a subset of the functions of a protein are known, and whether this protein has other functions is unknown. It is evident that protein function prediction tasks suffer from weak-label problem; thus protein function prediction with incomplete annotation matches well with the MIML with weak-label learning framework. In this paper, we have applied the state-of-the-art MIML with weak-label learning algorithm MIMLwel for predicting protein functions in two typical real-world electricigens organisms which have been widely used in microbial fuel cells (MFCs) researches. Our experimental results validate the effectiveness of MIMLwel algorithm in predicting protein functions with incomplete annotation.

Genome-Wide Protein Function Prediction through Multi-Instance Multi-Label Learning

Automated annotation of protein function is challenging. As the number of sequenced genomes rapidly grows, the vast majority of proteins can only be annotated computationally. Nature often brings several domains together to form multi-domain and multi-functional proteins with a vast number of possibilities, and each domain may fulfill its own function independently or in a concerted manner with its neighbors. Thus, it is evident that the protein function prediction problem is naturally and inherently Multi-Instance Multi-Label (MIML) learning tasks. Based on the state-of-the-art MIML algorithm MIMLNN, we propose a novel ensemble MIML learning framework EnMIMLNN and design three algorithms for this task by combining the advantage of three kinds of Hausdorff distance metrics. Experiments on seven real-world organisms covering the biological three-domain system, i.e., archaea, bacteria, and eukaryote, show that the EnMIMLNN algorithms are superior to most state-of-the-art MIML and Multi-Label learning algorithms.

Salt-bridge energetics in halophilic proteins

Halophilic proteins have greater abundance of acidic over basic and very low bulky hydrophobic residues. Classical electrostatic stabilization was suggested as the key determinant for halophilic adaptation of protein. However, contribution of specific electrostatic interactions (i.e. salt-bridges) to overall stability of halophilic proteins is yet to be understood. To understand this, we use Adaptive-Poison-Boltzmann-Solver Methods along with our home-built automation to workout net as well as associated component energy terms such as desolvation energy, bridge energy and background energy for 275 salt-bridges from 20 extremely halophilic proteins. We then perform extensive statistical analysis on general and energetic attributes on these salt-bridges. On average, 8 salt-bridges per 150 residues protein were observed which is almost twice than earlier report. Overall contributions of salt-bridges are −3.0 kcal mol⁻¹. Majority (78%) of salt-bridges in our dataset are stable and conserved in nature. Although, average contributions of component energy terms are equal, their individual details vary greatly from one another indicating their sensitivity to local micro-environment. Notably, 35% of salt-bridges in our database are buried and stable. Greater desolvation penalty of these buried salt-bridges are counteracted by stable network salt-bridges apart from favorable equal contributions of bridge and background terms. Recruitment of extensive network salt-bridges (46%) with a net contribution of −5.0 kcal mol⁻¹ per salt-bridge, seems to be a halophilic design wherein favorable average contribution of background term (−10 kcal mol⁻¹) exceeds than that of bridge term (−7 kcal mol⁻¹). Interiors of proteins from halophiles are seen to possess relatively higher abundance of charge and polar side chains than that of mesophiles which seems to be satisfied by cooperative network salt-bridges. Overall, our theoretical analyses provide insight into halophilic signature in its specific electrostatic interactions which we hope would help in protein engineering and bioinformatics studies.

Kinetic Analysis of Guanidine Hydrochloride Inactivation of β-Galactosidase in the Presence of Galactose

Inactivation of purified β-Galactosidase was done with GdnHCl in the absence and presence of varying [galactose] at 50°C and at pH 4.5. Lineweaver-Burk plots of initial velocity data, in the presence and absence of guanidine hydrochloride (GdnHCl) and galactose, were used to determine the relevant K_m and V_max values, with p-nitrophenyl β-D-galactopyranoside (pNPG) as substrate, S. Plots of ln([P]_∞ − [P]_t) against time in the presence of GdnHCl yielded the inactivation rate constant, A. Plots of A versus [S] at different galactose concentrations were straight lines that became increasingly less steep as the [galactose] increased, showing that A was dependent on [S]. Slopes and intercepts of the 1/[P]_∞ versus 1/[S] yielded k₊₀ and k'₊₀, the microscopic rate constants for the free enzyme and the enzyme-substrate complex, respectively. Plots of k₊₀ and k'₊₀ versus [galactose] showed that galactose protected the free enzyme as well as the enzyme-substrate complex (only at the lowest and highest [galactose]) against GdnHCl inactivation. In the absence of galactose, GdnHCl exhibited some degree of non-competitive inhibition. In the presence of GdnHCl, galactose exhibited competitive inhibition at the lower [galactose] of 5 mM which changed to non-competitive as the [galactose] increased. The implications of our findings are further discussed.

COMPACTNESS DETERMINES PROTEIN FOLDING TYPE

We have demonstrated here that protein compactness, which we define as the ratio of the accessible surface area of a protein to that of the ideal sphere of the same volume, is one of the factors determining the mechanism of protein folding. Proteins with multi-state kinetics, on average, are more compact (compactness is 1.49 ± 0.02 for proteins within the size range of 101–151 amino acid residues) than proteins with two-state kinetics (compactness is 1.59 ± 0.03 for proteins within the same size range of 101–151 amino acid residues). We have shown that compactness for homologous proteins can explain both the difference in folding rates and the difference in folding mechanisms.

Comparison of the structural basis for thermal stability between archaeal and bacterial proteins

In this study, the structural basis for thermal stability in archaeal and bacterial proteins was investigated. There were many common factors that confer resistance to high temperature in both archaeal and bacterial proteins. These factors include increases in the Lys content, the bends and blanks of secondary structure, the Glu content of salt bridge; decreases in the number of main-side chain hydrogen bond and exposed surface area, and changes in the bends and blanks of amino acids. Certainly, the utilization of charged amino acids to form salt bridges is a primary factor. In both heat-resistant archaeal and bacterial proteins, most Glu and Asp participate in the formation of salt bridges. Other factors may influence either archaeal or bacterial protein thermostability, which includes the more frequent occurrence of shorter 3(10)-helices and increased hydrophobicity in heat-resistant archaeal proteins. However, there were increases in average helix length, the Glu content in salt bridges, temperature factors and decreases in the number of main-side chain hydrogen bonds, uncharged-uncharged hydrogen bonds, hydrophobicity, and buried and exposed polar surface area in heat-resistant bacterial proteins. Evidently, there are few similarities and many disparities between the heat-resistant mechanisms of archaeal and bacterial proteins.

Backbone flexibility controls the activity and specificity of a protein-protein interface: specificity in snake venom metalloproteases

Protein-protein interfaces have crucial functions in many biological processes. The large interaction areas of such interfaces show complex interaction motifs. Even more challenging is the understanding of (multi)specificity in protein-protein binding. Many proteins can bind several partners to mediate their function. A perfect paradigm to study such multispecific protein-protein interfaces are snake venom metalloproteases (SVMPs). Inherently, they bind to a variety of basement membrane proteins of capillaries, hydrolyze them, and induce profuse bleeding. However, despite having a high sequence homology, some SVMPs show a strong hemorrhagic activity, while others are (almost) inactive. We present computer simulations indicating that the activity to induce hemorrhage, and thus the capability to bind the potential reaction partners, is related to the backbone flexibility in a certain surface region. A subtle interplay between flexibility and rigidity of two loops seems to be the prerequisite for the proteins to carry out their damaging function. Presumably, a significant alteration in the backbone dynamics makes the difference between SVMPs that induce hemorrhage and the inactive ones.

Protein structure classification based on conserved hydrophobic residues

Protein folding is frequently guided by local residue interactions that form clusters in the protein core. The interactions between residue clusters serve as potential nucleation sites in the folding process. Evidence postulates that the residue interactions are governed by the hydrophobic propensities that the residues possess. An array of hydrophobicity scales has been developed to determine the hydrophobic propensities of residues under different environmental conditions. In this work, we propose a graph-theory-based data mining framework to extract and isolate protein structural features that sustain invariance in evolutionary-related proteins, through the integrated analysis of five well-known hydrophobicity scales over the 3D structure of proteins. We hypothesize that proteins of the same homology contain conserved hydrophobic residues and exhibit analogous residue interaction patterns in the folded state. The results obtained demonstrate that discriminatory residue interaction patterns shared among proteins of the same family can be employed for both the structural and the functional annotation of proteins. We obtained on the average 90 percent accuracy in protein classification with a significantly small feature vector compared to previous results in the area. This work presents an elaborate study, as well as validation evidence, to illustrate the efficacy of the method and the correctness of results reported.

Analysis of protein contacts into Protein Units

Three-dimensional structures of proteins are the support of their biological functions. Their folds are maintained by inter-residue interactions which are one of the main focuses to understand the mechanisms of protein folding and stability. Furthermore, protein structures can be composed of single or multiple functional domains that can fold and function independently. Hence, dividing a protein into domains is useful for obtaining an accurate structure and function determination. In previous studies, we enlightened protein contact properties according to different definitions and developed a novel methodology named Protein Peeling. Within protein structures, Protein Peeling characterizes small successive compact units along the sequence called protein units (PUs). The cutting done by Protein Peeling maximizes the number of contacts within the PUs and minimizes the number of contacts between them. This method is so a relevant tool in the context of the protein folding research and particularly regarding the hierarchical model proposed by George Rose. Here, we accurately analyze the PUs at different levels of cutting, using a non-redundant protein databank. Distribution of PU sizes, number of PUs or their accessibility are screened to determine their common and different features. Moreover, we highlight the preferential amino acid interactions inside and between PUs. Our results show that PUs are clearly an intermediate level between secondary structures and protein structural domains.

Intra-molecular chaperone: the role of the N-terminal in conformational selection and kinetic control

The vast majority of the proteins in nature are under thermodynamic control, consistent with the universally accepted notion that proteins exist in their thermodynamically most stable state. Yet, recently a number of examples of proteins whose fold is under kinetic control have come to light. Their functions and environments vary. The first among these are some proteases, discovered in the early 1990s. There, an N-terminal proregion is self-cleaved after the protein folded, leaving the remainder of the chain in a kinetically trapped state. A related scenario was observed for microcin J25, an antibacterial peptide. This peptide presents a trapped covalently knotted conformation. The third and the most recently discovered case is the multidrug-resistant transporter protein, P-glycoprotein. There, a synonymous 'silent' mutation leads to ribosome stalling with a consequent altered kinetically trapped state. Here we argue that in all three examples, the N-terminal plays the role of an intra-molecular chaperone, that is, the N-terminal conformation selects among all competing local conformations of a downstream segment. By providing a pattern, the N-terminal chaperone segment assists the protein folding process. If the N-terminal is subsequently cleaved, the protein can be under kinetic control, since it is trapped in a thermodynamically less-stable state.

Architectures and functional coverage of protein-protein interfaces

The diverse range of cellular functions is performed by a limited number of protein folds existing in nature. One may similarly expect that cellular functional diversity would be covered by a limited number of protein-protein interface architectures. Here, we present 8205 interface clusters, each representing a unique interface architecture. This data set of protein-protein interfaces is analyzed and compared with older data sets. We observe that the number of both biological and crystal interfaces increases significantly compared to the number of Protein Data Bank entries. Furthermore, we find that the number of distinct interface architectures grows at a much faster rate than the number of folds and is yet to level off. We further analyze the growth trend of the functional coverage by constructing functional interaction networks from interfaces. The functional coverage is also found to steadily increase. Interestingly, we also observe that despite the diversity of interface architectures, some are more favorable and frequently used, and of particular interest, are the ones that are also preferred in single chains.

Kinetic Folding and Assembly Mechanisms Differ for Two Homologous Heptamers

Here we investigate the time-resolved folding and assembly mechanism of the heptameric co-chaperonin protein 10 (cpn10) in vitro. The structure of cpn10 is conserved throughout nature: seven beta-barrel subunits are non-covalently assembled through beta-strand pairings in an overall doughnut-like shape. Kinetic folding/assembly experiments of chemically denatured cpn10 from Homo sapiens (hmcpn10) and Aquifex aeolicus (Aacpn10) were monitored by far-UV circular dichroism and fluorescence. We find the processes to be complex, involving several kinetic steps, and to differ between the mesophilic and hyper-thermophilic proteins. The hmcpn10 molecules partition into two parallel pathways, one involving polypeptide folding before protein-protein assembly and another in which inter-protein interactions take place prior to folding. In contrast, the Aacpn10 molecules follow a single sequential path that includes initial monomer misfolding, relaxation to productive intermediates and, subsequently, final folding and heptamer assembly. An A. aeolicus variant lacking the unique C-terminal extension of Aacpn10 displays the same kinetic mechanism as Aacpn10, signifying that the tail is not responsible for the rapid misfolding step. This study demonstrates that molecular details can overrule similarity of native-state topology in defining apparent protein-biophysical properties.

Protein Peeling 2: a web server to convert protein structures into series of protein units

Protein Peeling 2 (PP2) is a web server for the automatic identification of protein units (PUs) given the 3D coordinates of a protein. PUs are an intermediate level of protein structure description between protein domains and secondary structures. It is a new tool to better understand and analyze the organization of protein structures. PP2 uses only the matrices of protein contact probabilities and cuts the protein structures optimally using Matthews' coefficient correlation. An index assesses the compactness quality of each PU. Results are given both textually and graphically using JMol and PyMol softwares. The server can be accessed from .

Hydrophobic collapse in late-stage folding (in silico) of bovine pancreatic trypsin inhibitor

Hydrophobic collapse is commonly considered as a process of significance for protein folding. Many models have been applied for description of this event. A model introducing an external force field mimicking the hydrophobic environment and simultaneously the driving force for the folding process is presented in this paper. Bovine pancreatic trypsin inhibitor (BPTI) was taken as a test protein. An early-stage folding (in silico) model presented elsewhere was used to create the starting structure for hydrophobic collapse. The resulting structure was energy-refined using molecular dynamics simulation in an explicit solvent. The similarity versus the crystal structure of BPTI is estimated using visual analysis, residue-residue contacts, phi, psi angle distributions, RMSD, accessible solvent area, radii of gyration and hydrodynamic radii. A program allowing creation of early-stage folding structural forms to be created for any protein is available from http://bioinformatics.cm-uj.krakow.pl/earlystage. The program for late-stage folding simulation is available on request.

Mapping the Distribution of Conformational Information Throughout a Protein Sequence

The three-dimensional structure of protein is encoded in the sequence, but many amino acid residues carry no essential conformational information, and the identity of those that are structure-determining is elusive. By circular permutation and terminal deletion, we produced and purified 25 Bacillus licheniformis beta-lactamase (ESBL) variants that lack 5-21 contiguous residues each, and collectively have 82% of the sequence and 92% of the non-local atom-atom contacts eliminated. Circular dichroism and size-exclusion chromatography showed that most of the variants form conformationally heterogeneous mixtures, but by measuring catalytic constants, we found that all populate, to a greater or lesser extent, conformations with the essential features of the native fold. This suggests that no segment of the ESBL sequence is essential to the structure as a whole, which is congruent with the notion that local information and modular organization can impart most of the tertiary fold specificity and cooperativity.

'Protein Peeling': an approach for splitting a 3D protein structure into compact fragments

MOTIVATION

The object of this study is to propose a new method to identify small compact units that compose protein three-dimensional structures. These fragments, called 'protein units (PU)', are a new level of description to well understand and analyze the organization of protein structures. The method only works from the contact probability matrix, i.e. the inter Calpha-distances translated into probabilities. It uses the principle of conventional hierarchical clustering, leading to a series of nested partitions of the 3D structure. Every step aims at dividing optimally a unit into 2 or 3 subunits according to a criterion called 'partition index' assessing the structural independence of the subunits newly defined. Moreover, an entropy-derived squared correlation R is used for assessing globally the protein structure dissection. The method is compared to other splitting algorithms and shows relevant performance.

AVAILABILITY

An Internet server with dedicated tools is available at http://www.ebgm.jussieu.fr/~gelly/

Recent developments in structural proteomics for protein structure determination

The major challenges in structural proteomics include identifying all the proteins on the genome-wide scale, determining their structure-function relationships, and outlining the precise three-dimensional structures of the proteins. Protein structures are typically determined by experimental approaches such as X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. However, the knowledge of three-dimensional space by these techniques is still limited. Thus, computational methods such as comparative and de novo approaches and molecular dynamic simulations are intensively used as alternative tools to predict the three-dimensional structures and dynamic behavior of proteins. This review summarizes recent developments in structural proteomics for protein structure determination; including instrumental methods such as X-ray crystallography and NMR spectroscopy, and computational methods such as comparative and de novo structure prediction and molecular dynamics simulations.

Potential folding-function interrelationship in proteins

The possibility is addressed that protein folding and function may be related via regions that are critical for both folding and function. This approach is based on the building blocks folding model that describes protein folding as binding events of conformationally fluctuating building blocks. Within these, we identify building block fragments that are critical for achieving the native fold. A library of such critical building blocks (CBBs) is constructed. Then, it is asked whether the functionally important residues fall in these CBB fragments. We find that for over two-thirds of the proteins in our library with available functional information, the catalytic or binding site residues lie within the CBB regions. From the evolutionary standpoint, a folding-function relationship is advantageous, since the need to guard against mutations is limited to one region. Furthermore, conformationally similar CBBs are found in globally unrelated proteins with different functions. Hence, substituting CBBs may lead to designed proteins with altered functions. We further find that the CBBs in our library are conformationally unstable.

Different roles of electrostatics in heat and in cold: adaptation by citrate synthase

Electrostatics plays a major role in heat adaptation by thermophilic proteins. Here we ask whether electrostatics similarly contributes to cold adaptation in psychrophilic proteins. We compare the sequences and structures of citrate synthases from the psychrophile Arthobacter Ds2-3R, from chicken, and from the hyperthermophile Pyrococcus furiosus. The three enzymes share similar packing, burial of nonpolar surface area, and main-chain hydrogen bonding. However, both psychrophilic and hyperthermophilic citrate synthases contain more charged residues, salt bridges, and salt-bridge networks than the mesophile. The electrostatic free-energy contributions toward protein stability by individual charged residues show greater variabilities in the psychrophilic citrate synthase than in the hyperthermophilic enzyme. The charged residues in the active-site regions of the psychrophile are more destabilizing than those in the active-site regions of the hyperthermophile. In the hyperthermophilic enzyme, salt bridges and their networks largely cluster in the active-site regions and at the dimer interface. In contrast, in the psychrophile, they are more dispersed throughout the structure. On average, salt bridges and their networks provide greater electrostatic stabilization to the thermophilic citrate synthase at 100 degrees C than to the psychrophilic enzyme at 0 degrees C. Electrostatics appears to play an important role in both heat and cold adaptation of citrate synthase. However, remarkably, the role may be different in the two types of enzyme: In the hyperthermophile, it may contribute to the integrity of both the protein dimer and the active site by possibly countering conformational disorder at high temperatures. On the other hand, in the psychrophile at low temperatures, electrostatics may contribute to enhance protein solvation and to ensure active-site flexibility.

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.

Analysis of the Information Structure of Protein Sequences: A New Method for Analyzing the Domain Organization of Proteins

The amino acid sequences of gamma-crystallin, Haloalkane Dehalogenase, Phthalate Dioxygenase, Porphobilinogen Deaminase and Myosin Regulatory Domain c-chain were analyzed for their information content. Sites of increased degree of information coordination between residues (IDIC-sites) were identified, and their organization was studied by means of analyzing the information structure of the protein sequences. Relationships between the structural units forming the spatial and informational structure of proteins were demonstrated. Associations of information-coordinated structural elements (IDIC-associations) were mapped onto compact structural domains found in the spatial structures of globular proteins. The proposed method of analyzing the information structure of protein sequences may find applications in the biotechnology and structural chemistry of proteins.

Atomic contact vectors in protein-protein recognition

The ability to analyze and compare protein-protein interactions on the structural level is critical to our understanding of various aspects of molecular recognition and the functional interplay of components of biochemical networks. In this study, we introduce atomic contact vectors (ACVs) as an intuitive way to represent the physico-chemical characteristics of a protein-protein interface as well as a way to compare interfaces to each other. We test the utility of ACVs in classification by using them to distinguish between homodimers and crystal contacts. Our results compare favorably with those reported by other authors. We then apply ACVs to mine the PDB for all known protein-protein complexes and separate transient recognition complexes from permanent oligomeric ones. Getting at the basis of this difference is important for our understanding of recognition and we achieved a success rate of 91% for distinguishing these two classes of complexes. Although accessible surface area of the interface is a major discriminating feature, we also show that there are distinct differences in the contact preferences between the two kinds of complexes. Illustrating the superiority of ACVs as a basic comparison measure over a sequence-based approach, we derive a general rule of thumb to determine whether two protein-protein interfaces are redundant. With this method, we arrive at a nonredundant set of 209 recognition complexes--the largest set reported so far.

Energy landscape and dynamics of the beta-hairpin G peptide and its isomers: Topology and sequences

We have investigated free energy landscape [MM/PBSA + normal modes entropy] of permutations in the G peptide (41-56) from the protein G B1 domain by studying six isomers corresponding to moving the hydrophobic cluster along the beta-strands (toward the turn: T1, AGEWTYDDKTFTVTET; T2, GEDTWDYATFTVTKTE; T3, GEDDWTYATFTVTKTE; toward the end: E1, WTYDDAGETKTFTVT; E2, WEYTGDDATKTETFTV; E3, WTYEGDDATKTETFTV). The free energy terms include molecular mechanics energy, Poisson-Boltzmann electrostatic solvation energy, surface area solvation energy, and conformational entropy estimated by using normal mode analysis. From the wild type to T1, then T3, and finally T2, we see a progressively changing energy landscape, toward a less stable beta-hairpin structure. Moving the hydrophobic cluster outside toward the end region causes a greater change in the energy landscape. alpha-Helical instead of a beta-hairpin structure was the most stable form for the E2 isomer. However, no matter how much the sequence changes, for all variants studied, ideal "native" beta-hairpin topologies remain as minima (regardless of whether global or local) in the energy landscape. In general, we find that the energy landscape is dependent on the hydrophobic cluster topology and on the sequence. Our present study indicates that the key is the relative conformational energies of the different conformations. Changes in the sequence strongly modulate the relative stabilities of topologically similar regions in the energy landscape, rather than redefine the topology space. This finding is consistent with a population redistribution in the process of protein folding. The limited variation of topological space, compared with the number of possible sequence changes, may relate to the observation that the number of known protein folds are far less than the sequential allowance.

Comparison of protein fragments identified by limited proteolysis and by computational cutting of proteins

Here we present a comparison between protein fragments produced by limited proteolysis and those identified by computational cutting based on the building block folding model. The principles upon which the two methods are based are different. Limited proteolysis of natively folded proteins occurs at flexible sites and never at the level of chain segments of regular secondary structure such as alpha-helices. Therefore, the targets for limited proteolysis are locally unfolded regions. In contrast, the computational cutting algorithm considers the compactness of the fragments, their nonpolar buried surface area, and their isolatedness, that is, the surface area which was buried prior to the cutting and becomes exposed subsequently. Despite the different criteria, there is an overall correspondence between sites or regions of limited proteolysis with those identified by computational cutting. The computational cutting method has been applied to several model proteins for which detailed limited proteolysis data are available, namely apomyoglobin, cytochrome c, ribonuclease A, alpha-lactalbumin, and thermolysin. As expected, more cuts are obtained computationally than experimentally and the agreement is better when a number of proteolytic enzymes are used. For example, cytochrome c is cleaved by thermolysin at 56-57, 45-46, and at 80-81, and by proteinase K at 48-49 and 50-51. Incubation of the noncovalent and native-like complex of cytochrome c fragments 1-56 and 57-104 with proteinase K yielded the gapped protein species 1-48/57-104 and finally 1-40/57-104. Computational cutting of cytochrome c reproduced the major experimental observations, with cuts at 47, 64-65 or 65-66 and 80-81 and an unstable 32-47 region not assigned to any building block. The next step, not addressed in this work, is to probe the ability of the generated fragments to fold independently. Since both the computational algorithm and limited proteolysis attempt to dissect the protein folding problem, the general agreement between the two procedures is gratifying. This consistency allows us to propose the use of limited proteolysis to produce protein fragments that can adopt an independent folding and, therefore, to study folding intermediates. The results of the present study appear to validate the building block folding model and are in line with the proposal that protein folding is a hierarchical process, where parts constituting local minima of energy fold first, with their subsequent association and mutual stabilization to finally yield the global fold.

Interdomain interactions in hinge-bending transitions

BACKGROUND

The mechanisms that allow or constrain protein movement have not been understood. Here we study interdomain interactions in proteins to investigate hinge-bending motions.

RESULTS

We find a limited number of salt bridges and hydrogen bonds at the interdomain interface, in both the "closed" and the "open" conformations. Consistently, analysis of 222 salt bridges in an independently selected database indicates that most salt bridges form within rather than between independently folding hydrophobic units. Calculations show that these interdomain salt bridges either destabilize or only marginally stabilize the closed conformation in most proteins. In contrast, the nonpolar buried surface area between the moving parts can be extensive in the closed conformations. However, when the nonpolar buried surface area is large, we find that at the interdomain interface in the open conformation it may be as large or larger than in the closed conformation. Hence, the energetic penalty of opening the closed conformation is overcome. Consistently, a large nonpolar surface area buried in the closed interdomain interface accompanies limited opening of the domains, yielding a larger interface.

CONCLUSIONS

Short-range electrostatic interactions are largely absent between moving domains. Interdomain nonpolar buried surface area may be large in the closed conformation, but it is largely offset by the area buried in the open conformation. In such cases the opening of the domains appears to be relatively small. This may allow prediction of the extent of domain opening. Such predictions may have implications for the shape and size of the binding pockets in drug/protein design.

Building blocks, hinge-bending motions and protein topology

Here we show that the locations of molecular hinges in protein structures fall between building block elements. Building blocks are fragments of the protein chain which constitute local minima. These elements fold first. In the next step they associate through a combinatorial assembly process. While chain-linked building blocks may be expected to trial-associate first, if unstable, alternate more stable associations will take place. Hence, we would expect that molecular hinges will be at such inter-building block locations, or at the less stable, unassigned regions. On the other hand, hinge-bending motions are well known to be critical for protein function. Hence, protein folding and protein function are evolutionarily related. Further, the pathways through which proteins attain their three dimensional folds are determined by protein topology. However, at the same time the locations of the hinges, and hinge-bending motions are also an outcome of protein topology. Thus, protein folding and function appear coupled, and relate to protein topology. Here we provide some results illustrating such a relationship.

Very empirical treatment of solvation and entropy: a force field derived from log Po/w

A non-covalent interaction force field model derived from the partition coefficient of 1-octanol/water solubility is described. This model, HINT for Hydropathic INTeractions, is shown to include, in very empirical and approximate terms, all components of biomolecular associations, including hydrogen bonding, Coulombic interactions, hydrophobic interactions, entropy and solvation/desolvation. Particular emphasis is placed on: (1) demonstrating the relationship between the total empirical HINT score and free energy of association, deltaGinteraction; (2) showing that the HINT hydrophobic-polar interaction sub-score represents the energy cost of desolvation upon binding for interacting biomolecules; and (3) a new methodology for treating constrained water molecules as discrete independent small ligands. An example calculation is reported for dihydrofolate reductase (DHFR) bound with methotrexate (MTX). In that case the observed very tight binding, deltaGinteraction < or = -13.6 kcal mol(-1), is largely due to ten hydrogen bonds between the ligand and enzyme with estimated strength ranging between -0.4 and -2.3 kcal mol(-1). Four water molecules bridging between DHFR and MTX contribute an additional -1.7 kcal mol(-1) stability to the complex. The HINT estimate of the cost of desolvation is +13.9 kcal mol(-1).

The building block folding model and the kinetics of protein folding

Here we show that qualitatively, the building blocks folding model accounts for three-state versus the two-state protein folding. Additionally, it is consistent with the faster versus slower folding rates of the two-state proteins. Specifically, we illustrate that the building blocks size, their mode of associations in the native structure, the number of ways they can combinatorially assemble, their population times and the way they are split in the iterative, step-by-step structural dissection which yields the anatomy trees, explain a broad range of folding rates. We further show that proteins with similar general topologies may have different folding pathways, and hence different folding rates. On the other hand, the effect of mutations resembles that of changes in conditions, shifting the population times and hence the energy landscapes. Hence, together with the secondary structure type and the extent of local versus non-local interactions, a coherent, consistent rationale for folding kinetics can be outlined, in agreement with experimental results. Given the native structure of a protein, these guidelines enable a qualitative prediction of the folding kinetics. We further describe these in the context of the protein folding energy landscape. Quantitatively, in principle, the diffusion-collision model for the building block association can be used. However, the folding rates of the building blocks and traps in their formation and association, need to be considered.

Fluctuations in ion pairs and their stabilities in proteins

This report investigates the effect of systemic protein conformational flexibility on the contribution of ion pairs to protein stability. Toward this goal, we use all NMR conformer ensembles in the Protein Data Bank (1) that contain at least 40 conformers, (2) whose functional form is monomeric, (3) that are nonredundant, and (4) that are large enough. We find 11 proteins adhering to these criteria. Within these proteins, we identify 22 ion pairs that are close enough to be classified as salt bridges. These are identified in the high-resolution crystal structures of the respective proteins or in the minimized average structures (if the crystal structures are unavailable) or, if both are unavailable, in the "most representative" conformer of each of the ensembles. We next calculate the electrostatic contribution of each such ion pair in each of the conformers in the ensembles. This results in a comprehensive study of 1,201 ion pairs, which allows us to look for consistent trends in their electrostatic contributions to protein stability in large sets of conformers. We find that the contributions of ion pairs vary considerably among the conformers of each protein. The vast majority of the ion pairs interconvert between being stabilizing and destabilizing to the structure at least once in the ensembles. These fluctuations reflect the variabilities in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other, and with respect to other charged groups in the remainder of the protein. The higher crystallographic B-factors for the respective side-chains are consistent with these fluctuations. The major conclusion from this study is that salt bridges observed in crystal structure may break, and new salt bridges may be formed. Hence, the overall stabilizing (or, destabilizing) contribution of an ion pair is conformer population dependent.

Protein folding and function: the N-terminal fragment in adenylate kinase

Three-dimensional protein folds range from simple to highly complex architectures. In complex folds, some building block fragments are more important for correct protein folding than others. Such fragments are typically buried in the protein core and mediate interactions between other fragments. Here we present an automated, surface area-based algorithm that is able to indicate which, among all local elements of the structure, is critical for the formation of the native fold, and apply it to structurally well-characterized proteins. In particular, we focus on adenylate kinase. The fragment containing the phosphate binding, P-loop (the "giant anion hole") flanked by a beta-strand and an alpha-helix near the N-terminus, is identified as a critical building block. This building block shows a high degree of sequence and structural conservation in all adenylate kinases. The results of our molecular dynamics simulations are consistent with this identification. In its absence, the protein flips to a stable, non-native state. In this misfolded conformation, the other local elements of the structure are in their native-like conformations; however, their association is non-native. Furthermore, this element is critically important for the function of the enzyme, coupling folding, and function.

A proposed structural model for amyloid fibril elongation: domain swapping forms an interdigitating beta-structure polymer

We propose a model illustrating how proteins, which differ in their overall sequences and structures, can form the propagating, twisted beta-sheet conformations, characteristic of amyloids. Some cases of amyloid formation can be explained through a "domain swapping" event, where the swapped segment is either a beta-hairpin or an unstable conformation which can partially unfold and assume a beta-hairpin structure. As in domain swapping, here the swapped beta-hairpin is at the edge of the structure, has few (if any) salt bridges and hydrogen bonds connecting it to the remainder of the structure and variable extents of buried non-polar surface areas. Additionally, in both cases the swapped piece constitutes a transient "building block" of the structure, with a high population time. Whereas in domain swapping the swapped fragment has been shown to be an alpha-helix, loop, strand or an entire domain, but so far not a beta-hairpin, despite the large number of cases in which it was already detected, here swapping may involve such a structural motif. We show how the swapping of beta-hairpins would form an interdigitated, twisted beta-sheet conformation, explaining the remarkable high stability of the protofibril in vitro. Such a swapping mechanism is attractive as it involves a universal mechanism in proteins, critical for their function, namely hinge-bending motions. Our proposal is consistent with structural superpositioning of mutational variants. While the overall r.m.s.d.s of the wild-type and mutants are small, the proposed hinge-bending region consistently shows larger deviations. These larger deviations illustrate that this region is more prone to respond to the mutational changes, regardless of their location in the sequence or in the structure. Nevertheless, above all, we stress that this proposition is hypothetical, since it is based on assumptions lacking definitive experimental support.

Molecular dynamics simulation of Escherichia coli dihydrofolate reductase and its protein fragments: relative stabilities in experiment and simulations

We have carried out molecular dynamics simulations of the native dihydrofolate reductase from Escherichia coli and several of its folded protein fragments at standard temperature. The simulations have shown fragments 1--36, 37--88, and 89--159 to be unstable, with a C(alpha)RMSD (C(alpha) root mean squared deviation) >5 A after 3.0 nsec of simulation. The unfolding of fragment 1--36 was immediate, whereas fragments 37--88 and 89--159 gradually unfolded because of the presence of the beta-sheet core structure. In the absence of residues 1--36, the two distinct domains comprising fragment 39--159 associated with each other, resulting in a stable conformation. This conformation retained most of its native structural elements. We have further simulated fragments derived from computational protein cutting. These were also found to be unstable, with the exception of fragment 104--159. In the absence of alpha(4), the loose loop region of residues 120--127 exhibited a beta-strand-like behavior, associating itself with the beta-sheet core of the protein fragment. The current study suggests that the folding of dihydrofolate reductase involves cooperative folding of distinct domains which otherwise would have been unstable as independent folded units in solution. Finally, the critical role of residues 1--36 in allowing the two distinct domains of fragment 104--159 to fold into the final native conformation is discussed.

A systematic study of the vibrational free energies of polypeptides in folded and random states

Molecular vibrations, especially low frequency motions, may be used as an indication of the rigidity or the flatness of the protein folding energy landscape. We have studied the vibrational properties of native folded as well as random coil structures of more than 60 polypeptides. The picture we obtain allows us to perceive how and why the energy landscape progressively rigidifies while still allowing potential flexibility. Compared with random coil structures, both alpha-helices and beta-hairpins are vibrationally more flexible. The vibrational properties of loop structures are similar to those of the corresponding random coil structures. Inclusion of an alpha-helix tends to rigidify peptides and so-called building blocks of the structure, whereas the addition of a beta-structure has less effect. When small building blocks coalesce to form larger domains, the protein rigidifies. However, some folded native conformations are still found to be vibrationally more flexible than random coil structures, for example, beta(2)-microglobulin and the SH3 domain. Vibrational free energy contributes significantly to the thermodynamics of protein folding and affects the distribution of the conformational substates. We found a weak correlation between the vibrational folding energy and the protein size, consistent with both previous experimental estimates and theoretical partition of the heat capacity change in protein folding.

Anatomy of protein structures: visualizing how a one-dimensional protein chain folds into a three-dimensional shape

Here, we depict the anatomy of protein structures in terms of the protein folding process. Via an iterative, top-down dissecting procedure, tertiary structures are spliced down to reveal their anatomy: first, to produce domains (defined by visual three-dimensional inspection criteria); then, hydrophobic folding units (HFU); and, at the end of a multilevel process, a set of building blocks. The resulting anatomy tree organization not only clearly depicts the organization of a one-dimensional polypeptide chain in three-dimensional space but also straightforwardly describes the most likely folding pathway(s). Comparison of the tree with the formation of the hydrophobic folding units through combinatorial assembly of the building blocks illustrates how the chain folds in a sequential or a complex folding pathway. Further, the tree points to the kinetics of the folding, whether the chain is a fast or a slow folder, and the probability of misfolding. Our ability to successfully dissect the protein into an anatomy tree illustrates that protein folding is a hierarchical process and further validates a building blocks protein folding model.

A rapid test for identification of autonomous folding units in proteins

The structure of a protein is dictated by a large number of weak interactions that cooperatively stabilize the native state. Usually, excised fragments smaller than a domain have little if any residual structure. When autonomous units of structure are found within domains, this challenges common assumptions about the cooperativity of protein structure. Such autonomous folding units (AFUs) are of wide interest and have applications in protein engineering and as simple model systems for studying the determinants of stability and specificity. A new method of identifying AFUs within proteins is presented here. The rapid autonomous fragment test (RAFT) identifies AFUs based on analysis of inter-residue contacts present in the three-dimensional structure of a protein. RAFT is fast enough to mine the entire PDB for AFUs and provide a library of potential small stable folds. We show that RAFT is able to predict whether a protein fragment will be structured if isolated from its parent domain.

Homology modelling of the DNA 5mC methyltransferase M.BssHII. Is permutation of functional subdomains common to all subfamilies of DNA methyltransferases?

This work presents a full tertiary model of the M.BssHII methyltransferase (MTase) complexed with substrate DNA and cofactor S-adenosyl-L-methionine, built by homology modelling based on previously solved complete structures of DNA MTases M.HaeIII and M. HhaI. M.BssHII and the template proteins show high sequence similarity, which indicates that they are evolutionary related. However, they are topologically different: M.BssHII is a circularly permuted variant of template MTases (Xu et al. Nucleic Acids Res 1997;25:3991). The model developed in this work will be a good starting point and valuable help in designing mutagenesis experiments to better understand the biological function of methyltransferases and the process of domain swapping.

Factors enhancing protein thermostability

Several sequence and structural factors have been proposed to contribute toward greater stability of thermophilic proteins. Here we present a statistical examination of structural and sequence parameters in representatives of 18 non-redundant families of thermophilic and mesophilic proteins. Our aim was to look for systematic differences among thermophilic and mesophilic proteins across the families. We observe that both thermophilic and mesophilic proteins have similar hydrophobicities, compactness, oligomeric states, polar and non-polar contribution to surface areas, main-chain and side-chain hydrogen bonds. Insertions/deletions and proline substitutions do not show consistent trends between the thermophilic and mesophilic members of the families. On the other hand, salt bridges and side chain-side chain hydrogen bonds increase in the majority of the thermophilic proteins. Additionally, comparisons of the sequences of the thermophile-mesophile homologous protein pairs indicate that Arg and Tyr are significantly more frequent, while Cys and Ser are less frequent in thermophilic proteins. Thermophiles both have a larger fraction of their residues in the alpha-helical conformation, and they avoid Pro in their alpha-helices to a greater extent than the mesophiles. These results indicate that thermostable proteins adapt dual strategies to withstand high temperatures. Our intention has been to explore factors contributing to the stability of proteins from thermophiles with respect to the melting temperatures (T(m)), the best descriptor of thermal stability. Unfortunately, T(m) values are available only for a few proteins in our high resolution dataset. Currently, this limits our ability to examine correlations in a meaningful way.

Electrostatic strengths of salt bridges in thermophilic and mesophilic glutamate dehydrogenase monomers

Here we seek to understand the higher frequency of occurrence of salt bridges in proteins from thermophiles as compared to their mesophile homologs. We focus on glutamate dehydrogenase, owing to the availability of high resolution thermophilic (from Pyrococcus furiosus) and mesophilic (from Clostridium symbiosum) protein structures, the large protein size and the large difference in melting temperatures. We investigate the location, statistics and electrostatic strengths of salt bridges and of their networks within corresponding monomers of the thermophilic and mesophilic enzymes. We find that many of the extra salt bridges which are present in the thermophilic glutamate dehydrogenase monomer but absent in the mesophilic enzyme, form around the active site of the protein. Furthermore, salt bridges in the thermostable glutamate dehydrogenase cluster within the hydrophobic folding units of the monomer, rather than between them. Computation of the electrostatic contribution of salt bridge energies by solving the Poisson equation in a continuum solvent medium, shows that the salt bridges in Pyrococcus furiosus glutamate dehydrogenase are highly stabilizing. In contrast, the salt bridges in the mesophilic Clostridium symbiosum glutamate dehydrogenase are only marginally stabilizing. This is largely the outcome of the difference in the protein environment around the salt bridges in the two proteins. The presence of a larger number of charges, and hence, of salt bridges contributes to an electrostatically more favorable protein energy term. Our results indicate that salt bridges and their networks may have an important role in resisting deformation/unfolding of the protein structure at high temperatures, particularly in critical regions such as around the active site.