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.

Nucleotide-sequence of 5S RNA from mycobacterium-smegmatis

32P labelled 5S RNA isolated fromMycobacterium smegmatis was digested withT 1 and pancreatic ribonucleases separately and fingerprinted by two dimensional high voltage electrophoresis on thin-layer DEAE-cellulose plates. The radioactive spots were sequenced and their molar yields were determined. The chain length of the 5S RNA was found to be 120. It showed resemblances to both prokaryotic and eukaryotic 5S RNAs.

Primary structure of belladonna mottle virus coat protein.

The coat protein of belladonna mottle virus (a tymovirus) was cleaved by trypsin and chymotrypsin, and the peptides were separated by high performance liquid chromatography using a combination of gel permeation, reverse phase, and ion pair chromatography. The peptides were sequenced manually using the 4-N, N-dimethylaminoazobenzene-4'-isothiocyanate/phenyl isothiocyanate double-coupling method. The chymotryptic peptides were aligned by overlapping sequences of tryptic peptides and by homology with another tymovirus, eggplant mosaic virus. The belladonna mottle virus is more closely related to eggplant mosaic virus than to turnip yellow mosaic virus, the type member of this group, as evident from the sequence homologies of 57 and 32%, respectively. The accumulation of basic residues at the amino terminus implicated in RNA-protein interactions in many spherical plant viruses was absent in all the three sequences. Interestingly, the amino-terminal region is the least conserved among the tymoviruses. The longest stretch of conserved sequence between belladonna mottle virus and eggplant mosaic virus was residues 34-44, whereas it was residues 96-102 in the case of belladonna mottle virus and turnip yellow mosaic virus. A tetrapeptide in the region (residues 154-157) was found to be common for all the three sequences. It is possible that these conserved regions (residues 34-44, 96-102, 154-157) are involved in either intersubunit or RNA-protein interactions.

Sequence specificity in spermine-induced structural changes in CG-oligomers

The role of spermine in inducing A-DNA conformation in deoxyoligonucleotides has been studied using CCGG and GGCC as model sequences. It has been found that while CCGG adopts an alternating B-DNA conformation in low salt solution at low temperature, addition of spermine to this medium induces a B --greater than A transition. In contrast, the A-DNA-like structure of GGCC in low salt solution at low temperature does not change under the influence of spermine. This suggests a sequence-dependent behaviour of spermine. Further these results suggest that the A-DNA conformation observed in the crystals of d(iCCGG) and d(GGCC)2 might have been due to the presence of spermine in the crystallization cocktail.

Characterization of a mercury-reducing Bacillus cereus strain isolated from the Pulicat Lake sediments, south east coast of India

Pulicat Lake sediments are often severely polluted with the toxic heavy metal mercury. Several mercury-resistant strains of Bacillus species were isolated from the sediments and all the isolates exhibited broad spectrum resistance (resistance to both organic and inorganic mercuric compounds). Plasmid curing assay showed that all the isolated Bacillus strains carry chromosomally borne mercury resistance. Polymerase chain reaction and southern hybridization analyses using merA and merB3 gene primers/probes showed that five of the isolated Bacillus strains carry sequences similar to known merA and merB3 genes. Results of multiple sequence alignment revealed 99% similarity with merA and merB3 of TnMERI1 (class II transposons). Other mercury resistant Bacillus species lacking homology to these genes were not able to volatilize mercuric chloride, indicating the presence of other modes of resistance to mercuric compounds.

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.

Characterization of a mercury-reducing Bacillus cereus strain isolated from the Pulicat Lake sediments, south east coast of India

Pulicat Lake sediments are often severely polluted with the toxic heavy metal mercury. Several mercury-resistant strains of Bacillus species were isolated from the sediments and all the isolates exhibited broad spectrum resistance (resistance to both organic and inorganic mercuric compounds). Plasmid curing assay showed that all the isolated Bacillus strains carry chromosomally borne mercury resistance. Polymerase chain reaction and southern hybridization analyses using merA and merB3 gene primers/probes showed that five of the isolated Bacillus strains carry sequences similar to known merA and merB3 genes. Results of multiple sequence alignment revealed 99% similarity with merA and merB3 of TnMERI1 (class II transposons). Other mercury resistant Bacillus species lacking homology to these genes were not able to volatilize mercuric chloride, indicating the presence of other modes of resistance to mercuric compounds.

Estradiol - 17β induces polyaromatic hydrocarbon-inducible cytochrome P-450 in chicken liver

A cDNA clone has been isolated from a chicken liver library prepared against messenger RNA isolated after chronic estradiol-17β treatment. The clone, pP-450 IA - 61, has an insert of 900nt and the sequence shows high homology to CYPIA2 subfamily from four other species. A single injection of estradiol-17β to immature chicken results in a striking induction of mRNA hybridizing to labeled pP-450IA - 61. The probe also hybridizes to mRNA induced by 3 — methylcholanthrene in chicken. These results offer direct proof for the similarity in the mode of action at the transcriptional level of polyaromatic hydrocarbons and estrogenic compounds.

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.

Src Homology 3-interacting Domain of Rv1917c of Mycobacterium tuberculosis Induces Selective Maturation of Human Dendritic Cells by Regulating PI3K-MAPK-NF-kappa B Signaling and Drives Th2 Immune Responses

Mycobacterium tuberculosis, an etiological agent of pulmonary tuberculosis, causes significant morbidity and mortality worldwide. Pathogenic mycobacteria survive in the host by subverting host innate immunity. Dendritic cells (DCs) are professional antigen-presenting cells that are vital for eliciting immune responses to infectious agents, including pathogenic mycobacteria. DCs orchestrate distinct Th responses based on the signals they receive. In this perspective, deciphering the interactions of the proline-glutamic acid/proline-proline-glutamic acid (PE/PPE) family of proteins of M. tuberculosis with DCs assumes significant pathophysiological attributes. In this study, we demonstrate that Rv1917c (PPE34), a representative member of the proline-proline-glutamic-major polymorphic tandem repeat family, interacts with TLR2 and triggers functional maturation of human DCs. Signaling perturbations implicated a critical role for integrated cross-talk among PI3K-MAPK and NF-kappa B signaling cascades in Rv1917c-induced maturation of DCs. However, this maturation of DCs was associated with a secretion of high amounts of anti-inflammatory cytokine IL-10, whereas Th1-polarizing cytokine IL-12 was not induced. Consistent with these results, Rv1917c-matured DCs favored secretion of IL-4, IL-5, and IL-10 from CD4(+) T cells and contributed to Th2-skewed cytokine balance ex vivo in healthy individuals and in patients with pulmonary tuberculosis. Interestingly, the Rv1917c-skewed Th2 immune response involved induced expression of cyclooxygenase-2 (COX-2) in DCs. Taken together, these results indicate that Rv1917c facilitates a shift in the ensuing immunity toward the Th2 phenotype and could aid in immune evasion by mycobacteria.

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.

A hypothesis on a specific sequence-dependent conformation of DNA and its relation to the binding of the lac-repressor protein

The structure and properties of the double-helical form of the alternating copolymer poly(dA-dT) are considered. Different lines of evidence are interpreted in terms of a structure in which every second phosphate-diester linkage has a conformation different from that of the normal B form. A rationale for this “alternating-B” structure is given which provides an explanation for the effects of chemical modifications of the T residues on the binding of the poly(dA-dT)· poly(dA-dT) to the lac repressor of Escherichia coli.