SIRIUS. An automated method for the analysis of the preferred packing arrangements between protein groups

Automated methods have been developed to determine the preferred packing arrangement between interacting protein groups. A suite of FORTRAN programs, SIRIUS, is described for calculating and analysing the geometries of interacting protein groups using crystallographically derived atomic co-ordinates. The programs involved in calculating the geometries search for interacting pairs of protein groups using a distance criterion, and then calculate the spatial disposition and orientation of the pair. The second set of programs is devoted to analysis. This involves calculating the observed and expected distributions of the angles and assessing the statistical significance of the difference between the two. A database of the geometries of the 400 combinations of side-chain to side-chain interaction has been created. The approach used in analysing the geometrical information is illustrated here with specific examples of interactions between side-chains, peptide groups and particular types of atom. At the side-chain level, an analysis of aromatic-amino interactions, and the interactions of peptide carbonyl groups with arginine residues is presented. At the atomic level the analyses include the spatial disposition of oxygen atoms around tyrosine residues, and the frequency and type of contact between carbon, nitrogen and oxygen atoms. This information is currently being applied to the modelling of protein interactions.     

Cation-pi interactions in protein-protein interfaces

Arginine is an abundant residue in protein-protein interfaces. The importance of this residue relates to the versatility of its side chain in intermolecular interactions. Different classes of protein-protein interfaces were surveyed for cation-pi interactions. Approximately half of the protein complexes and one-third of the homodimers analyzed were found to contain at least one intermolecular cation-pi pair. Interactions between arginine and tyrosine were found to be the most abundant. The electrostatic interaction energy was calculated to be approximately 3 kcal/mol, on average. A distance-based search of guanidinium:aromatic interactions was also performed using the Macromolecular Structure Database (MSD). This search revealed that half of the guanidinium:aromatic pairs pack in a coplanar manner. Furthermore, it was found that the cationic group of the cation-pi pair is frequently involved in intermolecular hydrogen bonds. In this manner the arginine side chain can participate in multiple interactions, providing a mechanism for inter-protein specificity. Thus, the cation-pi interaction is established as an important contributor to protein-protein interfaces.     

A combined computational and experimental investigation of the [2Fe–2S] cluster in biotin synthase

Biotin synthase was the first example of what is now regarded as a distinctive enzyme class within the radical S-adenosylmethionine superfamily, the members of which use Fe/S clusters as the sulphur source in radical sulphur insertion reactions. The crystal structure showed that this enzyme contains a [2Fe–2S] cluster with a highly unusual arginine ligand, besides three normal cysteine ligands. However, the crystal structure is at such a low resolution that neither the exact coordination mode nor the role of this exceptional ligand has been elucidated yet, although it has been shown that it is not essential for enzyme activity. We have used quantum refinement of the crystal structure and combined quantum mechanical and molecular mechanical calculations to explore possible coordination modes and their influences on cluster properties. The investigations show that the protonation state of the arginine ligand has little influence on cluster geometry, so even a positively charged guanidinium moiety would be in close proximity to the iron atom. Nevertheless, the crystallised enzyme most probably contains a deprotonated (neutral) arginine coordinating via the NH group. Furthermore, the Fe···Fe distance seems to be independent of the coordination mode and is in perfect agreement with distances in other structurally characterised [2Fe–2S] clusters. The exceptionally large Fe···Fe distance found in the crystal structure could not be reproduced.     

Salt bridges destabilize a leucine zipper designed for maximized ion pairing between helices

Interhelical salt bridges are common in leucine zippers and are thought to stabilize the coiled coil conformation. Here we present a detailed thermodynamic investigation of the designed, disulfide-linked leucine zipper AB(SS) whose high-resolution NMR structure shows six interhelical ion pairs between heptad positions g of one helix and e' of the other helix but no ion pairing within single helices. The average pK(a) value of the Glu side chain carboxyl groups of AB(SS) is slightly higher than the pK(a) of a freely accessible Glu in an unfolded peptide [Marti, D. N., Jelesarov, I., and Bosshard, H. R. (2000) Biochemistry 39, 12804-12818]. This indicates that the salt bridges are destabilizing, a prediction we now have confirmed by determining the pH +/- stability profile of AB(SS). Circular dichroism-monitored unfolding by urea and by heating and differential scanning calorimetry show that the coiled coil conformation is approximately 5 kJ/mol more stable when salt bridges are broken by protonation of the carboxyl side chains. Using guanidinium chloride as the denaturant, the increase in the free energy of unfolding on protonation of the carboxyl side chains is larger, approximately 17 kJ/mol. The discrepancy between urea and guanidinium chloride unfolding can be ascribed to the ionic nature of guanidinium chloride, which screens charge-charge interactions. This work demonstrates the difficulty of predicting the energetic contribution of salt bridges from structural data alone even in a case where the ion pairs are seen in high-resolution NMR structures. The reason is that the contribution to stability results from a fine balance between energetically favorable Coulombic attractions and unfavorable desolvation of charges and conformational constraints of the residues involved in ion pairing. The apparent discrepancy between the results presented here and mutational studies indicating stabilization by salt bridges is discussed and resolved. An explanation is proposed for why interhelical salt bridges are frequently found in natural coiled coils despite evidence that they do not directly contribute to stability.     

Netropsin interactions in the minor groove of d(GGCCAATTGG) studied by a combination of resolution enhancement and ab initio calculations

The structure of the complex between the minor groove binder netropsin and d(GGCCAATTGG) was determined via single-crystal X-ray techniques. The structure was refined to completion using refmac5.1.24, resulting in a residual R-factor of 20.0% (including 68 water molecules). Using crystal engineering and cryocooling techniques, the resolution could be enhanced to 1.75 A, resulting in an unambiguous determination of the drug conformation and orientation. As previously noticed, bifurcated hydrogen bonds are formed between the amide nitrogen atoms of the drug and the N3 and O2 atoms of A and T base pairs, respectively, clearly cataloging the structure to class I. As the bulky NH2 group on guanine was believed to prevent binding of the drug in the minor groove, the detailed nature of several of the amidinium and guanidinium end contacts were further investigated by ab initio quantum chemical methods.     

Database of non-canonical base pairs found in known RNA structures

Atomic resolution RNA structures are being published at an increasing rate. It is common to find a modest number of non-canonical base pairs in these structures in addition to the usual Watson-Crick pairs. This database summarizes the occurrence of these rare base pairs in accordance with standard nomenclature. The database, http://prion.bchs.uh.edu/, contains information such as sequence context, sugar pucker conformation, anti / syn base conformations, chemical shift, p K (a)values, melting temperature and free energy. Of the 29 anticipated pairs with two or more hydrogen bonds, 20 have been encountered to date. In addition, four unexpected pairs with two hydrogen bonds have been reported bringing the total to 24. Single hydrogen bond versions of five of the expected geometries have been encountered among the single hydrogen bond interactions. In addition, 18 different types of base triplets have been encountered, each of which involves three to six hydrogen bonds. The vast majority of the rare base pairs are antiparallel with the bases in the anti configuration relative to the ribose. The most common are the GU wobble, the Sheared GA pair, the Reverse Hoogsteen pair and the GA imino pair.     

Nonsequence-specific arginine interactions in the nucleosome core particle

We have analyzed the relative orientation of basic amino acid side chains towards DNA in the nucleosome core particle. The electric field created by DNA phosphates has no apparent preferential orientation: no favored orientation of the arginine guanidinium group is found. Arginine may be either directly hydrogen bonded to a phosphate oxygen or stabilized in the minor groove by van der Waals contacts and the local negative electric field. On the other hand, the phosphate oxygen atoms hydrogen bonded to arginines are always found close to the plane defined by the guanidinium group. Thus it can be concluded that the interactions of arginine are strongly directional, those of phosphate are not. We also find that a highly charged fragment of histone H2B, which is placed between two DNA turns, has a very variable conformation. An increase in protein positive charge density apparently allows multiple nonspecific protein conformations when interacting with DNA.     

Spectrophotometric and steady-state kinetic analysis of the biosynthetic arginine decarboxylase of Yersinia pestis utilizing arginine analogues as inhibitors and alternative substrates

The PLP-dependent, biosynthetic arginine decarboxylase (ADC) of Yersinia pestis was investigated using steady-state kinetics employing structural analogues of arginine as both alternative substrates and competitive inhibitors. The inhibitor analysis indicates that binding of the carboxyl and guanidinium groups of the substrate, l-arginine, provides essentially all of the free energy change realized upon substrate binding in the ground state. Furthermore, recognition of the guanidinium group is primarily responsible for substrate specificity. Comparison of the steady-state parameters for a series of alternative substrates that contained chemically modified guanidinium moieties provides evidence of a role for induced fit in ADC catalysis. ADC was also characterized by UV/vis and fluorescence spectrophotometry in the presence or absence of a number of arginine analogues. The enzyme complexes formed served as models for the adsorption complex and the external aldimine complex of the enzyme with the substrate.     

Hydration of transfer RNA molecules: a crystallographic study

Four crystal structures of transfer RNA molecules were refined at 3 A resolution with the inclusion of the solvent molecules found in the difference maps: yeast tRNA-phe in the orthorhombic form, yeast tRNA-phe in the monoclinic form and yeast tRNA-asp in the A and B forms. Over 100 solvent molecules were located in each tRNA crystal. Several hydration schemes are found repeatedly in the 4 crystals. The tertiary interactions in the corner of the L-shaped molecule attract numerous solvent molecules which bridge the ribose hydroxyl O(2') atoms, base exocyclic atoms and phosphate anionic oxygen atoms. Conservation of bases leads to conservative localized hydration patterns. Several solvent molecules are found stabilizing unusual base pairs like the G-U pairs and those involving the pseudouridine base. Water bridges between the O(2') and the exocyclic atom O2 of pyrimidines or the N3 atom of purines are common. Water bridges occur frequently between successive anionic oxygen atoms of each strand as well as between N7 or other exocyclic atoms of successive bases in the major groove. Magnesium ions or spermine molecules are found to bind in the major groove of tRNA helices without specific interactions.     

Relationship between ion pair geometries and electrostatic strengths in proteins

The electrostatic free energy contribution of an ion pair in a protein depends on two factors, geometrical orientation of the side-chain charged groups with respect to each other and the structural context of the ion pair in the protein. Conformers in NMR ensembles enable studies of the relationship between geometry and electrostatic strengths of ion pairs, because the protein structural contexts are highly similar across different conformers. We have studied this relationship using a dataset of 22 unique ion pairs in 14 NMR conformer ensembles for 11 nonhomologous proteins. In different NMR conformers, the ion pairs are classified as salt bridges, nitrogen-oxygen (N-O) bridges and longer-range ion pairs on the basis of geometrical criteria. In salt bridges, centroids of the side-chain charged groups and at least a pair of side-chain nitrogen and oxygen atoms of the ion-pairing residues are within a 4 A distance. In N-O bridges, at least a pair of the side-chain nitrogen and oxygen atoms of the ion-pairing residues are within 4 A distance, but the distance between the side-chain charged group centroids is greater than 4 A. In the longer-range ion pairs, the side-chain charged group centroids as well as the side-chain nitrogen and oxygen atoms are more than 4 A apart. Continuum electrostatic calculations indicate that most of the ion pairs have stabilizing electrostatic contributions when their side-chain charged group centroids are within 5 A distance. Hence, most (approximately 92%) of the salt bridges and a majority (68%) of the N-O bridges are stabilizing. Most (approximately 89%) of the destabilizing ion pairs are the longer-range ion pairs. In the NMR conformer ensembles, the electrostatic interaction between side-chain charged groups of the ion-pairing residues is the strongest for salt bridges, considerably weaker for N-O bridges, and the weakest for longer-range ion pairs. These results suggest empirical rules for stabilizing electrostatic interactions in proteins.     

The environment of amide groups in protein–ligand complexes: H-bonds and beyond

A comprehensive structural analysis of interactions involving amide NH and C=O groups in protein-ligand complexes has been performed based on 3,275 published crystal structures (resolution < or =2.5 A). Most of the amide C=O and NH groups at the protein-ligand interface are highly buried within the binding site and involved in H-bonds with corresponding counter-groups. Small percentages of C=O and NH groups are solvated or embedded in hydrophobic environments. In particular, C=O groups show a higher propensity to be solvated or embedded in a hydrophobic environment than NH groups do. A small percentage of carbonyl groups is involved in weak hydrogen bonds with CH. Cases of dipolar interactions, involving carbonyl oxygen and electrophilic carbon atoms, such as amide, amidinium, guanidium groups, are also identified. A higher percentage of NH are in contact with aromatic carbons, interacting either through hydrogen bonds (preferably with the NH group pointing towards a ring carbon atom) or through stacking between amide plane and ring plane. Comprehensive studies such as the present one are thought to be important for future improvements in the molecular design area, in particular for the development of new scoring functions. [Figure: see text].     

Geometry of guanidinium groups in arginines

The restraints in common usage today have been obtained based on small molecule X‐ray crystal structures available 25 years ago and recent reports have shown that the values of bond lengths and valence angles can be, in fact, significantly different from those stored in libraries, for example for the peptide bond or the histidine ring geometry. We showed that almost 50% of outliers found in protein validation reports released in the Protein Data Bank on 23 March 2016 come from geometry of guanidine groups in arginines. Therefore, structures of small molecules and atomic resolution protein crystal structures have been used to derive new target values for the geometry of this group. The most significant difference was found for NE‐CZ‐NH1 and NE‐CZ‐NH2 angles, showing that the guanidinium group is not symmetric. The NE‐CZ‐NH1 angle is larger, 121.5(10)?, than NE‐CZ‐NH2, 119.2(10)?, due to the repulsive interaction between NH1 and CD1 atom.     

Mechanism of charge compensation and impairment of co-operative functions in haemoglobin Tacoma (Arg B12(30)beta leads to Ser).

A difference Fourier synthesis of deoxyhaemoglobin Tacoma minus deoxyhaemoglobin A at 3.5 Å resolution has been calculated. The map shows a large negative peak due to the removal of the guanidinium group of Arg B12(30)β, surrounded by positive and negative peaks indicative of some atoms moving towards, and others away from, the vacated site. Among the latter, the most important is the carboxylate of Glu B8(26)β which is hydrogen-bonded to the guanidinium group of the arginine in haemoglobin A, but swings round its α-β carbon bond towards the imidazoles of histidines G18(116) and 19(117)β in haemoglobin Tacoma. This movement would raise the pK values of the histidines, so that their positive charges compensated for the loss of the arginine. This may explain why haemoglobin Tacoma has the same electrophoretic mobility as haemoglobin A. It is shown that haemoglobin Tacoma has a lower oxygen equilibrium constant K_T and a larger allosteric constant L than haemoglobin A. The lowering of K_T may be due to the loosening of the T structure in haemoglobin Tacoma consequent upon the removal of the hydrogen bonds made by the guanidinium group of Arg B12(30)β at the α₁β₁ contact. Their removal also accounts for the decreased stability of haemoglobin Tacoma. We cannot yet explain its diminished Bohr effect, nor the increase in L.     

Conserved interactions of the active carboxyls in pepsin-like enzymes and retroviral proteases

In addition to previous studies, 30 crystal structures of retroviral proteases corresponding to the highest resolution were inspected to analyze the interactions of the active carboxyl with surroundings groups. The outer oxygen of the active carboxyl in retroviral enzymes form contacts only with the water molecule participating in catalysis. This is an important difference between retroviral proteases and pepsin-like enzymes, which form a net of hydrogen bonds of these outer oxygen with residues neighboring the catalytic site in 3D structures. At the same time, it was found that in all aspartic proteases the inner oxygen of the active carboxyl are also involved in the chain of interactions through peptide groups Thr-Gly adjacent to the active residues. Polarization of these peptide groups influences the donor-acceptor properties of the active carboxyl. The found chain of interactions is more extensive in retroviral than in pepsin-like proteases; however, its main part is conserved for the whole class of these enzymes. Some implications of the role of these interactions are discussed.     

2.0 Å resolution structure of a ternary complex of pig muscle phosphoglycerate kinase containing 3-phospho-D-glycerate and the nucleotide Mn adenylylimidodiphosphate

The crystal structure of a ternary complex of pig muscle phosphoglycerate kinase (PGK) containing 3-phosphoglycerate (3-PG) and manganese adenylylimidodiphosphate (Mn AMP-PNP) has been determined and refined at 2.0 A resolution. The complex differs from the true substrate ternary complex only in the presence of an imido- rather than an oxylink between β- and γ-phosphates of the bound nucleotide. The 3-PG is bound in a similar manner to that observed in binary complexes. The nucleotide is bound in a similar manner to Mg ADP except that the metal ion is coordinated by all three α-, β-, and γ-phosphates, but not by the protein. The γ-phosphate, which is transferred in the reaction, is not bound by the protein. One further characteristic of the ternary complex is that Arg-38 moves to a position where its guanidinium group makes a triple interaction with the N-terminal domain, the C-terminal domain, and the 1-carboxyl group of the bound 3-PG. Although a hinge-bending conformation change is seen in the ternary complex, it is no larger than that observed in the 3-PG binary complex. To reduce that distance between two bound substrates to a value consistent with the direct in-line transfer known to occur in PGK, we modeled the closure of a pronounced cleft in the protein structure situated between the bound substrates. This closure suggested a mechanism of catalysis that involves the “capture” of the γ-phosphate by Arg-38 and the N-terminus of helix-14, which has a conserved Gly-Gly-Gly phosphate binding motif. We propose that nucleophilic attack by the 1-carboxyl group of the 3-PG on the γ-phosphorus follows the capture of the γ-phosphate, leading to a pentacoordinate transition state that may be stabilized by hydrogen bonds donated by the NH groups in the N-terminus of helix 14 and the guanidinium group of Arg-38. During the course of the reaction the metal ion is proposed to migrate to a position coordinating the α- and β-phosphates and the carboxyl group of Asp-374. The mechanism is consistent with the structural information from binary and ternary substrate complexes and much solution data, and gives a major catalytic role to Arg-38, as indicated by site-directed mutagenesis.     

SuperStar: improved knowledge-based interaction fields for protein binding sites

SuperStar is an empirical method for identifying interaction sites in proteins, based entirely on experimental information about non-bonded interactions occurring in small-molecule crystal structures, taken from the IsoStar database. We describe recent modifications and additions to SuperStar, validating the results on a test set of 122 X-ray structures of protein-ligand complexes. In this validation, propensity maps are generated for all the binding sites of these proteins, using four different probes: a charged NH(+)(3) nitrogen atom, a carbonyl oxygen atom, a hydroxyl oxygen atom and a methyl carbon atom. Next, the maps are compared with the experimentally observed positions of ligand atoms of these types. A peak-searching algorithm is introduced that highlights potential interaction hot spots. For the three hydrogen-bonding probes - NH(+)(3) nitrogen atom, carbonyl oxygen atom and hydroxyl oxygen atom - the average distance from the ligand atom to the nearest SuperStar peak is 1.0-1.2 A (0.8-1.0 A for solvent-inaccessible ligand atoms). For the methyl carbon atom probe, this distance is about 1.5 A, probably because interactions to methyl groups are much less directional.The most important addition to SuperStar is the enabling of propensity maps around metal centres - Ca(2+), Mg(2+) and Zn(2+) - in protein binding sites. The results are validated on a test set of 24 protein-ligand complexes that have a metal ion in their binding site. Coordination geometries are derived automatically, using only the protein atoms that coordinate to the metal ion. The correct coordination geometry is derived in approximately 75 % of the cases. If the derived geometry is assumed during the SuperStar calculation, the average distance from a ligand atom coordinating to the metal ion to the nearest peak in the propensity map for an oxygen probe is 0.87(7) A. If the correct coordination geometry is imposed, this distance reduces to 0.59(7)A. This indicates that the SuperStar predictions around metal-binding sites are at least as good as those around other protein groups. Using clustering techniques, a non-redundant set of probes is selected from the set of probes available in the IsoStar database. The performance in SuperStar of all these probes is tested on the test set of protein-ligand complexes. With the exception of the "ether oxygen" probe and the "any NH(+)" probe, all new probes perform as well as the four probes introduced first.     

Sequence specificity, statistical potentials, and three-dimensional structure prediction with self-correcting distance geometry calculations of beta-sheet formation in proteins

A statistical analysis of a representative data set of 169 known protein structures was used to analyze the specificity of residue interactions between spatial neighboring strands in beta-sheets. Pairwise potentials were derived from the frequency of residue pairs in nearest contact, second nearest and third nearest contacts across neighboring beta-strands compared to the expected frequency of residue pairs in a random model. A pseudo-energy function based on these statistical pairwise potentials recognized native beta-sheets among possible alternative pairings. The native pairing was found within the three lowest energies in 73% of the cases in the training data set and in 63% of beta-sheets in a test data set of 67 proteins, which were not part of the training set. The energy function was also used to detect tripeptides, which occur frequently in beta-sheets of native proteins. The majority of native partners of tripeptides were distributed in a low energy range. Self-correcting distance geometry (SECODG) calculations using distance constraints sets derived from possible low energy pairing of beta-strands uniquely identified the native pairing of the beta-sheet in pancreatic trypsin inhibitor (BPTI). These results will be useful for predicting the structure of proteins from their amino acid sequence as well as for the design of proteins containing beta-sheets.     

Pyruvate phosphate dikinase: sequence of the histidyl peptide, the pyrophosphoryl and phosphoryl carrier

Pyruvate phosphate dikinase contains a pivotal histidyl residue which functions to mediate the transfer of phosphoryl moieties during the reaction catalyzed by the enzyme. The tryptic peptide which contains this essential histidyl residue has been isolated by a two-step procedure originally developed by Wang and co-workers [Wang, T., Jurasek, L., & Bridger, W. A. (1972) Biochemistry 11, 2067]. This peptide has been sequenced by the manual dansyl-Edman procedure and is shown to be NH2-Gly-Gly-Met-Thr-Ser-His-Ala-Ala-Val-Val-Ala-Arg-CO2H. There is no readily interpretable homology between this peptide and other phosphorylated histidyl peptides previously isolated from other enzymes. By use of Chou & Fasman [Chou, P. Y., & Fasman, G. D. (1974) Biochemistry 13, 222], it is predicted that the sequence contains an alpha helix from the methionine residue through to the carboxyl terminal arginine residue.     

T1-T1 interactions occur in ER membranes while nascent Kv peptides are still attached to ribosomes

For voltage-gated K+ channels (Kv), it is not clear at which stage during biosynthesis in the endoplasmic reticulum (ER) oligomerization occurs, specifically whether it can begin while nascent peptide chains of individual subunits are still attached to ribosomes. Kv channels possess a T1-recognition domain in the NH2-terminus, which confers subfamily specificity for intersubunit assembly and forms a tetramer. Using pairs of cysteines engineered into the T1-T1 interface and cross-linking methods, we show that specific residues in the T1-T1 interface of different Kv1.3 subunits come into close proximity in the ER, both in microsomal membranes and in Xenopus oocytes. Furthermore, using translocation intermediates containing pairs of engineered cysteines in the T1 interface, we demonstrate that specific residues in the folded T1 domain interface can approach within 2 A of each other and form tetramers while the nascent Kv1.3 peptides are still attached to ribosomes and have translocated across the membrane. ER membranes are required for this interaction, and T1-T1 interactions occur inter-polysomally. Thus, folding of the T1 domain and intersubunit interaction may represent the first assembly event in channel formation.