Structural assessment of single amino acid mutations: application to TP53 function.

Single amino acid substitution is the type of protein alteration most related to human diseases. Current studies seek primarily to distinguish neutral mutations from harmful ones. Very few methods offer an explanation of the final prediction result in terms of the probable structural or functional effect on the protein. In this study, we describe the use of three novel parameters to identify experimentally-verified critical residues of the TP53 protein (p53). The first two parameters make use of a surface clustering method to calculate the protein surface area of highly conserved regions or regions with high nonlocal atomic interaction energy (ANOLEA) score. These parameters help identify important functional regions on the surface of a protein. The last parameter involves the use of a new method for pseudobinding free-energy estimation to specifically probe the importance of residue side-chains to the stability of protein fold. A decision tree was designed to optimally combine these three parameters. The result was compared to the functional data stored in the International Agency for Research on Cancer (IARC) TP53 mutation database. The final prediction achieved a prediction accuracy of 70% and a Matthews correlation coefficient of 0.45. It also showed a high specificity of 91.8%. Mutations in the 85 correctly identified important residues represented 81.7% of the total mutations recorded in the database. In addition, the method was able to correctly assign a probable functional or structural role to the residues. Such information could be critical for the interpretation and prediction of the effect of missense mutations, as it not only provided the fundamental explanation of the observed effect, but also helped design the most appropriate laboratory experiment to verify the prediction results.

Effect of environmental conditions on the N uptake route of soil microorganisms and adaption of a method to determine amino acid oxidase in soil

Soil microorganisms have evolved two possible mechanisms for their uptake of organic N: the direct route and the mobilization-immobilization-turnover (MIT) route. In the direct route, simple organic molecules are taken up via various mechanisms directly into the cell. In the MIT route, the deamination occurs outside the cell and all N is mineralized to NH4+ before assimilation. A better understanding of the mechanisms controlling the different uptake routes of soil microorganisms under different environmental conditions is crucial for understanding mineralization processes of organic material in soil. For the first experiment we incubated soil samples from the long term trial in Bad Lauchstädt with corn residues with different C to N ratios and inorganic N for 21 days at 20 °C. Under the assumption that all added amino acids were taken up or mineralized, the direct uptake route was more important in soil amended with corn residues with a wide C to N ratio. After 21 days of incubation the direct uptake of added amino acids increased in the order addition of corn residue with a: “C to N ratio of 40 & (NH4)2SO4 and no addition (control)” (69% and 68%, respectively) < “C to N ratio of 20” (73%) < “C to N ratio of 40” (95%). In all treatments the proportion of the added amino acids that were mineralized increased with time, indicating that the MIT route became more important over time. To investigate the effects of soil depth on the N uptake route of soil microorganisms (experiment II), soil samples in two soil depths (0-5 cm; 30-40 cm) were incubated with corn residues with different C to N ratios and inorganic N for 21 days at 20 °C and 60% (WHC). The addition of corn residue resulted in a marked increase of protease activity in both depths due to the induction from the added substrate. Addition of corn residue with a wide C to N ratio resulted in a significantly greater part of the direct uptake (97% and 94%) than without the addition of residues (85% and 80%) or addition of residue with a small C to N ratio (90% and 84%) or inorganic N (91% and 79% in the surface soil and subsoil, respectively), suggesting that under conditions of sufficient mineralizable N (C to N ratio of 20) or increased concentrations of NH4+, the enzyme system involved in the direct uptake is slightly repressed. Substrate additions resulted in an initially significantly higher increase of the direct uptake in the surface soil than in the subsoil. As a large proportion of the organic N input into soil is in form of proteinaceous material, the deamination of amino acids is a key reaction of the MIT route. Therefore the enzyme amino acid oxidase contribute to the extracellular N mineralization in soil. The objective of experiment III was to adapt a method to determine amino acid oxidase in soil. The detection via synthetic fluorescent Lucifer Yellow derivatives of the amino acid lysine is possible in soil. However, it was not possible to find the substrate concentration at which the reaction rate is independent of substrate concentration and therefore we were not able to develop a valid soil enzyme assay.

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favourability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB) versus Thr(HB)-Asn(nHB). Significantly (p <= 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi(1) rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi(1) conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi(1) and chi(2) conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favourability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http:// www.rubic.rdg.ac.uk/betapairprefsparallel/). (c) 2005 Elsevier Ltd. All rights reserved.

A single amino acid in the HA of pH1N1 2009 influenza virus affects cell tropism in human airway epithelium, but not transmission in ferrets

The first pandemic of the 21(st) century, pandemic H1N1 2009 (pH1N1 2009), emerged from a swine-origin source. Although human infections with swine-origin influenza have been reported previously, none went on to cause a pandemic or indeed any sustained human transmission. In previous pandemics, specific residues in the receptor binding site of the haemagglutinin (HA) protein of influenza have been associated with the ability of the virus to transmit between humans. In the present study we investigated the effect of residue 227 in HA on cell tropism and transmission of pH1N1 2009. In pH1N1 2009 and recent seasonal H1N1 viruses this residue is glutamic acid, whereas in swine influenza it is alanine. Using human airway epithelium, we show a differential cell tropism of pH1N1 2009 compared to pH1N1 2009 E227A and swine influenza suggesting this residue may alter the sialic acid conformer binding preference of the HA. Furthermore, both pH1N1 2009 E227A and swine influenza multi-cycle viral growth was found to be attenuated in comparison to pH1N1 2009 in human airway epithelium. However this altered tropism and viral growth in human airway epithelium did not abrogate respiratory droplet transmission of pH1N1 2009 E227A in ferrets. Thus, acquisition of E at residue 227 was not solely responsible for the ability of pH1N1 2009 to transmit between humans.

Rapid protein domain assignment from amino acid sequence using predicted secondary structure

The elucidation of the domain content of a given protein sequence in the absence of determined structure or significant sequence homology to known domains is an important problem in structural biology. Here we address how successfully the delineation of continuous domains can be accomplished in the absence of sequence homology using simple baseline methods, an existing prediction algorithm (Domain Guess by Size), and a newly developed method (DomSSEA). The study was undertaken with a view to measuring the usefulness of these prediction methods in terms of their application to fully automatic domain assignment. Thus, the sensitivity of each domain assignment method was measured by calculating the number of correctly assigned top scoring predictions. We have implemented a new continuous domain identification method using the alignment of predicted secondary structures of target sequences against observed secondary structures of chains with known domain boundaries as assigned by Class Architecture Topology Homology (CATH). Taking top predictions only, the success rate of the method in correctly assigning domain number to the representative chain set is 73.3%. The top prediction for domain number and location of domain boundaries was correct for 24% of the multidomain set (±20 residues). These results have been put into context in relation to the results obtained from the other prediction methods assessed

Amino acid sequence of a myotoxic Lys49-phospholipase A(2) homologue from the venom of Cerrophidion (Bothrops) godmani

The complete amino acid sequence of myotoxin II (godMT-II), a myotoxic phospholipase A( 2 )(PLA(2)) homologue from the venom of the Central American crotaline snake Cerrophidion (Bothrops) godmani, was determined by direct protein sequencing methods. GodMT-II is a class II PLA, showing a Lys instead of Asp at position 49. An additional substitution in the calcium binding loop region (Asn instead of Tyr at position 28) suggests the lack of enzymatic activity observed in this toxin is due to loss of its ability to bind the co-factor Ca2+, since the residues involved in forming the catalytic network of PLA(2)s (His-48, Tyr-52 and Asp-99) an conserved in godMT-II. This myotoxin shows highest sequence homology with other Lys-49 PLA(2)s from Bothrops, Agkistrodon and Trimeresurus species, suggesting that they constitute a conserved family of proteins, yet in contrast presents lower homology with Bothrops asper myotoxin III, a catalytically-active PLA(2). The C-terminal region of godMT-II, which is rich in cationic and hydrophobic residues, shares high sequence homology to the corresponding region in the myotoxin II from B. asper, which has been proposed to play an important role in the Ca2+-independent membrane damaging activity. (C) 1998 Elsevier B.V. B.V. All rights reserved.

Amino acid sequence and crystal structure of BaP1, a metalloproteinase from Bothrops asper snake venom that exerts multiple tissue-damaging activities

BaP1 is a 22.7-kD P-I-type zinc-dependent metalloproteinase isolated from the venom of the snake Bothrops asper, a medically relevant species in Central America. This enzyme exerts multiple tissue-damaging activities, including hemorrhage, myonecrosis, dermonecrosis, blistering, and edema. BaP1 is a single chain of 202 amino acids that shows highest sequence identity with metalloproteinases isolated front the venoms of snakes of the subfamily Crotalinae. It has six Cys residues involved in three disulfide bridges (Cys 117-Cys 197, Cys 159-Cys 181, Cys 157-Cys 164). It has the consensus sequence H(142)E(143)XXH(146)XXGXXH(152), as well as the sequence C164I165M166, which characterize the metzincin superfamily of metalloproteinases. The active-site cleft separates a major subdomain (residues 1-152), comprising four a-helices and a five-stranded beta-sheet, from the minor subdomain, which is formed by a single a-helix and several loops. The catalytic zinc ion is coordinated by the N-epsilon2 nitrogen atoms of His 142, His 146, and His 152, in addition to a solvent water molecule, which in turn is bound to Glu 143. Several conserved residues contribute to the formation of the hydrophobic pocket, and Met 166 serves as a hydrophobic base for the active-site groups. Sequence and structural comparisons of hemorrhagic and nonhemorrhagic P-I metalloproteinases from snake venoms revealed differences in several regions. In particular, the loop comprising residues 153 to 176 has marked structural differences between metalloproteinases with very different hemorrhagic activities. Because this region lies in close proximity to the active-site microenvironment, it may influence the interaction of these enzymes with physiologically relevant substrates in the extracellular matrix.

AMINO-ACID-SEQUENCE OF TSTX-V, AN ALPHA-TOXIN FROM TITYUS-SERRULATUS SCORPION-VENOM, AND ITS EFFECT ON K+ PERMEABILITY OF BETA-CELLS FROM ISOLATED RAT ISLETS OF LANGERHANS

Highly purified Tityustoxin V (TsTX-V), an alpha-toxin isolated from the venom of the Brazilian scorpion Tityus serrulatus, was obtained by ion exchange chromatography on carboxymethylcellulose-52. It was shown to be homogeneous by reverse phase high performance liquid chromatography, N-terminal sequencing (first 39 residues) of the reduced and alkylated protein and by polyacrylamide gel electrophoresis in the presence of sodium dodecylsulfate and tricine. Following enzymatic digestion, the complete amino acid sequence (64 residues) was determined. The sequence showed higher homology with the toxins from the venoms of the North African than with those of the North and South American scorpions. Using the rate of Rb-86(+) release from depolarized rat pancreatic beta-cells as a measure of K+ permeability changes, TsTX-V (5.6 mu g/ml) was found to increase by 2.0-2.4-fold the rate of marker outflow in the presence of 8.3 mM glucose. This effect was persistent and slowly reversible, showing similarity to that induced by 100 mu-M veratridine, an agent that increases the open period of Na+ channels, delaying their inactivation. It is suggested that, by extending the depolarized period, TsTX-V indirectly affects beta-cell voltage-dependent K+ channels, thus increasing K+ permeability.

Identification of amino acid substitutions with compensational effects in the attachment protein of canine distemper virus.

The hemagglutinin (H) gene of canine distemper virus (CDV) encodes the receptor-binding protein. This protein, together with the fusion (F) protein, is pivotal for infectivity since it contributes to the fusion of the viral envelope with the host cell membrane. Of the two receptors currently known for CDV (nectin-4 and the signaling lymphocyte activation molecule [SLAM]), SLAM is considered the most relevant for host susceptibility. To investigate how evolution might have impacted the host-CDV interaction, we examined the functional properties of a series of missense single nucleotide polymorphisms (SNPs) naturally accumulating within the H-gene sequences during the transition between two distinct but related strains. The two strains, a wild-type strain and a consensus strain, were part of a single continental outbreak in European wildlife and occurred in distinct geographical areas 2 years apart. The deduced amino acid sequence of the two H genes differed at 5 residues. A panel of mutants carrying all the combinations of the SNPs was obtained by site-directed mutagenesis. The selected mutant, wild type, and consensus H proteins were functionally evaluated according to their surface expression, SLAM binding, fusion protein interaction, and cell fusion efficiencies. The results highlight that the most detrimental functional effects are associated with specific sets of SNPs. Strikingly, an efficient compensational system driven by additional SNPs appears to come into play, virtually neutralizing the negative functional effects. This system seems to contribute to the maintenance of the tightly regulated function of the H-gene-encoded attachment protein. Importance: To investigate how evolution might have impacted the host-canine distemper virus (CDV) interaction, we examined the functional properties of naturally occurring single nucleotide polymorphisms (SNPs) in the hemagglutinin gene of two related but distinct strains of CDV. The hemagglutinin gene encodes the attachment protein, which is pivotal for infection. Our results show that few SNPs have a relevant detrimental impact and they generally appear in specific combinations (molecular signatures). These drastic negative changes are neutralized by compensatory mutations, which contribute to maintenance of an overall constant bioactivity of the attachment protein. This compensational mechanism might reflect the reaction of the CDV machinery to the changes occurring in the virus following antigenic variations critical for virulence.

Structure and function in rhodopsin: Rhodopsin mutants with a neutral amino acid at E134 have a partially activated conformation in the dark state*

The Glu-134–Arg-135 residues in rhodopsin, located near the cytoplasmic end of the C helix, are involved in G protein binding, or activation, or both. Furthermore, the charge-neutralizing mutation Glu-134 to Gln-134 produces hyperactivity in the activated state and produces constitutive activity in opsin. The Glu/Asp-Arg charge pair is highly conserved in equivalent positions in other G protein-coupled receptors. To investigate the structural consequences of charge-neutralizing mutations at Glu-134 and Arg-135 in rhodopsin, single spin-labeled side chains were introduced at sites in the cytoplasmic domains of helices C (140), E (227), F (250), or G (316) to serve as “molecular sensors” of the local helix bundle conformation. In each of the spin-labeled rhodopsins, a Gln substitution was introduced at either Glu-134 or Arg-135, and the electron paramagnetic resonance spectrum of the spin label was used to monitor the structural response of the helix bundle. The results indicate that a Gln substitution at Glu-134 induces a photoactivated conformation around helices C and G even in the dark state, an observation of potential relevance to the hyperactivity and constitutive activity of the mutant. In contrast, little change is induced in helix F, which has been shown to undergo a dominant motion upon photoactivation. This result implies that the multiple helix motions accompanying photoactivation are not strongly coupled and can be induced to take place independently. Gln substitution at Arg-135 produces only minor structural changes in the dark- or light-activated conformation, suggesting that this residue is not a determinant of structure in the regions investigated, although it may be functionally important.

Two amino acid substitutions convert a guanylyl cyclase, RetGC-1, into an adenylyl cyclase

Guanylyl cyclases (GCs) and adenylyl cyclases (ACs) have fundamental roles in a wide range of cellular processes. Whereas GCs use GTP as a substrate to form cGMP, ACs catalyze the analogous conversion of ATP to cAMP. Previously, a model based on the structure of adenylate cyclase was used to predict the structure of the nucleotide-binding pocket of a membrane guanylyl cyclase, RetGC-1. Based on this model, we replaced specific amino acids in the guanine-binding pocket of GC with their counterparts from AC. A change of two amino acids, E925K together with C995D, is sufficient to completely alter the nucleotide specificity from GTP to ATP. These experiments strongly validate the AC-derived RetGC-1 structural model and functionally confirm the role of these residues in nucleotide discrimination.

CKAAPs DB: a conserved key amino acid positions database

The Conserved Key Amino Acid Positions DataBase (CKAAPs DB) provides access to an analysis of structurally similar proteins with dissimilar sequences where key residues within a common fold are identified. The derivation and significance of CKAAPs starting from pairwise structure alignments is described fully in Reddy et al. [Reddy,B.V.B., Li,W.W., Shindyalov,I.N. and Bourne,P.E. (2000) Proteins, in press]. The CKAAPs identified from this theoretical analysis are provided to experimentalists and theoreticians for potential use in protein engineering and modeling. It has been suggested that CKAAPs may be crucial features for protein folding, structural stability and function. Over 170 substructures, as defined by the Combinatorial Extension (CE) database, which are found in approximately 3000 representative polypeptide chains have been analyzed and are available in the CKAAPs DB. CKAAPs DB also provides CKAAPs of the representative set of proteins derived from the CE and FSSP databases. Thus the database contains over 5000 representative poly­peptide chains, covering all known structures in the PDB. A web interface to a relational database permits fast retrieval of structure-sequence alignments, CKAAPs and associated statistics. Users may query by PDB ID, protein name, function and Enzyme Classification number. Users may also submit protein alignments of their own to obtain CKAAPs. An interface to display CKAAPs on each structure from a web browser is also being implemented. CKAAPs DB is maintained by the San Diego Supercomputer Center and accessible at the URL http://ckaaps.sdsc.edu.

Global properties of the mapping between local amino acid sequence and local structure in proteins.

Local protein structure prediction efforts have consistently failed to exceed approximately 70% accuracy. We characterize the degeneracy of the mapping from local sequence to local structure responsible for this failure by investigating the extent to which similar sequence segments found in different proteins adopt similar three-dimensional structures. Sequence segments 3-15 residues in length from 154 different protein families are partitioned into neighborhoods containing segments with similar sequences using cluster analysis. The consistency of the sequence-to-structure mapping is assessed by comparing the local structures adopted by sequence segments in the same neighborhood in proteins of known structure. In the 154 families, 45% and 28% of the positions occur in neighborhoods in which one and two local structures predominate, respectively. The sequence patterns that characterize the neighborhoods in the first class probably include virtually all of the short sequence motifs in proteins that consistently occur in a particular local structure. These patterns, many of which occur in transitions between secondary structural elements, are an interesting combination of previously studied and novel motifs. The identification of sequence patterns that consistently occur in one or a small number of local structures in proteins should contribute to the prediction of protein structure from sequence.

A single amino acid in gamma-aminobutyric acid rho 1 receptors affects competitive and noncompetitive components of picrotoxin inhibition.

A class of bicuculline-insensitive gamma-aminobutyric acid (GABA) receptors, GABAC, has been identified in retina. Several lines of evidence indicate that GABAC receptors are formed partially or wholly of GABA rho subunits. These receptors generate a Cl- current in response to GABA but differ from GABAA receptors in a number of ways. Picrotoxin, widely accepted as a noncompetitive antagonist of GABAA receptors, displays competitive and noncompetitive antagonism of GABAC receptors in perch and bovine retina and GABA rho 1 receptors expressed in Xenopus oocytes. The aim of this study was to identify the molecular basis of the two components of picrotoxin inhibition of GABA rho 1 receptors. By using a domain-swapping and mutagenesis strategy, a difference in picrotoxin sensitivity between rho 1 and rho 2 receptors was localized to a single amino acid in the putative second transmembrane domain. Substitution of this amino acid with residues found in the analogous position in highly picrotoxin-sensitive glycine alpha and GABAA subunits increased the sensitivity of rho 1 mutants 10- to 500-fold. Importantly, the competitive component of picrotoxin inhibition of the rho 1 mutant receptors was almost eliminated. These findings demonstrate that an amino acid in the putative channel domain of GABA rho 1 receptors influences picrotoxin sensitivity and mediates agonist binding by an allosteric mechanism.