Crystal structure of the EF-hand parvalbumin at atomic resolution (0.91 A) and at low temperature (100 K). Evidence for conformational multistates within the hydrophobic core.

Several crystal structures of parvalbumin (Parv), a typical EF-hand protein, have been reported so far for different species with the best resolution achieving 1.5 A. Using a crystal grown under microgravity conditions, cryotechniques (100 K), and synchrotron radiation, it has now been possible to determine the crystal structure of the fully Ca2+-loaded form of pike (component pI 4.10) Parv.Ca2 at atomic resolution (0.91 A). The availability of such a high quality structure offers the opportunity to contribute to the definition of the validation tools useful for the refinement of protein crystal structures determined to lower resolution. Besides a better definition of most of the elements in the protein three-dimensional structure than in previous studies, the high accuracy thus achieved allows the detection of well-defined alternate conformations, which are observed for 16 residues out of 107 in total. Among them, six occupy an internal position within the hydrophobic core and converge toward two small buried cavities with a total volume of about 60 A3. There is no indication of any water molecule present in these cavities. It is probable that at temperatures of physiological conditions there is a dynamic interconversion between these alternate conformations in an energy-barrier dependent manner. Such motions for which the amplitudes are provided by the present study will be associated with a time-dependent remodeling of the void internal space as part of a slow dynamics regime (millisecond timescales) of the parvalbumin molecule. The relevance of such internal dynamics to function is discussed.     

Comparison of solution and crystal structures of maize nonspecific lipid transfer protein: A model for a potential in vivo lipid carrier protein

The three-dimensional solution structure of maize nonspecific lipid transfer protein (nsLTP) obtained by nuclear magnetic resonance (NMR) is compared to the X-ray structure. Although both structures are very similar, some local structural differences are observed in the first and the fourth helices and in several side-chain conformations. These discrepancies arise partly from intermolecular contacts in the crystal lattice. The main characteristic of nsLTP structures is the presence of an internal hydrophobic cavity whose volume was found to vary from 237 to 513 ų without major variations in the 15 solution structures. Comparison of crystal and NMR structures shows the existence of another small hollow at the periphery of the protein containing a water molecule in the X-ray structure, which could play an important structural role. A model of the complexed form of maize nsLTP by α-lysopalmitoylphosphatidylcholine was built by docking the lipid inside the protein cavity of the NMR structure. The main structural feature is a hydrogen bond found also in the X-ray structure of the complex maize nsLTP/palmitate between the hydroxyl of Tyr81 and the carbonyl of the lipid. Comparison of 12 primary sequences of nsLTPs emphasizes that all residues delineating the cavities calculated on solution and X-ray structures are conserved, which suggests that this large cavity is a common feature of all compared plant nsLTPs. Furthermore several conserved basic residues seem to be involved in the stabilization of the protein architecture. Proteins 31:160–171, 1998. © 1998 Wiley-Liss, Inc.     

Complexes of porcine odorant binding protein with odorant molecules belonging to different chemical classes

Porcine odorant binding protein (pOBP) is a monomer of 157 amino acid residues, purified in abundance from pig nasal mucosa. In contrast to the observation on lipocalins as retinol binding protein (RBP), major urinary protein (MUP) or bovine odorant binding protein (bOBP), no naturally occurring ligand was found in the beta-barrel cavity of pOBP. Porcine OBP was therefore chosen as a simple model for structure/function studies with odorant molecules. In competition experiments with tritiated pyrazine, the affinity of pOBP towards several odorant molecules belonging to different chemical classes has been found to be of the micromolar order, with a 1:1 stoichiometry. The X-ray structures of pOBP complexed to these molecules were determined at resolution between 2.15 and 1.4 A. As expected, the electron density of the odorant molecules was observed into the hydrophobic beta-barrel of the lipocalin. Inside this cavity, very few specific interactions were established between the odorant molecule and the amino acid side-chains, which did not undergo significant conformational change. The high B-factors observed for the odorant molecules as well as the existence of alternative conformations reveal a non-specific mode of binding of the odorant molecules in the cavity.     

Effect of xenon binding to a hydrophobic cavity on the proton pumping cycle in bacteriorhodopsin

To understand the functional role of apolar cavities in bacteriorhodopsin, a light-driven proton pump found in Halobacterium salinarum, we investigated the crystal structure in pressurized xenon or krypton. Diffraction data from the P622 crystal showed that one Xe or Kr atom binds to a preexisting hydrophobic cavity buried between helices C and D, located at the same depth from the membrane surface as Asp96, a key residue in the proton uptake pathway. The occupation fraction of Xe or Kr was calculated as approximately 0.32 at a pressure of 1 MPa. In the unphotolyzed state, the binding of Xe or Kr caused no large deformation of the cavity. However, the proton pumping cycle was greatly perturbed when an aqueous suspension of purple membrane was pressurized with xenon gas; that is, the decay of the M state was accelerated significantly (~5 times at full occupancy), while the decay of an equilibrium state of N and O was slightly decelerated. A similar but much smaller perturbation in the reaction kinetics was observed upon pressurization with krypton gas. In a glycerol/water mixture, xenon-induced acceleration of M decay became less significant in proportion to the water activity. Together with the structure of the xenon-bound protein, these observations suggest that xenon binding helps water molecules permeate into apolar cavities in the proton uptake pathway, thereby accelerating the water-mediated proton transfer from Asp96 to the Schiff base.     

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.     

Automated electron-density sampling reveals widespread conformational polymorphism in proteins

Although proteins populate large structural ensembles, X-ray diffraction data are traditionally interpreted using a single model. To search for evidence of alternate conformers, we developed a program, Ringer, which systematically samples electron density around the dihedral angles of protein side chains. In a diverse set of 402 structures, Ringer identified weak, nonrandom electron-density features that suggest of the presence of hidden, lowly populated conformations for >18% of uniquely modeled residues. Although these peaks occur at electron-density levels traditionally regarded as noise, statistically significant (P < 10⁻⁵) enrichment of peaks at successive rotameric χ angles validates the assignment of these features as unmodeled conformations. Weak electron density corresponding to alternate rotamers also was detected in an accurate electron density map free of model bias. Ringer analysis of the high-resolution structures of free and peptide-bound calmodulin identified shifts in ensembles and connected the alternate conformations to ligand recognition. These results show that the signal in high-resolution electron density maps extends below the traditional 1 σ cutoff, and crystalline proteins are more polymorphic than current crystallographic models. Ringer provides an objective, systematic method to identify previously undiscovered alternate conformations that can mediate protein folding and function.     

Homology modelling of the Frankia nitrogenase iron protein

The NifH protein contains an iron-sulfur cluster performing different functions during nitrogen fixation. Frankia is an actinomycete, entering into symbiotic association with a number of dicotyledonous plants and fixing nitrogen. The structure of the Frankia NifH protein was determined using homology modelling technique. Metal binding sites and functionally important regions of the protein were analyzed. Thiol ligands and active sites help in protein functioning and conformations. Structurally important nests were recognized. Clefts and cavities contain biologically important residues. Site-directed mutagenesis results reveal that mutations in functional residues hamper nitrogen fixation. The structure is rigid with an accessible surface for solvents. The structure is reliable offering insights into the 3D structural framework as well as structure-function relation of NifH protein.     

A molecular dynamics study of thermodynamic and structural aspects of the hydration of cavities in proteins.

The structure and activity of a protein molecule are strongly influenced by the extent of hydration of its cavities. This is, in turn, related to the free energy change on transfer of a water molecule from bulk solvent into a cavity. Such free energy changes have been calculated for two cavities in a sulfate-binding protein. One of these cavities contains a crystallographically observed water molecule while the other does not. Thermodynamic integration and perturbation methods were used to calculate free energies of hydration for each of the cavities from molecular dynamics simulations of two separate events: the removal of a water molecule from pure water, and the introduction of a water molecule into each protein cavity. From the simulations for the pure water system, the excess chemical potential of water was computed to be -6.4 +/- 0.4 kcal/mol, in accord with experiment and with other recent theoretical calculations. For the protein cavity containing an experimentally observed water molecule, the free energy change on hydrating it with one water molecule was calculated as -10.0 +/- 1.3 kcal/mol, indicating the high probability that this cavity is occupied by a water molecule. By contrast, for the cavity in which no water molecules were experimentally observed, the free energy change on hydrating it with one water molecule was calculated as 0.2 +/- 1.5 kcal/mol, indicating its low occupancy by water. The agreement of these results with experiment suggests that thermodynamic simulation methods may become useful for the prediction and analysis of internal hydration in proteins.     

Sub-atomic resolution crystal structure of cholesterol oxidase: what atomic resolution crystallography reveals about enzyme mechanism and the role of the FAD cofactor in redox activity

The crystal structure of cholesterol oxidase, a 56kDa flavoenzyme was anisotropically refined to 0.95A resolution. The final crystallographic R-factor and R(free) value is 11.0% and 13.2%, respectively. The quality of the electron density maps has enabled modeling of alternate conformations for 83 residues in the enzyme, many of which are located in the active site. The additional observed structural features were not apparent in the previous high-resolution structure (1.5A resolution) and have enabled the identification of a narrow tunnel leading directly to the isoalloxazine portion of the FAD prosthetic group. The hydrophobic nature of this narrow tunnel suggests it is the pathway for molecular oxygen to access the isoalloxazine group for the oxidative half reaction. Resolving the alternate conformations in the active site residues provides a model for the dynamics of substrate binding and a potential oxidation triggered gating mechanism involving access to the hydrophobic tunnel. This structure reveals that the NE2 atom of the active site histidine residue, H447, critical to the redox activity of this flavin oxidase, acts as a hydrogen bond donor rather than as hydrogen acceptor. The atomic resolution structure of cholesterol oxidase has revealed the presence of hydrogen atoms, dynamic aspects of the protein and how side-chain conformations are correlated with novel structural features such as the oxygen tunnel. This new structural information has provided us with the opportunity to re-analyze the roles played by specific residues in the mechanism of the enzyme.     

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.     

Ligand migration and protein fluctuations in myoglobin mutant L29W

We have determined eight X-ray structures of myoglobin mutant L29W at various experimental conditions. In addition, infrared spectroscopic experiments are presented, which are discussed in the light of the X-ray structures. Two distinct conformations of the CO-ligated protein were identified, giving rise to two stretching bands of heme-bound CO. If L29W MbCO crystals are illuminated around 180 K, a deoxy species is formed. The CO molecules migrate to the proximal side of the heme and remain trapped in the so-called Xe1 cavity upon temperature decrease to 105 K. The structure of this photoproduct is almost identical to the equilibrium high-temperature deoxy Mb structure. If the temperature is cycled to increasingly higher values, CO recombination is observed. Three intermediate structures have been determined during the rebinding process. Efficient recombination occurs only above 180 K, the characteristic temperature for the onset of protein dynamics. Rebinding is remarkably slow because bulky residues His64 and Trp29 block important migration pathways of the CO molecule.     

Quantized Water Clusters Around Apolar Molecules

Abstract

In order to understand the mechanism of gas hydrate kinetics and to explore the existance of other new cavities in the hydrate structure, we have used Molecular Dynamics (MD) simulation to study a system comprising two Lennard-Jones particles and 214 water molecules. Equilibrium structure and properties of twelve cases have been investigated. Our findings were as follows: ? Apolar molecules promote spherical liquid water clusters in a hydrate-like labile cavity.

? The size of the cavity and the coordination number is dependent upon the size of the apolar molecule.

? The coordination number of water molecules is quantized in jumps of four.

? Similarities are observed between the labile cavities and cavities in solid hydrates and in other chemical structures such as Buckminsterfullerene.

? Such a simulation procedure suggests the possibility of other clusters which may exist in yet-to-be-found hydrates. A separate question involves whether such suggested cavities can be combined with other cavities into a space-filling crystal.

     

CryoEM structure at 9A resolution of an adenovirus vector targeted to hematopoietic cells

We report a sub-nanometer resolution cryo-electron microscopy (cryoEM) structural analysis of an adenoviral vector, Ad35F, comprised of an adenovirus type 5 (Ad5) capsid pseudo-typed with an Ad35 fiber. This vector transduces human hematopoietic cells via association of its fiber protein with CD46, a member of the complement regulatory protein family. Major advances in data acquisition and image processing allowed a significant improvement in resolution compared to earlier structures. Analysis of the cryoEM density was enhanced by docking the crystal structures of both the hexon and penton base capsid proteins. CryoEM density was observed for hexon residues missing from the crystal structure that include hypervariable regions and the epitope of a neutralizing monoclonal antibody. Within the penton base, density was observed for the integrin-binding RGD loop missing from the crystal structure and for the flexible beta ribbon of the variable loop on the side of the penton base. The Ad35 fiber is flexible, consistent with the sequence insert in the third beta-spiral repeat. On the inner capsid surface density is revealed at the base of the hexons and below the penton base. A revised model is presented for protein IX within the virion. Well-defined density was assigned to a conserved domain in the N terminus of protein IX required for incorporation into the virion. For the C-terminal domain of protein IX two alternate conformations are proposed, either binding on the capsid surface or extending away from the capsid. This model is consistent with the tolerance of the C terminus for inserted ligands and its potential use in vector retargeting. This structural study increases our knowledge of Ad capsid assembly, antibody neutralization mechanisms, and may aid further improvements in gene delivery to important human cell types.     

Local structure of cytochrome c from horse heart in solution. Conformational analysis using data of two-dimensional nuclear Overhauser effect spectroscopy]

Using the earlier suggested method the calculation of the backbone conformations of horse heart cytochrome c in oxidized (ferricytochrome c) and reduced (ferrocytochrome c) states has been performed by the two-dimensional nuclear Overhauser effect spectroscopy data. For both protein forms the secondary structure elements have been revealed and the conformations of the irregular polypeptide chain segments have been analysed. The similarity of the secondary structures of ferri- and ferrocytochrome c in solution was established from the comparison of their conformations. Small differences between the conformations of two molecule forms are shown to be localized within the polypeptide chain fragments situated in the spatial structure near the heme crevice. The comparison of the dihedral phi and psi angles in the calculated conformations of horse cytochrome C with the corresponding characteristics of X-ray structures of tuna ferri- and ferrocytochrome c made for the oxidized and reduced protein forms using the quantitative criteria testifies the similarity of their conformations in solution and crystal. In is shown that the conformational changes of the separate amino acid residues which take place as the result of the "solution-to-crystal" transition occur on the surface fragments of protein globule and do not lead to essential alterations of the secondary molecule structure.     

A procedure for detection and quantitation of cavity volumes proteins. Application to measure the strength of the hydrophobic driving force in protein folding

Accurate identification of cavities is important in the study of protein structure, stability, design, and ligand binding. Identification and quantitation of cavities is a nontrivial problem because most cavities are connected to the protein exterior. We describe a computational procedure for quantitating cavity volumes and apply this to derive an estimate of the hydrophobic driving force in protein folding. A grid-based Monte Carlo procedure is used to position water molecules on the surface of a protein. A Voronoi procedure is used to identify and quantitate empty space within the solvated protein. Additional cavities not detected by other existing procedures can be identified. Most of these are close to surface concavities. Residue volumes for both the interior and the surface residues as well as cavity volumes are in good agreement with volumes calculated from fully hydrated protein structures obtained from molecular dynamic simulations. We show that the loss of stability because of cavity-creating mutations correlates better with cavity volumes determined by this procedure than with cavity volumes determined by other methods. Available structural and thermodynamic data for a number of cavity-containing mutants were analyzed to obtain estimates of 26.1 cal x mol(-1) x A(-3) and 18.5 cal x mol(-1) x A(-2) for the relative contributions of cavity formation and the hydrophobic effect to the observed stability changes. The present estimate for the hydrophobic driving force is at the lower end of estimates derived from model compound studies and considerably lower than previous estimates of approximately 50 cal x mol(-1) x A(-2) derived from protein mutational data. In the absence of structural rearrangement, on average, deletion of a single methylene group is expected to result in losses in stability of 0.41 and 0.70 kcal x mol(-1) resulting from decrease in hydrophobicity and packing, respectively.     

Crystal structure of the trigonal form of bovine beta-lactoglobulin and of its complex with retinol at 2.5 A resolution

The structure of the trigonal crystal form of bovine beta-lactoglobulin has been determined by X-ray diffraction methods. An electron density map, calculated with phases obtained by the multiple isomorphous replacement method, served as a starting point for alternate cycles of model building and restrained least-squares refinement. The model of the molecule fitted to the initial Fourier map was the one built for the orthorhombic crystal form of beta-lactoglobulin, solved at 2.8 A resolution (1 A = 0.1 nm). The final R factor for 1456 atoms (1276 non-hydrogen protein atoms and 180 solvent atoms) is 0.22, including 5245 reflections from 6.0 to 2.5 A. The molecule shows significant differences in the two crystal forms mentioned, mainly due to different packing. In the trigonal form, the species crystallized does not appear to be dimeric, but a linear polymer with tight intermolecular contacts. A difference electron density map between the complex of beta-lactoglobulin with retinol and the native protein shows no significant peaks in the cavity which, in the similar retinol-binding protein, binds the chromophore. Instead, differences are found at a surface pocket, which is limited almost completely by hydrophobic residues.     

Crystal structure of cytoglobin: the fourth globin type discovered in man displays heme hexa-coordination

Cytoglobin is a recently discovered hemeprotein belonging to the globin superfamily together with hemoglobin, myoglobin and neuroglobin. Although distributed in almost all human tissues, cytoglobin has not been ascribed a specific function. Human cytoglobin is composed of 190 amino acid residues. Sequence alignments show that a protein core region (about 150 residues) is structurally related to hemoglobin and myoglobin, being complemented by about 20 extra residues both on the N and C termini. In the absence of exogenous ligands (e.g. O2), the cytoglobin distal HisE7 residue is coordinated to the heme Fe atom, thus decreasing the ligand affinity. The crystal structure of human cytoglobin (2.1 A resolution, 21.3% R-factor) highlights a three-over-three alpha-helical globin fold, covering residues 18-171; the 1-17 N-terminal, and the 172-190 C-terminal residue segments are disordered in both molecules of the crystal asymmetric unit. Heme hexa-coordination is evident in one of the two cytoglobin chains, whereas alternate conformation for the heme distal region, achieving partial heme penta-coordination, is observed in the other. Human cytoglobin displays a large apolar protein matrix cavity, next to the heme, not related to the myoglobin cavities recognized as temporary ligand docking stations. The cavity, which may provide a heme ligand diffusion pathway, is connected to the external space through a narrow tunnel nestled between the globin G and H helices.     

The volume of cavities in proteins and virus capsids

An improved algorithm for the calculation of the volume of internal cavities within protein structures and virus capsids as well as the volumes occupied by single amino acid residues were presented. The geometrical approach was based on atomic van der Waals radii. The results obtained with two sets of the radii, those proposed by Pauling and those determined by Tsai et al were compared. The main improvement compared with our previous approach is a more elaborate treatment of the regions at the very boundary of the cavities, which yields a more accurate volume estimate. The cavity volume of a number of Plant Pathogenesis‐Related proteins of class 10 (PR‐10) were reevaluated and the volumes and other geometrical parameters for about 400 capsids of icosahedral viruses were reported. Using the same approach the volumes of amino acid residues in polypeptides as mean values averaged over multiple conformations of the side chain were also estimated. Proteins 2016; 84:1275–1286. © 2016 Wiley Periodicals, Inc.     

Stereochemistry of the interaction between methionine sulfur and the protein core.

The stereochemical features of the interaction between the sulfur atom of methionine residues and surrounding atoms are examined on a large set of known protein crystal structures. It appears that the minimum energy conformations observed in small molecule crystals are not observed within the protein core. This suggests that these interactions are either of little intensity, though they might contribute to regulate the protein physiological behavior, or physicochemically different from their counterpart in small molecule crystals.