A Simple Definition of Structural Regions in Proteins and Its Use in Analyzing Interface Evolution

Analysis of proteins commonly requires the partition of their structure into regions such as the surface, interior, or interface. Despite the frequent use of such categorization, no consensus definition seems to exist. This study thus aims at providing a definition that is general, is simple to implement, and yields new biological insights. This analysis relies on 397, 196, and 701 protein structures from Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens, respectively, and the conclusions are consistent across all three species. A threshold of 25% relative accessible surface area best segregates amino acids at the interior and at the surface. This value is further used to extend the core-rim model of protein-protein interfaces and to introduce a third region called support. Interface core, rim, and support regions contain similar numbers of residues on average, but core residues contribute over two-thirds of the contact surface. The amino acid composition of each region remains similar across different organisms and interface types. The interface core composition is intermediate between the surface and the interior, but the compositions of the support and the rim are virtually identical with those of the interior and the surface, respectively. The support and rim could thus “preexist” in proteins, and evolving a new interaction could require mutations to form an interface core only. Using the interface regions defined, it is shown through simulations that only two substitutions are necessary to shift the average composition of a  1000-Ų surface patch involving ∼ 28 residues to that of an equivalent interface. This analysis and conclusions will help understand the notion of promiscuity in protein-protein interaction networks.     

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.     

The Subunit Interfaces of Weakly Associated Homodimeric Proteins

We analyzed subunit interfaces in 315 homodimers with an X-ray structure in the Protein Data Bank, validated by checking the literature for data that indicate that the proteins are dimeric in solution and that, in the case of the “weak” dimers, the homodimer is in equilibrium with the monomer. The interfaces of the 42 weak dimers, which are smaller by a factor of 2.4 on average than in the remainder of the set, are comparable in size with antibody-antigen or protease-inhibitor interfaces. Nevertheless, they are more hydrophobic than in the average transient protein-protein complex and similar in amino acid composition to the other homodimer interfaces. The mean numbers of interface hydrogen bonds and hydration water molecules per unit area are also similar in homodimers and transient complexes. Parameters related to the atomic packing suggest that many of the weak dimer interfaces are loosely packed, and we suggest that this contributes to their low stability. To evaluate the evolutionary selection pressure on interface residues, we calculated the Shannon entropy of homologous amino acid sequences at 60% sequence identity. In 93% of the homodimers, the interface residues are better conserved than the residues on the protein surface. The weak dimers display the same high degree of interface conservation as other homodimers, but their homologs may be heterodimers as well as homodimers. Their interfaces may be good models in terms of their size, composition, and evolutionary conservation for the labile subunit contacts that allow protein assemblies to share and exchange components, allosteric proteins to undergo quaternary structure transitions, and molecular machines to operate in the cell.     

A dissection of the protein-protein interfaces in icosahedral virus capsids

We selected 49 icosahedral virus capsids whose crystal structures are reported in the Protein Data Bank. They belong to the T=1, T=3, pseudo T=3 and other lattice types. We identified in them 779 unique interfaces between pairs of subunits, all repeated by icosahedral symmetry. We analyzed the geometric and physical chemical properties of these interfaces and compared with interfaces in protein-protein complexes and homodimeric proteins, and with crystal packing contacts. The capsids contain one to 16 subunits implicated in three to 66 unique interfaces. Each subunit loses 40-60% of its accessible surface in contacts with an average of 8.5 neighbors. Many of the interfaces are very large with a buried surface area (BSA) that can exceed 10,000 A(2), yet 39% are small with a BSA<800 A(2) comparable to crystal packing contacts. Pairwise capsid interfaces overlap, so that one-third of the residues are part of more than one interface. Those with a BSA>800 A(2) resemble homodimer interfaces in their chemical composition. Relative to the protein surface, they are non-polar, enriched in aliphatic residues and depleted of charged residues, but not of neutral polar residues. They contain one H-bond per about 200 A(2) BSA. Small capsid interfaces (BSA<800 A(2)) are only slightly more polar. They have a similar amino acid composition, but they bury fewer atoms and contain fewer H-bonds for their size. Geometric parameters that estimate the quality of the atomic packing suggest that the small capsid interfaces are loosely packed like crystal packing contacts, whereas the larger interfaces are close-packed as in protein-protein complexes and homodimers. We discuss implications of these findings on the mechanism of capsid assembly, assuming that the larger interfaces form first to yield stable oligomeric species (capsomeres), and that medium-size interfaces allow the stepwise addition of capsomeres to build larger intermediates.     

Comparative studies on ion-pair energetic,distribution among three domains of life: Archaea,eubacteria, and eukarya

Salt-bridges play a unique role in the structural and functional stability of proteins, especially under harsh environments. How these salt-bridges contribute to the overall thermodynamic stability of protein structure and function across different domains of life is elusive still date. To address the issue, statistical analyses on the energies of salt-bridges, involved in proteins' structure and function, are performed across three domains of life, that is, archaea, eubacteria, and eukarya. Results show that although the majority of salt-bridges are stable and conserved, yet the stability of archaeal proteins (∆∆G_net = −5.06 ± 3.8) is much more than that of eubacteria (∆∆G_net = −3.7 ± 2.9) and eukarya (∆∆G_net = −3.54 ± 3.1). Unlike earlier study with archaea, in eukarya and eubacteria, not all buried salt-bridge in our dataset are stable. Buried salt-bridges play surprising role in protein stability, whose variations are clearly observed among these domains. Greater desolvation penalty of buried salt-bridges is compensated by stable network of salt-bridges apart from equal contribution of bridge and background energy terms. On the basis proteins' secondary structure, topology, and evolution, our observation shows that salt-bridges when present closer to each other in sequence tend to form a greater number. Overall, our comparative study provides insight into the role of specific electrostatic interactions in proteins from different domains of life, which we hope, would be useful for protein engineering and bioinformatics study.     

Contribution of electrostatic interactions, compactness and quaternary structure to protein thermostability: lessons from structural genomics of Thermotoga maritima

Studies of the structural basis of protein thermostability have produced a confusing picture. Small sets of proteins have been analyzed from a variety of thermophilic species, suggesting different structural features as responsible for protein thermostability. Taking advantage of the recent advances in structural genomics, we have compiled a relatively large protein structure dataset, which was constructed very carefully and selectively; that is, the dataset contains only experimentally determined structures of proteins from one specific organism, the hyperthermophilic bacterium Thermotoga maritima, and those of close homologs from mesophilic bacteria. In contrast to the conclusions of previous studies, our analyses show that oligomerization order, hydrogen bonds, and secondary structure play minor roles in adaptation to hyperthermophily in bacteria. On the other hand, the data exhibit very significant increases in the density of salt-bridges and in compactness for proteins from T.maritima. The latter effect can be measured by contact order or solvent accessibility, and network analysis shows a specific increase in highly connected residues in this thermophile. These features account for changes in 96% of the protein pairs studied. Our results provide a clear picture of protein thermostability in one species, and a framework for future studies of thermal adaptation.     

Insights into protein-protein interfaces using a Bayesian network prediction method

Identifying the interface between two interacting proteins provides important clues to the function of a protein, and is becoming increasing relevant to drug discovery. Here, surface patch analysis was combined with a Bayesian network to predict protein-protein binding sites with a success rate of 82% on a benchmark dataset of 180 proteins, improving by 6% on previous work and well above the 36% that would be achieved by a random method. A comparable success rate was achieved even when evolutionary information was missing, a further improvement on our previous method which was unable to handle incomplete data automatically. In a case study of the Mog1p family, we showed that our Bayesian network method can aid the prediction of previously uncharacterised binding sites and provide important clues to protein function. On Mog1p itself a putative binding site involved in the SLN1-SKN7 signal transduction pathway was detected, as was a Ran binding site, previously characterized solely by conservation studies, even though our automated method operated without using homologous proteins. On the remaining members of the family (two structural genomics targets, and a protein involved in the photosystem II complex in higher plants) we identified novel binding sites with little correspondence to those on Mog1p. These results suggest that members of the Mog1p family bind to different proteins and probably have different functions despite sharing the same overall fold. We also demonstrated the applicability of our method to drug discovery efforts by successfully locating a number of binding sites involved in the protein-protein interaction network of papilloma virus infection. In a separate study, we attempted to distinguish between the two types of binding site, obligate and non-obligate, within our dataset using a second Bayesian network. This proved difficult although some separation was achieved on the basis of patch size, electrostatic potential and conservation. Such was the similarity between the two interacting patch types, we were able to use obligate binding site properties to predict the location of non-obligate binding sites and vice versa.     

Dissociation of Abeta(16-22) amyloid fibrils probed by molecular dynamics

The mechanisms of deposition and dissociation are implicated in the assembly of amyloid fibrils. To investigate the kinetics of unbinding of Abeta(16-22) monomers from preformed fibrils, we use molecular dynamics (MD) simulations and the structures for Abeta(16-22) amyloid fibrils. Consistent with experimental studies, the dissociation of Abeta(16-22) peptides involves two main stages, locked and docked, after which peptides unbind. The lifetime of the locked state, in which a peptide retains fibril-like structure and interactions, extends up to 0.5 micros under normal physiological conditions. Upon cooperative rupture of all fibril-like hydrogen bonds (HBs) with the fibril, a peptide enters a docked state. This state is populated by disordered random coil conformations and its lifetime ranges from approximately 10 to 200 ns. The docked state is stabilized by hydrophobic side chain interactions, while the contribution from HBs is small. Our simulations also suggest that the peptides located on fibril edges may form stable beta-strand conformations distinct from the fibril "bulk". We propose that such edge peptides can act as fibril caps, which impede fibril elongation. Our results indicate that the interactions between unbinding peptides constitute the molecular basis for cooperativity of peptide dissociation. The kinetics of fibril growth is reconstructed from unbinding assuming the reversibility of deposition/dissociation pathways. The relation of in silica dissociation kinetics to experimental observations is discussed.     

The Membrane Insertase Oxa1 Is Required for Efficient Import of Carrier Proteins into Mitochondria

Oxa1 serves as a protein insertase of the mitochondrial inner membrane that is evolutionary related to the bacterial YidC insertase. Its activity is critical for membrane integration of mitochondrial translation products and conservatively sorted inner membrane proteins after their passage through the matrix. All Oxa1 substrates identified thus far have bacterial homologs and are of endosymbiotic origin. Here, we show that Oxa1 is critical for the biogenesis of members of the mitochondrial carrier proteins. Deletion mutants lacking Oxa1 show reduced steady‐state levels and activities of the mitochondrial ATP/ADP carrier protein Aac2. To reduce the risk of indirect effects, we generated a novel temperature-sensitive oxa1 mutant that allows rapid depletion of a mutated Oxa1 variant in situ by mitochondrial proteolysis. Oxa1-depleted mitochondria isolated from this mutant still contain normal levels of the membrane potential and of respiratory chain complexes. Nevertheless, in vitro import experiments showed severely reduced import rates of Aac2 and other members of the carrier family, whereas the import of matrix proteins was unaffected. From this, we conclude that Oxa1 is directly or indirectly required for efficient biogenesis of carrier proteins. This was unexpected, since carrier proteins are inserted into the inner membrane from the intermembrane space side and lack bacterial homologs. Our observations suggest that the function of Oxa1 is relevant not only for the biogenesis of conserved mitochondrial components such as respiratory chain complexes or ABC transporters but also for mitochondria-specific membrane proteins of eukaryotic origin.     

Membrane-SPINE: an improved method to identify protein-protein interaction partners of membrane proteins in vivo

Membrane proteins are crucial for many essential cellular processes. As membrane proteins function in complexes, methods to detect and to characterize membrane protein-protein interactions are undoubtedly required. Therefore, we developed the "Membrane-Strep-tagged protein interaction experiment" (Membrane-SPINE) that combines the specific purification of a Strep-tagged membrane protein with the reversible fixation of protein complexes by formaldehyde cross-linking. In combination with MS analysis, we suggest Membrane-SPINE as a powerful tool to identify unknown interaction partners of membrane proteins in vivo.     

Comparison of X‐ray and NMR structures: Is there a systematic difference in residue contacts between X‐ray‐ and NMR‐resolved protein structures?

We have compared structures of 78 proteins determined by both NMR and X-ray methods. It is shown that X-ray and NMR structures of the same protein have more differences than various X-ray structures obtained for the protein, and even more than various NMR structures of the protein. X-ray and NMR structures of 18 of these 78 proteins have obvious large-scale structural differences that seem to reflect a difference of crystal and solution structures. The other 60 pairs of structures have only small-scale differences comparable with differences between various X-ray or various NMR structures of a protein; we have analyzed these structures more attentively. One of the main differences between NMR and X-ray structures concerns the number of contacts per residue: (1) NMR structures presented in PDB have more contacts than X-ray structures at distances below 3.0 A and 4.5-6.5 A, and fewer contacts at distances of 3.0-4.5 A and 6.5-8.0 A; (2) this difference in the number of contacts is greater for internal residues than for external ones, and it is larger for beta-containing proteins than for all-alpha proteins. Another significant difference is that the main-chain hydrogen bonds identified in X-ray and NMR structures often differ. Their correlation is 69% only. However, analogous difference is found for refined and rerefined NMR structures, allowing us to suggest that the observed difference in interresidue contacts of X-ray and NMR structures of the same proteins is due mainly to a difference in mathematical treatment of experimental results.     

Glycosphingolipid behaviour in complex membranes

Glycosphingolipids, due to their tendency to form laterally separated liquid-ordered phases, possess a high potential for the creation of order in biological membranes. The formation of glycosphingolipid-rich domains within the membrane has profound consequences on the membrane organization at different levels, and on the conformational and biological properties of membrane-associated proteins and multimolecular protein complexes. In this review, we will discuss 1) how glycosphingolipids influence the lateral organization of biological membranes; 2) how glycosphingolipids influence the function of membrane-associated proteins.     

Distribution, lateral mobility and function of membrane proteins incorporated into giant unilamellar vesicles

GUVs have been widely used for studies on lipid mobility, membrane dynamics and lipid domain (raft) formation, using single molecule techniques like fluorescence correlation spectroscopy. Reports on membrane protein dynamics in these types of model membranes are by far less advanced due to the difficulty of incorporating proteins into GUVs in a functional state. We have used sucrose to prevent four distinct membrane protein(s) (complexes) from inactivating during the dehydration step of the GUV-formation process. The amount of sucrose was optimized such that the proteins retained 100% biological activity, and many proteo-GUVs were obtained. Although GUVs could be formed by hydration of lipid mixtures composed of neutral and anionic lipids, an alternate current electric field was required for GUV formation from neutral lipids. Distribution, lateral mobility, and function of an ATP-binding cassette transport system, an ion-linked transporter, and a mechanosensitive channel in GUVs were determined by confocal imaging, fluorescence correlation spectroscopy, patch-clamp measurements, and biochemical techniques. In addition, we show that sucrose slows down the lateral mobility of fluorescent lipid analogs, possibly due to hydrogen-bonding with the lipid headgroups, leading to larger complexes with reduced mobility.     

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.     

Protein complexes of the Escherichia coli cell envelope

Protein complexes are an intrinsic aspect of life in the membrane. Knowing which proteins are assembled in these complexes is therefore essential to understanding protein function(s). Unfortunately, recent high throughput protein interaction studies have failed to deliver any significant information on proteins embedded in the membrane, and many membrane protein complexes remain ill defined. In this study, we have optimized the blue native-PAGE technique for the study of membrane protein complexes in the inner and outer membranes of Escherichia coli. In combination with second dimension SDS-PAGE and mass spectrometry, we have been able to identify 43 distinct protein complexes. In addition to a number of well characterized complexes, we have identified known and orphan proteins in novel oligomeric states. For two orphan proteins, YhcB and YjdB, our findings enable a tentative functional assignment. We propose that YhcB is a hitherto unidentified additional subunit of the cytochrome bd oxidase and that YjdB, which co-localizes with the ZipA protein, is involved in cell division. Our reference two-dimensional blue native-SDS-polyacrylamide gels will facilitate future studies of the assembly and composition of E. coli membrane protein complexes during different growth conditions and in different mutant backgrounds.     

Probing protein conformation with a minimal photochemical reagent

3H-diazirine (3H-DZN), a photoreactive gas similar in size to water, was used to probe the topography of the surface and inner space of proteins. On photolysis 3H-DZN generates 3H-methylene carbene, which reacts unselectively with its molecular cage, inserting even into C-H bonds. Labeling of bovine alpha-lactalbumin (alpha-LA, MW: 14,200) with 1 mM (3)H-DZN yielded 0.0041 mol CH2/mol of protein, in agreement with the expectation for an unspecific surface-labeling phenomenon. The cooperative urea-induced unfolding of alpha-LA, as monitored by the extent of 3H-methylene labeling, agrees with that measured by circular dichroism spectroscopy in the far and near ultraviolet regions. At 8 M urea, the unfolded state U was labeled 25-30% more than the native state N primarily because of the increase in the accessible surface area (ASA) of the protein occurring upon unfolding. However, this result lies below the approximately 100% increment expected from theoretical estimates of ASA of state U. Among other factors, most likely the existence of a residual structure in U, that involves helices H2 and H4 of the alpha subdomain, might account for this fact, as shown by a comparative analysis of peptide labeling patterns of N and U samples. In this paper, we demonstrate the usefulness of the 3H-methylene labeling method to monitor conformational transitions and map solvent accessibility along the polypeptide sequence, thus opening the possibility of outlining structural features of nonnative states (i.e., denatured states, molten globule). We anticipate that this technique also would help to identify ligand binding and oligomerization sites in proteins.     

Structure of a hyperthermophilic archaeal homing endonuclease, I-Tsp061I: contribution of cross-domain polar networks to thermostability

A novel LAGLIDADG-type homing endonuclease (HEase), I-Tsp061I, from the hyperthermophilic archaeon Thermoproteus sp. IC-061 16 S rRNA gene (rDNA) intron was characterized with respect to its structure, catalytic properties and thermostability. It was found that I-Tsp061I is a HEase isoschizomer of the previously described I-PogI and exhibits the highest thermostability among the known LAGLIDADG-type HEases. Determination of the crystal structure of I-Tsp061I at 2.1 A resolution using the multiple isomorphous replacement and anomalous scattering method revealed that the overall fold is similar to that of other known LAGLIDADG-type HEases, despite little sequence similarity between I-Tsp061I and those HEases. However, I-Tsp061I contains important cross-domain polar networks, unlike its mesophilic counterparts. Notably, the polar network Tyr6-Asp104-His180-107O-HOH12-104O-Asn177 exists across the two packed alpha-helices containing both the LAGLIDADG catalytic motif and the GxxxG hydrophobic helix bundle motif. Another important structural feature is the salt-bridge network Asp29-Arg31-Glu182 across N and C-terminal domain interface, which appears to contribute to the stability of the domain/domain packing. On the basis of these structural analyses and extensive mutational studies, we conclude that such cross-domain polar networks play key roles in stabilizing the catalytic center and domain packing, and underlie the hyperthermostability of I-Tsp061I.     

Image correlation spectroscopy of randomly distributed disks

Image correlation spectroscopy (ICS) has been widely used to quantify spatiotemporal distributions of fluorescently labelled cell membrane proteins and receptors. When the membrane proteins are randomly distributed, ICS may be used to estimate protein densities, provided the proteins behave as point-like objects. At high protein area fraction, however, even randomly placed proteins cannot obey Poisson statistics, because of excluded area. The difficulty can arise if the protein effective area is quite large, or if proteins form large complexes or aggregate into clusters. In these cases, there is a need to determine the correct form of the intensity correlation function for hard disks in two dimensions, including the excluded area effects. We present an approximate but highly accurate algorithm for the computation of this correlation function. The correlation function was verified using test images of randomly distributed hard disks of uniform intensity convolved with the microscope point spread function. This algorithm can be readily modified to compute exact intensity correlation functions for any probe geometry, interaction potential, and fluorophore distribution; we show how to apply it to describe a random distribution of large proteins labeled with a single fluorophore.     

Multiple roles of the vesicular-SNARE TI-VAMP in post-Golgi and endosomal trafficking

SNARE (Soluble N-ethylmaleimide-sensitive factor attachment protein receptor) proteins are the core machinery of membrane fusion. Vesicular SNAREs (v-SNAREs) interact with their target SNAREs (t-SNAREs) to form SNARE complexes which mediate membrane fusion. Here we review the basic properties and functions of the v-SNARE TI-VAMP/VAMP7 (Tetanus neurotoxin insensitive-vesicle associated membrane protein). TI-VAMP interacts with its t-SNARE partners, particularly plasmalemmal syntaxins, to mediate membrane fusion and with several regulatory proteins especially via its amino-terminal regulatory Longin domain. Partners include AP-3, Hrb/(Human immunodeficiency virus Rev binding) protein, and Varp (Vps9 domain and ankyrin repeats containing protein) and regulate TI-VAMP’s function and targeting. TI-VAMP is involved both in secretory and endocytic pathways which mediate neurite outgrowth and synaptic transmission, plasma membrane remodeling and lysosomal secretion.