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.

Structural studies of model peptides containing beta-, gamma- and delta-amino acids

The crystal structures of five model peptides Piv-Pro-Gly-NHMe (1), Piv-Pro-beta Gly-NHMe (2), Piv-Pro-beta Gly-OMe (3), Piv-Pro-delta Ava-OMe (4) and Boc-Pro-gamma Abu-OH (5) are described (Piv:pivaloyl; NHMe: N-methylamide; beta Gly:beta-glycine; OMe:O-methyl ester; delta Ava:delta-aminovaleric acid; gamma Abu:gamma-aminobutyric acid). A comparison of the structures of peptides 1 and 2 illustrates the dramatic consequences upon backbone homologation in short sequences. 1 adopts a type II beta-turn conformation in the solid state, while in 2, the molecule adopts an open conformation with the beta-residue being fully extended. Piv-Pro-beta Gly-OMe (3), which differs from 2 by replacement of the C-terminal NH group by an O-atom, adopts an almost identical molecular conformation and packing arrangement in the solid state. In peptide 4, the observed conformation resembles that determined for 2 and 3, with the delta Ava residue being fully extended. In peptide 5, the molecule undergoes a chain reversal, revealing a beta-turn mimetic structure stabilized by a C-H center dot center dot center dot O hydrogen bond.

Ferroelectricity in Bi26-xMxO40-delta (M = Al and Ga) with the gamma-Bi2O3 structure

We report on the dielectric proper-ties of bismuth aluminate and gallate with Bi:AI(Ga) ratio of 1: 1 and 12:1 prepared at high temperature and ambient pressure. These compounds crystallize in a noncentrosymmetric body-centered cubic structure (space group 123) with a similar to 10.18 angstrom rather than in the perovskite structure.This cubic phase is related to the gamma-Bi2O3 structure which has the actual chemical formula Bi-24(3+) (Bi3+Bi5+)O40-delta. In the aluminates and gallates studied by us, the Al and Ga ions are distributed over the 24f and 2a sites. These compounds exibit ferroclectric hysteresis at room temperature with a weak polarization. (c) 2006 Elsevier Ltd. All rights reserved.

Computer modelling approach to study the modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein

The modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein (ABP) have been studied by energy minimization using the low resolution (2.4 A) X-ray data of the protein. These studies suggest that these sugars preferentially bind in the alpha-form to ABP, unlike L-arabinose where both alpha- and beta-anomers bind almost equally. The best modes of binding of alpha- and beta-anomers of D-galactose and D-fucose differ slightly in the nature of the possible hydrogen bonds with the protein. The residues Arg 151 and Asn 232 of ABP from bidentate hydrogen bonds with both L-arabinose and D-galactose, but not with D-fucose or D-glucose. However in the case of L-arabinose, Arg 151 forms hydrogen bonds with the hydroxyl group at the C-4 atom and the ring oxygen, whereas in case of D-galactose it forms bonds with the hydroxyl groups at the C-4 and C-6 atoms of the pyranose ring. The calculated conformational energies also predict that D-galactose is a better inhibitor than D-fucose and D-glucose, in agreement with kinetic studies. The weak inhibitor D-glucose binds preferentially to one domain of ABP leading to the formation of a weaker complex. Thus these studies provide information about the most probable binding modes of these sugars and also provide a theoretical explanation for the observed differences in their binding affinities.

Preparation and characterization of monoclonal antibodies to nucleocapsid protein N and H glycoprotein of rinderpest virus

Fifteen stable mouse spleen cell myeloma hybrids (hybridomas) producing monoclonal antibodies to rinderpest virus proteins were produced. The specificity of these monoclonal antibodies was established by radioimmunoprecipitation followed by polyacrylamide gel analysis and immunofluorescence. Nine antibodies were specific for the surface glycoprotein H. All the nine clones showed inhibition of haemagglutination by measles virus. The antibodies from two clones (A7D2 and B2F6) neutralise infectious virus. Six clones produce antibodies reacting with the nucleocapsid protein N. Three antigenic sites designated I–III, with sites I and II partially overlapping, were topographically mapped on the H molecule by competitive binding assay. Similarly, two antigenic sites I and II were delineated on the N protein. The monoclonal antibodies were used to study the antigenic relationships of H and N proteins of rinderpest virus, measles virus and canine distemper virus.

Computer simulation of protein—carbohydrate complexes: application to arabinose-binding protein and pea lectin

The CCEM method (Contact Criteria and Energy Minimisation) has been developed and applied to study protein-carbohydrate interactions. The method uses available X-ray data even on the native protein at low resolution (above 2.4 Å) to generate realistic models of a variety of proteins with various ligands.The two examples discussed in this paper are arabinose-binding protein (ABP) and pea lectin. The X-ray crystal structure data reported on ABP-β-l-arabinose complex at 2.8, 2.4 and 1.7 Å resolution differ drastically in predicting the nature of the interactions between the protein and ligand. It is shown that, using the data at 2.4 Å resolution, the CCEM method generates complexes which are as good as the higher (1.7 Å) resolution data. The CCEM method predicts some of the important hydrogen bonds between the ligand and the protein which are missing in the interpretation of the X-ray data at 2.4 Å resolution. The theoretically predicted hydrogen bonds are in good agreement with those reported at 1.7 Å resolution. Pea lectin has been solved only in the native form at 3 Å resolution. Application of the CCEM method also enables us to generate complexes of pea lectin with methyl-α-d-glucopyranoside and methyl-2,3-dimethyl-α-d-glucopyranoside which explain well the available experimental data in solution.

cAMP-regulated Protein Lysine Acetylases in Mycobacteria

Cyclic AMP synthesized by Mycobacterium tuberculosis has been shown to play a role in pathogenesis. However, the high levels of intracellularcAMP found in both pathogenic and nonpathogenic mycobacteria suggest that additional and important biological processes are regulated by characterization of novel cAMP-binding proteins in M. smegmatis and M. tuberculosis (MSMEG_5458 and Rv0998, respectively) that contain a cyclic nucleotide binding domain fused to a domain that shows similarity to the GNAT family of acetyltransferases. We detect protein lysine acetylation in mycobacteria and identify a universal stress protein (USP) as a substrate of MSMEG_5458. Acetylation of a lysine residue in USP is regulated by cAMP, and using a strain deleted for MSMEG_5458, we show that USP is indeed an in vivo substrate for MSMEG_5458. The Rv0998 protein shows a strict cAMP-dependent acetylation of USP, despite a lower affinity for cAMP than MSMEG_5458. Thus, this report not only represents the first demonstration of protein lysine acetylation in mycobacteria but also describes a unique functional interplay between a cyclic nucleotide binding domain and a protein acetyltransferase.

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.

Dynamic changes in mitogen-activated protein kinase (MAPK) activities in the corpus luteum of the bonnet monkey (Macaca radiata) during development, induced luteolysis, and simulated early pregnancy: A role for p38 MAPK in the regulation of luteal function

Changes in MAPK activities were examined in the corpus luteum (CL) during luteolysis and pregnancy, employing GnRH antagonist (Cetrorelix)-induced luteolysis, stages of CL, and hCG treatment to mimic early pregnancy as model systems in the bonnet monkey. We hypothesized that MAPKs could serve to phosphorylate critical phosphoproteins to regulate luteal function. Analysis of several indices for structural (caspase-3 activity and DNA fragmentation) and functional (progesterone and steroidogenic acute regulatory protein expression) changes in the CL revealed that the decreased luteal function observed during Cetrorelix treatment and late luteal phase was associated with increased caspase-3 activity and DNA fragmentation. As expected, human chorionic gonadotropin treatment dramatically increased luteal function, but the indices for structural changes were only partially attenuated. All three MAPKs appeared to be constitutively active in the mid-luteal-phase CL, and activities of ERK-1/2 and p38-MAPK (p38), but not Jun N-terminal kinase (JNK)-1/2, decreased significantly (P < 0.05) within 12 - 24 h after Cetrorelix treatment. During the late luteal phase, in contrast to decreased ERK-1/2 and p38 activities, JNK-1/2 activities increased significantly (P < 0.05). Although human chorionic gonadotropin treatment increased ERK-1/2 and p38 activities, it decreased JNK-1/2 activities. The activation status of p38 was correlated with the phosphorylation status of an upstream activator, MAPK kinase-3/6 and the expression of MAPK activated protein kinase-3, a downstream target. Intraluteal administration of p38 kinase inhibitor (SB203580), but not MAPK kinase-1/2 inhibitor (PD98059), decreased the luteal function. Together, these data suggest an important role for p38 in the regulation of CL function in primates.

Transcriptional analysis of persistent Chlamydia pneumoniae infection in vitro

Chlamydia pneumoniae can cause acute respiratory infections including pneumonia. Repeated and persistent Chlamydia infections occur and persistent C. pneumoniae infection may have a role in the pathogenesis of atherosclerosis and coronary heart disease and may also contribute to the development of chronic inflammatory lung diseases like chronic obstructive pulmonary disease (COPD) and asthma. In this thesis in vitro models for persistent C. pneumonia infection were established in epithelial and monocyte/macrophage cell lines. Expression of host cell genes in the persistent C. pneumoniae infection model of epithelial cells was studied by microarray and RT-PCR. In the monocyte/macrophage infection model expression of selected C. pneumoniae genes were studied by RT-PCR and immunofluorescence microscopy. Chlamydia is able to modulate host cell gene expression and apoptosis of host cells, which may assist Chlamydia to evade the host cells' immune responses. This, in turn, may lead to extended survival of the organism inside epithelial cells and promote the development of persistent infection. To simulate persistent C. pneumoniae infection in vivo, we set up a persistent infection model exposing the HL cell cultures to IFN-gamma. When HL cell cultures were treated with moderate concentration of IFN-gamma, the replication of C. pneumoniae DNA was unaffected while differentiation into infectious elementary bodies (EB) was strongly inhibited. By transmission electron microscopy small atypical inclusions were identified in IFN-gamma treated cultures. No second cycle of infection was observed in cells exposed to IFN-gamma , whereas C. pneumoniae was able to undergo a second cycle of infection in unexposed HL cells. Although monocytic cells can naturally restrict chlamydial growth, IFN-gamma further reduced production of infectious C. pneumoniae in Mono Mac 6 cells. Under both studied conditions no second cycle of infection could be detected in monocytic cell line suggesting persistent infection in these cells. As a step toward understanding the role of host genes in the development and pathogenesis of persistent C. pneumoniae infection, modulation of host cell gene expression during IFN-gamma induced persistent infection was examined and compared to that seen during active C. pneumoniae infection or IFN-gamma treatment. Total RNA was collected at 6 to 150 h after infection of an epithelial cell line (HL) and analyzed by a cDNA array (available at that time) representing approximately 4000 human transcripts. In initial analysis 250 of the 4000 genes were identified as differentially expressed upon active and persistent chlamydial infection and IFN-gamma treatment. In persistent infection more potent up-regulation of many genes was observed in IFN-gamma induced persistent infection than in active infection or in IFN-gamma treated cell cultures. Also sustained up-regulation was observed for some genes. In addition, we could identify nine host cell genes whose transcription was specifically altered during the IFN-gamma induced persistent C. pneumoniae infection. Strongest up-regulation in persistent infection in relation to controls was identified for insulin like growth factor binding protein 6, interferon-stimulated protein 15 kDa, cyclin D1 and interleukin 7 receptor. These results suggest that during persistent infection, C. pneumoniae reprograms the host transcriptional machinery regulating a variety of cellular processes including adhesion, cell cycle regulation, growth and inflammatory response, all of which may play important roles in the pathogenesis of persistent C. pneumoniae infection. C. pneumoniae DNA can be detected in peripheral blood mononuclear cells indicating that the bacterium can also infect monocytic cells in vivo and thereby monocytes can assist the spread of infection from the lungs to other anatomical sites. Persistent infection established at these sites could promote inflammation and enhance pathology. Thus, the mononuclear cells are in a strategic position in the development of persistent infection. To investigate the intracellular replication and fate of C. pneumoniae in mononuclear cells we analyzed the transcription of 11 C. pneumoniae genes in Mono Mac 6 cells during infection by real time RT-PCR. Our results suggest that the transcriptional profile of the studied genes in monocytes is different from that seen in epithelial cells and that IFN-gamma has a less significant effect on C. pneumoniae transcription in monocytes. Furthermore, our study shows that type III secretion system (T3SS) related genes are transcribed and that Chlamydia possesses a functional T3SS during infection in monocytes. Since C. pneumoniae infection in monocytes has been implicated to have reduced antibiotic susceptibility, this creates opportunities for novel therapeutics targeting T3SS in the management of chlamydial infection in monocytes.

Glioblastoma-Specific Protein Interaction Network Identifies PP1A and CSK21 as Connecting Molecules between Cell Cycle-Associated Genes

Glioblastoma (GBM; grade IV astrocytoma) is a very aggressive form of brain cancer with a poor survival and few qualified predictive markers. This study integrates experimentally validated genes that showed specific upregulation in GBM along with their protein-protein interaction information. A system level analysis was used to construct GBM-specific network. Computation of topological parameters of networks showed scale-free pattern and hierarchical organization. From the large network involving 1,447 proteins, we synthesized subnetworks and annotated them with highly enriched biological processes. A careful dissection of the functional modules, important nodes, and their connections identified two novel intermediary molecules CSK21 and protein phosphatase 1 alpha (PP1A) connecting the two subnetworks CDC2-PTEN-TOP2A-CAV1-P53 and CDC2-CAV1-RB-P53-PTEN, respectively. Real-time quantitative reverse transcription-PCR analysis revealed CSK21 to be moderately upregulated and PP1A to be overexpressed by 20-fold in GBM tumor samples. Immunohistochemical staining revealed nuclear expression of PP1A only in GBM samples. Thus, CSK21 and PP1A, whose functions are intimately associated with cell cycle regulation, might play key role in gliomagenesis. Cancer Res; 70(16); 6437-47. (C)2010 AACR.