Computer modelling studies on the mechanism of action of ribonuclease T1

The mechanism of action of ribonuclease (RNase) T1 is still a matter of considerable debate as the results of x-ray, 2-D nmr and site-directed mutagenesis studies disagree regarding the role of the catalytically important residues. Hence computer modelling studies were carried out by energy minimisation of the complexes of RNase T1 and some of its mutants (His40Ala, His40Lys, and Glu58Ala) with the substrate guanyl cytosine (GpC), and of native RNase T1 with the reaction intermediate guanosine 2',3'-cyclic phosphate (G greater than p). The puckering of the guanosine ribose moiety in the minimum energy conformer of the RNase T1-GpC (substrate) complex was found to be O4'-endo and not C3'-endo as in the RNase T1-3'-guanylic acid (inhibitor/product) complex. A possible scheme for the mechanism of action of RNase T1 has been proposed on the basis of the arrangement of the catalytically important amino acid residues His40, Glu58, Arg77, and His92 around the guanosine ribose and the phosphate moiety in the RNase T1-GpC and RNase T1-G greater than p complexes. In this scheme, Glu58 serves as the general base group and His92 as the general acid group in the transphosphorylation step. His40 may be essential for stabilising the negatively charged phosphate moiety in the enzyme-transition state complex.

Nucleotide sequence of the 3terminal region of belladonna mottle virus (renamed Physalis mottle virus) RNA and analysis of the evolutionary relationships among tymoviral coat proteins

The 3' terminal 1255 nt sequence of Physalis mottle virus (PhMV) genomic RNA has been determined from a set of overlapping cDNA clones. The open reading frame (ORF) at the 3' terminus corresponds to the amino acid sequence of the coat protein (CP) determined earlier except for the absence of the dipeptide, Lys-Leu, at position 110-111. In addiition, the sequence upstream of the CP gene contains the message coding for 178 amino acid residues of the C-terminus of the putative replicase protein (RP). The sequence downstream of the CP gene contains an untranslated region whose terminal 80 nucleotides can be folded into a characteristic tRNA-like structure. A phylogenetic tree constructed after aligning separately the sequence of the CP, the replicase protein (RP) and the tRNA-like structure determined in this study with the corresponding sequences of other tymoviruses shows that PhMV wrongly named belladonna mottle virus [BDMV(I)] is a separate tymovirus and not another strain of BDMV(E) as originally envisaged. The phylogenetic tree in all the three cases is identical showing that any subset of genomic sequence of sufficient length can be used for establishing evolutionary relationships among tymoviruses.

Psi-angle dependence of geminal proton–carbonyl-13C coupling constants in peptides. A theoretical investigation†

Theoretical calculations of the geminal carbonyl-13C- proton coupling constant, 2J(C′H), in α-amino acids have been carried out using Dirac Vector model and Penney-Dirac bond order formulations. The results indicate that the couplings are dependent on the backbone torsion angle psi (ψ) of the amino acid residues in peptides. The meagre available experimental data seem to support the theoretical findings.

A Deterministic Optimization Approach to Protein Sequence Design Using Continuous Models

Determining the sequence of amino acid residues in a heteropolymer chain of a protein with a given conformation is a discrete combinatorial problem that is not generally amenable for gradient-based continuous optimization algorithms. In this paper we present a new approach to this problem using continuous models. In this modeling, continuous "state functions" are proposed to designate the type of each residue in the chain. Such a continuous model helps define a continuous sequence space in which a chosen criterion is optimized to find the most appropriate sequence. Searching a continuous sequence space using a deterministic optimization algorithm makes it possible to find the optimal sequences with much less computation than many other approaches. The computational efficiency of this method is further improved by combining it with a graph spectral method, which explicitly takes into account the topology of the desired conformation and also helps make the combined method more robust. The continuous modeling used here appears to have additional advantages in mimicking the folding pathways and in creating the energy landscapes that help find sequences with high stability and kinetic accessibility. To illustrate the new approach, a widely used simplifying assumption is made by considering only two types of residues: hydrophobic (H) and polar (P). Self-avoiding compact lattice models are used to validate the method with known results in the literature and data that can be practically obtained by exhaustive enumeration on a desktop computer. We also present examples of sequence design for the HP models of some real proteins, which are solved in less than five minutes on a single-processor desktop computer Some open issues and future extensions are noted.

Interaction Energy Based Protein Structure Networks

The three dimensional structure of a protein is formed and maintained by the noncovalent interactions among the amino acid residues of the polypeptide chain These interactions can be represented collectively in the form of a network So far such networks have been investigated by considering the connections based on distances between the amino acid residues Here we present a method of constructing the structure network based on interaction energies among the amino acid residues in the protein We have investigated the properties of such protein energy based networks (PENs) and have shown correlations to protein structural features such as the clusters of residues involved in stability formation of secondary and super secondary structural units Further we demonstrate that the analysis of PENs in terms of parameters such as hubs and shortest paths can provide a variety of biologically important information such as the residues crucial for stabilizing the folded units and the paths of communication between distal residues in the protein Finally the energy regimes for different levels of stabilization in the protein structure have clearly emerged from the PEN analysis

Antigenic determinants on chicken riboflavin carrier protein. A study with monoclonal antibodies

Monoclonal antibodies raised against chicken egg white riboflavin carrier protein were classified into seven categories each recognizing a distinct epitope. Of these, six were directed against conformation dependent epitopes and one to a sequential epitope. The roles of lysine residues and the post-translationally attached phosphate and oligosaccharide moieties in the antigenicity of riboflavin carrier protein recognized by the monoclonal antibodies were investigated. The binding region of three monoclonal antibodies could be located within the 87–219 amino acid sequence of the protein and one antibody among these recognized a sequence of 182–204 amino acid residues. All the monoclonal antibodies were able to recognize riboflavin carrier proteins present in the sera of pregnant rats, cows and humans indicating that the epitopes to which they are directed are conserved through evolution from chicken to the human.

Computer modeling studies on the binding of 2',5'-linked dinucleoside phosphates to ribonuclease T1-influence of subsite interactions on the substrate specificity

The modes of binding of Gp(2',5')A, Gp(2',5')C, Gp(2',5')G and Gp(2',5')U to RNase T1 have been determined by computer modelling studies. All these dinucleoside phosphates assume extended conformations in the active site leading to better interactions with the enzyme. The 5'-terminal guanine of all these ligands is placed in the primary base binding site of the enzyme in an orientation similar to that of 2'-GMP in the RNase T1-2'-GMP complex. The 2'-terminal purines are placed close to the hydrophobic pocket formed by the residues Gly71, Ser72, Pro73 and Gly74 which occur in a loop region. However, the orientation of the 2'-terminal pyrimidines is different from that of 2'-terminal purines. This perhaps explains the higher binding affinity of the 2',5'-linked guanine dinucleoside phosphates with 2'-terminal purines than those with 2'-terminal pyrimidines. A comparison of the binding of the guanine dinucleoside phosphates with 2',5'- and 3',5'-linkages suggests significant differences in the ribose pucker and hydrogen bonding interactions between the catalytic residues and the bound nucleoside phosphate implying that 2',5'-linked dinucleoside phosphates may not be the ideal ligands to probe the role of the catalytic amino acid residues. A change in the amino acid sequence in the surface loop region formed by the residues Gly71 to Gly74 drastically affects the conformation of the base binding subsite, and this may account for the inactivity of the enzyme with altered sequence i.e., with Pro, Gly and Ser at positions 71 to 73 respectively. These results thus suggest that in addition to recognition and catalytic sites, interactions at the loop regions which constitute the subsite for base binding are also crucial in determining the substrate specificity.

Identification of tropomyosin as the major shrimp allergen and characterization of its IgE-binding epitopes

The major heat-stable shrimp allergen (designated as Sa-II), capable of provoking IgE-mediated immediate type hypersensitivity reactions after the ingestion of cooked shrimp, has been shown to be a 34-kDa heat- stable protein containing 300 amino acid residues. Here, we report that a comparison of amino acid sequences of different peptides generated by proteolysis of Sa-II revealed an 86% homology with tropomyosin from Drosophila melanogaster, suggesting that Sa-II could be the shrimp muscle protein tropomyosin. To establish that Sa-II is indeed tropomyosin, the latter was isolated from uncooked shrimp (Penaeus indicus) and its physicochemical and immunochemical properties were compared with those of Sa-II. Both tropomyosin and Sa-II had the same molecular mass and focused in the isoelectric pH range of 4.8 to 5.4. In the presence of 6 M urea, the mobility of both Sa-II and shrimp tropomyosin shifted to give an apparent molecular mass of 50 kDa, which is a characteristic property of tropomyosins. Shrimp tropomyosin bound to specific IgE antibodies in the sera of shrimp-sensitive patients as assessed by competitive ELISA inhibition and Western blot analysis. Tryptic maps of both Sa-II and tropomyosin as obtained by reverse phase HPLC were superimposable. Dot-blot and competitive ELISA inhibition using sera of shrimp-sensitive patients revealed that antigenic as well as allergenic activities were associated with two peptide fractions. These IgE-binding tryptic peptides were purified and sequenced. Mouse anti-anti-idiotypic antibodies raised against Sa-II specific human idiotypic antibodies recognized not only tropomyosin but also the two allergenic peptides, thus suggesting that these peptides represent the major IgE binding epitopes of tropomyosin. A comparison of the amino acid sequence of shrimp tropomyosin in the region of IgE binding epitopes (residues 50-66 and 153-161) with the corresponding regions of tropomyosins from different vertebrates confirmed lack of allergenic cross-reactivity between tropomyosins from phylogenetically distinct species.

The complete primary structure of a unique mannose/glucose-specific lectin from field bean (Dolichos lab lab)

The complete amino acid sequence of two non identical subunits of the glucose/mannose-specific lectin from Dolichos lab lab (field bean) has been determined by sequential Edman analyses of the intact subunits and peptides derived by enzymatic and chemical cleavage. Peptides were purified by reverse phase high performance liquid chromatography and ion pair chromatography. The D. lab lab lectin is a glycoprotein having two polypeptide chains of 132 and 105 amino acid residues. The amino acid sequence of the D. Lab lab lectin is compared with the various lectins of the family Leguminosae. The D. lab lab lectin is the only species of the tribe Phaseoleae that contains two nonidentical subunits of almost equal size and that shows a specificity to glucose/ mannose. The lectin shows a greater homology to the glucose/mannose specific lectins, especially concanavalin A. The unique subunit architecture of the D. lab lab lectin indicates the presence of new post translational cleavage sites.

The primary structure of sheep liver cytosolic serine hydroxymethyltransferase and an analysis of the evolutionary relationships among serine hydroxymethyltransferases

The complete amino-acid sequence of sheep liver cytosolic serine hydroxymethyltransferase was determined from an analysis of tryptic, chymotryptic, CNBr and hydroxylamine peptides. Each subunit of sheep liver serine hydroxymethyltransferase consisted of 483 amino-acid residues. A comparison of this sequence with 8 other serine hydroxymethyltransferases revealed that a possible gene duplication event could have occurred after the divergence of animals and fungi. This analysis also showed independent duplication of SHMT genes in Neurospora crassa. At the secondary structural level, all the serine hydroxymethyltransferases belong to the alpha/beta category of proteins. The predicted secondary structure of sheep liver serine hydroxymethyltransferase was similar to that of the observed structure of tryptophan synthase, another pyridoxal 5'-phosphate containing enzyme, suggesting that sheep liver serine hydroxymethyltransferase might have a similar pyridoxal 5'-phosphate binding domain. In addition, a conserved glycine rich region, G L Q G G P, was identified in all the serine hydroxymethyltransferases and could be important in pyridoxal 5'-phosphate binding. A comparison of the cytosolic serine hydroxymethyltransferases from rabbit and sheep liver with other proteins sequenced from both these sources showed that serine hydroxymethyltransferase was a highly conserved protein. It was slightly less conserved than cytochrome c but better conserved than myoglobin, both of which are well known evolutionary markers. C67 and C203 were specifically protected by pyridoxal 5'-phosphate against modification with [C-14]iodoacetic acid, while C247 and C261 were buried in the native serine hydroxymethyltransferase. However, the cysteines are not conserved among the various serine hydroxymethyltransferases. The exact role of the cysteines in the reaction catalyzed by serine hydroxymethyltransferase remains to be elucidated.

Synthesis and Antioxidant Activity of Peptide-Based Ebselen Analogues

A series of di- and tripeptide-based ebselen analogues has been synthesized. The compounds were characterized by H-1, C-13, and Se-77 NMR spectroscopy and mass spectral techniques. The glutathione peroxidase (GPx)-like antioxidant activity has been studied by using H2O2, tert-butyl hydroperoxide (tBuOOH), and cumene hydroperoxide (Cum-OOH) as substrates, and glutathione (GSH) as a co-substrate. Although all the peptide-based compounds have a selenazole ring similar to that of ebselen, the GPx activity of these compounds highly depends on the nature of the peptide moiety attached to the nitrogen atom of the selenazole ring. It was observed that the introduction of a phenylalanine (Phe) amino acid residue in the N-terminal reduces the activity in all three peroxide systems. On the other hand, the introduction of aliphatic amino acid residues such as valine (Val) significantly enhances the GPx activity of the ebselen analogues. The difference in the catalytic activity of dipeptide-based ebselen derivatives can be ascribed mainly to the change in the reactivity of these compounds toward GSH and peroxide. Although the presence of the Val-Ala-CO2Me moiety facilitates the formation of a catalytically active selenol species, the reaction of ebselen analogues that has a Phe-Ile-CO2Me residue with GSH does not generate the corresponding selenol. To understand the antioxidant activity of the peptide-based ebselen analogues in the absence of GSH, these compounds were studied for their ability to inhibit peroxynitrite (PN)-mediated nitration of bovine serum albumin (BSA) and oxidation of dihydrorhodamine 123. In contrast to the GPx activity, the PN-scavenging activity of the Phe-based peptide analogues was found to be comparable to that of the Val-based compounds. However, the introduction of an additional Phe residue to the ebselen analogue that had a Val-Ala dipeptide significantly reduced the potency of the parent compound in PN-mediated nitration.

Primary structure of a cytotoxin-like basic protein from Naja naja naja (Indian Cobra) venom

The complete amino acid sequence of a cytotoxin-like basic protein (CLBP) from the venom of Naja naja naja (Indian Cobra) was determined by manual degradation using a 4-dimethylaminoazobenzene-4'-isothiocyanate double-coupling method. Peptide fragments obtained by chemical cleavage with cyanogen bromide and enzymic cleavages with trypsin and Staphylococcus aureus proteases for sequence analysis were purified by reversed-phase chromatography. The total number of amino acid residues was 61, with leucine as the C-terminal residue. (C) Munksgaard 1995.

RPS: Repeats in Protein Sequences

Repeats are two or more contiguous segments of amino acid residues that are believed to have arisen as a result of intragenic duplication, recombination and mutation events. These repeats can be utilized for protein structure prediction and can provide insights into the protein evolution and phylogenetic relationship. Therefore, to aid structural biologists and phylogeneticists in their research, a computing resource (a web server and a database), Repeats in Protein Sequences (RPS), has been created. Using RPS, users can obtain useful information regarding identical, similar and distant repeats (of varying lengths) in protein sequences. In addition, users can check the frequency of occurrence of the repeats in sequence databases such as the Genome Database, PIR and SWISS-PROT and among the protein sequences available in the Protein Data Bank archive. Furthermore, users can view the three-dimensional structure of the repeats using the Java visualization plug-in Jmol. The proposed computing resource can be accessed over the World Wide Web at http://bioserver1.physics.iisc.ernet.in/rps/.

Dynamics of water at the interface of a small protein, enterotoxin

Fully atomistic molecular dynamics simulations have been carried out to investigate the correlation of biological activity with dynamics of water molecules in an aqueous protein solution of the toxic domain of enterotoxin (PDB ID: 1ETN). This is a small protein of 13 amino acid residues. Our study of this water soluble protein clearly reveals that water dynamics slows down in the hydration layer. Despite this general slowing down, water molecules in the vicinity of the second beta turn of this protein exhibit faster dynamics than those near other regions of the protein. Since this beta turn is believed to play a critical role in the receptor binding of this protein, the faster dynamics of water near the beta turn m ay have biological significance. The collective orientational dynamics of the water molecules in the protein solution exhibits a characteristic long time component of 27 ps, which agrees well with dielectric relaxation experiments.

Helix and hairpin nucleation in short peptides using centrally positioned conformationally constrained dipeptide segments

The effect of incorporation of a centrally positioned Ac(6)c-Xxx segment where Xxx = (L)Val/(D)Val into a host oligopeptide composed of L-amino acid residues has been investigated. Studies of four designed octapeptides Boc-Leu-Phe-Val-Ac(6)c-Xxx-Leu-Phe-Val-OMe (Xxx = (D)Val 1, (L)Val 2) Boc-Leu-Val-Val-Ac(6)c-Xxx-Leu-Val-Val-OMe (Xxx = (D)Val 3, (L)Val 4) are reported. Diagnostic nuclear Overhouse effects characteristic of hairpin conformations are observed for Xxx = (D)Val peptides (1 and 3) while continuous helical conformation characterized by sequential NiH <-> Ni+1H NOEs are favored for Xxx = (L)Val peptides (2 and 4) in methanol solutions. Temperature co-efficient of NH chemical shifts are in agreement with distinctly different conformational preferences upon changing the configuration of the residue at position 5. Crystal structures of peptides 2 and 4 (Xxx = (L)Val) establish helical conformations in the solid state, in agreement with the structures deduced from NMR data. The results support the design principle that centrally positioned type I beta-turns may be used to nucleate helices in short peptides, while type I' beta-turns can facilitate folding into beta-hairpins.