An amino acid as a cofactor for a catalytic polynucleotide

Natural ribozymes require metal ion cofactors that aid both in structural folding and in chemical catalysis. In contrast, many protein enzymes produce dramatic rate enhancements using only the chemical groups that are supplied by their constituent amino acids. This fact is widely viewed as the most important feature that makes protein a superior polymer for the construction of biological catalysts. Herein we report the in vitro selection of a catalytic DNA that uses histidine as an active component for an RNA cleavage reaction. An optimized deoxyribozyme from this selection requires l-histidine or a closely related analog to catalyze RNA phosphoester cleavage, producing a rate enhancement of ≈1-million-fold over the rate of substrate cleavage in the absence of enzyme. Kinetic analysis indicates that a DNA–histidine complex may perform a reaction that is analogous to the first step of the proposed catalytic mechanism of RNase A, in which the imidazole group of histidine serves as a general base catalyst. Similarly, ribozymes of the “RNA world” may have used amino acids and other small organic cofactors to expand their otherwise limited catalytic potential.

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.

Twenty-first aminoacyl-tRNA synthetase–suppressor tRNA pairs for possible use in site-specific incorporation of amino acid analogues into proteins in eukaryotes and in eubacteria

Two critical requirements for developing methods for the site-specific incorporation of amino acid analogues into proteins in vivo are (i) a suppressor tRNA that is not aminoacylated by any of the endogenous aminoacyl-tRNA synthetases (aaRSs) and (ii) an aminoacyl-tRNA synthetase that aminoacylates the suppressor tRNA but no other tRNA in the cell. Here we describe two such aaRS–suppressor tRNA pairs, one for use in the yeast Saccharomyces cerevisiae and another for use in Escherichia coli. The “21st synthetase–tRNA pairs” include E. coli glutaminyl-tRNA synthetase (GlnRS) along with an amber suppressor derived from human initiator tRNA, for use in yeast, and mutants of the yeast tyrosyl-tRNA synthetase (TyrRS) along with an amber suppressor derived from E. coli initiator tRNA, for use in E. coli. The suppressor tRNAs are aminoacylated in vivo only in the presence of the heterologous aaRSs, and the aminoacylated tRNAs function efficiently in suppression of amber codons. Plasmids carrying the E. coli GlnRS gene can be stably maintained in yeast. However, plasmids carrying the yeast TyrRS gene could not be stably maintained in E. coli. This lack of stability is most likely due to the fact that the wild-type yeast TyrRS misaminoacylates the E. coli proline tRNA. By using error-prone PCR, we have isolated and characterized three mutants of yeast TyrRS, which can be stably expressed in E. coli. These mutants still aminoacylate the suppressor tRNA essentially quantitatively in vivo but show increased discrimination in vitro for the suppressor tRNA over the E. coli proline tRNA by factors of 2.2- to 6.8-fold.

Identification of thermophilic species by the amino acid compositions deduced from their genomes

The global amino acid compositions as deduced from the complete genomic sequences of six thermophilic archaea, two thermophilic bacteria, 17 mesophilic bacteria and two eukaryotic species were analysed by hierarchical clustering and principal components analysis. Both methods showed an influence of several factors on amino acid composition. Although GC content has a dominant effect, thermophilic species can be identified by their global amino acid compositions alone. This study presents a careful statistical analysis of factors that affect amino acid composition and also yielded specific features of the average amino acid composition of thermophilic species. Moreover, we introduce the first example of a ‘compositional tree’ of species that takes into account not only homologous proteins, but also proteins unique to particular species. We expect this simple yet novel approach to be a useful additional tool for the study of phylogeny at the genome level.

The interacting domains of three MutL heterodimers in man: hMLH1 interacts with 36 homologous amino acid residues within hMLH3, hPMS1 and hPMS2

In human cells, hMLH1, hMLH3, hPMS1 and hPMS2 are four recognised and distinctive homologues of MutL, an essential component of the bacterial DNA mismatch repair (MMR) system. The hMLH1 protein forms three different heterodimers with one of the other MutL homologues. As a first step towards functional analysis of these molecules, we determined the interacting domains of each heterodimer and tried to understand their common features. Using a yeast two-hybrid assay, we show that these MutL homologues can form heterodimers by interacting with the same amino acid residues of hMLH1, residues 492–742. In contrast, three hMLH1 partners, hMLH3, hPMS1 and hPMS2 contain the 36 homologous amino acid residues that interact strongly with hMLH1. Contrary to the previous studies, these homologous residues reside at the N-terminal regions of three subdomains conserved in MutL homologues in many species. Interestingly, these residues in hPMS2 and hMLH3 may form coiled-coil structures as predicted by the MULTICOIL program. Furthermore, we show that there is competition for the interacting domain in hMLH1 among the three other MutL homologues. Therefore, the quantitative balance of these three MutL heterodimers may be important in their functions.

Control of gating mode by a single amino acid residue in transmembrane segment IS3 of the N-type Ca2+ channel

N-type Ca2+ channels can be inhibited by neurotransmitter-induced release of G protein βγ subunits. Two isoforms of Cav2.2 α1 subunits of N-type calcium channels from rat brain (Cav2.2a and Cav2.2b; initially termed rbB-I and rbB-II) have different functional properties. Unmodulated Cav2.2b channels are in an easily activated “willing” (W) state with fast activation kinetics and no prepulse facilitation. Activating G proteins shifts Cav2.2b channels to a difficult to activate “reluctant” (R) state with slow activation kinetics; they can be returned to the W state by strong depolarization resulting in prepulse facilitation. This contrasts with Cav2.2a channels, which are tonically in the R state and exhibit strong prepulse facilitation. Activating or inhibiting G proteins has no effect. Thus, the R state of Cav2.2a and its reversal by prepulse facilitation are intrinsic to the channel and independent of G protein modulation. Mutating G177 in segment IS3 of Cav2.2b to E as in Cav2.2a converts Cav2.2b tonically to the R state, insensitive to further G protein modulation. The converse substitution in Cav2.2a, E177G, converts it to the W state and restores G protein modulation. We propose that negatively charged E177 in IS3 interacts with a positive charge in the IS4 voltage sensor when the channel is closed and produces the R state of Cav2.2a by a voltage sensor-trapping mechanism. G protein βγ subunits may produce reluctant channels by a similar molecular mechanism.

An essential amino acid residue in the protein translocation channel revealed by targeted random mutagenesis of SecY

The SecY/Sec61α family of membrane proteins are the central subunits of the putative protein translocation channel. We introduced random mutations into a segment of Escherichia coli SecY within its cytoplasmic domain 5, which was shown previously to be important for the SecA-dependent translocation activity. Mutations were classified into those retaining function and those gaining a dominant-interfering ability caused by a loss of function. These analyses showed that Arg-357, Pro-358, Gly-359, and Thr-362 are functionally important; Arg-357, conserved in almost all organisms, was identified as an indispensable residue.

Direct evidence for modified solvent structure within the hydration shell of a hydrophobic amino acid.

Neutron scattering experiments are used to determine scattering profiles for aqueous solutions of hydrophobic and hydrophilic amino acid analogs. Solutions of hydrophobic solutes show a shift in the main diffraction peak to smaller angle as compared with pure water, whereas solutions of hydrophilic solutes do not. The same difference for solutions of hydrophobic and hydrophilic side chains is also predicted by molecular dynamics simulations. The neutron scattering curves of aqueous solutions of hydrophobic amino acids at room temperature are qualitatively similar to differences between the liquid molecular structure functions measured for ambient and supercooled water. The nonpolar solute-induced expansion of water structure reported here is also complementary to recent neutron experiments where compression of aqueous solvent structure has been observed at high salt concentration.

Mapping rat megalin: the second cluster of ligand binding repeats contains a 46-amino acid pathogenic epitope involved in the formation of immune deposits in Heymann nephritis.

Megalin (gp330), an epithelial endocytic receptor, is a major target antigen of Heymann nephritis (HN), an autoimmune disease in rats. To elucidate the mechanisms of HN, we have mapped a pathogenic epitope in megalin that binds anti-megalin antibodies. We focused our attention on four clusters of cysteine-rich, low density lipoprotein receptor (LDLR) ligand binding repeats in the extracellular domain of megalin because they represent putative ligand binding regions and therefore would be expected to be exposed in vivo and to be able to bind circulating antibodies. Rat megalin cDNA fragments I through IV encoding the first through fourth clusters of ligand-binding repeats, respectively, were expressed in a baculovirus system. All four expression products were detected by immunoblotting with two antisera capable of inducing passive HN (pHN). When antibodies eluted from glomeruli of rats with pHN were used for immunoblotting, only the expression product encoded by fragment II was detected. This indicates that the second cluster of LDLR ligand binding repeats is directly involved in binding anti-megalin antibodies and in the induction of pHN. To narrow the major epitope in this domain, fragment II was used to prepare proteins sequentially truncated from the C- and N-terminal ends by in vitro translation. Analysis of the truncated translation products by immunoprecipitation with anti-megalin IgG revealed that the fifth ligand-binding repeat (amino acids 1160-1205) contains the major epitope recognized. This suggests that a 46-amino acid sequence in the second cluster of LDLR ligand binding repeats contains a major pathogenic epitope that plays a key role in pHN. Identification of this epitope will facilitate studies on the pathogenesis of HN.

Anti-peptide aptamers recognize amino acid sequence and bind a protein epitope.

In vitro selection of nucleic acid binding species (aptamers) is superficially similar to the immune response. Both processes produce biopolymers that can recognize targets with high affinity and specificity. While antibodies are known to recognize the sequence and conformation of protein surface features (epitopes), very little is known about the precise interactions between aptamers and their epitopes. Therefore, aptamers that could recognize a particular epitope, a peptide fragment of human immunodeficiency virus type I Rev, were selected from a random sequence RNA pool. Several of the selected RNAs could bind the free peptide more tightly than a natural RNA ligand, the Rev-binding element. In accord with the hypothesis that protein and nucleic acid binding cusps are functionally similar, interactions between aptamers and the peptide target could be disrupted by sequence substitutions. Moreover, the aptamers appeared to be able to bind peptides with different solution conformations, implying an induced fit mechanism for binding. Just as anti-peptide antibodies can sometimes recognize the corresponding epitope when presented in a protein, the anti-peptide aptamers were found to specifically bind to Rev.

Crystal structure of D-amino acid oxidase: a case of active site mirror-image convergent evolution with flavocytochrome b2.

D-amino acid oxidase is the prototype of the FAD-dependent oxidases. It catalyses the oxidation of D-amino acids to the corresponding alpha-ketoacids. The reducing equivalents are transferred to molecular oxygen with production of hydrogen peroxide. We have solved the crystal structure of the complex of D-amino acid oxidase with benzoate, a competitive inhibitor of the substrate, by single isomorphous replacement and eightfold averaging. Each monomer is formed by two domains with an overall topology similar to that of p-hydroxybenzoate hydroxylase. The benzoate molecule lays parallel to the flavin ring and is held in position by a salt bridge with Arg-283. Analysis of the active site shows that no side chains are properly positioned to act as the postulated base required for the catalytic carboanion mechanism. On the contrary, the benzoate binding mode suggests a direct transfer of the substrate alpha-hydrogen to the flavin during the enzyme reductive half-reaction.The active site Of D-amino acid oxidase exhibits a striking similarity with that of flavocytochrome b2, a structurally unrelated FMN-dependent flavoenzyme. The active site groups (if these two enzymes are in fact superimposable once the mirror-image of the flavocytochrome b2 active site is generated with respect to the flavin plane. Therefore, the catalytic sites of D-amino acid oxidase and flavocytochrome b2 appear to have converged to a highly similar but enantiomeric architecture in order to catalvze similar reactions (oxidation of alpha-amino acids or alpha-hydroxy acids), although with opposite stereochemistry.

Glycerol-induced development of catalytically active conformation of Crotalus adamanteus L-amino acid oxidase in vitro.

The reconstitutable apoprotein of Crotalus adamanteus L-amino acid oxidase was prepared using hydrophobic interaction chromatography. After reconstitution with flavin adenine dinucleotide, the resulting protein was inactive, with a perturbed conformation of the flavin binding site. Subsequently, a series of cosolvent-dependent compact intermediates was identified. The nearly complete activation of the reconstituted apoprotein and the restoration of its native flavin binding site was achieved in the presence of 50% glycerol. We provide evidence that in addition to a merely stabilizing effect of glycerol on native proteins, glycerol can also have a restorative effect on their compact equilibrium intermediates, and we suggest the hydrophobic effect as a dominating force in this in vitro-assisted restorative process.

Interactions between tRNA identity nucleotides and their recognition sites in glutaminyl-tRNA synthetase determine the cognate amino acid affinity of the enzyme.

Sequence-specific interactions between aminoacyl-tRNA synthetases and their cognate tRNAs both ensure accurate RNA recognition and prevent the binding of noncognate substrates. Here we show for Escherichia coli glutaminyl-tRNA synthetase (GlnRS; EC 6.1.1.18) that the accuracy of tRNA recognition also determines the efficiency of cognate amino acid recognition. Steady-state kinetics revealed that interactions between tRNA identity nucleotides and their recognition sites in the enzyme modulate the amino acid affinity of GlnRS. Perturbation of any of the protein-RNA interactions through mutation of either component led to considerable changes in glutamine affinity with the most marked effects seen at the discriminator base, the 10:25 base pair, and the anticodon. Reexamination of the identity set of tRNA(Gln) in the light of these results indicates that its constituents can be differentiated based upon biochemical function and their contribution to the apparent Gibbs' free energy of tRNA binding. Interactions with the acceptor stem act as strong determinants of tRNA specificity, with the discriminator base positioning the 3' end. The 10:25 base pair and U35 are apparently the major binding sites to GlnRS, with G36 contributing both to binding and recognition. Furthermore, we show that E. coli tryptophanyl-tRNA synthetase also displays tRNA-dependent changes in tryptophan affinity when charging a noncognate tRNA. The ability of tRNA to optimize amino acid recognition reveals a novel mechanism for maintaining translational fidelity and also provides a strong basis for the coevolution of tRNAs and their cognate synthetases.

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.