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.

Characterization of human symptomatic rotavirus isolates MP409 and MP480 having ‘long’ RNA electropherotype and subgroup I specificity, highly related to the P6[1],G8 type bovine rotavirus A5, from Mysore, India

In an epidemiological study of symptomatic human rotaviruses in Mysore, India during 1993 and 1994, isolates MP409 and MP480 were isolated from two children suffering from severe, acute dehydrating diarrhea. Both isolates exhibited 'long' RNA pattern and subgroup I specificity suggesting the likelihood of their animal origin. Both isolates did not react with monoclonal antibodies (MAbs) specific for serotypes G1 to G6 as well as CIO. To determine the genetic origin of these isolates, complete nucleotide sequences of genes encoding the outer capsid proteins VP4 and VP7, nonstructural proteins NSP1 and NSP3 and viral enterotoxin protein NSP4 from MP409 and partial sequences of genes from MP480 were determined. Comparison of the 5' and 3' terminal sequences of 250 nucleotides revealed complete identity of the gene sequences in both strains suggesting that MP409 and MP480 are two different isolates of a single strain. Comparison of the nucleotide and deduced amino acid sequences of VP4, VP7, NSP1 and NSP3 of MP409 with published sequences of strains belonging to different serotypes revealed that both outer capsid proteins VP4 and VP7 and NSP1 are highly related to the respective proteins from the P6[1], G8 type bovine rotavirus A5 isolated from a calf with diarrhoea in Thailand and that the NSP3 is highly homologous to that of bovine rotaviruses. The NSP 1 protein showed greatest sequence identity with NSP4s belonging to the KUN genetic group to which NSP4s from human G2 type strains and bovine rotaviruses belong. MP409 and MP480 likely signify interspecies transmission of P6[1], G8 type strains from cattle to humans and represent the first P6[1] type rotaviruses isolated in humans. These and our previous studies on the asymptomatic neonatal strain I321 are of evolutionary and epidemiological significance in the context of close association of majority of the Indian population with cattle.

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.

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.

Interrelations between Dry Eye Syndrome and Tear Fluid Phospholipid Transfer Protein

The simplified model of human tear fluid (TF) is a three-layered structure composed of a homogenous gel-like layer of hydrated mucins, an aqueous phase, and a lipid-rich outermost layer found in the tear-air interface. It is assumed that amphiphilic phospholipids are found adjacent to the aqueous-mucin layer and externally to this a layer composed of non-polar lipids face the tear-air interface. The lipid layer prevents evaporation of the TF and protects the eye, but excess accumulation of lipids may lead to drying of the corneal epithelium. Thus the lipid layer must be controlled and maintained by some molecular mechanisms. In the circulation, phospholipid transfer protein (PLTP) and cholesteryl ester transfer protein (CETP) mediate lipid transfers. The aim of this thesis was to investigate the presence and molecular mechanisms of lipid transfer proteins in human TF. The purpose was also to study the role of these proteins in the development of dry eye syndrome (DES). The presence of TF PLTP and CETP was studied by western blotting and mass spectrometry. The concentration of these proteins was determined by ELISA. The activities of the enzymes were determined by specific lipid transfer assays. To study the molecular mechanisms involved in PLTP mediated lipid transfer Langmuir monolayers and asymmetrical flow field-flow fractionation (AsFlFFF) was used. Ocular tissue samples were stained with monoclonal antibodies against PLTP to study the secretion route of PLTP. Heparin-Sepharose affinity chromatography was used for PLTP pull-down experiments and co-eluted proteins were identified with MALDI-TOF mass spectrometry or Western blot analysis. To study whether PLTP plays any functional role in TF PLTP-deficient mice were examined. The activity of PLTP was also studied in dry eye patients. PLTP is a component of normal human TF, whereas CETP is not. TF PLTP concentration was about 2-fold higher than that in human plasma. Inactivation of PLTP by heat treatment or immunoinhibition abolished the phospholipid transfer activity in tear fluid. PLTP was found to be secreted from lacrimal glands. PLTP seems to be surface active and is capable of accepting lipid molecules without the presence of lipid-protein complexes. The active movement of radioactively labeled lipids and high activity form of PLTP to acceptor particles suggested a shuttle model of PLTP-mediated lipid transfer. In this model, PLTP physically transports lipids between the donor and acceptor. Protein-protein interaction assays revealed ocular mucins as PLTP interaction partners in TF. In mice with a full deficiency of functional PLTP enhanced corneal epithelial damage, increased corneal permeability to carboxyfluorescein, and decreased corneal epithelial occludin expression was demonstrated. Increased tear fluid PLTP activity was observed among human DES patients. These results together suggest a scavenger property of TF PLTP: if the corneal epithelium is contaminated by hydrophobic material, PLTP could remove them and transport them to the superficial layer of the TF or, alternatively, transport them through the naso-lacrimal duct. Thus, PLTP might play an integral role in tear lipid trafficking and in the protection of the corneal epithelium. The increased PLTP activity in human DES patients suggests an ocular surface protective role for this lipid transfer protein.