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.     

Structure of the YibK methyltransferase from Haemophilus influenzae (HI0766): a cofactor bound at a site formed by a knot

Despite the suitability of various lattice geometries for coarse-grained modeling of proteins, the actual packing geometry of residues in folded structures has remained largely unexplored. A strong tendency to assume a regular packing geometry is shown here by optimally reorienting and superimposing clusters of neighboring residues from databank structures examined on a coarse-grained (single-site-per-residue) scale. The orientation function (or order parameter) of the examined coordination clusters with respect to fcc lattice directions is found to be 0.82. The observed geometry, which may be termed an incomplete distorted face-centered cubic (fcc) packing, is apparently favored by the drive to maximize packing density, in a fashion analogous to the way identical spheres pack densely and follow fcc geometry. About 2/3 of all residues obey this packing geometry, while the remainder occupy other context-dependent positions. The preferred coordination directions show relatively small variations over the various amino acid types, consistent with uniform residue viewpoint. Both the extremes of solvent-exposed and completely buried residue neighborhoods approximate the same generic packing, the only difference being in the numbers (and not the orientations) of coordination sites that are occupied (or left void for solvent occupancy). We observe the prevalence of a rather uniform (tight) residue packing density throughout the structure, including even the residues packed near solvent-exposed regions. The observed orientation distribution reveals an underlying, intrinsic orientation lattice for proteins.     

Incomplete protein packing as a selectivity filter in drug design

The conservation of structure across paralog proteins promotes alternative protein-ligand associations often leading to side effects in drug-based inhibition. However, sticky packing defects are typically not conserved across paralogs, making them suitable targets to reduce drug toxicity. This observation enables a strategy for the design of highly specific inhibitors involving ligands that wrap nonconserved packing defects. The selectivity of these inhibitors is evidenced in affinity assays on a cancer-related pharmacokinome: a powerful inhibitor is redesigned by using the wrapping technology to enhance its selectivity and affinity for a target kinase. In this way, the packing defects of a soluble protein may be used as selectivity filters for drug design.     

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.      

Cavities and packing defects in the structural dynamics of myoglobin

Small globular proteins contain internal cavities and packing defects that reduce thermodynamic stability but seem to play a role in controlling function by defining pathways for the diffusion of the ligand/substrate to the active site. In the case of myoglobin (Mb), a prototype for structure–function relationship studies, the photosensitivity of the adduct of the reduced protein with CO, O₂ and NO allows events related to the migration of the ligand through the matrix to be followed. The crystal structures of intermediate states of wild-type (wt) and mutant Mbs show the photolysed CO to be located either in the distal heme pocket (primary docking site) or in one of two alternative cavities (secondary docking sites) corresponding to packing defects accessible to an atom of xenon. These results convey the general picture that pre-existing internal cavities are involved in controlling the dynamics and reactivity of the reactions of Mb with O₂ and other ligands, including NO.     

Crystal packing effects on protein loops

The effects of crystal packing on protein loop structures are examined by (1) a comparison of loops in proteins that have been crystallized in alternate packing arrangements, and (2) theoretical prediction of loops both with and without the inclusion of the crystal environment. Results show that in a minority of cases, loop geometries are dependent on crystal packing effects. Explicit representation of the crystal environment in a loop prediction algorithm can be used to model these effects and to reconstruct the structures, and relative energies, of a loop in alternative packing environments. By comparing prediction results with and without the inclusion of the crystal environment, the loop prediction algorithm can further be used to identify cases in which a crystal structure does not represent the most stable state of a loop in solution. We anticipate that this capability has implications for structural biology.     

The interpretation of protein structures: total volume, group volume distributions and packing density

Following the general procedure of Bernal &; Finney (1967) using Voronoi polyhedra, volumes occupied by all the atoms, or groups of atoms, in lysozyme and ribonuclease S have been estimated from the atomic co-ordinates provided by crystal structure studies. The average packing density for the interior of both proteins is close to 0.75, which is in the middle of the range found for crystals of most small organic molecules. For all atom types the mean packing densities fall between 0.7 and 0.8 with standard deviations between ± 0.1 and ± 0.2. It is suggested that simple geometrical packing considerations may provide useful criteria in guiding and evaluating trial structures in theoretical studies of protein folding, especially the association of distant parts of a peptide chain.Packing densities averaged over a relatively small number of atoms (5 to 15) appear to vary substantially in different parts of the same protein. Low densities representing packing defects may permit relatively easy motions, for example in an active site. Surrounding areas of high density may serve as relatively incompressible regions which transmit or correlate motions over considerable distances.The total volume of each of these two proteins as derived from the co-ordinate list appears to be larger than that estimated from the partial specific volume by 7 to 10%. If this volume difference is attributed to a change in the packing of water in a monolayer surrounding the protein, it would correspond to an average decrease, relative to bulk water, of 1 to 2 ų per water molecule in this monolayer. Such a change is about half of the decrease that occurs on the melting of ice.     

Cavity defects in the procapsid of bacteriophage P22 and the mechanism of capsid maturation

Bacteriophage P22 belongs to a family of double-stranded DNA viruses that share common morphogenetic features like DNA packaging into a procapsid precursor and maturation. Maturation involves cooperative expansion of the procapsid shell with concomitant lattice stabilization. The expansion is thought to be mediated by movement of two coat protein domains around a hinge. The metastable conformation of subunit within the procapsid lattice is considered to constitute a late folding intermediate. In order to understand the mechanism of expansion it is necessary to characterize the interactions stabilizing procapsid and mature capsid lattices, respectively. We employ pressure dissociation to compare subunit packing within the procapsid and expanded lattice. Procapsid shells contain larger cavities than the expanded shells, presumably due to polypeptide packing defects. These defects contribute to the metastable nature of the procapsid lattice and are cured during expansion. Improved packing contributes to the increased stability of the expanded shell. Comparison of two temperature-sensitive folding (tsf) mutants of coat protein (T294I and W48Q) with wild-type coat revealed that both mutations markedly destabilized the procapsid shell and yet had little effect on relative stability of the monomeric subunit. Thus, the regions affected by these packing defects constitute subunit interfaces of the procapsid shell. The larger activation volume of pressure dissociation observed for both T294I and W48Q indicates that the decreased stability of these particles is due to increase of cavity defects. These defects in the procapsid lattice are cured upon expansion suggesting that the intersubunit contacts affected by tsf mutations are absent or rearranged in the mature shell. The energetics of the in vitro expansion reaction also suggests that entropic stabilization contributes to the large free energy barrier for expansion.     

Reducing PEX13 expression ameliorates physiological defects of late-acting peroxin mutants

Proteins are targeted to the peroxisome matrix via processes that are mechanistically distinct from those used by other organelles. Protein entry into peroxisomes requires peroxin (PEX) proteins, including early-acting receptor (e.g. PEX5) and docking peroxins (e.g. PEX13 and PEX14) and late-acting PEX5-recycling peroxins (e.g. PEX4 and PEX6). We examined genetic interactions among Arabidopsis peroxin mutants and found that the weak pex13-1 allele had deleterious effects when combined with pex5-1 and pex14-2, which are defective in early-acting peroxins, as shown by reduced matrix protein import and enhanced physiological defects. In contrast, combining pex13-1 with pex4-1 or pex6-1, which are defective in late-acting peroxins, unexpectedly ameliorated mutant growth defects. Matrix protein import remained impaired in pex4-1 pex13-1 and pex6-1 pex13-1, suggesting that the partial suppression of pex4-1 and pex6-1 physiological defects by a weak pex13 allele may result from restoring the balance between import and export of PEX5 or other proteins that are retrotranslocated from the peroxisome with the assistance of PEX4 and PEX6. Our results suggest that symptoms caused by pex mutants defective in late-acting peroxins may result not only from defects in matrix protein import but also from inefficient removal of PEX5 from the peroxisomal membrane following cargo delivery.     

Receptor antagonist and selective agonist derivatives of mouse interleukin-2.

Mouse interleukin-2 (mIL-2) proteins with substitutions at two residues (D34 and Q141) that interact specifically with different signalling subunits (respectively, beta and gamma) of the IL-2 receptor (IL-2R) were examined using several in vitro cellular assays. Proteins with specific substitutions at both residues were partial agonists and their maximal responses varied widely in different IL-2-responsive cell types. Two of these cell types had comparable numbers of IL-2R and similar affinities for wild-type mIL-2 and mutant mIL-2 proteins. However, the more responsive cell type had 'spare' IL-2R. Various mIL-2 proteins with substitutions at Q141 had modest defects in IL-2R-binding and were potent antagonists of native mIL-2 action. Proteins with bulky or basic substitutions at residue D34 were weak antagonists due to severely reduced IL-2 binding and their reduced binding paralleled their defects in IL-2R activation. Our results suggest that interaction of mIL-2 with IL-2R beta is more important for binding than activation and that the converse holds for mIL-2 interaction with IL-2R gamma. Also genetic manipulation of the interaction of IL-2 with IL-2R beta and IL-2R gamma has led to the discovery of potentially useful IL-2 antagonists and selective agonists.     

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.     

Trends in rates of multiple vascular disruption defects, Atlanta, 1968-1989: is there evidence of a cocaine teratogenic epidemic?

Research suggests that, perhaps through mechanisms initiated by vasoconstriction and leading to vessel thrombosis or embolism, cocaine causes vascular disruption defects, and that frequent cocaine use during early pregnancy could disrupt multiple organ systems in the fetus. We hypothesized that if cocaine is an important cause of multiple vascular disruption defects, a rising prevalence of cocaine use by mothers during pregnancy should be accompanied by rising rates of these defects in their offspring. Using data from the Metropolitan Atlanta Congenital Defects Program, we identified all infants born in Atlanta from 1968 through 1989 who had nonsyndromic, provisional vascular disruption defects affecting more than one organ system: 61 infants (78%) had gastrointestinal and genitourinary defects, 7 (9%) had gastrointestinal and abdominal wall defects, 2 (3%) had gastrointestinal and limb reduction defects, 2 (3%) had limb reduction and abdominal wall defects, 2 (3%) had central nervous system and gastrointestinal defects, 2 (3%) had genitourinary and limb reduction defects, 1 (1%) had genitourinary and abdominal wall defects, and 1 (1%) had central nervous system and genitourinary defects. The prevalence of Atlanta infants with more than one vascular disruption defect is 0.13 per 1,000 live births. Chi-square analysis for trends showed no increase in prevalence during the study period. Our data are from one of the first population-based studies in which trends for defects potentially caused by maternal cocaine use are examined; the results of our study show no significant change in the prevalence of multiple vascular disruption defects over time.(ABSTRACT TRUNCATED AT 250 WORDS)     

Calculation of protein volumes: An alternative to the Voronoi procedure

All the methods that have been applied to assessing local volume occupation or packing in proteins have particular defects. For example, the Voronoi method used by Richards (method A) and by Finney misallocates both non-bonded and covalent contacts in a geometrically rigorous, though chemically inconsistent, manner. Richards' method B, in which covalent and non-bonded contacts are partitioned in chemically sensible ways, is unfortunately not completely rigorous, in that every polyhedron vertex has associated with it a small vertex-error polyhedron, which is not allocated to any atom.We present here a generalization of the Voronoi method that is particularly suited to multicomponent assemblies such as proteins. This radical plane partitioning of volume is completely rigorous; it gives rise to no vertex error, and handles the more numerous non-bonded contacts realistically. Its application to RNAase-S is described and the results compared with both Voronoi's method and Richards' method B. A particular advantage of both the radical plane and Richards' methods is a relative insensitivity to the treatment of the surface, a problem that has plagued other approaches to describing packing in proteins. Although the radical plane is seen to misallocate volume chemically between covalently-bonded neighbours, this problem vanishes when groups of atoms in side-chain residues are considered.     

An AFM study of solid-phase bilayers of unsaturated PC lipids and the lateral distribution of the transmembrane model peptide WALP23 in these bilayers

An altered lipid packing can have a large influence on the properties of the membrane and the lateral distribution of proteins and/or peptides that are associated with the bilayer. Here, it is shown by contact-mode atomic force microscopy that the surface topography of solid-phase bilayers of PC lipids with an unsaturated cis bond in their acyl chains shows surfaces with a large number of line-type packing defects, in contrast to the much smoother surfaces observed for saturated PC lipids. Di-n:1-PC (n = 20, 22, 24) and (16:0,18:1)-PC (POPC) were used. Next, the influence of an altered lipid environment on the lateral distribution of the single α-helical model peptide WALP23 was studied by incorporating the peptide in the bilayers of di-n:1-PC (n = 20, 22, 24) and (16:0,18:1)-PC unsaturated lipids. The presence of WALP23 leads to an increase in the number of packing defects but does not lead to the formation of the striated domains that were previously observed in bilayers of saturated PC lipids and WALP. This is ascribed to the less efficient lateral lipid packing of the unsaturated lipids, while the increase in packing defects is probably an indirect effect of the peptide. Finally, the fact that an altered lipid packing affects the distribution of WALP23 is also confirmed in an additional experiment where the solvent TFE (2,2,2-trifluorethanol) is added to bilayers of di-16:0-PC/WALP23. At 3.5 vol% TFE, the previous striated ordering of the peptide is abolished and replaced by loose lines.     

The role of ion-lipid interactions and lipid packing in transient defects caused by phenolic compounds

The transient disruption of membranes for the passive permeation of ions or small molecules is a complex process relevant to understanding physiological processes and biotechnology applications. Phenolic compounds are widely studied for their antioxidant and antimicrobial properties, and some of these activities are based on the interactions of the phenolic compound with membranes. Ions are ubiquitous in cells and are known to alter the structure of phospholipid bilayers. Yet, ion-lipid interactions are usually ignored when studying the membrane-altering properties of phenolic compounds. This study aims to assess the role of Ca²⁺ ions on the membrane-disrupting activity of two phenolic acids and to highlight the role of local changes in lipid packing in forming transient defects or pores. Results from tethered bilayer lipid membrane electrical impedance spectroscopy experiments showed that Ca²⁺ significantly reduces membrane disruption by caffeic acid methyl ester and caffeic acid. As phenolic acids are known metal chelators, we used UV-vis and fluorescence spectroscopy to exclude the possibility that Ca²⁺ interferes with membrane disruption by binding to the phenolic compound and subsequently preventing membrane binding. Molecular dynamics simulations showed that Ca²⁺ but not caffeic acid methyl ester or caffeic acid increases lipid packing in POPC bilayers. The combined data confirm that Ca²⁺ reduces the membrane-disrupting activity of the phenolic compounds, and that Ca²⁺-induced changes to lipid packing govern this effect. We discuss our data in the context of ion-induced pores and transient defects and how lipid packing affects membrane disruption by small molecules.     

Feature-similarity protein classifier as a ligand engineering tool

Kinases have been often targeted in drug therapy aimed at blocking signaling pathways. However, the conservation of protein structure across homologs often leads to uncontrolled cross-reactivity. On the other hand, sticky packing defects in proteins are typically not conserved across homologs, making them ligand-anchoring sites potentially important to enhance selectivity. Thus, we introduce a hierarchical clustering of PDB-reported kinases according to packing differences. This kinome partitioning is highly correlated with proximity relations arising from the pharmacological profiling of kinases. A variable packing sensitivity is observed for individual drugs, with highly promiscuous ligands being the most insensitive to packing differences. Our classifier enables a strategy to design selective inhibitors.     

Mechanism of membrane curvature sensing by amphipathic helix containing proteins

There are several examples of membrane-associated protein domains that target curved membranes. This behavior is believed to have functional significance in a number of essential pathways, such as clathrin-mediated endocytosis, which involve dramatic membrane remodeling and require the recruitment of various cofactors at different stages of the process. This work is motivated in part by recent experiments that demonstrated that the amphipathic N-terminal helix of endophilin (H0) targets curved membranes by binding to hydrophobic lipid bilayer packing defects which increase in number with increasing membrane curvature. Here we use state-of-the-art atomistic simulation to explore the packing defect structure of curved membranes, and the effect of this structure on the folding of H0. We find that not only are packing defects increased in number with increasing membrane curvature, but also that their size distribution depends nontrivially on the curvature, falling off exponentially with a decay constant that depends on the curvature, and crucially that even on highly curved membranes defects large enough to accommodate the hydrophobic face of H0 are never observed. We furthermore find that a percolation model for the defects explains the defect size distribution, which implies that larger defects are formed by coalescence of noninteracting smaller defects. We also use the recently developed metadynamics algorithm to study in detail the effect of such defects on H0 folding. It is found that the comparatively larger defects found on a convex membrane promote H0 folding by several kcal/mol, while the smaller defects found on flat and concave membrane surfaces inhibit folding by kinetically trapping the peptide. Together, these observations suggest H0 folding is a cooperative process in which the folding peptide changes the defect structure relative to an unperturbed membrane.     

Role of secondary structure in protein-phospholipid surface interactions: reconstitution and denaturation of apolipoprotein C-I:DMPC complexes

Binding of protein to a phospholipid surface is commonly mediated by amphipathic alpha-helices. To understand the role of alpha-helical structure in protein-lipid interactions, we used discoidal lipoproteins reconstituted from dimyristoylphosphatidylcholine (DMPC) and human apolipoprotein C-I (apoC-I, 6 kDa) or its mutants containing single Pro substitutions along the sequence and differing in their alpha-helical content in solution (0-48%) and on DMPC (40-75%). Thermal denaturation revealed that lipoprotein stability correlates weakly with the protein helix content: proteins with higher alpha-helical content on DMPC may form more stable complexes. Lipoprotein reconstitution upon cooling from the heat-denatured state and DMPC clearance studies revealed that protein secondary structure in solution and on DMPC correlates strongly with the maximal temperature of lipoprotein reconstitution: more helical proteins can reconstitute lipoproteins at higher temperatures. Interestingly, at Tc = 24 degrees C of the DMPC gel-to-liquid crystal transition, the clearance rate is independent of the protein helical content. Consequently, if the packing defects at the phospholipid surface are readily available (e.g., at the lipid phase boundary), insertion of protein into these defects is independent of the secondary structure in solution. However, if hydrophobic defects are limited, protein binding and insertion are aided by other surface-bound proteins and depend on their helical propensity: the larger the propensity, the faster the binding and the broader its temperature range. This positive cooperativity in binding of alpha-helices to phospholipid surface, which may result from direct and/or lipid-mediated protein-protein interactions, may be important for lipoprotein metabolism and for protein-membrane binding.