The ultrahigh resolution crystal structure of ribonuclease A containing an isoaspartyl residue: hydration and sterochemical analysis

Revisiting macromolecular hydration with HullRadSAS

Hydration of biological macromolecules is important for their stability and function. Historically, attempts have been made to describe the degree of macromolecular hydration using a single parameter over a narrow range of values. Here, we describe a method to calculate two types of hydration: surface shell water and entrained water. A consideration of these two types of hydration helps to explain the "hydration problem" in hydrodynamics. The combination of these two types of hydration allows accurate calculation of hydrodynamic volume and related macromolecular properties such as sedimentation and diffusion coefficients, intrinsic viscosities, and the concentration-dependent non-ideality identified with sedimentation velocity experiments.

Deamidation in Moxetumomab Pasudotox Leading to Conformational Change and Immunotoxin Activity Loss

Asparagine (Asn) deamidation is a common posttranslational modification in which Asn is converted to aspartic acid or isoaspartic acid. By introducing a negative charge, deamidation could potentially impact the binding interface and biological activities of protein therapeutics. We identified a deamidation variant in moxetumomab pasudotox, an immunotoxin Fv fusion protein drug derived from a 38-kDa truncated Pseudomonas exotoxin A (PE38) for the treatment of hairy-cell leukemia. Although the deamidation site, Asn-358, was outside of the binding interface, the modification had a significant impact on the biological activity of moxetumomab pasudotox. Surprisingly, the variant eluted earlier than its unmodified form on anion exchange chromatography, which often leads to the conclusion that it has a higher positive charge. Here we describe the characterization of the deamidation variant with differential scanning calorimetry and hydrogen-deuterium exchange mass spectrometry, which revealed that the Asn-358 deamidation caused the conformational changes in the catalytic domain of the PE38 region. These results provide an explanation for why the deamidation affected the biological activity of moxetumomab pasudotox and suggest the approach that can be used for process control to ensure product quality and process consistency.

Structural and biochemical basis of the formation of isoaspartate in the complementarity-determining region of antibody 64M-5 Fab

The formation of the isoaspartate (isoAsp) is one of spontaneous degradation processes of proteins, affecting their stability and activity. Here, we report for the first time the crystal structures of an antibody Fab that contains isoAsp in the complementarity-determining region (CDR), along with biochemical studies to detect isoAsp. By comparing the elution profiles of cation-exchange chromatography, it was clarified that the antibody 64M-5 Fab is converted from the normal form to isoAsp form spontaneously and time-dependently under physiological conditions. The isoAsp residue was identified with tryptic peptide mapping, N-terminal sequencing, and the protein isoaspartyl methyltransferase assay. Based on the fluorescence quenching method, the isoAsp form of 64M-5 Fab shows a one order of magnitude lower binding constant for its dinucleotide ligand dT(6-4)T than the normal form. According to the structure of the isoAsp form, the conformation of CDR L1 is changed from the normal form to isoAsp form; the loss of hydrogen bonds involving the Asn28L side-chain, and structural conversion of the β-turn from type I to type II'. The formation of isoAsp leads to a large displacement of the side chain of His27dL, and decreased electrostatic interactions with the phosphate group of dT(6-4)T. Such structural changes should be responsible for the lower affinity of the isoAsp form for dT(6-4)T than the normal form. These findings may provide insight into neurodegenerative diseases (NDDs) and related diseases caused by misfolded proteins.

Deamidation of Amino Acids on the Surface of Adeno-Associated Virus Capsids Leads to Charge Heterogeneity and Altered Vector Function

Post-translational modification of the adeno-associated virus capsids is a poorly understood factor in the development of these viral vectors into pharmaceutical products. Here we report the extensive capsid deamidation of adeno-associated virus serotype 8 and seven other diverse adeno-associated virus serotypes, with supporting evidence from structural, biochemical, and mass spectrometry approaches. The extent of deamidation at each site depended on the vector's age and multiple primary-sequence and three-dimensional structural factors. However, the extent of deamidation was largely independent of the vector recovery and purification conditions. We demonstrate the potential for deamidation to impact transduction activity and, moreover, correlate an early time point loss in vector activity to rapidly progressing spontaneous deamidation at several adeno-associated virus 8 asparagines. We explore mutational strategies that stabilize side-chain amides, improving vector transduction and reducing the lot-to-lot molecular variability that presents a key concern in biologics manufacturing. This study illuminates a previously unknown aspect of adeno-associated virus capsid heterogeneity and highlights its importance in the development of these vectors for gene therapy.

Extensive deamidation of RNase A inhibits its oligomerization through 3D domain swapping

Bovine pancreatic ribonuclease A (RNase A) is the monomeric prototype of the so-called secretory 'pancreatic-type' RNase super-family. Like the naturally domain-swapped dimeric bovine seminal variant, BS-RNase, and its glycosylated RNase B isoform, RNase A forms N- and C-terminal 3D domain-swapped oligomers after lyophilization from acid solutions, or if subjected to thermal denaturation at high protein concentration. All mentioned RNases can undergo deamidation at Asn67, forming Asp or isoAsp derivatives that modify the protein net charge and consequently its enzymatic activity. In addition, deamidation slightly affects RNase B self-association through the 3D domain swapping (3D-DS) mechanism. We report here the influence of extensive deamidation on RNase A tendency to oligomerize through 3D-DS. In particular, deamidation of Asn67 alone slightly decreases the propensity of the protein to oligomerize, with the Asp derivative being less affected than the isoAsp one. Contrarily, the additional Asp and/or isoAsp conversion of residues other than N67 almost nullifies RNase A oligomerization capability. In addition, Gln deamidation, although less kinetically favorable, may affect RNase A self-association. Using 2D and 3D NMR we identified the Asn/Gln residues most prone to undergo deamidation. Together with CD spectroscopy, NMR also indicates that poly-deamidated RNase A generally maintains its native tertiary structure. Again, we investigated in silico the effect of the residues undergoing deamidation on RNase A dimers structures. Finally, the effect of deamidation on RNase A oligomerization is discussed in comparison with studies on deamidation-prone proteins involved in amyloid formation.

Bond distances in polypeptide backbones depend on the local conformation

By combining quantum-mechanical analysis of small model peptides and statistical surveys of high-resolution protein structures, a systematic conformational dependence of bond lengths in polypeptide backbones has been unveiled which involves both the peptide bond (C—O and C—N) and those bonds centred on the Cα atom. All of these bond lengths indeed display a systematic variability in the ψ angle according to both calculations and surveys of protein structures. The overall agreement between the computed and the statistical data suggests that these trends are essentially driven by local effects. The dependence of Cα distances on ψ is governed by interactions between the σ system of the Cα moiety and the C—O π system of the peptide bond. Maximum and minimum values for each bond distance are found for conformations with the specific bond perpendicular and parallel to the adjacent CONH peptide plane, respectively. On the other hand, the variability of the C—O and C—N distances is related to the strength of the interactions between the lone pair of the N atom and the C—O π* system, which is modulated by the ψ angle. The C—O and C—N distances are related but their trends are not strictly connected to peptide-bond planarity, although a correlation amongst all of these parameters is expected on the basis of the classical resonance model.

Structural effects of protein aging: terminal marking by deamidation in human triosephosphate isomerase

Deamidation, the loss of the ammonium group of asparagine and glutamine to form aspartic and glutamic acid, is one of the most commonly occurring post-translational modifications in proteins. Since deamidation rates are encoded in the protein structure, it has been proposed that they can serve as molecular clocks for the timing of biological processes such as protein turnover, development and aging. Despite the importance of this process, there is a lack of detailed structural information explaining the effects of deamidation on the structure of proteins. Here, we studied the effects of deamidation on human triosephosphate isomerase (HsTIM), an enzyme for which deamidation of N15 and N71 has been long recognized as the signal for terminal marking of the protein. Deamidation was mimicked by site directed mutagenesis; thus, three mutants of HsTIM (N15D, N71D and N15D/N71D) were characterized. The results show that the N71D mutant resembles, structurally and functionally, the wild type enzyme. In contrast, the N15D mutant displays all the detrimental effects related to deamidation. The N15D/N71D mutant shows only minor additional effects when compared with the N15D mutation, supporting that deamidation of N71 induces negligible effects. The crystal structures show that, in contrast to the N71D mutant, where minimal alterations are observed, the N15D mutation forms new interactions that perturb the structure of loop 1 and loop 3, both critical components of the catalytic site and the interface of HsTIM. Based on a phylogenetic analysis of TIM sequences, we propose the conservation of this mechanism for mammalian TIMs.

Assignments of RNase A by ADAPT-NMR and enhancer

We report here backbone (1)H and (15)N assignments for ribonuclease A obtained by using ADAPT-NMR, a fully-automated approach for combined data collection, spectral analysis and resonance assignment. ADAPT-NMR was able to assign 98% of the resonances with 93% agreement with traditional data collection and assignment. Further refinement of the automated results with ADAPT-NMR enhancer led to complete (100%) assignments with 96% agreement with assignments by the traditional approach.

Interplay between peptide bond geometrical parameters in nonglobular structural contexts

Several investigations performed in the last two decades have unveiled that geometrical parameters of protein backbone show a remarkable variability. Although these studies have provided interesting insights into one of the basic aspects of protein structure, they have been conducted on globular and water-soluble proteins. We report here a detailed analysis of backbone geometrical parameters in nonglobular proteins/peptides. We considered membrane proteins and two distinct fibrous systems (amyloid-forming and collagen-like peptides). Present data show that in these systems the local conformation plays a major role in dictating the amplitude of the bond angle N-C(α)-C and the propensity of the peptide bond to adopt planar/nonplanar states. Since the trends detected here are in line with the concept of the mutual influence of local geometry and conformation previously established for globular and water-soluble proteins, our analysis demonstrates that the interplay of backbone geometrical parameters is an intrinsic and general property of protein/peptide structures that is preserved also in nonglobular contexts. For amyloid-forming peptides significant distortions of the N-C(α)-C bond angle, indicative of sterical hidden strain, may occur in correspondence with side chain interdigitation. The correlation between the dihedral angles Δω/ψ in collagen-like models may have interesting implications for triple helix stability.

On the reproducibility of protein crystal structures: five atomic resolution structures of trypsin

Structural studies of proteins usually rely on a model obtained from one crystal. By investigating the details of this model, crystallographers seek to obtain insight into the function of the macromolecule. It is therefore important to know which details of a protein structure are reproducible or to what extent they might differ. To address this question, the high-resolution structures of five crystals of bovine trypsin obtained under analogous conditions were compared. Global parameters and structural details were investigated. All of the models were of similar quality and the pairwise merged intensities had large correlation coefficients. The C(α) and backbone atoms of the structures superposed very well. The occupancy of ligands in regions of low thermal motion was reproducible, whereas solvent molecules containing heavier atoms (such as sulfur) or those located on the surface could differ significantly. The coordination lengths of the calcium ion were conserved. A large proportion of the multiple conformations refined to similar occupancies and the residues adopted similar orientations. More than three quarters of the water-molecule sites were conserved within 0.5 Å and more than one third were conserved within 0.1 Å. An investigation of the protonation states of histidine residues and carboxylate moieties was consistent for all of the models. Radiation-damage effects to disulfide bridges were observed for the same residues and to similar extents. Main-chain bond lengths and angles averaged to similar values and were in agreement with the Engh and Huber targets. Other features, such as peptide flips and the double conformation of the inhibitor molecule, were also reproducible in all of the trypsin structures. Therefore, many details are similar in models obtained from different crystals. However, several features of residues or ligands located in flexible parts of the macromolecule may vary significantly, such as side-chain orientations and the occupancies of certain fragments.

A novel interdomain interface in crystallins: structural characterization of the βγ-crystallin from Geodia cydonium at 0.99 Å resolution

The βγ-crystallin superfamily includes highly diverse proteins belonging to all of the kingdoms of life. Based on structural topology, these proteins are considered to be evolutionarily related to the long-lived βγ-crystallins that constitute the vertebrate eye lens. This study reports the crystallographic structure at 0.99 Å resolution of the two-domain βγ-crystallin (geodin) from the sponge Geodia cydonium. This is the most ancient member of the βγ-crystallin superfamily in metazoans. The X-ray structure shows that the geodin domains adopt the typical βγ-crystallin fold with a paired Greek-key motif, thus confirming the hypothesis that the crystallin-type scaffold used in the evolution of bacteria and moulds was recruited very early in metazoans. As a significant new structural feature, the sponge protein possesses a unique interdomain interface made up by pairing between the second motif of the first domain and the first motif of the second domain. The atomic resolution also allowed a detailed analysis of the calcium-binding site of the protein.

Carbon-oxygen hydrogen bonding in biological structure and function

Carbon-oxygen (CH···O) hydrogen bonding represents an unusual category of molecular interactions first documented in biological structures over 4 decades ago. Although CH···O hydrogen bonding has remained generally underappreciated in the biochemical literature, studies over the last 15 years have begun to yield direct evidence of these interactions in biological systems. In this minireview, we provide a historical context of biological CH···O hydrogen bonding and summarize some major advancements from experimental studies over the past several years that have elucidated the importance, prevalence, and functions of these interactions. In particular, we examine the impact of CH···O bonds on protein and nucleic acid structure, molecular recognition, and enzyme catalysis and conclude by exploring overarching themes and unresolved questions regarding unconventional interactions in biomolecular structure.

Significant reduction in errors associated with nonbonded contacts in protein crystal structures: automated all-atom refinement with PrimeX

All-atom models are essential for many applications in molecular modeling and computational chemistry. Nonbonded atomic contacts much closer than the sum of the van der Waals radii of the two atoms (clashes) are commonly observed in such models derived from protein crystal structures. A set of 94 recently deposited protein structures in the resolution range 1.5-2.8 Å were analyzed for clashes by the addition of all H atoms to the models followed by optimization and energy minimization of the positions of just these H atoms. The results were compared with the same set of structures after automated all-atom refinement with PrimeX and with nonbonded contacts in protein crystal structures at a resolution equal to or better than 0.9 Å. The additional PrimeX refinement produced structures with reasonable summary geometric statistics and similar R(free) values to the original structures. The frequency of clashes at less than 0.8 times the sum of van der Waals radii was reduced over fourfold compared with that found in the original structures, to a level approaching that found in the ultrahigh-resolution structures. Moreover, severe clashes at less than or equal to 0.7 times the sum of atomic radii were reduced 15-fold. All-atom refinement with PrimeX produced improved crystal structure models with respect to nonbonded contacts and yielded changes in structural details that dramatically impacted on the interpretation of some protein-ligand interactions.

Chain termini cross-talk in the swapping process of bovine pancreatic ribonuclease

3D domain swapping is the process by which two or more protein molecules exchange part of their structure to form intertwined dimers or higher oligomers. Bovine pancreatic ribonuclease (RNase A) is able to swap the N-terminal α-helix (residues 1-13) and/or the C-terminal β-strand (residues 116-124), thus forming a variety of oligomers, including two different dimers. Cis-trans isomerization of the Asn113-Pro114 peptide group was observed when the protein formed the C-terminal swapped dimer. To study the effect of the substitution of Pro114 on the swapping process of RNase A, we have prepared and characterized the P114A monomeric and dimeric variants of the enzyme. In contrast with previous reports, the crystal structure and NMR data on the monomer reveals a mixed cis-trans conformation for the Asn113-Ala114 peptide group, whereas the X-ray structure of the C-terminal swapped dimer of the variant is very close to that of the corresponding dimer of RNase A. The mutation at the C-terminus affects the capability of the N-terminal α-helix to swap and the stability of both dimeric forms. The present results underscore the importance of the hydration shell in determining the cross-talk between the chain termini in the swapping process of RNase A.

Automatic α-helix identification in Patterson maps

α-Helices are peculiar atomic arrangements characterizing protein structures. Their occurrence can be used within crystallographic methods as minimal a priori information to drive the phasing process towards solution. Recently, brute-force methods have been developed which search for all possible positions of α-helices in the crystal cell by molecular replacement and explore all of them systematically. Knowing the α-helix orientations in advance would be a great advantage for this kind of approach. For this purpose, a fully automatic procedure to find α-helix orientations within the Patterson map has been developed. The method is based on Fourier techniques specifically addressed to the identification of helical shapes and operating on Patterson maps described in spherical coordinates. It supplies a list of candidate orientations, which are then refined by using a figure of merit based on a rotation function calculated for a template polyalanine helix oriented along the current direction. The orientation search algorithm has been optimized to work at 3 Å resolution, while the candidates are refined against all measured reflections. The procedure has been applied to a large number of protein test structures, showing an overall efficiency of 77% in finding α-helix orientations, which decreases to 48% on limiting the number of candidate solutions (to 13 on average). The information obtained may be used in many aspects in the framework of molecular-replacement phasing, as well as to constrain the generation of models in computational modelling programs. The procedure will be accessible through the next release of IL MILIONE and could be decisive in the solution of new unknown structures.

Peptide bond distortions from planarity: new insights from quantum mechanical calculations and peptide/protein crystal structures

By combining quantum-mechanical analysis and statistical survey of peptide/protein structure databases we here report a thorough investigation of the conformational dependence of the geometry of peptide bond, the basic element of protein structures. Different peptide model systems have been studied by an integrated quantum mechanical approach, employing DFT, MP2 and CCSD(T) calculations, both in aqueous solution and in the gas phase. Also in absence of inter-residue interactions, small distortions from the planarity are more a rule than an exception, and they are mainly determined by the backbone ψ dihedral angle. These indications are fully corroborated by a statistical survey of accurate protein/peptide structures. Orbital analysis shows that orbital interactions between the σ system of C(α) substituents and the π system of the amide bond are crucial for the modulation of peptide bond distortions. Our study thus indicates that, although long-range inter-molecular interactions can obviously affect the peptide planarity, their influence is statistically averaged. Therefore, the variability of peptide bond geometry in proteins is remarkably reproduced by extremely simplified systems since local factors are the main driving force of these observed trends. The implications of the present findings for protein structure determination, validation and prediction are also discussed.

On the complexity of Engh and Huber refinement restraints: the angle τ as example

The angle τ (backbone N—C^α—C) is the most contested Engh and Huber refinement target parameter. It is shown that this parameter is ‘correct’ as a PDB-wide average, but can be improved by taking into account residue types, secondary structures and many other aspects of our knowledge of the biophysical relations between residue type and protein structure.

The Engh and Huber parameters for bond lengths and bond angles have been used uncontested in macromolecular structure refinement from 1991 until very recently, despite critical discussion of their ubiquitous validity by many authors. An extensive analysis of the backbone angle τ (N—C^α—C) illustrates that the Engh and Huber parameters can indeed be improved and a recent study [Tronrud et al. (2010 ▶), Acta Cryst. D66, 834–842] confirms these ideas. However, the present study of τ shows that improving the Engh and Huber parameters will be considerably more complex than simply making the parameters a function of the backbone ϕ, ψ angles. Many other aspects, such as the cooperativity of hydrogen bonds, the bending of secondary-structure elements and a series of biophysical aspects of the 20 amino-acid types, will also need to be taken into account. Different sets of Engh and Huber parameters will be needed for conceptually different refinement programs.

Structural changes induced by the deamidation and isomerization of asparagine revealed by the crystal structure of Ustilago sphaerogena ribonuclease U2B

Under physiological conditions, the deamidation and isomerization of asparagine to isoaspartate (isoAsp) proceeds nonenzymatically via succinimide. Although a large number of proteins have been reported to contain isoAsp, information concerning the three-dimensional structure of proteins containing isoaspartate is still limited. We have crystallized isoAsp containing Ustilago sphaerogena ribonuclease U2B, and determined the crystal structure at 1.32 Å resolution. The structure revealed that the formation of isoAsp32 induces a single turn unfolding of the α-helix from Asp29 to Asp34, and the region from Asp29 to Arg35 forms a U-shaped loop structure. The electron density map shows that isoAsp32 retained the L-configuration at the C(α) atom. IsoAsp32 is in gauche conformation about a C(α)--C(β) bond, and the polypeptide chain bends by ∼90° at isoAsp32. IsoAsp32 protrudes from the surface of the protein, and the abnormal β-peptide bond in the main-chain and α-carboxylate in the side-chain is fully exposed. The structure suggests that the deamidation of the Asn and the isoAsp formation in proteins could confer immunogenicity.

Variation in quadrupole couplings of alpha deuterons in ubiquitin suggests the presence of C(alpha)-H(alpha)...O=C hydrogen bonds

Nuclear quadrupolar couplings are sensitive probes of hydrogen bonding. Experimental quadrupolar coupling constants of alpha deuterons (D(alpha) QCC) are reported for the residues of human ubiquitin that do not experience large-amplitude internal dynamics on the pico- to nanosecond time scale. Two different methods for D(alpha) QCC estimation are employed: (i) direct estimation of D(alpha) QCC values from R(1) and R(2) (2)H D(alpha) rates using the dynamics parameters (S(C(alpha)-H(alpha))(2)) derived from 1 micros molecular dynamics simulations as well as from (13)C(alpha) relaxation measurements and (ii) indirect measurements via scalar relaxation of the second kind that affects (13)C(alpha) relaxation rates in (13)C(alpha)-D(alpha) spin systems. A relatively large variability of D(alpha) QCC values is produced by both methods. The average value of 170.6 +/- 3 kHz is derived from the combined data set, with D(alpha) QCC values ranging from 159.2 to 177.2 kHz. The set of lowest quadrupolar couplings in all data sets corresponds to the residues that are likely to form weak C(alpha)-H(alpha)...O=C hydrogen bonds as predicted from the analysis of short H(alpha)...O distances in three-dimensional structures of ubiquitin. These D(alpha) nuclei show up to 10 kHz reduction in their QCC values, which is in agreement with earlier solid-state NMR measurements in alpha deuterons of glycine. A statistically significant correlation is observed between the QCC values of alpha-deuterons and the inverse cube of C(alpha)-H(alpha)...O=C distances in ubiquitin.

X-ray crystallographic studies of RNase A variants engineered at the most destabilizing positions of the main hydrophobic core: further insight into protein stability

To investigate the structural origin of decreased pressure and temperature stability, the crystal structure of bovine pancreatic ribonuclease A variants V47A, V54A, V57A, I81A, I106A, and V108A was solved at 1.4-2.0 A resolution and compared with the structure of wild-type protein. The introduced mutations had only minor influence on the global structure of ribonuclease A. The structural changes had individual character that depends on the localization of mutated residue, however, they seemed to expand from mutation site to the rest of the structure. Several different parameters have been evaluated to find correlation with decrease of free energy of unfolding DeltaDeltaG(T), and the most significant correlation was found for main cavity volume change. Analysis of the difference distance matrices revealed that the ribonuclease A molecule is organized into five relatively rigid subdomains with individual response to mutation. This behavior consistent with results of unfolding experiments is an intrinsic feature of ribonuclease A that might be surviving remnants of folding intermediates and reflects the dynamic nature of the molecule.

Protein repair L-isoaspartyl methyltransferase 1 is involved in both seed longevity and germination vigor in Arabidopsis

The formation of abnormal amino acid residues is a major source of spontaneous age-related protein damage in cells. The protein l-isoaspartyl methyltransferase (PIMT) combats protein misfolding resulting from l-isoaspartyl formation by catalyzing the conversion of abnormal l-isoaspartyl residues to their normal l-aspartyl forms. In this way, the PIMT repair enzyme system contributes to longevity and survival in bacterial and animal kingdoms. Despite the discovery of PIMT activity in plants two decades ago, the role of this enzyme during plant stress adaptation and in seed longevity remains undefined. In this work, we have isolated Arabidopsis thaliana lines exhibiting altered expression of PIMT1, one of the two genes encoding the PIMT enzyme in Arabidopsis. PIMT1 overaccumulation reduced the accumulation of l-isoaspartyl residues in seed proteins and increased both seed longevity and germination vigor. Conversely, reduced PIMT1 accumulation was associated with an increase in the accumulation of l-isoaspartyl residues in the proteome of freshly harvested dry mature seeds, thus leading to heightened sensitivity to aging treatments and loss of seed vigor under stressful germination conditions. These data implicate PIMT1 as a major endogenous factor that limits abnormal l-isoaspartyl accumulation in seed proteins, thereby improving seed traits such as longevity and vigor. The PIMT repair pathway likely works in concert with other anti-aging pathways to actively eliminate deleterious protein products, thus enabling successful seedling establishment and strengthening plant proliferation in natural environments.

Structural basis for distinctive recognition of fibrinogen gammaC peptide by the platelet integrin alphaIIbbeta3

Hemostasis and thrombosis (blood clotting) involve fibrinogen binding to integrin α_IIbβ₃ on platelets, resulting in platelet aggregation. α_vβ₃ binds fibrinogen via an Arg-Asp-Gly (RGD) motif in fibrinogen's α subunit. α_IIbβ₃ also binds to fibrinogen; however, it does so via an unstructured RGD-lacking C-terminal region of the γ subunit (γC peptide). These distinct modes of fibrinogen binding enable α_IIbβ₃ and α_vβ₃ to function cooperatively in hemostasis. In this study, crystal structures reveal the integrin α_IIbβ₃–γC peptide interface, and, for comparison, integrin α_IIbβ₃ bound to a lamprey γC primordial RGD motif. Compared with RGD, the GAKQAGDV motif in γC adopts a different backbone configuration and binds over a more extended region. The integrin metal ion–dependent adhesion site (MIDAS) Mg²⁺ ion binds the γC Asp side chain. The adjacent to MIDAS (ADMIDAS) Ca²⁺ ion binds the γC C terminus, revealing a contribution for ADMIDAS in ligand binding. Structural data from this natively disordered γC peptide enhances our understanding of the involvement of γC peptide and integrin α_IIbβ₃ in hemostasis and thrombosis.

Algorithm for backrub motions in protein design

Motivation: The Backrub is a small but kinematically efficient side-chain-coupled local backbone motion frequently observed in atomic-resolution crystal structures of proteins. A backrub shifts the C_α–C_β orientation of a given side-chain by rigid-body dipeptide rotation plus smaller individual rotations of the two peptides, with virtually no change in the rest of the protein. Backrubs can therefore provide a biophysically realistic model of local backbone flexibility for structure-based protein design. Previously, however, backrub motions were applied via manual interactive model-building, so their incorporation into a protein design algorithm (a simultaneous search over mutation and backbone/side-chain conformation space) was infeasible.

Results: We present a combinatorial search algorithm for protein design that incorporates an automated procedure for local backbone flexibility via backrub motions. We further derive a dead-end elimination (DEE)-based criterion for pruning candidate rotamers that, in contrast to previous DEE algorithms, is provably accurate with backrub motions. Our backrub-based algorithm successfully predicts alternate side-chain conformations from ≤0.9 Å resolution structures, confirming the suitability of the automated backrub procedure. Finally, the application of our algorithm to redesign two different proteins is shown to identify a large number of lower-energy conformations and mutation sequences that would have been ignored by a rigid-backbone model.

Availability: Contact authors for source code.

Contact:brd+ismb08@cs.duke.edu

Protein structural variation in computational models and crystallographic data

Normal mode analysis offers an efficient way of modeling the conformational flexibility of protein structures. We use anisotropic displacement parameters from crystallography to test the quality of prediction of both the magnitude and directionality of conformational flexibility. Normal modes from four simple elastic network model potentials and from the CHARMM force field are calculated for a data set of 83 diverse, ultrahigh-resolution crystal structures. While all five potentials provide good predictions of the magnitude of flexibility, all-atom potentials have a clear edge at prediction of directionality, and the CHARMM potential has the highest prediction quality. The low-frequency modes from different potentials are similar, but those computed from the CHARMM potential show the greatest difference from the elastic network models. The comprehensive evaluation demonstrates the costs and benefits of using normal mode potentials of varying complexity.

Statistical and conformational analysis of the electron density of protein side chains

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.

The backrub motion: how protein backbone shrugs when a sidechain dances

Surprisingly, the frozen structures from ultra-high-resolution protein crystallography reveal a prevalent, but subtle, mode of local backbone motion coupled to much larger, two-state changes of sidechain conformation. This "backrub" motion provides an influential and common type of local plasticity in protein backbone. Concerted reorientation of two adjacent peptides swings the central sidechain perpendicular to the chain direction, changing accessible sidechain conformations while leaving flanking structure undisturbed. Alternate conformations in sub-1 angstroms crystal structures show backrub motions for two-thirds of the significant Cbeta shifts and 3% of the total residues in these proteins (126/3882), accompanied by two-state changes in sidechain rotamer. The Backrub modeling tool is effective in crystallographic rebuilding. For homology modeling or protein redesign, backrubs can provide realistic, small perturbations to rigid backbones. For large sidechain changes in protein dynamics or for single mutations, backrubs allow backbone accommodation while maintaining H bonds and ideal geometry.

Optimization of the GB/SA solvation model for predicting the structure of surface loops in proteins

Implicit solvation models are commonly optimized with respect to experimental data or Poisson-Boltzmann (PB) results obtained for small molecules, where the force field is sometimes not considered. In previous studies, we have developed an optimization procedure for cyclic peptides and surface loops in proteins based on the entire system studied and the specific force field used. Thus, the loop has been modeled by the simplified solvation function E(tot) = E(FF) (epsilon = 2r) + Sigma(i) sigma(i)A(i), where E(FF) (epsilon = nr) is the AMBER force field energy with a distance-dependent dielectric function, epsilon = nr, A(i) is the solvent accessible surface area of atom i, and sigma(i) is its atomic solvation parameter. During the optimization process, the loop is free to move while the protein template is held fixed in its X-ray structure. To improve on the results of this model, in the present work we apply our optimization procedure to the physically more rigorous solvation model, the generalized Born with surface area (GB/SA) (together with the all-atom AMBER force field) as suggested by Still and co-workers (J. Phys. Chem. A 1997, 101, 3005). The six parameters of the GB/SA model, namely, P(1)-P(5) and the surface area parameter, sigma (programmed in the TINKER package) are reoptimized for a "training" group of nine loops, and a best-fit set is defined from the individual sets of optimized parameters. The best-fit set and Still's original set of parameters (where Lys, Arg, His, Glu, and Asp are charged or neutralized) were applied to the training group as well as to a "test" group of seven loops, and the energy gaps and the corresponding RMSD values were calculated. These GB/SA results based on the three sets of parameters have been found to be comparable; surprisingly, however, they are somewhat inferior (e.g, of larger energy gaps) to those obtained previously from the simplified model described above. We discuss recent results for loops obtained by other solvation models and potential directions for future studies.

Structure of murine angiogenin: features of the substrate- and cell-binding regions and prospects for inhibitor-binding studies

Angiogenin is an unusual member of the pancreatic ribonuclease superfamily that induces blood-vessel formation and is a promising anticancer target. The three-dimensional structure of murine angiogenin (mAng) has been determined by X-ray crystallography. Two structures are presented: one is a complex with sulfate ions (1.5 Angstroms resolution) and the other a complex with phosphate ions (1.6 Angstroms resolution). Residues forming the putative B(1), P(1) and B(2) subsites occupy positions similar to their hAng counterparts and are likely to play similar roles. The anions occupy the P(1) subsite, sulfate binding conventionally and phosphate adopting two orientations, one of which is novel. The B(1) subsite is obstructed by Glu116 and Phe119, with the latter assuming a less invasive position than its hAng counterpart. Hydrophobic interactions between the C-terminal segment and the main body of the protein are more extensive than in hAng and may underly the lower enzymatic activity of the murine protein. Elsewhere, the structure of the H3-B2 loop supports the view that hAng Asn61 interacts directly with cell-surface molecules and does not merely stabilize adjacent regions of the hAng structure. mAng crystals appear to offer small-molecule inhibitors a clear route to the active site and may even withstand a reorientation of the C-terminal segment that provides access to the cryptic B(1) subsite. These features represent considerable advantages over crystalline hAng and bAng.

How thick is the layer of thermal volume surrounding the protein?

Investigation on the volume properties of protein hydration layers is reported. Presented results are based on combination of Monte Carlo modeling and available experimental data. Six globular proteins with known data are chosen for analysis. Analyzing the model and the experimental results we found that water molecules bound to proteins by hydrogen bond are preferentially located at the places with local depressions on the protein surface. Consequently, the hydration level is not strictly proportional to the area of charged and polar surfaces, but also depends on the shape of the molecular surface. The thickness of the thermal volume layer as calculated in the framework of the scaled particle theory is 0.6-0.65 A for chosen proteins. The obtained value is significantly lower than that presented for proteins in earlier papers (where proportionality between the hydration level and the area of charged and polar surfaces was assumed), but is close to the value published for small solute molecules. Discussion including the influence of protein size and the thermal motion of the surface is presented.

An Investigation of Weak CH···O Hydrogen Bonds in Maltose Anomers by a Combination of Calculation and Experimental Solid-State NMR Spectroscopy

Two-dimensional (1)H-(13)C MAS-J-HMQC solid-state NMR spectra of the two anomeric forms of maltose at natural abundance are presented. The experimental (1)H chemical shifts of the CH and CH(2) protons are assigned using first-principles chemical shift calculations that employ a plane-wave pseudopotential approach. Further calculations show that the calculated change in the (1)H chemical shift when comparing the full crystal and an isolated molecule is a quantitative measure of intermolecular C-H...O weak hydrogen bonding. Notably, a clear correlation between a large chemical shift change (up to 2 ppm) and both a short H...O distance (<2.7 A) and a CHO bond angle greater than 130 degrees is observed, thus showing that directionality is important in C-H...O hydrogen bonding.

The Importance of Dynamic Effects on the Enzyme Activity

Onconase (ONC), a member of the RNase A superfamily extracted from oocytes of Rana pipiens, is an effective cancer killer. It is currently used in treatment of various forms of cancer. ONC antitumor properties depend on its ribonucleolytic activity that is low in comparison with other members of the superfamily. The most damaging side effect from Onconase treatment is renal toxicity, which seems to be caused by the unusual stability of the enzyme. Therefore, mutants with reduced thermal stability and/or increased catalytic activity may have significant implications for human cancer chemotherapy. In this context, we have determined the crystal structures of two Onconase mutants (M23L-ONC and C87S,des103-104-ONC) and performed molecular dynamic simulations of ONC and C87S,des103-104-ONC with the aim of explaining on structural grounds the modifications of the activity and thermal stability of the mutants. The results also provide the molecular bases to explain the lower catalytic activity of Onconase compared with RNase A and the unusually high thermal stability of the amphibian enzyme.

Insight into Ribonuclease A Domain Swapping by Molecular Dynamics Unfolding Simulations

Bovine pancreatic ribonuclease (RNase A) deserves a special place among the numerous proteins that form oligomers by three-dimensional domain swapping. In fact, under destabilizing conditions and at high protein concentrations, it can swap two different domains, the N-terminal alpha-helix or the C-terminal beta-strand, leading to dimers with different quaternary structures. With the change in the unfolding conditions, the relative abundance of the two dimers varies, and the prevalence of one dimer over the other is inverted. To investigate the dynamic behavior of the termini, four independent 10 ns high-temperature molecular dynamics simulations of RNase A were carried out at two different pH values in an attempt to reproduce the experimental conditions of neutral and very low pH that favor the formation of the N- and C-terminal domain-swapped dimers, respectively. In agreement with experimental data, under mild unfolding conditions, a partial or complete opening of the N-terminal arm is observed, whereas the dislocation of the C-terminus away from the core of the structure occurs only during the low-pH simulations. Furthermore, the picture emerging from this study indicates that the same protein can have different pathways for domain swapping. Indeed, in RNase A the C-terminal swapping requires a substantial unfolding of the monomers, whereas the N-terminal swapping can occur through only partial unfolding.

Protein Functional Cycle Viewed at Atomic Resolution: Conformational Change and Mobility in Nitrophorin 4 as a Function of pH and NO Binding,

The blood-sucking insect Rhodnius prolixus uses nitrophorin 4, a heme protein, to deliver nitric oxide (NO) to a victim, causing vasodilation and improved feeding. Binding of NO occurs at a ferric heme and is modulated by pH. NO binding at lower pH induces a large conformational change involving loops A-B and G-H that leads to distal pocket desolvation and protection of the nitrosyl heme complex. We have determined the crystal structures of Rhodnius nitrophorin 4 to ultrahigh resolution in four functional states: +/-NO at pH = 7.4 and +/-NO at pH = 5.6. The structure with NO at pH 7.4 (1.08 A) is newly determined while the other complexes have been modeled to resolutions much greater than previously reported (1.0-0.85 A). The ultrahigh resolution allowed us to resolve multiple conformers in binding-site loops, leading to a detailed description of the dynamics involved with storing NO in the insect salivary gland at low pH, and releasing NO in response to the increased pH of a victim's tissue. Strikingly, features for both the "open" and "closed" conformers exist under all conditions, suggesting that the flexible loops can transition with relative ease between conformational states. Yet, release of NO from rNP4 is much slower than found for other ferric heme proteins. The structures suggest that highly mobile loops can limit diffusion of diatomic molecules into and out of a protein cavity, a result with implications for the role of protein dynamics in function.

Protein hydration dynamics in solution: a critical survey

The properties of water in biological systems have been studied for well over a century by a wide range of physical techniques, but progress has been slow and erratic. Protein hydration--the perturbation of water structure and dynamics by the protein surface--has been a particularly rich source of controversy and confusion. Our aim here is to critically examine central concepts in the description of protein hydration, and to assess the experimental basis for the current view of protein hydration, with the focus on dynamic aspects. Recent oxygen-17 magnetic relaxation dispersion (MRD) experiments have shown that the vast majority of water molecules in the protein hydration layer suffer a mere twofold dynamic retardation compared with bulk water. The high mobility of hydration water ensures that all thermally activated processes at the protein-water interface, such as binding, recognition and catalysis, can proceed at high rates. The MRD-derived picture of a highly mobile hydration layer is consistent with recent molecular dynamics simulations, but is incompatible with results deduced from intermolecular nuclear Overhauser effect spectroscopy, dielectric relaxation and fluorescence spectroscopy. It is also inconsistent with the common view of hydration effects on protein hydrodynamics. Here, we show how these discrepancies can be resolved.

A second protein L-isoaspartyl methyltransferase gene in Arabidopsis produces two transcripts whose products are sequestered in the nucleus

The spontaneous and deleterious conversion of l-asparaginyl and l-aspartyl protein residues to l-iso-Asp or d-Asp occurs as proteins age and is accelerated under stressful conditions. Arabidopsis (Arabidopsis L. Heynh.) contains two genes (At3g48330 and At5g50240) encoding protein-l-isoaspartate methyltransferase (EC 2.1.1.77; PIMT), an enzyme capable of correcting this damage. The gene located on chromosome 5 (PIMT2) produces two proteins differing by three amino acids through alternative 3' splice site selection in the first intron. Recombinant protein from both splicing variants has PIMT activity. Subcellular localization using cell fractionation followed by immunoblot detection, as well as confocal visualization of PIMT:GFP fusions, demonstrated that PIMT1 is cytosolic while a canonical nuclear localization signal, present in PIMT2psi and the shorter PIMT2omega, is functional. Multiplex reverse transcription-PCR was used to establish PIMT1 and PIMT2 transcript presence and abundance, relative to beta-TUBULIN, in various tissues and under a variety of stresses imposed on seeds and seedlings. PIMT1 transcript is constitutively present but can increase, along with PIMT2, in developing seeds presumably in response to increasing endogenous abscisic acid (ABA). Transcript from PIMT2 also increases in establishing seedlings due to exogenous ABA and applied stress presumably through an ABA-dependent pathway. Furthermore, cleaved amplified polymorphic sequences from PIMT2 amplicons determined that ABA preferentially enhances the production of PIMT2omega transcript in leaves and possibly in tissues other than germinating seeds.

The occurrence of CH…O hydrogen bonds in α-helices and helix termini in globular proteins

A comprehensive database analysis of C--H...O hydrogen bonds in 3124 alpha-helices and their corresponding helix termini has been carried out from a nonredundant data set of high-resolution globular protein structures resolved at better than 2.0 A in order to investigate their role in the helix, the important protein secondary structural element. The possible occurrence of 5 --> 1 C--H...O hydrogen bond between the ith residue CH group and (i - 4)th residue C==O with C...O < or = 3.8 A is studied, considering as potential donors the main-chain Calpha and the side-chain carbon atoms Cbeta, Cgamma, Cdelta and Cepsilon. Similar analysis has been carried out for 4 --> 1 C--H...O hydrogen bonds, since the C--H...O hydrogen bonds found in helices are predominantly of type 5 --> 1 or 4 --> 1. A total of 17,367 (9310 of type 5 --> 1 and 8057 of type 4 --> 1) C--H...O hydrogen bonds are found to satisfy the selected criteria. The average stereochemical parameters for the data set suggest that the observed C--H...O hydrogen bonds are attractive interactions. Our analysis reveals that the Cgamma and Cbeta hydrogen atom(s) are frequently involved in such hydrogen bonds. A marked preference is noticed for aliphatic beta-branched residue Ile to participate in 5 --> 1 C--H...O hydrogen bonds involving methylene Cgamma 1 atom as donor in alpha-helices. This may be an enthalpic compensation for the greater loss of side-chain conformational entropy for beta-branched amino acids due to the constraint on side-chain torsion angle, namely, chi1, when they occur in helices. The preference of amino acids for 4 --> 1 C--H...O hydrogen bonds is found to be more for Asp, Cys, and for aromatic residues Trp, Phe, and His. Interestingly, overall propensity for C--H...O hydrogen bonds shows that a majority of the helix favoring residues such as Met, Glu, Arg, Lys, Leu, and Gln, which also have large side-chains, prefer to be involved in such types of weak attractive interactions in helices. The amino acid side-chains that participate in C--H...O interactions are found to shield the acceptor carbonyl oxygen atom from the solvent. In addition, C--H...O hydrogen bonds are present along with helix stabilizing salt bridges. A novel helix terminating interaction motif, X-Gly with Gly at C(cap) position having 5 --> 1 Calpha--H...O, and a chain reversal structural motif having 1 --> 5 Calpha-H...O have been identified and discussed. Our analysis highlights that a multitude of local C--H...O hydrogen bonds formed by a variety of amino acid side-chains and Calpha hydrogen atoms occur in helices and more so at the helix termini. It may be surmised that the main-chain Calpha and the side-chain CH that participate in C--H...O hydrogen bonds collectively augment the cohesive energy and thereby contribute together with the classical N--H...O hydrogen bonds and other interactions to the overall stability of helix and therefore of proteins.

Biomolecular cryocrystallography: structural changes during flash-cooling

To minimize radiation damage, crystal structures of biological macromolecules are usually determined after rapid cooling to cryogenic temperatures, some 150-200 K below the normal physiological range. The biological relevance of such structures relies on the assumption that flash-cooling is sufficiently fast to kinetically trap the macromolecule and associated solvent in a room-temperature equilibrium state. To test this assumption, we use a two-state model to calculate the structural changes expected during rapid cooling of a typical protein crystal. The analysis indicates that many degrees of freedom in a flash-cooled protein crystal are quenched at temperatures near 200 K, where local conformational and association equilibria may be strongly shifted toward low-enthalpy states. Such cryoartifacts should be most important for strongly solvent-coupled processes, such as hydration of nonpolar cavities and surface regions, conformational switching of solvent-exposed side chains, and weak ligand binding. The dynamic quenching that emerges from the model considered here can also rationalize the glass transition associated with the atomic fluctuations in the protein.

Inelastic Incoherent Neutron Scattering Measurements of Intact Cells and Tissues and Detection of Interfacial Water

We have previously used inelastic incoherent neutron scattering spectroscopy to investigate the properties of aqueous suspensions of biomolecules as a function of hydration. These experiments led to the identification of signals corresponding to interfacial (hydration) water at low water content. A prediction from these studies was that in the crowded environment inside living cells, a significant proportion of the water would be interfacial, with profound implications for biological function. Here we describe the first inelastic incoherent neutron scattering spectroscopy studies of living cells and tissues. We find that the interfacial water signal is similar to that observed for water interacting with purified biomolecules and other solutes, i.e., it is strongly perturbed in the librational and translational intermolecular optical regions of the spectrum at 20-150 meV. The ratio of interfacial water compared to total water in cells (approximately 30%) is in line with previous experimental data for hydration water and calculations based on simple assumptions.

Population shift vs induced fit: The case of bovine seminal ribonuclease swapping dimer

Bovine seminal ribonuclease (BS-RNase) is a unique member of the pancreatic-like ribonuclease superfamily. This enzyme exists as two conformational isomers with distinctive biological properties. The structure of the major isomer is characterized by the swapping of the N-terminal segment (MxM BS-RNase). In this article, the crystal structures of the ligand-free MxM BS-RNase and its complex with 2'-deoxycitidylyl(3',5')-2'-deoxyadenosine derived from isomorphous crystals have been refined. Interestingly, the comparison between this novel ligand-free form and the previously published sulfate-bound structure reveals significant differences. In particular, the ligand-free MxM BS-RNase is closer to the structure of MxM BS-RNase productive complexes than to the sulfate-bound form. These results reveal that MxM BS-RNase presents a remarkable flexibility, despite the structural constraints of the interchain disulfide bridges and the swapping of the N-terminal helices. These findings have important implications to the ligand binding mechanism of MxM BS-RNase. Indeed, a population shift rather than a substrate-induced conformational transition may occur in the MxM BS-RNase ligand binding process.

Direct Observation of Cα−Hα···OC Hydrogen Bonds in Proteins by Interresidue h3JCαC‘ Scalar Couplings

The role of C-H...O hydrogen bonds in the stabilization of biomolecules is increasingly being recognized from the evidence of close C-H...O contacts in crystal structures. However, relatively little is known about their strength. Here, we report the observation of NMR scalar couplings (h3JCalphaC') between the two carbons on each side of Calpha-Halpha...O=C H-bonds in proteins. These couplings give direct evidence of the correlation of the electronic wave functions in the donor and acceptor groups of Calpha-Halpha...O=C H-bonds. A long-range H(NCO)CA experiment or a selective long-range H(NCA)CO experiment was used for the detection of h3JCalphaC' correlations in the beta-sheet regions of the immunoglobulin binding domain of protein G. In total, six such correlations were detectable. These correspond to half of the Calpha-Halpha...O=C H-bonds of protein G with Halpha...O distances shorter than 2.5 A. The h3JCalphaC' couplings range from 0.2 to 0.3 Hz and are in good agreement with predicted average values based on DFT/FPT calculations. An anticorrelation is observed with the size of h3JNC' coupling constants across N-HN...O=C H-bonds, which share the same acceptor carbonyl oxygen.

Glycosylation and specific deamidation of ribonuclease B affect the formation of three-dimensional domain-swapped oligomers

RNase A oligomerizes via the three-dimensional domain-swapping mechanism to form a variety of oligomers, including two dimers. One, called the N-dimer, forms by swapping of the N termini of the protein; the other, called the C-dimer, forms by swapping of the C termini. RNase B is identical in protein sequence and conformation to RNase A, but its Asn34 bears an oligosaccharide chain that might affect oligomerization. The ability of RNase B to oligomerize under two sets of conditions has been examined. The amount of oligomers formed via lyophilization was somewhat lower for RNase B than RNase A, and RNase B oligomerized more rapidly in 40% ethanol solution at high temperature than RNase A. The ratio of the N-dimer to C-dimer formed increased with the size of the carbohydrate chain under both sets of conditions. These results suggest that the oligosaccharide chain either favors productive collisions or stabilizes the oligomers, especially the N-dimer. Endoglycosidase H treatment of RNase B partially restored RNase A-like oligomerization. Derivatives of RNase A conjugated at the amine groups to polyethylene glycol chains showed a greatly reduced capacity for oligomerization, suggesting that oligomerization can be impeded sterically. Commercial preparations of RNase B eluted as two main peaks by cation exchange chromatography. Using chromatography, mass spectroscopy, and two-dimensional NMR, the major peak was identified as RNase B selectively deamidated at Asn67. This deamidated protein showed a >4 degrees C drop in thermal stability, disruption of the native structure of residues 67-69, and a decreased ability to oligomerize compared with unmodified RNase B.

WaterScore: a novel method for distinguishing between bound and displaceable water molecules in the crystal structure of the binding site of protein-ligand complexes

We have performed a multivariate logistic regression analysis to establish a statistical correlation between the structural properties of water molecules in the binding site of a free protein crystal structure, with the probability of observing the water molecules in the same location in the crystal structure of the ligand-complexed form. The temperature B-factor, the solvent-contact surface area, the total hydrogen bond energy and the number of protein-water contacts were found to discriminate between bound and displaceable water molecules in the best regression functions obtained. These functions may be used to identify those bound water molecules that should be included in structure-based drug design and ligand docking algorithms. FIGURE The binding site ( thin sticks) of penicillopepsin (3app) with its crystallographically determined water molecules ( spheres) and superimposed ligand (in thick sticks, from complexed structure 1ppk). Water molecules sterically displaced by the ligand upon complexation are shown in cyan. Bound water molecules are shown in blue. Displaced water molecules are shown in yellow. Water molecules removed from the analysis due to a lack of hydrogen bonds to the protein are shown in white. WaterScore correctly predicted waters in blue as Probability=1 to remain bound and waters in yellow as Probability<1x10(-20) to remain bound.

Ultrahigh resolution structure of a class A beta-lactamase: on the mechanism and specificity of the extended-spectrum SHV-2 enzyme

Bacterial beta-lactamases hydrolyze beta-lactam antibiotics such as penicillins and cephalosporins. The TEM-type class A beta-lactamase SHV-2 is a natural variant that exhibits activity against third-generation cephalosporins normally resistant to hydrolysis by class A enzymes. SHV-2 contains a single Gly238Ser change relative to the wild-type enzyme SHV-1. Crystallographic refinement of a model including hydrogen atoms gave R and R(free) of 12.4% and 15.0% for data to 0.91 A resolution. The hydrogen atom on the O(gamma) atom of the reactive Ser70 is clearly seen for the first time, bridging to the water molecule activated by Glu166. Though hydrogen atoms on the nearby Lys73 are not seen, this observation of the Ser70 hydrogen atom and the hydrogen bonding pattern around Lys73 indicate that Lys73 is protonated. These findings support a role for the Glu166-water couple, rather than Lys73, as the general base in the deprotonation of Ser70 in the acylation process of class A beta-lactamases. Overlay of SHV-2 with SHV-1 shows a significant 1-3 A displacement in the 238-242 beta-strand-turn segment, making the beta-lactam binding site more open to newer cephalosporins with large C7 substituents and thereby expanding the substrate spectrum of the variant enzyme. The OH group of the buried Ser238 side-chain hydrogen bonds to the main-chain CO of Asn170 on the Omega loop, that is unaltered in position relative to SHV-1. This structural role for Ser238 in protein-protein binding makes less likely its hydrogen bonding to oximino cephalosporins such as cefotaxime or ceftazidime.