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.     

Cavities and Atomic Packing in Protein Structures and Interfaces

A comparative analysis of cavities enclosed in a tertiary structure of proteins and interfaces formed by the interaction of two protein subunits in obligate and non-obligate categories (represented by homodimeric molecules and heterocomplexes, respectively) is presented. The total volume of cavities increases with the size of the protein (or the interface), though the exact relationship may vary in different cases. Likewise, for individual cavities also there is quantitative dependence of the volume on the number of atoms (or residues) lining the cavity. The larger cavities tend to be less spherical, solvated, and the interfaces are enriched in these. On average 15 ų of cavity volume is found to accommodate single water, with another 40–45 ų needed for each additional solvent molecule. Polar atoms/residues have a higher propensity to line solvated cavities. Relative to the frequency of occurrence in the whole structure (or interface), residues in β-strands are found more often lining the cavities, and those in turn and loop the least. Any depression in one chain not complemented by a protrusion in the other results in a cavity in the protein–protein interface. Through the use of the Voronoi volume, the packing of residues involved in protein–protein interaction has been compared to that in the protein interior. For a comparable number of atoms the interface has about twice the number of cavities relative to the tertiary structure.     

Cavities in protein–DNA and protein–RNA interfaces

An analysis of cavities present in protein–DNA and protein–RNA complexes is presented. In terms of the number of cavities and their total volume, the interfaces formed in these complexes are akin to those in transient protein–protein heterocomplexes. With homodimeric proteins protein–DNA interfaces may contain cavities involving both the protein subunits and DNA, and these are more than twice as large as cavities involving a single protein subunit and DNA. A parameter, cavity index, measuring the degree of surface complementarity, indicates that the packing of atoms in protein–protein/DNA/RNA is very similar, but it is about two times less efficient in the permanent interfaces formed between subunits in homodimers. As within the tertiary structure and protein–protein interfaces, protein–DNA interfaces have a higher inclination to be lined by β-sheet residues; from the DNA side, base atoms, in particular those in minor grooves, have a higher tendency to be located in cavities. The larger cavities tend to be less spherical and solvated. A small fraction of water molecules are found to mediate hydrogen-bond interactions with both the components, suggesting their primary role is to fill in the void left due to the local non-complementary nature of the surface patches.     

A review about nothing: Are apolar cavities in proteins really empty?

Cavities within proteins that are strictly apolar typically appear to be empty. It has been suggested, however, that water molecules may be present within such cavities but are too disordered to be seen in conventional crystallographic analyses. In contrast, it is argued here that solvent mobility will be limited by the size of the cavity and for this reason high‐occupancy solvent in cavities of typical volume should be readily detectable using X‐ray crystallography. Recent experimental studies of cavity hydration are reviewed. Such studies are consistent with theoretical predictions that it is energetically unfavorable to have a single water molecule in an apolar cavity. As apolar cavities become larger, a point is reached where it is favorable to have the cavity occupied by a cluster of mutually H‐bonded water molecules. The exact size of such a cavity in a protein is yet to be verified.     

Filling small, empty protein cavities: structural and energetic consequences

Most proteins contain small cavities that can be filled by replacing cavity-lining residues by larger ones. Since shortening mutations in hydrophobic cores tend to destabilize proteins, it is expected that cavity-filling mutations may conversely increase protein stability. We have filled three small cavities in apoflavodoxin and determined by NMR and equilibrium unfolding analysis their impact in protein structure and stability. The smallest cavity (14 A3) has been filled, at two different positions, with a variety of residues and, in all cases, the mutant proteins are locally unfolded, their structure and energetics resembling those of an equilibrium intermediate of the thermal unfolding of the wild-type protein. In contrast, two slightly larger cavities of 20 A3 and 21 A3 have been filled with Val to Ile or Val to Leu mutations and the mutants preserve both the native fold and the equilibrium unfolding mechanism. From the known relationship, observed in shortening mutations, between stability changes and the differential hydrophobicity of the exchanged residues and the volume of the cavities, the filling of these apoflavodoxin cavities is expected to stabilize the protein by approximately 1.5 kcal mol(-1). However, both urea and thermal denaturation analysis reveal much more modest stabilizations, ranging from 0.0 kcal mol(-1) to 0.6 kcal mol(-1), which reflects that the accommodation of single extra methyl groups in small cavities requires some rearrangement, necessarily destabilizing, that lowers the expected theoretical stabilization. As the size of these cavities is representative of that of the typical small, empty cavities found in most proteins, it seems unlikely that filling this type of cavities will give rise to large stabilizations.     

Protein-ligand dynamics. A 96 picosecond simulation of a myoglobin-xenon complex

A 96 picosecond dynamics trajectory of myoglobin with five xenon-probe ligands in internal cavities is examined to study the effect of protein motions on ligand motion and internal cavity fluctuations. Average structural and energetic properties indicate that the simulation is well behaved. The average protein volume is similar to the volume of the X-ray model and the main-chain atom root-mean-square deviation between the X-ray model and the average dynamical structure is 1.25 A. The protein volume oscillates 3 to 4% around the volume of the X-ray structure. These fluctuations lead to changes in the internal free volume and in the size, shape and location of atom-sized cavity features. Transient cavities produced in the simulation have a crucial role in the movement of two of the ligands. One of the ligands escapes to the protein surface, whilst a second ligand travels through the protein interior. Complex gating processes involving several protein residues are responsible for producing the necessary pores through which the ligand passes between transient cavities or packing defects.     

Size versus polarizability in protein-ligand interactions: binding of noble gases within engineered cavities in phage T4 lysozyme

To investigate the relative importance of size and polarizability in ligand binding within proteins, we have determined the crystal structures of pseudo wild-type and cavity-containing mutant phage T4 lysozymes in the presence of argon, krypton, and xenon. These proteins provide a representative sample of predominantly apolar cavities of varying size and shape. Even though the volumes of these cavities range up to the equivalent of five xenon atoms, the noble gases bind preferentially at highly localized sites that appear to be defined by constrictions in the walls of the cavities, coupled with the relatively large radii of the noble gases. The cavities within pseudo wild-type and L121A lysozymes each bind only a single atom of noble gas, while the cavities within mutants L133A and F153A have two independent binding sites, and the L99A cavity has three interacting sites. The binding of noble gases within two double mutants was studied to characterize the additivity of binding at such sites. In general, when a cavity in a protein is created by a "large-to-small" substitution, the surrounding residues relax somewhat to reduce the volume of the cavity. The binding of xenon and, to a lesser degree, krypton and argon, tend to expand the volume of the cavity and to return it closer to what it would have been had no relaxation occurred. In nearly all cases, the extent of binding of the noble gases follows the trend xenon>krypton>argon. Pressure titrations of the L99A mutant have confirmed that the crystallographic occupancies accurately reflect fractional saturation of the binding sites. The trend in noble gas affinity can be understood in terms of the effects of size and polarizability on the intermolecular potential. The plasticity of the protein matrix permits repulsion due to increased ligand size to be more than compensated for by attraction due to increased ligand polarizability. These results have implications for the mechanism of general anesthesia, the migration of small ligands within proteins, the detection of water molecules within apolar cavities and the determination of crystallographic phases.     

A simple semiempirical model for the effect of molecular confinement upon the rate of protein folding

A simple two-state model for the dependence of the rate of folding of a polypeptide confined within a spherical cavity upon the size of the cavity relative to that of the polypeptide is presented. A general prediction of the model is that decreasing the size of the cavity will increase the rate of refolding until the cavity becomes only slightly larger than the native state of the protein, at which point a further decrease in cavity size decreases the rate of refolding. The model qualitatively accounts for the behavior of several previously published simulations of folding within a cavity, as well as recently reported experimental measurements of the relative rate of refolding of each of five proteins encapsulated within wild-type and mutant GroEL-GroES complexes that have been engineered to provide internal cavities with similar surface composition and varying volume.     

Lattice models, packing density, and Boltzmann-like distribution of cavities in proteins

A model reproducing the experimental Boltzmann-like distribution of empty cavity sizes in proteins is introduced. Proteins are represented by lattices of different dimensionalities, corresponding to different numbers of nearest neighbor contacts. Small cavities emerge and join into larger ones in a random process that can be related to random mutations. Simulations of cavity creation are performed under the constraint of a limiting total packing density. Cavities sufficiently large (20 A(3) or more), that they might accommodate at least one additional methyl group produced by a mutation, are counted and compared to the distribution of cavities according to their sizes from protein statistics. The distributions calculated with this very simple model within a realistic range of packing densities are in good agreement with the empirical cavity distribution. The results suggest that the Boltzmann-like distribution of cavities in proteins might be affected by a mechanism controlled by limiting packing density and maximum allowed protein destabilization. This supports an earlier suggestion that the agreement between the free energies of cavity formation from the mutational experiments and from the statistics of the empty cavity distribution in X-ray protein structures is nonfortuitous. A possible relation of the suggested model to the Boltzmann hypothesis is discussed.     

Plant Proteins Are Smaller Because They Are Encoded by Fewer Exons than Animal Proteins

Protein size is an important biochemical feature since longer proteins can harbor more domains and therefore can display more biological functionalities than shorter proteins. We found remarkable differences in protein length, exon structure, and domain count among different phylo-genetic lineages. While eukaryotic proteins have an average size of 472 amino acid residues (aa), average protein sizes in plant genomes are smaller than those of animals and fungi. Proteins unique to plants are ?81 aa shorter than plant proteins conserved among other eukaryotic lineages. The smaller average size of plant proteins could neither be explained by endosymbiosis nor subcellular compartmentation nor exon size, but rather due to exon number. Metazoan proteins are encoded on average by ?10 exons of small size [?176 nucleotides (nt)]. Streptophyta have on average only ?5.7 exons of medium size (?230 nt). Multicellular species code for large proteins by increasing the exon number, while most unicellular organisms employ rather larger exons (>400 nt). Among sub-cellular compartments, membrane proteins are the largest (?520 aa), whereas the smallest proteins correspond to the gene ontology group of ribosome (?240 aa). Plant genes are encoded by half the number of exons and also contain fewer domains than animal proteins on average. Interestingly, endosymbiotic proteins that migrated to the plant nucleus became larger than their cyanobacterial orthologs. We thus conclude that plants have proteins larger than bacteria but smaller than animals or fungi. Compared to the average of eukaryotic species, plants have ?34%more but ?20%smal-ler proteins. This suggests that photosynthetic organisms are unique and deserve therefore special attention with regard to the evolutionary forces acting on their genomes and proteomes.     

Projected Availability of Natural Cavities for Wood Ducks in Southern Illinois

Abstract: Wood ducks (Aix sponsa) and other species use tree cavities in forested wetlands and adjacent upland forests for nest sites and cover. The availability of tree cavities suitable for nesting is important to the population dynamics of hole-nesting species, but there is little quantitative information on how forest succession and maturation affect densities of suitable nest sites in eastern deciduous forests. Several studies have measured availability of tree cavities for nesting wood ducks, but data on cavity formation and persistence rates are needed to model changes in cavity abundance. We measured abundance and persistence of tree cavities suitable for nesting wood ducks in southern Illinois, USA, during 1993-2002. We simulated changes in abundance of nest cavities in the Mississippi River floodplain and adjacent upland forests using estimates of tree cavity densities by tree-diameter size classes and 10-year cavity persistence rates by tree species. Cavities were disproportionately common in the largest size classes, but tree species varied in their propensity to form cavities. Beech (Fagus grandifolia; 0.41 cavities/tree) and sycamore (Plantanus occidentalis; 0.50 cavities/tree) were prolific cavity producers, whereas a small proportion (0.05 cavities/tree) of cottonwoods (Populus deltoides) contained cavities. Kaplan-Meier estimates of annual and 10-year cavity persistence averaged 0.95 and 0.64, respectively. Cavity persistence also differed among species (P = 0.02): cottonwoods had the lowest (0.54) and sycamores had the highest (0.89) 10-year tree cavity persistence rates. Tree fall (50.0%), cavity floor deterioration (37.5%), and narrowing of the cavity entrance (12.5%) were the most prevalent causes of tree cavity loss. Forest stand projections indicated that cavity abundance will increase up to 34% over recent levels during the first 10 years and by 44% after 50 years. Most of this increase will be contributed by tree species that are not commonly used by wood ducks, but cavities will increase in oaks (Quercus spp.) and beeches as the forest matures into cavity-bearing size classes. Sycamores will steadily contribute cavities, but cottonwood is predicted to provide fewer cavities due to low survival of cavity-bearing size classes. Our results suggest that availability of nest and den sites for cavity-dependent wildlife will increase as eastern deciduous forests mature over the next half century. Cost-effectiveness of artificial nest box programs should be reevaluated in light of projected changes in tree cavity availability as deciduous forests mature in the eastern United States.     

Modeling the evolution of protein domain architectures using maximum parsimony

Domains are basic evolutionary units of proteins and most proteins have more than one domain. Advances in domain modeling and collection are making it possible to annotate a large fraction of known protein sequences by a linear ordering of their domains, yielding their architecture. Protein domain architectures link evolutionarily related proteins and underscore their shared functions. Here, we attempt to better understand this association by identifying the evolutionary pathways by which extant architectures may have evolved. We propose a model of evolution in which architectures arise through rearrangements of inferred precursor architectures and acquisition of new domains. These pathways are ranked using a parsimony principle, whereby scenarios requiring the fewest number of independent recombination events, namely fission and fusion operations, are assumed to be more likely. Using a data set of domain architectures present in 159 proteomes that represent all three major branches of the tree of life allows us to estimate the history of over 85% of all architectures in the sequence database. We find that the distribution of rearrangement classes is robust with respect to alternative parsimony rules for inferring the presence of precursor architectures in ancestral species. Analyzing the most parsimonious pathways, we find 87% of architectures to gain complexity over time through simple changes, among which fusion events account for 5.6 times as many architectures as fission. Our results may be used to compute domain architecture similarities, for example, based on the number of historical recombination events separating them. Domain architecture "neighbors" identified in this way may lead to new insights about the evolution of protein function.     

Horizontal gene transfers as metagenomic gene duplications

While it is well accepted that horizontal gene transfer plays an important role in the evolution and the diversification of prokaryotic genomes, many questions remain open regarding its functional mechanisms of action and its interplay with the extant genome. This study addresses the relationship between proteome innovation by horizontal gene transfer and genome content in Proteobacteria. We characterize the transferred genes, focusing on the protein domain compositions and their relationships with the existing protein domain superfamilies in the genome. In agreement with previous observations, we find that the protein domain architectures of horizontally transferred genes are significantly shorter than the genomic average. Furthermore, protein domains that are more common in the total pool of genomes appear to have a proportionally higher chance to be transferred. This suggests that transfer events behave as if they were drawn randomly from a cross-genomic community gene pool, much like gene duplicates are drawn from a genomic gene pool. Finally, horizontally transferred genes carry domains of exogenous families less frequently for larger genomes, although they might do it more than expected by chance.     

Supra-domains: evolutionary units larger than single protein domains

Domains are the evolutionary units that comprise proteins, and most proteins are built from more than one domain. Domains can be shuffled by recombination to create proteins with new arrangements of domains. Using structural domain assignments, we examined the combinations of domains in the proteins of 131 completely sequenced organisms. We found two-domain and three-domain combinations that recur in different protein contexts with different partner domains. The domains within these combinations have a particular functional and spatial relationship. These units are larger than individual domains and we term them "supra-domains". Amongst the supra-domains, we identified some 1400 (1203 two-domain and 166 three-domain) combinations that are statistically significantly over-represented relative to the occurrence and versatility of the individual component domains. Over one-third of all structurally assigned multi-domain proteins contain these over-represented supra-domains. This means that investigation of the structural and functional relationships of the domains forming these popular combinations would be particularly useful for an understanding of multi-domain protein function and evolution as well as for genome annotation. These and other supra-domains were analysed for their versatility, duplication, their distribution across the three kingdoms of life and their functional classes. By examining the three-dimensional structures of several examples of supra-domains in different biological processes, we identify two basic types of spatial relationships between the component domains: the combined function of the two domains is such that either the geometry of the two domains is crucial and there is a tight constraint on the interface, or the precise orientation of the domains is less important and they are spatially separate. Frequently, the role of the supra-domain becomes clear only once the three-dimensional structure is known. Since this is the case for only a quarter of the supra-domains, we provide a list of the most important unknown supra-domains as potential targets for structural genomics projects.     

A guest molecule-host cavity fitting algorithm to mine PDB for small molecule targets

Inhaled anesthetic molecule occupancy of a protein internal cavity depends in part on the volumes of the guest molecule and the host site. Current algorithms to determine volume and surface area of cavities in proteins whose structures have been determined and cataloged make no allowance for shape or small degrees of shape adjustment to accommodate a guest. We developed an algorithm to determine spheroid dimensions matching cavity volume and surface area and applied it to screen the cavities of 6,658 nonredundant structures stored in the Protein Data Bank (PDB) for potential targets of halothane (2-bromo-2-chloro-1,1,1-trifluoroethane). Our algorithm determined sizes of prolate and oblate spheroids matching dimensions of each cavity found. If those spheroids could accommodate halothane (radius 2.91 A) as a guest, we determined the packing coefficient. 394,766 total cavities were identified. Of 58,681 cavities satisfying the fit criteria for halothane, 11,902 cavities had packing coefficients in the range of 0.46-0.64. This represents 20.3% of cavities large enough to hold halothane, 3.0% of all cavities processed, and found in 2,432 protein structures. Our algorithm incorporates shape dependence to screen guest-host relationships for potential small molecule occupancy of protein cavities. Proteins with large numbers of such cavities are more likely to be functionally altered by halothane.     

Clutch size relative to tree cavity size in Northern Flickers

We analysed clutch size versus nest size in 153 broods of the Northern Flicker Colaptes auratus , a woodpecker using natural cavities in British Columbia, Canada. Larger volume cavities were less susceptible to predation and cavity size was positively associated with the age and body size of males and with the body condition of female parents. Although clutches varied between 4 and 11 eggs, and the floor area of cavities varied about 5-fold, we found no relationship between clutch size and floor area or cavity volume. To see if there were fitness consequences to clutch size relative to nest size, we examined hatching success and nestling mortality in flicker broods. Hatching success was not related to cavity size, but crowding slightly reduced nestling survival even when clutch size was controlled statistically. However, there was no effect of cavity size on the total number of nestlings fledged. Newly excavated flicker cavities were smaller than reused cavities suggesting a cost to excavation. This cost, coupled with the minimal fitness consequences of overcrowding, may explain why flickers do not adjust clutch size to cavity size.     

Analysis of void volumes in proteins and application to stability of the p53 tumour suppressor protein

We have developed a new method for the analysis of voids in proteins (defined as empty cavities not accessible to solvent). This method combines analysis of individual discrete voids with analysis of packing quality. While these are different aspects of the same effect, they have traditionally been analysed using different approaches. The method has been applied to the calculation of total void volume and maximum void size in a non-redundant set of protein domains and has been used to examine correlations between thermal stability and void size. The tumour-suppressor protein p53 has then been compared with the non-redundant data set to determine whether its low thermal stability results from poor packing. We found that p53 has average packing, but the detrimental effects of some previously unexplained mutations to p53 observed in cancer can be explained by the creation of unusually large voids.     

Binding cavities and druggability of intrinsically disordered proteins

To assess the potential of intrinsically disordered proteins (IDPs) as drug design targets, we have analyzed the ligand-binding cavities of two datasets of IDPs (containing 37 and 16 entries, respectively) and compared their properties with those of conventional ordered (folded) proteins. IDPs were predicted to possess more binding cavity than ordered proteins at similar length, supporting the proposed advantage of IDPs economizing genome and protein resources. The cavity number has a wide distribution within each conformation ensemble for IDPs. The geometries of the cavities of IDPs differ from the cavities of ordered proteins, for example, the cavities of IDPs have larger surface areas and volumes, and are more likely to be composed of a single segment. The druggability of the cavities was examined, and the average druggable probability is estimated to be 9% for IDPs, which is almost twice that for ordered proteins (5%). Some IDPs with druggable cavities that are associated with diseases are listed. The optimism versus obstacles for drug design for IDPs is also briefly discussed.     

Dynamics alignment: comparison of protein dynamics in the SCOP database

A novel methodology for comparison of protein dynamics is presented. Protein dynamics is calculated using the Gaussian network model and the modes of motion are globally aligned using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. The alignment is fast and can be used to analyze large sets of proteins. The methodology is applied to the four major classes of the SCOP database: "all alpha proteins," "all beta proteins," "alpha and beta proteins," and "alpha/beta proteins". We show that different domains may have similar global dynamics. In addition, we report that the dynamics of "all alpha proteins" domains are less specific to structural variations within a given fold or superfamily compared with the other classes. We report that domain pairs with the most similar and the least similar global dynamics tend to be of similar length. The significance of the methodology is that it suggests a new and efficient way of mapping between the global structural features of protein families/subfamilies and their encoded dynamics.     

Internal cavities and buried waters in globular proteins

A fast algorithm that detects internal cavities in proteins and predicts the positions of buried water molecules is described. The cavities are characterized in terms of volume, surface area, polarity, and the presence of bound waters. The algorithm is applied to 12 proteins whose structures are known to high resolution and successfully predicts the locations of over 80% of internal water molecules. Most proteins are found to have a number of internal cavities ranging in volume from 10 to 180 A3. Some of these cavities contain water and some do not, with the probability of containing a buried water increasing with cavity size. However, many large cavities are found to be empty (i.e., they do not contain a crystallographically determined water). For multidomain proteins over half of the total cavity volume is at the interdomain interface. Possible implications for the energetics of cavity formation and for the functional role of internal cavities are discussed.