Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In α-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of α-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 Å. The sch:mch interactions involving ionizable side chains that belong either to β-strands or to the central part of α-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.     

Structural analysis of protein‐ligand interactions: the binding of endogenous compounds and of synthetic drugs

The large number of macromolecular structures deposited with the Protein Data Bank (PDB) describing complexes between proteins and either physiological compounds or synthetic drugs made it possible a systematic analysis of the interactions occurring between proteins and their ligands. In this work, the binding pockets of about 4000 PDB protein‐ligand complexes were investigated and amino acid and interaction types were analyzed. The residues observed with lowest frequency in protein sequences, Trp, His, Met, Tyr, and Phe, turned out to be the most abundant in binding pockets. Significant differences between drug‐like and physiological compounds were found. On average, physiological compounds establish with respect to drugs about twice as many hydrogen bonds with protein atoms, whereas drugs rely more on hydrophobic interactions to establish target selectivity. The large number of PDB structures describing homologous proteins in complex with the same ligand made it possible to analyze the conservation of binding pocket residues among homologous protein structures bound to the same ligand, showing that Gly, Glu, Arg, Asp, His, and Thr are more conserved than other amino acids. Also in the cases in which the same ligand is bound to unrelated proteins, the binding pockets showed significant conservation in the residue types. In this case, the probability of co‐occurrence of the same amino acid type in the binding pockets could be up to thirteen times higher than that expected on a random basis. The trends identified in this study may provide an useful guideline in the process of drug design and lead optimization. Copyright © 2014 John Wiley & Sons, Ltd.     

Suggestions for "safe" residue substitutions in site-directed mutagenesis

The conserved topological structure observed in various molecular families such as globins or cytochromes c allows structural equivalencing of residues in every homologous structure and defines in a coherent way a global alignment in each sequence family. A search was performed for equivalent residue pairs in various topological families that were buried in protein cores or exposed at the protein surface and that had mutated but maintained similar unmutated environments. Amino acid residues with atoms in contact with the mutated residue pairs defined the environment. Matrices of preferred amino acid exchanges were then constructed and preferred or avoided amino acid substitutions deduced. Given the conserved atomic neighborhoods, such natural in vivo substitutions are subject to similar constrains as point mutations performed in site-directed mutagenesis experiments. The exchange matrices should provide guidelines for "safe" amino acid substitutions least likely to disturb the protein structure, either locally or in its overall folding pathway, and most likely to allow probing the structural and functional significance of the substituted site.     

A persulfurated cysteine promotes active site reactivity in Azotobacter vinelandii Rhodanese

Active site reactivity and specificity of RhdA, a thiosulfate:cyanide sulfurtransferase (rhodanese) from Azotobacter vinelandii, have been investigated through ligand binding, site-directed mutagenesis, and X-ray crystallographic techniques, in a combined approach. In native RhdA the active site Cys230 is found persulfurated; fluorescence and sulfurtransferase activity measurements show that phosphate anions interact with Cys230 persulfide sulfur atom and modulate activity. Crystallographic analyses confirm that phosphate and hypophosphite anions react with native RhdA, removing the persulfide sulfur atom from the active site pocket. Considering that RhdA and the catalytic subunit of Cdc25 phosphatases share a common three-dimensional fold as well as active site Cys (catalytic) and Arg residues, two RhdA mutants carrying a single amino acid insertion at the active site loop were designed and their phosphatase activity tested. The crystallographic and functional results reported here show that specific sulfurtransferase or phosphatase activities are strictly related to precise tailoring of the catalytic loop structure in RhdA and Cdc25 phosphatase, respectively.     

The crystal structure of a sulfurtransferase from Azotobacter vinelandii highlights the evolutionary relationship between the rhodanese and phosphatase enzyme families

Rhodanese is an ubiquitous enzyme that in vitro catalyses the transfer of a sulfur atom from suitable donors to nucleophilic acceptors by way of a double displacement mechanism. During the catalytic process the enzyme cycles between a sulfur-free and a persulfide-containing form, via formation of a persulfide linkage to a catalytic Cys residue. In the nitrogen-fixing bacteria Azotobacter vinelandii the rhdA gene has been identified and the encoded protein functionally characterized as a rhodanese. The crystal structure of the A. vinelandii rhodanese has been determined and refined at 1.8 A resolution in the sulfur-free and persulfide-containing forms. Conservation of the overall three-dimensional fold of bovine rhodanese is observed, with substantial modifications of the protein structure in the proximity of the catalytic residue Cys230. Remarkably, the native enzyme is found as the Cys230-persulfide form; in the sulfur-free state the catalytic Cys residue adopts two alternate conformations, reflected by perturbation of the neighboring active-site residues, which is associated with a partly reversible loss of thiosulfate:cyanide sulfurtransferase activity. The catalytic mechanism of A. vinelandii rhodanese relies primarily on the main-chain conformation of the 230 to 235 active-site loop and on a surrounding strong positive electrostatic field. Substrate recognition is based on residues which are entirely different in the prokaryotic and eukaryotic enzymes. The active-site loop of A. vinelandii rhodanese displays striking structural similarity to the active-site loop of the similarly folded catalytic domain of dual specific phosphatase Cdc25, suggesting a common evolutionary origin of the two enzyme families.     

ENVIRON: a software package to compare protein three-dimensional structures with homologous sequences using local structural motifs

This work presents a method to compare local clusters of interactingresidues as observed in a known three-dimensional protein structurewith corresponding clusters inferred from homologous proteinsequences, assuming conserved protein folding. For this purposethe local environment of a selected residue in a known proteinstructure is defined as the ensemble of amino acids in contactwith it in the folded state. Using a multiple sequence alignmentto identify corresponding residues in homologous proteins, adetailed comparison can be performed between the local environmentof a selected amino acid in the template protein structure andthe expected local environments at the sets of equivalent residues,derived from the aligned protein sequences. The comparison makesit possible to detect conserved local features such as hydrogenbonding or complementarity in residue substitution. A globalmeasure of environmental similarity is also defined, to searchfor conserved amino acid clusters subject to functional or structural constraints. The proposed approach is useful for investigatingprotein function as well as for site-directed mutagenesis experiments,where appropriate amino acid substitutions can be suggestedby observing naturally occurring protein variants.