Oncogenic mutations and packing defects in protein structure

Oncogenic mutations in expressed proteins are of primary interest to understand tumor formation but their structural consequences bearing on protein function are not clearly understood. In this contribution I report on two illustrative examples, p21ras and p57, revealing that such mutations have an effect on specific structural deficiencies in the packing of the protein structure, i. e., on backbone hydrogen bonds insufficiently shielded from water attack. These structural deficiencies in the wild type are typically "corrected intermolecularly" by protein complexation or protein-ligand association. However, in the oncogenic mutants, these binding signals are partially or completely suppressed: the mutated residues properly wrap or desolvate the hydrogen bonds intramolecularly. Thus, the interactivity of the proteins becomes impaired: their binding affinity decreases sharply, as there is no thermodynamic benefit from removing water surrounding properly desolvated hydrogen bonds. The results, specialized for p21ras and p53, reveal how oncogenic mutations determine a hindrance to GAP-induced hydrolysis (p21) and decrease binding affinity for DNA (p53). Furthermore, the oncogenic potential of mutations in residues not directly engaged in the interface electrostatics is assessed. The results suggest that a high sensitivity of structural defects to genetic accident might be a necessary condition to establish the existence of a proto-oncogene, an angle that merits a systematic study.     

Increasing the thermostability of staphylococcal nuclease: implications for the origin of protein thermostability

Seven hyper-stable multiple mutants have been constructed in staphylococcal nuclease by various combinations of eight different stabilizing single mutants. The stabilities of these multiple mutants determined by guanidine hydrochloride denaturation were 3.4 to 5.6 kcal/mol higher than that of the wild-type. Their thermal denaturation midpoint temperatures were 12.6 to 22.9 deg. C higher than that of the wild-type. These are among the greatest increases in protein stability and thermal denaturation midpoint temperature relative to the wild-type yet attained. There has been great interest in understanding how proteins found in thermophilic organisms are stabilized. One frequently cited theory is that the packing of hydrophobic side-chains is improved in the cores of proteins isolated from thermophiles when compared to proteins from mesophiles. The crystal structures of four single and five multiple stabilizing mutants of staphylococcal nuclease were solved to high resolution. No large overall structural change was found, with most changes localized around the sites of mutation. Rearrangements were observed in the packing of side-chains in the major hydrophobic core, although none of the mutations was in the core. It is surprising that detailed structural analysis showed that packing had improved, with the volume of the mutant protein's hydrophobic cores decreasing as protein stability increased. Further, the number of van der Waals interactions in the entire protein showed an experimentally significant increase correlated with increasing stability. These results indicate that optimization of packing follows as a natural consequence of increased protein thermostability and that good packing is not necessarily the proximate cause of high stability. Another popular theory is that thermostable proteins have more electrostatic and hydrogen bonding interactions and these are responsible for the high stabilities. The mutants here show that increased numbers of electrostatic and hydrogen bonding interactions are not obligatory for large increases in protein stability.     

Roles for holes: are cavities in proteins mere packing defects?

Atomic packing in proteins is not optimized, most structures containing internal cavities, which have been identified by molecular modelling and characterized experimentally. Cavities seem to play a role in assisting conformational changes between domains or subunit interfaces. Comparison between homologous proteins from thermophiles and mesophiles indicates that optimizing packing enhances stabilization at the expense of flexibility. For proteins which interact with small ligands or substrates, cavities seem to play a role in controlling binding and catalysis, rather than being mere "packing defects". We believe that a more complete analysis on the localization, conservation and role of cavities in protein structures (by modelling and site-directed mutagenesis), will reveal that rather than being randomly distributed, they are located in key positions to allow structural dynamics and thereby functional control.     

Analysis of pH-dependent elements in proteins: geometry and properties of pairs of hydrogen-bonded carboxylic acid side-chains

A rather frequent but so far little discussed observation is that pairs of carboxylic acid side-chains in proteins can share a proton in a hydrogen bond. In the present article, quantum chemical calculations of simple model systems for carboxyl-carboxylate interactions are compared with structural observations from proteins. A detailed structural analysis of the proteins deposited in the PDB revealed that, in a subset of proteins sharing less than 90% sequence identity, 19% (314) contain at least one pair of carboxylic acids with their side-chain oxygen atoms within hydrogen-bonding distance. As the distance between those interacting oxygen atoms is frequently very short ( approximately 2.55 A), many of these carboxylic acids are suggested to share a proton in a strong hydrogen bond. When situated in an appropriate structural environment (low dielectric constant), some might even form a low barrier hydrogen bond. The quantum chemical studies show that the most frequent geometric features of carboxyl-carboxylate pairs found in proteins, and no or symmetric ligation, are also the most stable arrangements at low dielectric constants, and they also suggest at medium and low pH a higher stability than for isosteric amide-carboxylate pairs. The presence of these pairs in 119 different enzymes found in the BRENDA database is set in relation to their properties and functions. This analysis shows that pH optima of enzymes with carboxyl-carboxylate pairs are shifted to lower than average values, whereas temperature optima seem to be increased. The described structural principles can be used as guidelines for rational protein design (e.g., in order to improve pH or temperature stability).     

Oncogenic Mutations and Packing Defects in Protein Structure

Abstract

Oncogenic mutations in expressed proteins are of primary interest to understand tumor formation but their structural consequences bearing on protein function are not clearly understood. In this contribution I report on two illustrative examples, p21ras and p57, revealing that such mutations have an effect on specific structural deficiencies in the packing of the protein structure, i. e., on backbone hydrogen bonds insufficiently shielded from water attack. These structural deficiencies in the wild type are typically “corrected intermolecu- larly” by protein complexation or protein-ligand association. However, in the oncogenic mutants, these binding signals are partially or completely suppressed: the mutated residues properly wrap or desolvate the hydrogen bonds intramolecularly. Thus, the interactivity of the proteins becomes impaired: their binding affinity decreases sharply, as there is no thermodynamic benefit from removing water surrounding properly desolvated hydrogen bonds. The results, specialized for p21ras and p53, reveal how oncogenic mutations determine a hindrance to GAP-induced hydrolysis (p21) and decrease binding affinity for DNA (p53). Furthermore, the oncogenic potential of mutations in residues not directly engaged in the interface electrostatics is assessed. The results suggest that a high sensitivity of structural defects to genetic accident might be a necessary condition to establish the existence of a proto-oncogene, an angle that merits a systematic study.     

Lipid packing defects and membrane charge control RAB GTPase recruitment

Specific intracellular localization of RAB GTPases has been reported to be dependent on protein factors, but the contribution of the membrane physicochemical properties to this process has been poorly described. Here, we show that three RAB proteins (RAB1/RAB5/RAB6) preferentially bind in vitro to disordered and curved membranes, and that this feature is uniquely dependent on their prenyl group. Our results imply that the addition of a prenyl group confers to RAB proteins, and most probably also to other prenylated proteins, the ability to sense lipid packing defects induced by unsaturated conical‐shaped lipids and curvature. Consistently, RAB recruitment increases with the amount of lipid packing defects, further indicating that these defects drive RAB membrane targeting. Membrane binding of RAB35 is also modulated by lipid packing defects but primarily dependent on negatively charged lipids. Our results suggest that a balance between hydrophobic insertion of the prenyl group into lipid packing defects and electrostatic interactions of the RAB C‐terminal region with charged membranes tunes the specific intracellular localization of RAB proteins.      

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.     

Hydration of protein-protein interfaces

We present an analysis of the water molecules immobilized at the protein-protein interfaces of 115 homodimeric proteins and 46 protein-protein complexes, and compare them with 173 large crystal packing interfaces representing nonspecific interactions. With an average of 15 waters per 1000 A2 of interface area, the crystal packing interfaces are more hydrated than the specific interfaces of homodimers and complexes, which have 10-11 waters per 1000 A2, reflecting the more hydrophilic composition of crystal packing interfaces. Very different patterns of hydration are observed: Water molecules may form a ring around interfaces that remain "dry," or they may permeate "wet" interfaces. A majority of the specific interfaces are dry and most of the crystal packing interfaces are wet, but counterexamples exist in both categories. Water molecules at interfaces form hydrogen bonds with protein groups, with a preference for the main-chain carbonyl and the charged side-chains of Glu, Asp, and Arg. These interactions are essentially the same in specific and nonspecific interfaces, and very similar to those observed elsewhere on the protein surface. Water-mediated polar interactions are as abundant at the interfaces as direct protein-protein hydrogen bonds, and they may contribute to the stability of the assembly.     

Backbone cluster identification in proteins by a graph theoretical method

A graph theoretical algorithm has been developed to identify backbone clusters of residues in proteins. The identified clusters show protein sites with the highest degree of interactions. An adjacency matrix is constructed from the non-bonded connectivity information in proteins. The diagonalization of such a matrix yields eigenvalues and eigenvectors, which contain the information on clusters. In graph theory, distinct clusters can be obtained from the second lowest eigenvector components of the matrix. However, in an interconnected graph, all the points appear as one single cluster. We have developed a method of identifying highly interacting centers (clusters) in proteins by truncating the vector components of high eigenvalues. This paper presents in detail the method adopted for identifying backbone clusters and the application of the algorithm to families of proteins like RNase-A and globin. The objective of this study was to show the efficiency of the algorithm as well as to detect conserved or similar backbone packing regions in a particular protein family. Three clusters in topologically similar regions in the case of the RNase-A family and three clusters around the porphyrin ring in the globin family were observed. The predicted clusters are consistent with the features of the family of proteins such as the topology and packing density. The method can be applied to problems such as identification of domains and recognition of structural similarities in proteins.     

Atomic size packing defects in proteins

The three-dimensional refined high resolution structures of 20 proteins were examined for the presence of packing defects of atomic size or larger. Of the proteins examined, 12 had no such packing defects, 6 proteins had just 1 packing defect, and 2 proteins had 2 or 3 packing defects. These results confirm earlier studies on smaller samples of proteins which demonstrated that proteins are well packed. The atoms that surround the packing defects are almost always hydrophobic (carbon or sulfur). This study also tabulated the number of internal waters in each protein, which varied from 0 to 28.     

Formation of an infinite beta-sheet arrangement dominates the crystallization behavior of lambda-type antibody light chains

The packing interactions in crystals of human lambda-type antibody light chain dimers have been reviewed. These homologous proteins are composed of individually specific variable domains, but all have very similar constant domain sequences. The proteins do not emulate each other in their overall crystallization behavior: each attains an individually characteristic space group or unit cell dimensions. However, each of these protein crystals has one unit cell dimension in common, 72.4(+/- 0.2) A. Examination of the protein packing in these crystals reveals that the common cell dimension is a consequence of a packing arrangement of their constant domains, which is conserved in all three crystals. In this striking arrangement, beta-sheets of adjacent constant domains are placed in juxta-position to form an "infinite chain". Although this constant domain packing pattern is rigorously conserved, the variable domain packing arrangements in each of these crystals are different. The conservation of the "infinite" beta-sheet pattern suggests that the constant domain interactions dominate the thermodynamic energy of lattice formation, probably through a combination of specific hydrogen bond formations and by a decrease in the solvent-accessible surface. A single amino acid substitution prohibits this characteristic interneighbor hydrogen bond pattern in the homologous kappa-type light chains. This may explain the observation that very few kappa-type light chains have been crystallized.     

Cavities and packing at protein interfaces.

An analysis of internal packing defects or "cavities" (both empty and water-containing) within protein structures has been undertaken and includes 3 cavity classes: within domains, between domains, and between protein subunits. We confirm several basic features common to all cavity types but also find a number of new characteristics, including those that distinguish the classes. The total cavity volume remains only a small fraction of the total protein volume and yet increases with protein size. Water-filled "cavities" possess a more polar surface and are typically larger. Their constituent waters are necessary to satisfy the local hydrogen bonding potential. Cavity-surrounding atoms are observed to be, on average, less flexible than their environments. Intersubunit and interdomain cavities are on average larger than the intradomain cavities, occupy a larger fraction of their resident surfaces, and are more frequently water-filled. We observe increased cavity volume at domain-domain interfaces involved with shear type domain motions. The significance of interfacial cavities upon subunit and domain shape complementarity and the protein docking problem, as well as in their structural and functional role in oligomeric proteins, will be discussed. The results concerning cavity size, polarity, solvation, general abundance, and residue type constituency should provide useful guidelines for protein modeling and design.     

Buried waters and internal cavities in monomeric proteins

We have analyzed the buried water molecules and internal cavities in a set of 75 high-resolution, nonhomologous, monomeric protein structures. The number of hydrogen bonds formed between each water molecule and the protein varies from 0 to 4, with 3 being most common. Nearly half of the water molecules are found in pairs or larger clusters. Approximately 90% are shown to be associated with large cavities within the protein, as determined by a novel program, PRO_ACT. The total volume of a protein's large cavities is proportional to its molecular weight and is not dependent on structural class. The largest cavities in proteins are generally elongated rather than globular. There are many more empty cavities than hydrated cavities. The likelihood of a cavity being occupied by a water molecule increases with cavity size and the number of available hydrogen bond partners, with each additional partner typically stabilizing the occupied state by 0.6 kcal/mol.     

beta-Trefoil fold. Patterns of structure and sequence in the Kunitz inhibitors interleukins-1 beta and 1 alpha and fibroblast growth factors.

Previous crystallographic analyses of the Kunitz inhibitors from soybean. Erythrina caffra and wheat, the interleukins-1 beta and 1 alpha and the acidic and basic fibroblast growth factors have shown that they contain a most unusual fold. It is formed by six two-stranded hairpins. Three of these form a barrel structure and the other three are in a triangular array that caps the barrel. The arrangement of the secondary structures gives the molecules a pseudo 3-fold axis. Although the different proteins have very similar structures, many of their sequences have no significant similarities overall. The structural determinants of this fold are described and discussed in this paper. The barrels in the different proteins have the same geometrical features: six strands tilted at 56 degrees to the barrel axis; a barrel diameter of 16 A, and the beta-sheet hydrogen bonded so that it is staggered with a shear number of 12. These features fit McLachlan's equations for ideal barrels formed by beta-sheets. The wide diameter of the barrels is filled by layers of residues that, while not identical in the different proteins, are, in almost all cases, large. The structure of the triangular array of hairpins is determined by the coiling of the strands and the packing of hairpin residues against each other and against residues from the interior of the barrel. The major sequence requirements of this fold are large or medium hydrophobic residues at 18 buried sites. In the different structures the total volume of these residues is 3000 (+/- 120) A. The polyhedron model of protein architecture is used to demonstrate that the main, and in particular the symmetrical, features of this fold arise from the ideal and equal packing of six hairpins, modified only slightly to form hydrogen bonds between the hairpins.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Calculation of standard atomic volumes for RNA and comparison with proteins: RNA is packed more tightly

Traditionally, for biomolecular packing calculations research has focused on proteins. Besides proteins, RNA is the other large biomolecule that has tertiary structure interactions and complex packing. No one has yet quantitatively investigated RNA packing or compared its packing to that of proteins because, until recently, there were no large RNA structures. Here we address this question in detail, using Voronoi volume calculations on a set of high-resolution RNA crystal structures. We do a careful parameterization, taking into account many factors such as atomic radii, crystal packing, structural complexity, solvent, and associated protein to obtain a self-consistent, universal set of volumes that can be applied to both RNA and protein. We report this set of volumes, which we call the NucProt parameter set. Our measured values are consistent across the many different RNA structures and packing environments. When common atom types are compared between proteins and RNA, nine of 12 types show that RNA has a smaller volume and packs more tightly than protein, suggesting that close-packing may be as important for the folding of RNAs as it is for proteins. Moreover, calculated partial specific volumes show that RNA bases pack more densely than corresponding aromatic residues from proteins. Finally, we find that RNA bases have similar packing volumes to DNA bases, despite the absence of tertiary contacts in DNA. Programs, parameter sets and raw data are available online at.     

A study of protein structure changes during hydration by means of diffuse X-ray scattering. II. Fourier transform analysis of X-ray scattering data

Radial distribution functions were deduced by Fourier transform analysis of angular dependences of diffuse x-ray scattering intensities for the following proteins with different hydration degree: water-soluble a-protein myoglobin, water-soluble alpha+beta protein lysozyme, and transmembrane proteins of photosynthetic reaction centers from purple bacteria Rhodobacter sphaeroides and Blastochlorii viridis. The results of Fourier analysis of x-ray scattering intensities give the quantitative characteristics of the mechanisms underlying the influence of water on the formation of biomacromolecules. Water, on the one hand, weakens the intraglobular hydrogen bond net, loosens the protein structure, and increases the internal conformational dynamics. Concurrently water arranges the stability and ordering of the macromolecule. A sharp correlation is observed between the shift of the "first" peak of radial distribution functions, the weakening of the intraglobular hydrogen bond net, the increase in intraglobular mobility, and the appearance of functional activity in macromolecules. The behavior of the "first" peak is similar to that observed in transmembrane protein of reaction center and water-soluble proteins. The "first" peak for transmembrane protein of reaction center reaches its maximum value much faster (at smaller hydration degrees) than for water-soluble proteins. The fast transfer of reaction center protein to its native state during hydration is due to the fact that the dehydrated conformation of reaction center protein is very close to the native one. From a comparison of the radial distribution functions for water, water-soluble proteins and transmembrane proteins, one may conclude that water has the lowest packing density and the lowest order; water-soluble proteins have a larger packing density and are more ordered than water, and transmembrane proteins have the highest degree of packing density and ordering.     

Protein-DNA interactions: structural, thermodynamic and clustering patterns of conserved residues in DNA-binding proteins

Amino acid residues, which play important roles in protein function, are often conserved. Here, we analyze thermodynamic and structural data of protein-DNA interactions to explore a relationship between free energy, sequence conservation and structural cooperativity. We observe that the most stabilizing residues or putative hotspots are those which occur as clusters of conserved residues. The higher packing density of the clusters and available experimental thermodynamic data of mutations suggest cooperativity between conserved residues in the clusters. Conserved singlets contribute to the stability of protein-DNA complexes to a lesser extent. We also analyze structural features of conserved residues and their clusters and examine their role in identifying DNA-binding sites. We show that about half of the observed conserved residue clusters are in the interface with the DNA, which could be identified from their amino acid composition; whereas the remaining clusters are at the protein-protein or protein-ligand interface, or embedded in the structural scaffolds. In protein-protein interfaces, conserved residues are highly correlated with experimental residue hotspots, contributing dominantly and often cooperatively to the stability of protein-protein complexes. Overall, the conservation patterns of the stabilizing residues in DNA-binding proteins also highlight the significance of clustering as compared to single residue conservation.     

PPiClust: efficient clustering of 3D protein-protein interaction interfaces

The biological mechanisms through which proteins interact with one another are best revealed by studying the structural interfaces between interacting proteins. Protein-protein interfaces can be extracted from three-dimensional (3D) structural data of protein complexes and then clustered to derive biological insights. However, conventional protein interface clustering methods lack computational scalability and statistical support. In this work, we present a new method named "PPiClust" to systematically encode, cluster, and analyze similar 3D interface patterns in protein complexes efficiently. Experimental results showed that our method is effective in discovering visually consistent and statistically significant clusters of interfaces, and at the same time sufficiently time-efficient to be performed on a single computer. The interface clusters are also useful for uncovering the structural basis of protein interactions. Analysis of the resulting interface clusters revealed groups of structurally diverse proteins having similar interface patterns. We also found, in some of the interface clusters, the presence of well-known linear binding motifs which were noncontiguous in the primary sequences. These results suggest that PPiClust can discover not only statistically significant, but also biologically significant, protein interface clusters from protein complex structural data.     

A comparative analysis of interhelical polar interactions of various alpha-helix packings in proteins

One hundred twenty globular proteins and forty five "leucine zippers" representing all types of packing of long alpha-helices were studied in terms of revealing and comparing their interhelical hydrogen and salt bonds. Many previous studies of "leucine zippers" and their analogs showed that interhelical interactions between polar groups could impart specificity to packing of an alpha-helix. The current comparison demonstrated that basically, globular proteins and "leucine zippers" had similar interhelical polar interactions with presumably a similar structural role. However, depending on packing of alpha-helices, the networks of interhelical polar bonds were shown to be distinct and determined both by physicochemical properties of involved amino acid residues and by the relative positions of hydrophobic and hydrophilic residues on the surface of alpha-helices. The revealed distinction is probably crucial for selecting the unique packing of an alpha-helix.