Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Diversity of DNA methyltransferases that recognize asymmetric target sequences

DNA methyltransferases (MTases) are a group of enzymes that catalyze the methyl group transfer from S-adenosyl-L-methionine in a sequence-specific manner. Orthodox Type II DNA MTases usually recognize palindromic DNA sequences and add a methyl group to the target base (either adenine or cytosine) on both strands. However, there are a number of MTases that recognize asymmetric target sequences and differ in their subunit organization. In a bacterial cell, after each round of replication, the substrate for any MTase is hemimethylated DNA, and it therefore needs only a single methylation event to restore the fully methylated state. This is in consistent with the fact that most of the DNA MTases studied exist as monomers in solution. Multiple lines of evidence suggest that some DNA MTases function as dimers. Further, functional analysis of many restriction-modification systems showed the presence of more than one or fused MTase genes. It was proposed that presence of two MTases responsible for the recognition and methylation of asymmetric sequences would protect the nascent strands generated during DNA replication from cognate restriction endonuclease. In this review, MTases recognizing asymmetric sequences have been grouped into different subgroups based on their unique properties. Detailed characterization of these unusual MTases would help in better understanding of their specific biological roles and mechanisms of action. The rapid progress made by the genome sequencing of bacteria and archaea may accelerate the identification and study of species- and strain-specific MTases of host-adapted bacteria and their roles in pathogenic mechanisms.

Comparison of promoter-specific events during transcription initiation in mycobacteria

DNA protein interactions that occur during transcription initiation play an important role in regulating gene expression. To initiate transcription, RNA polymerase (RNAP) binds to promoters in a sequence-specific fashion. This is followed by a series of steps governed by the equilibrium binding and kinetic rate constants, which in turn determine the overall efficiency of the transcription process. We present here the first detailed kinetic analysis of promoter RNAP interactions during transcription initiation in the sigma(A)-dependent promoters P-rrnAPCL1, P-rrnB and P-gyr of Mycobacterium smegmatis. The promoters show comparable equilibrium binding affinity but differ significantly in open complex formation, kinetics of isomerization and promoter clearance. Furthermore, the two rrn promoters exhibit varied kinetic properties during transcription initiation and appear to be subjected to different modes of regulation. In addition to distinct kinetic patterns, each one of the housekeeping promoters studied has its own rate-limiting step in the initiation pathway, indicating the differences in their regulation.

Crystallization and preliminary X-ray studies of a galactose-specific lectin from the seeds of bitter gourd (Momordica charantia)

A galactose-specific lectin from the seeds of bitter gourd (Momordica charantia) is a four-chain type II ribosome-inactivating protein (RIP) resulting from covalent association through a disulfide bridge between two identical copies of a two-chain unit. The available structural information on such four-chain RIPs is meagre. The bitter gourd lectin was therefore crystallized for structural investigation and the crystals have been characterized. It is anticipated that the structure of the orthorhombic crystals will be analysed using molecular replacement by taking advantage of its sequence, and presumably structural, homology to normal two-chain type II RIPs.

Sequence specificity of Z-DNA formation in oligonucleotides

The sequence specific requirement for B----Z transition in solution was examined in d(CGTGCGCACG), d(CGTACGTACG), d(ACGTACGT) in presence of various Z-inducing factors. Conformational studies show that inspite of the alternating nature of purines and pyrimidines, the aforementioned sequences do not undergo B----Z transition under the influence of NaCl, hexamine cobalt chloride and ethanol. A comparison with the crystal structures of an assorted array of purine and pyrimidine sequences show that the sequence requirement for B----Z transition is much more stringent in solution as compared to the solid state. The disruptive influence of AT base pairs in B to Z transition is discussed.

Generation of a Manganese Specific Restriction Endonuclease with Nicking Activity

A typical feature of type II restriction endonucleases (REases) is their obligate sequence specificity and requirement for Mg2+ during catalysis. R.KpnI is an exception. Unlike most other type II REases, the active site of this enzyme can accommodate Mg2+, Mn2+, Ca2+, or Zn2+ and cleave DNA. The enzyme belongs to the HNH superfamily of nucleases and is characterized by the presence of a beta beta alpha-Me finger motif. Residues D148, H149, and Q175 together form the HNH active site and are essential for Mg2+ binding and catalysis. The unique ability of the enzyme to cleave DNA in the presence of different metal ions is exploited to generate mutants that are specific to one particular metal ion. We describe the generation of a Mn2+-dependent sequence specific endonuclease, defective in DNA cleavage with Mg2+ and other divalent metal ions. In the engineered mutant, only Mn2+ is selectively bound at the active site, imparting Mn2+-mediated cleavage. The mutant is impaired in concerted double-stranded DNA cleavage, leading to accumulation of nicked intermediates. The nicking activity of the mutant enzyme is further enhanced by altered reaction conditions. The active site fluidity of R Eases allowing flexible accommodation of catalytic cofactors thus forms a basis for engineering selective metal ion-dependent REase additionally possessing nicking activity.

Interaction of Sesbania Mosaic Virus Movement Protein with VPg and P10: Implication to Specificity of Genome Recognition

Sesbania mosaic virus (SeMV) is a single strand positive-sense RNA plant virus that belongs to the genus Sobemovirus. The mechanism of cell-to-cell movement in sobemoviruses has not been well studied. With a view to identify the viral encoded ancillary proteins of SeMV that may assist in cell-to-cell movement of the virus, all the proteins encoded by SeMV genome were cloned into yeast Matchmaker system 3 and interaction studies were performed. Two proteins namely, viral protein genome linked (VPg) and a 10-kDa protein (P10) c v gft encoded by OFR 2a, were identified as possible interacting partners in addition to the viral coat protein (CP). Further characterization of these interactions revealed that the movement protein (MP) recognizes cognate RNA through interaction with VPg, which is covalently linked to the 59 end of the RNA. Analysis of the deletion mutants delineated the domains of MP involved in the interaction with VPg and P10. This study implicates for the first time that VPg might play an important role in specific recognition of viral genome by MP in SeMV and shed light on the possible role of P10 in the viral movement.

Analysis of the effects of amino-acid-sequence on the structure of proteins

The conformation of amino acid side chains as observed in well-determined structures of globular proteins has earlier been extensively investigated. In contrast, the structural features of the polypeptide backbone that result from the occurrence of specific amino acids along the polypeptide have not been analysed. In this article, we present the statistically significant features in the backbone geometry that appear to be a consequence of the occurrence of rotamers of different amino acid side chains by analysing 102 well-refined structures that form a random collection of proteins. It is found that the persistence of helical segments around each residue is influenced by the residue type. Several residues exert asymmetrical influence between the carboxyl and amino terminal polypeptide segments. The degree to which secondary structures depart from an average geometry also appears to depend on residue type. These departures are correlated to the corresponding Chou and Fasman parameters of amino acid residues. The frequency distribution of the side chain rotamers is influenced by polypeptide secondary structure. In turn, the rotamer conformation of side chain affects the extension of the secondary structure of the backbone. The strongest correlation is found between the occurrence of g+ conformation and helix propagation on the carboxyl side of many residues.

Post-transcriptional Repair of a Split Heat Shock Protein 90 Gene by mRNA trans-Splicing

Heat shock protein 90 participates in diverse biological processes ranging from protein folding, cell cycle, signal transduction and development to evolution in all eukaryotes. It is also critically involved in regulating growth of protozoa such as Dictyostelium discoideum, Leishmania donovani, Plasmodium falciparum, Trypanosoma cruzi, and Trypanosoma evansi. Selective inhibition of Hsp90 has also been explored as an intervention strategy against important human diseases such as cancer, malaria, or trypanosomiasis. Giardia lamblia, a simple protozoan parasite of humans and animals, is an important cause of diarrheal disease with significant morbidity and some mortality in tropical countries. Here we show that the G. lamblia cytosolic hsp90 ( glhsp90) is split in two similar sized fragments located 777 kb apart on the same scaffold. Intrigued by this unique arrangement, which appears to be specific for the Giardiinae, we have investigated the biosynthesis of GlHsp90. We used genome sequencing to confirm the split nature of the giardial hsp90. However, a specific antibody raised against the peptide detected a product with a mass of about 80 kDa, suggesting a post-transcriptional rescue of the genomic defect. We show evidence for the joining of the two independent Hsp90 transcripts in-trans to one long mature mRNA presumably by RNA splicing. The splicing junction carries hallmarks of classical cis-spliced introns, suggesting that the regular cis-splicing machinery may be sufficient for repair of the open reading frame. A complementary 26-nt sequence in the ``intron'' regions adjacent to the splice sites may assist in positioning the two pre-mRNAs for processing. This is the first example of post-transcriptional rescue of a split gene by trans-splicing.

Cloning and sequence of a cDNA encoding a novel hybrid prolme-rich protein associated with cytokinin-induced haustoria formation in Cuscuta reflexa

A complete cDNA encoding a novel hybrid Pro-rich protein (HyPRP) was identified by differentially screening 3x10(4) recombinant plaques of a Cuscuta reflexa cytokinin-induced haustorial cDNA library constructed in lambda gt10. The nucleotide (nt) sequence consists of: (i) a 424-bp 5'-non coding region having five start codons (ATGs) and three upstream open reading frames (uORFs); (ii) an ORF of 987 bp with coding potential for a 329-amino-acid (aa) protein of M(r), 35203 with a hydrophobic N-terminal region including a stretch of nine consecutive Phe followed by a Pro-rich sequence and a Cys-rich hydrophobic C terminus; and (iii) a 178-bp 3'-UTR (untranslated region). Comparison of the predicted aa sequence with the NBRF and SWISSPROT databases and with a recent report of an embryo-specific protein of maize [Jose-Estanyol et al., Plant Cell 4 (1992) 413-423] showed it to be similar to the class of HyPRPs encoded by genes preferentially expressed in young tomato fruits, maize embryos and in vitro-cultured carrot embryos. Northern analysis revealed an approx. 1.8-kb mRNA of this gene expressed in the subapical region of the C. reflexa vine which exhibited maximum sensitivity to cytokinin in haustorial induction.

Interaction of ecoP15I DNA methyltransferase with oligonucleotides containing the asymmetric sequence 5?-CAGCAG-3

EcoP15I DNA methyltransferase (Mtase) recognizes the asymmeteric sequence CAGCAG and catalyzes the transfer of a methyl group from S-adenosyl-L-methionine to the second adenine residue. We have investigated the DNA binding properties of EcoP15I DNA Mtase using gel mobility shift assays. EcoP15I DNA Mtase binds approximately threefold more tightly to DNA containing its recognition sequence, CAGCAG, than to non-specific sequences in the absence or presence of cofactors. Interestingly, in the presence of ATP the discrimination between specific and non-specific sequences increases significantly. These results suggest for the first time a role for ATP in DNA recognition by type III restriction-modification enzymes. In addition, we have shown that bromodeoxyuridine-containing oligonucleotides form complexes with EcoP15I DNA Mtase that are crosslinked upon irradiation. More importantly, we have shown that the crosslink site is at the site of DNA binding, since it can be suppressed by an excess of unmodified oligonucleotide. EcoP15I DNA Mtase exhibited Michaelis-Menten kinetics with both unmodified and bromodeoxyuridine-substituted DNA, with a higher specificity constant for the latter. Furthermore, gel mobility shift assays showed that proteolyzed EcoP15I DNA Mtase formed a specific complex with DNA, which had similar mobility as the native protein-DNA complex. Taken together these results form the basis fora detailed structure-function analysis of EcoP15I DNA Mtase.

Testis-specific histone (H1t) is not phosphorylated and has a weak interaction with chromatin

Soluble chromatin was prepared from rat testes after a brief micrococcal nuclease digestion. After adsorption onto hydroxylapatite at low ionic strength, the histone Hl subtypes were eluted with a shallow salt gradient of 0.3 M NaCl to 0.7 M NaCl. Histone Hlt was eluted at 0.4 MNaCl, while histones H1a and Hlc were eluted at 0.43 M NaCl and 0.45 M respectively. The extreme divergence of the amino acid sequence of the C-terminal half of histone Hlt, the major DNA binding domain of histone Hl, from that of the somatic consensus sequence may contribute to the weaker interaction of histone Hlt with the rat testis chromatin. Further, histone Hlt was not phosphorylated in vivo in contrast to histone Hla and Hlc, as is evident from the observation that histone Hlt lacks the SPKK motif recognized by the CDC-2kinase or the RR/KXS motif recognized by protein kinase A.

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.

Analysis of a conformation-specific epitope of the alpha subunit of human chorionic gonadotropin: Study using monoclonal antibody probes

Monoclonal antibodies (MAbs) have been used extensively for identification of sequence-specific epitopes using either the ELISA or/and IRMA methods, However, attempts to use MAbs for identification of conformation-specific epitopes have been very few as they are considered very labile. We have investigated the stability of conformation-specific epitopes of human chorionic gonadotropin (hCG) using a quantitative solid-phase radioimmnunoassay (SPRIA) technique. Several epitopes are stable to mild modification (chemical and proteolytic) conditions, and epitopes show differential stability for these modifications. Based on these observations, a monoclonal antibody (MAb 16) for an a-subunit-specific epitope of hCG has been used to monitor changes at the epitopic site (identified as epitope 16) on modification of hCG, using SPRIA with immobilized MAb 16. Modifications of amino groups, hydroxyl group of tyrosine as well as carboxyl group of Asp/Glu all bring about sufficient changes in the epitope integrity. Peptide bond hydrolysis at lysine residues damages the epitope, but not at arginine residues, Hydrolysis at tyrosine does not affect the epitope, though modification of the side-chain of tyrosine inactivates the epitope. Destruction of the epitope occurs on reduction of the disulphide bonds. Partial retention of the epitope activity is seen on modification of carboxyl or the epsilon-amino groups of lysine. Based on these results four to six amino acids have been identified to be at the epitopic site, and the data suggest that two peptide segments are brought together by the disulphide bond Cys10-Cys60 to form the epitope.

Developmental analysis of a female-specific 16S rRNA gene from mycetome-associated endosymbionts of a mealybug, Planococcus lilacinus

A clone showing female-specific expression was identified from an embryonic cDNA library of a mealybug, Planococcus lilacinus, In Southern blots this clone (P7) showed hybridization to genomic DNA of females, but not to that of males, However, P7 showed no hybridization to nuclei of either sex, raising the possibility that it was extrachromosomal in origin, In sectioned adult females P7 hybridized to an abdominal organ called the mycetome. The mycetome is formed by mycetocytes, which are polyploid cells originating from the polar bodies and cleavage nuclei that harbour maternally transmitted, intracellular symbionts. Electron microscopy confirmed the presence of symbionts within the mycetocytes, Sequence analysis showed that P7 is a 16S rRNA gene, confirming its prokaryotic origin, P7 transcripts are localized to one pole in young embryos but are found in the pole as well as in the germ band during later stages of development, P7 expression is detectable in young embryos of both sexes but the absence of P7 in third instar and adult males suggests that this gene, and hence the endosymbionts, are subject to sex-specific elimination. Copyright (C) 1997 Elsevier Science Ltd.