Structural origins of pH‐dependent chemical shifts in the B1 domain of protein G

We report chemical shifts for H^N, N, and C′ nuclei in the His‐tagged B1 domain of protein G (GB1) over a range of pH values from pH 2.0 to 9.0, which fit well to standard pH‐dependent equations. We also report a 1.2 Å resolution crystal structure of GB1 at pH 3.0. Comparison of this crystal structure with published crystal structures at higher pHs provides details of the structural changes in GB1 associated with protonation of the carboxylate groups, in particular a conformational change in the C‐terminus of the protein at low pH. An additional change described recently is not seen in the crystal structure because of crystal contacts. We show that the pH‐dependent changes in chemical shifts can be almost entirely understood based on structural changes, thereby providing insight into the relationship between structure and chemical shift. In particular, we describe through‐bond effects extending up to five bonds, affecting N and C′ but not H^N; through‐space effects of carboxylates, which fit well to a simple electric field model; and effects due to conformational change, which have a similar magnitude to many of the direct effects. Finally, we discuss cooperative effects, demonstrating a lack of cooperative unfolding in the helix, and the existence of a β‐sheet “iceberg” extending over three of the four strands. This study therefore extends the application of chemical shifts to understanding protein structure. Proteins 2010; © 2010 Wiley‐Liss, Inc.     

Accurate geometries for “Mountain pass” regions of the Ramachandran plot using quantum chemical calculations

Unusual local arrangements of protein in Ramachandran space are not well represented by standard geometry tools used in either protein structure refinement using simple harmonic geometry restraints or in protein simulations using molecular mechanics force fields. In contrast, quantum chemical computations using small poly‐peptide molecular models can predict accurate geometries for any well‐defined backbone Ramachandran orientation. For conformations along transition regions—ϕ from −60 to 60°—a very good agreement with representative high‐resolution experimental X‐ray (≤1.5 Å) protein structures is obtained for both backbone C⁻¹‐N‐C_α angle and the nonbonded O⁻¹…C distance, while “standard geometry” leads to the “clashing” of O…C atoms and Amber FF99SB predicts distances too large by about 0.15 Å. These results confirm that quantum chemistry computations add valuable support for detailed analysis of local structural arrangements in proteins, providing improved or missing data for less understood high‐energy or unusual regions.     

Structures of the OmpF porin crystallized in the presence of foscholine‐12

The endogenous Escherichia coli porin OmpF was crystallized as an accidental by‐product of our efforts to express, purify, and crystallize the E. coli integral membrane protein KdpD in the presence of foscholine‐12 (FC12). FC12 is widely used in membrane protein studies, but no crystal structure of a protein that was both purified and crystallized with this detergent has been reported in the Protein Data Bank. Crystallization screening for KdpD yielded two different crystals of contaminating protein OmpF. Here, we report two OmpF structures, the first membrane protein crystal structures for which extraction, purification, and crystallization were done exclusively with FC12. The first structure was refined in space group P21 with cell parameters a = 136.7 Å, b = 210.5 Å, c = 137 Å, and β = 100.5°, and the resolution of 3.8 Å. The second structure was solved at the resolution of 4.4 Å and was refined in the P321 space group, with unit cell parameters a = 215.5 Å, b = 215.5 Å, c = 137.5 Å, and γ = 120°. Both crystal forms show novel crystal packing, in which the building block is a tetrahedral arrangement of four trimers. Additionally, we discuss the use of FC12 for membrane protein crystallization and structure determination, as well as the problem of the OmpF contamination for membrane proteins overexpressed in E. coli.     

Homology‐based hydrogen bond information improves crystallographic structures in the PDB

The Protein Data Bank (PDB) is the global archive for structural information on macromolecules, and a popular resource for researchers, teachers, and students, amassing more than one million unique users each year. Crystallographic structure models in the PDB (more than 100,000 entries) are optimized against the crystal diffraction data and geometrical restraints. This process of crystallographic refinement typically ignored hydrogen bond (H‐bond) distances as a source of information. However, H‐bond restraints can improve structures at low resolution where diffraction data are limited. To improve low‐resolution structure refinement, we present methods for deriving H‐bond information either globally from well‐refined high‐resolution structures from the PDB‐REDO databank, or specifically from on‐the‐fly constructed sets of homologous high‐resolution structures. Refinement incorporating HOmology DErived Restraints (HODER), improves geometrical quality and the fit to the diffraction data for many low‐resolution structures. To make these improvements readily available to the general public, we applied our new algorithms to all crystallographic structures in the PDB: using massively parallel computing, we constructed a new instance of the PDB‐REDO databank ( https://pdb-redo.eu ). This resource is useful for researchers to gain insight on individual structures, on specific protein families (as we demonstrate with examples), and on general features of protein structure using data mining approaches on a uniformly treated dataset.     

Evaluation at atomic resolution of the role of strain in destabilizing the temperature‐sensitive T4 lysozyme mutant Arg 96 → His

Mutant R96H is a classic temperature‐sensitive mutant of bacteriophage T4 lysozyme. It was in fact the first variant of the protein to be characterized structurally. Subsequently, it has been studied extensively by a variety of experimental and computational techniques, but the reasons for the loss of stability of the mutant protein remain controversial. In the crystallographic refinement of the mutant structure at 1.9 Å resolution one of the bond angles at the site of substitution appeared to be distorted by about 11^°, and it was suggested that this steric strain was one of the major factors in destabilizing the mutant. Different computationally‐derived models of the mutant structure, however, did not show such distortion. To determine the geometry at the site of mutation more reliably, we have extended the resolution of the data and refined the wildtype (WT) and mutant structures to be better than 1.1 Å resolution. The high‐resolution refinement of the structure of R96H does not support the bond angle distortion seen in the 1.9 Å structure determination. At the same time, it does confirm other manifestations of strain seen previously including an unusual rotameric state for His96 with distorted hydrogen bonding. The rotamer strain has been estimated as about 0.8 kcal/mol, which is about 25% of the overall reduction in stability of the mutant. Because of concern that contacts from a neighboring molecule in the crystal might influence the geometry at the site of mutation we also constructed and analyzed supplemental mutant structures in which this crystal contact was eliminated. High‐resolution refinement shows that the crystal contacts have essentially no effect on the conformation of Arg96 in WT or on His96 in the R96H mutant.     

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein‐ligand crystal structures? This value is important biophysically, for better understanding the protein‐ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh‐resolution crystal structures that can be refined without restraints.     

On interpretation of protein X‐ray structures: Planarity of the peptide unit

Pauling's mastery of peptide stereochemistry—based on small molecule crystal structures and the theory of chemical bonding—led to his realization that the peptide unit is planar and then to the Pauling–Corey–Branson model of the α‐helix. Similarly, contemporary protein structure refinement is based on experimentally determined diffraction data together with stereochemical restraints. However, even an X‐ray structure at ultra‐high resolution is still an under‐determined model in which the linkage among refinement parameters is complex. Consequently, restrictions imposed on any given parameter can affect the entire structure. Here, we examine recent studies of high resolution protein X‐ray structures, where substantial distortions of the peptide plane are found to be commonplace. Planarity is assessed by the ω‐angle, a dihedral angle determined by the peptide bond (C? N) and its flanking covalent neighbors; for an ideally planar trans peptide, ω = 180°. By using a freely available refinement package, Phenix [Afonine et al. (2012) Acta Cryst. D, 68:352–367], we demonstrate that tightening default restrictions on the ω‐angle can significantly reduce apparent deviations from peptide unit planarity without consequent reduction in reported evaluation metrics (e.g., R‐factors). To be clear, our result does not show that substantial non‐planarity is absent, only that an equivalent alternative model is possible. Resolving this disparity will ultimately require improved understanding of the deformation energy. Meanwhile, we urge inclusion of ω‐angle statistics in new structure reports in order to focus critical attention on the usual practice of assigning default values to ω‐angle constraints during structure refinement. Proteins 2015; 83:1687–1692. © 2015 Wiley Periodicals, Inc.     

On the use of structural templates for high‐resolution docking

Reliable high‐resolution prediction of protein complex structures starting from the free monomers is a considerable challenge toward large‐scale mapping of the structural details of protein‐protein interactions. The current major bottleneck is to model the conformational changes of the monomer backbone upon binding. We evaluate the use of homolog structures as source for conformational diversity, within the framework of RosettaDock—a leading high‐resolution docking protocol. We find that the use of homolog templates can improve significantly the modeling of a complex structure, including known difficult cases. Several conformational changes however are not sampled by any of the templates, indicating the need for additional sources of conformational variability. Interestingly, the successful homolog templates are not restricted to a confined range of sequence identity, highlighting the importance of the backbone conformation rather than the sequence. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Rapid increase of near atomic resolution virus capsid structures determined by cryo-electron microscopy

The recent technological advances in electron microscopes, detectors, as well as image processing and reconstruction software have brought single particle cryo-electron microscopy (cryo-EM) into prominence for determining structures of bio-molecules at near atomic resolution. This has been particularly true for virus capsids, ribosomes, and other large assemblies, which have been the ideal specimens for structural studies by cryo-EM approaches. An analysis of time series metadata of virus structures on the methods of structure determination, resolution of the structures, and size of the virus particles revealed a rapid increase in the virus structures determined by cryo-EM at near atomic resolution since 2010. In addition, the data highlight the median resolution (～3.0?Å) and size (～310.0?Å in diameter) of the virus particles determined by X-ray crystallography while no such limits exist for cryo-EM structures, which have a median diameter of 508?Å. Notably, cryo-EM virus structures in the last four years have a median resolution of 3.9?Å. Taken together with minimal sample requirements, not needing diffraction quality crystals, and being able to achieve similar resolutions of the crystal structures makes cryo-EM the method of choice for current and future virus capsid structure determinations.     

The macro domain as fusion tag for carrier‐driven crystallization

Obtaining well‐ordered crystals remains a significant challenge in protein X‐ray crystallography. Carrier‐driven crystallization can facilitate crystal formation and structure solution of difficult target proteins. We obtained crystals of the small and highly flexible SPX domain from the yeast vacuolar transporter chaperone 4 (Vtc4) when fused to a C‐terminal, non‐cleavable macro tag derived from human histone macroH2A1.1. Initial crystals diffracted to 3.3 Å resolution. Reductive protein methylation of the fusion protein yielded a new crystal form diffracting to 2.1 Å. The structures were solved by molecular replacement, using isolated macro domain structures as search models. Our findings suggest that macro domain tags can be employed in recombinant protein expression in E. coli, and in carrier‐driven crystallization.     

X‐ray vs. NMR structures as templates for computational protein design

Certain protein‐design calculations involve using an experimentally determined high‐resolution structure as a template to identify new sequences that can adopt the same fold. This approach has led to the successful design of many novel, well‐folded, native‐like proteins. Although any atomic‐resolution structure can serve as a template in such calculations, most successful designs have used high‐resolution crystal structures. Because there are many proteins for which crystal structures are not available, it is of interest whether nuclear magnetic resonance (NMR) templates are also appropriate. We have analyzed differences between using X‐ray and NMR templates in side‐chain repacking and design calculations. We assembled a database of 29 proteins for which both a high‐resolution X‐ray structure and an ensemble of NMR structures are available. Using these pairs, we compared the rotamericity, χ₁‐angle recovery, and native‐sequence recovery of X‐ray and NMR templates. We carried out design using RosettaDesign on both types of templates, and compared the energies and packing qualities of the resulting structures. Overall, the X‐ray structures were better templates for use with Rosetta. However, for ～20% of proteins, a member of the reported NMR ensemble gave rise to designs with similar properties. Re‐evaluating RosettaDesign structures with other energy functions indicated much smaller differences between the two types of templates. Ultimately, experiments are required to confirm the utility of particular X‐ray and NMR templates. But our data suggest that the lack of a high‐resolution X‐ray structure should not preclude attempts at computational design if an NMR ensemble is available. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Relaxation data in NMR structure determination: model calculations for the lysozyme-Gd3+ complex

The effect of including paramagnetic relaxation data as additional restraints in the determination of protein tertiary structures from NMR data has been explored by a systematic series of model calculations. The system used for testing the method was the 2.0 A resolution tetragonal crystal structure of hen egg white lysozyme (129 amino acid residues) and structures were generated using a version of the hybrid "distance geometry-dynamic simulated annealing" procedure. A limited set of 769 NOEs was used as restraints in all the calculations; the strengths of these were categorized into three classes on the basis of distances observed in the crystal structure. The values of 50 phi angles were also restrained on the basis of amide-alpha coupling constants calculated from the X-ray structure. Five sets of 12 structures were determined using differing sets of paramagnetic relaxation data as restraints additional to those involving the NOE and coupling constant data. The paramagnetic relaxation data were modeled on the basis of the distances of defined protons from the crystallographic binding site of Gd3+ in lysozyme. Analysis of the results showed that the relaxation data significantly improved the correspondence between the set of generated structures and the crystal structure, and that the more well defined the relaxation data, the more significant the improvement in the quality of the structures. The results suggest that the inclusion of paramagnetic relaxation restraints could be of significant value for the experimental determination of protein structures from NMR data.     

Solvent structure in crystals of trypsin determined by X-ray and neutron diffraction.

The solvent structure in orthorhombic crystals of bovine trypsin has been independently determined by X-ray diffraction to 1.35 A resolution and by neutron diffraction to 2.1 A resolution. A consensus model of the water molecule positions was obtained using oxygen positions identified in the electron density map determined by X-ray diffraction, which were verified by comparison to D2O-H2O difference neutron scattering density. Six of 184 water molecules in the X-ray structure, all with B-factors greater than 50 A2, were found to be spurious after comparison with neutron results. Roughly two-thirds of the water of hydration expected from thermodynamic data for proteins was localized by neutron diffraction; approximately one-half of the water of hydration was located by X-ray diffraction. Polar regions of the protein are well hydrated, and significant D2O-H2O difference density is seen for a small number of water molecules in a second shell of hydration. Hydrogen bond lengths and angles calculated from unconstrained refinement of water positions are distributed about values typically seen in small molecule structures. Solvent models found in seven other bovine trypsin and trypsinogen and rat trypsin structures determined by X-ray diffraction were compared. Internal water molecules are well conserved in all trypsin structures including anionic rat trypsin, which is 65% homologous to bovine trypsin. Of the 22 conserved waters in trypsin, 19 were also found in trypsinogen, suggesting that they are located in regions of the apoprotein that are structurally conserved in the transition to the mature protein. Seven waters were displaced upon activation of trypsinogen. Water structure at crystal contacts is not generally conserved in different crystal forms. Three groups of integral structural water molecules are highly conserved in all solvent structures, including a spline of water molecules inserted between two beta-strands, which may resemble an intermediate in the formation of beta sheets during the folding of a protein.     

Structural and functional analyses of catalytic domain of GH10 xylanase from Thermoanaerobacterium saccharolyticum JW/SL‐YS485

Xylanases are capable of decomposing xylans, the major components in plant cell wall, and releasing the constituent sugars for further applications. Because xylanase is widely used in various manufacturing processes, high specific activity, and thermostability are desirable. Here, the wild‐type and mutant (E146A and E251A) catalytic domain of xylanase from Thermoanaerobacterium saccharolyticum JW/SL‐YS485 (TsXylA) were expressed in Escherichia coli and purified subsequently. The recombinant protein showed optimal temperature and pH of 75°C and 6.5, respectively, and it remained fully active even after heat treatment at 75°C for 1 h. Furthermore, the crystal structures of apo‐form wild‐type TsXylA and the xylobiose‐, xylotriose‐, and xylotetraose‐bound E146A and E251A mutants were solved by X‐ray diffraction to high resolution (1.32–1.66 Å). The protein forms a classic (β/α)₈ folding of typical GH10 xylanases. The ligands in substrate‐binding groove as well as the interactions between sugars and active‐site residues were clearly elucidated by analyzing the complex structures. According to the structural analyses, TsXylA utilizes a double displacement catalytic machinery to carry out the enzymatic reactions. In conclusion, TsXylA is effective under industrially favored conditions, and our findings provide fundamental knowledge which may contribute to further enhancement of the enzyme performance through molecular engineering. Proteins 2013; 81:1256–1265. © 2013 Wiley Periodicals, Inc.     

BSP-SLIM: a blind low-resolution ligand-protein docking approach using predicted protein structures

We developed BSP‐SLIM, a new method for ligand–protein blind docking using low‐resolution protein structures. For a given sequence, protein structures are first predicted by I‐TASSER; putative ligand binding sites are transferred from holo‐template structures which are analogous to the I‐TASSER models; ligand–protein docking conformations are then constructed by shape and chemical match of ligand with the negative image of binding pockets. BSP‐SLIM was tested on 71 ligand–protein complexes from the Astex diverse set where the protein structures were predicted by I‐TASSER with an average RMSD 2.92 Å on the binding residues. Using I‐TASSER models, the median ligand RMSD of BSP‐SLIM docking is 3.99 Å which is 5.94 Å lower than that by AutoDock; the median binding‐site error by BSP‐SLIM is 1.77 Å which is 6.23 Å lower than that by AutoDock and 3.43 Å lower than that by LIGSITE^CSC. Compared to the models using crystal protein structures, the median ligand RMSD by BSP‐SLIM using I‐TASSER models increases by 0.87 Å, while that by AutoDock increases by 8.41 Å; the median binding‐site error by BSP‐SLIM increase by 0.69Å while that by AutoDock and LIGSITE^CSC increases by 7.31 Å and 1.41 Å, respectively. As case studies, BSP‐SLIM was used in virtual screening for six target proteins, which prioritized actives of 25% and 50% in the top 9.2% and 17% of the library on average, respectively. These results demonstrate the usefulness of the template‐based coarse‐grained algorithms in the low‐resolution ligand–protein docking and drug‐screening. An on‐line BSP‐SLIM server is freely available at http://zhanglab.ccmb.med.umich.edu/BSP‐SLIM . Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Crystal structure of Bacillus subtilis SPP1 phage gp23.1, a putative chaperone

SPP1 is a siphophage infecting the gram‐positive bacterium Bacillus subtilis. The SPP1 tail electron microscopy (EM) reconstruction revealed that it is mainly constituted by conserved structural proteins such as the major tail proteins (gp17.1), the tape measure protein (gp18), the Distal tail protein (Dit, gp19.1), and the Tail associated lysin (gp21). A group of five small genes (22–24.1) follows in the genome but it remains to be elucidated whether their protein products belong or not to the tail. Noteworthy, an unassigned EM density accounting for ～245 kDa is present at the distal end of the SPP1 tail‐tip. We report here the gp23.1 crystal structure at 1.6 Å resolution, a protein that lacks sequence identity to any known protein. We found that gp23.1 forms a hexamer both in the crystal lattice and in solution as revealed by light scattering measurements. The gp23.1 hexamer does not fit well in the unassigned SPP1 tail‐tip EM density and we hypothesize that this protein might act as a chaperone.     

Open and shut: Crystal structures of the dodecylmaltoside solubilized mechanosensitive channel of small conductance from Escherichia coli and Helicobacter pylori at 4.4 Å and 4.1 Å resolutions

The mechanosensitive channel of small conductance (MscS) contributes to the survival of bacteria during osmotic downshock by transiently opening large diameter pores for the efflux of cellular contents before the membrane ruptures. Two crystal structures of the Escherichia coli MscS are currently available, the wild type protein in a nonconducting state at 3.7 Å resolution (Bass et al., Science 2002; 298:1582–1587) and the Ala106Val variant in an open state at 3.45 Å resolution (Wang et al., Science 2008; 321:1179–1183). Both structures used protein solubilized in the detergent fos‐choline‐14. We report here crystal structures of MscS from E. coli and Helicobacter pylori solubilized in the detergent β‐dodecylmaltoside at resolutions of 4.4 and 4.2 Å, respectively. While the cytoplasmic domains are unchanged in these structures, distinct conformations of the transmembrane domains are observed. Intriguingly, β‐dodecylmaltoside solubilized wild type E. coli MscS adopts the open state structure of A106V E. coli MscS, while H. pylori MscS resembles the nonconducting state structure observed for fos‐choline‐14 solubilized E. coli MscS. These results highlight the sensitivity of membrane protein conformational equilibria to variations in detergent, crystallization conditions, and protein sequence.     

Water polygons in high‐resolution protein crystal structures

We have analyzed the interstitial water (ISW) structures in 1500 protein crystal structures deposited in the Protein Data Bank that have greater than 1.5 Å resolution with less than 90% sequence similarity with each other. We observed varieties of polygonal water structures composed of three to eight water molecules. These polygons may represent the time‐ and space‐averaged structures of “stable” water oligomers present in liquid water, and their presence as well as relative population may be relevant in understanding physical properties of liquid water at a given temperature. On an average, 13% of ISWs are localized enough to be visible by X‐ray diffraction. Of those, averages of 78% are water molecules in the first water layer on the protein surface. Of the localized ISWs beyond the first layer, almost half of them form water polygons such as trigons, tetragons, as well as expected pentagons, hexagons, higher polygons, partial dodecahedrons, and disordered networks. Most of the octagons and nanogons are formed by fusion of smaller polygons. The trigons are most commonly observed. We suggest that our observation provides an experimental basis for including these water polygon structures in correlating and predicting various water properties in liquid state.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.