Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to a-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the ``structurally variable'' regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of `variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.

Primary structure of sesbania mosaic virus coat protein: its implications to the assembly and architecture of the virus

Sesbania mosaic virus (SMV) is a plant virus that infects Sesbania grandiflora plants in Andhra Pradesh, India. The amino acid sequence of the coat protein of SMV was determined using purified peptides generated by cleavage with trypsin, chymotrypsin, V8 protease and clostripain. The 230 residues so far determined were compared to the corresponding residues of southern bean mosaic virus (SBMV), the type member of sobemoviruses. The overall identity between the sequences is 61.7%. The amino terminal 64 residues, which constitute an independent domain (R-domain) known to interact with RNA, are conserved to a lower extent (52.5%). Comparison of the positively charged residues in this domain suggests that the RNA-protein interactions are considerably weaker in SMV. The residues that constitute the major domain of the coat protein, the surface domain (S-domain, residues 65-260), are better conserved (66.5%). The positively charged residues of this domain that face the nucleic acid are well conserved. The longest conserved stretch of residues (131-142) corresponds to the loop involved in intersubunit interactions between subunits related by the quasi 3-fold symmetry. A unique cation binding site located on the quasi 3-fold axis contributes to the stability of SMV. These differences are reflected in the increased stability of the SMV coat protein and its ability to be reconstituted with RNA at pH 7.5. A major epitope was identified using monoclonal antibodies to SMV in the segment 201-223 which contains an exposed helix in the capsid structure. This region is highly conserved between SMV and SBMV (70%) suggesting that it could represent the site of an important function such as vector recognition.

Characterization of aerosol black carbon over a tropical semi-arid region of Anantapur, India

Black carbon (BC) aerosol mass concentrations measured using an aethalometer at Anantapur, a semi-arid tropical station in the southern part of peninsular India, from August 2006 to July 2007 are analyzed. Seasonal and diurnal variations of BC in relation to changes in the regional meteorological conditions have been studied along with the mass fraction of BC to the total aerosol mass concentration (M-t) and fine particle mass (FPM) concentration in different months. The data collected during the study period shows that the annual average BC mass concentration at Anantapur is 1.97 +/- 0.12 mu g m(-3). Seasonal variations of BC aerosol mass concentration showed high during the dry (winter and summer) seasons and low during the post-monsoon followed by the monsoon seasons. Diurnal variations of BC aerosols attain a gradual build up in BC concentration from morning and a sharp peak occurs between 07:00 and 09:00 h almost an hour after local sunrise and a broad nocturnal peak from 19:00 to 21:00 h with a minimum in noon hours. The ratio of BC to the fine particle mass concentration was high during the dry season and low during the monsoon season. The regression analysis between BC mass concentration and wind speed indicates that, with increase in wind speeds the BC mass concentrations would decrease and vice-versa. Aerosol BC mass concentration shows a significant positive correlation with total mass concentration (M-t) and aerosol optical depth (ACID, tau(p)) at 500 nm. (C) 2010 Elsevier B.V. All rights reserved.

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.

A novel method for mapping assembled epitopes in batches: Identification of three epitopes at the receptor binding region of human chorionic gonadotropin

Identification of conformation-specific epitopes of hCG beta has been done using a simple batch method, Chemically or enzymatically-modified hCG beta has been prepared in a batch and the effect of modifications on the integrity of different epitope regions has been investigated in a quantitative manner using monoclonal antibodies (MAbs) immobilized on plastic tubes from culture supernatants. Based on the extent of damage done to different regions by different modifications, three conformation-specific epitopes of hCG beta have been identified. The method has been shown to have important advantages over the existing methods on many considerations, Using this approach, these epitopes have been shown to be at/near the receptor-binding region.

Electron transfer reaction in the marcus inverted region: role of high frequency vibrational modes

A theoretical study of the dynamics of photo-electron transfer reactions in the Marcus inverted regime is presented. This study is motivated partly by the recent proposal of Barbara et al. (J. Phys. Chem. 96, 3728, 1991) that a minimal model of an electron transfer reaction should consist of a polar solvent mode (X), a low frequency vibrational mode (Q) and one high frequency mode (q). Interplay between these modes may be responsible for the crossover observed in the dynamics from a solvent controlled to a vibrational controlled electron transfer. The following results have been obtained. (i) In the case of slowly relaxing solvents, the proximity of the point of excitation to an effective sink on the excited surface is critical in determining the decay of the reactant population. This is because the Franck-Condon overlap between the reactant ground and the product excited states decreases rapidly with increase in the quantum number of the product vibrational state. (ii) Non-exponential solvation dynamics has an important effect in determining the rates of electron transfer. Especially, a biphasic solvation and a large coupling between the reactant and the product states both may be needed to explain the experimental results. ©1996 American Institute of Physics

Functional analysis of conserved motifs in EcoP15I DNA methyltransferase

EcoP15I DNA methyltransferase recognizes the sequence 5'-CAGCAG-3' and transfers a methyl group to N-6 of the second adenine residue in the recognition sequence. All N-6 adenine methyltransferases contain two highly conserved sequences, FxGxG (motif I), postulated to form part of the S-adenosyl-L-methionine binding site and (D/N/S)PP(Y/F) (motif IV) involved in catalysis. We have altered the second glycine residue in motif I to arginine and serine, and substituted tyrosine in motif IV with tryptophan in EcoP15I DNA methyltransferase, using site-directed mutagenesis. The mutant enzymes were overexpressed, purified and characterized by biochemical methods. The mutations in motif I completely abolished AdoMet binding but left target DNA recognition unaltered. Although the mutation in motif IV resulted in loss of enzyme activity, we observed enhanced crosslinking of S-adenosyl-L-methionine and DNA. This implies that DNA and AdoMet binding sites are close to motif IV. Taken together, these results reinforce the importance of motif I in AdoMet binding and motif IV in catalysis. Additionally, limited proteolysis and UV crosslinking experiments with EcoP15I DNA methyltransferase imply that DNA binds in a cleft formed by two domains in the protein. Methylation protection analysis provides evidence for the fact that EcoP15I DNA MTase makes contacts in the major groove of its substrate DNA. Interestingly, hypermethylation of the guanine residue next to the target adenine residue indicates that the protein probably flips out the target adenine residue. (C) 1996 Academic Press Limited

Comparative Nucleotide and Amino Acid Sequence Analysis of the Sequence-Specific RNA-Binding Rotavirus Nonstructural Protein NSP3

NSP3, an acidic nonstructural protein, encoded by gene 7 has been implicated as the key player in the assembly of the 11 viral plus-strand RNAs into the early replication intermediates during rotavirus morphogenesis. To date, the sequence or NSP3 from only three animal rotaviruses (SA11, SA114F, and bovine UK) has been determined and that from a human strain has not been reported. To determine the genetic diversity among gene 7 alleles from group A rotaviruses, the nucleotide sequence of the NSP3 gene from 13 strains belonging to nine different G serotypes, from both humans and animals, has been determined. Based on the amino acid sequence identity as well as phylogenetic analysis, NSP3 from group A rotaviruses falls into three evolutionarily related groups, i.e., the SA11 group, the Wa group, and the S2 group. The SA 11/SA114F gene appears to have a distant ancestral origin from that of the others and codes for a polypeptide of 315 amino acids (aa) in length. NSP3 from all other group A rotaviruses is only 313 aa in length because of a 2-amino-acid deletion near the carboxy-terminus, While the SA114F gene has the longest 3' untranslated region (UTR) of 132 nucleotides, that from other strains suffered deletions of varying lengths at two positions downstream of the translational termination codon. In spite of the divergence of the nucleotide (nt) sequence in the protein coding region, a stretch of about 80 nt in the 3' UTR is highly conserved in the NSP3 gene from all the strains. This conserved sequence in the 3' UTR might play an important role in the regulation of expression of the NSP3 gene. (C) 1995 Academic Press, Inc.

Quantum and non-Markovian effects in the electron transfer reaction dynamics in the Marcus inverted region

Electron transfer reactions in large molecules may often be coupled to both the polar solvent modes and the intramolecular vibrational modes of the molecule. This can give rise to a complex dynamics which may in some systems, like betaine, be controlled more by vibrational rather than by solvent effects. Additionally, a significant contribution from an ultrafast relaxation component in the solvation dynamics may enhance the complexity. To explain the wide range of behavior that has been observed experimentally, Barbara et al. recently proposed that a model of an electron transfer reaction should minimally consist of a low-frequency classical solvent mode (X), a low-frequency vibrational mode (Q), and a high-frequency quantum mode (q) (J. Phys. Chem. 1991, 96, 3728). In the present work, a theoretical study of this model is described. This study generalizes earlier work by including the biphasic solvent response and the dynamics of the low-frequency vibrational mode in the presence of a delocalized, extended reaction zone. A novel Green's function technique has been developed which allowed us to study the non-Markovian dynamics on a multidimensional surface. The contributions from the high-frequency vibrational mode and the ultrafast component in the non-Markovian solvent dynamics are found to be primarily responsible for the dramatic increase in charge transfer rate over the prediction of the classical theories that neglect both these factors. These, along with a large coupling between the reactant and the product states, may combine to render the electron transfer rate both very large and constant over a wide range of solvent relaxation rates. A study on the free energy gap dependence of the electron transfer rate reveals that the rates are sensitive to changes in the quantum frequency particularly when the free energy gap is very large.

Status Survey of the Malabar Pied Hornbill in the Dandeli Region, Northern Western Ghats, India

The Malabar Pied Hornbill, Anthracoceros coronatus, is a near threatened species, endemic to the tropical deciduous forests of central and southern India and Sri Lanka. The Dandeli region in Karnataka (India) is believed to be the last stronghold of this species in the Western Ghats biodiversity hotspot. Being a rapidly developing area with a growing human population, the threats to this species and their habitat are mounting, especially due to a large number of hydroelectric projects and habitat fragmentation caused by paper and plywood industries. This study evaluated the change in population status of the Malabar Pied Hornbill over a 23 year period and defined priorities for the long term conservation and monitoring of hornbills in Dandeli. Encounter rates of hornbills were also analysed in relation to the density and species richness of trees and fruiting trees, basal area, canopy cover and distance from river. Hornbill encounters were not significantly different compared to the earlier study carried out by Reddy in 1988, but were significantly different across the five sites in the current study. Higher numbers of hornbills were encountered closer to the river, but these results were only marginally significant. The mean numbers of hornbills recorded at the two roost sites identified in Dandeli were 26 +/- 4.47 (n=16 counts) and 31.78 +/- 3.53 (n=14 counts) respectively. The study also helped build local awareness about the species, train local Forest Department staff in monitoring hornbills and develop a management plan for its conservation.

Implications of the recent high statistics determination of the pion electromagnetic form factor in the timelike region

The recently evaluated two-pion contribution to the muon g - 2 and the phase of the pion electromagnetic form factor in the elastic region, known from pi pi scattering by Fermi-Watson theorem, are exploited by analytic techniques for finding correlations between the coefficients of the Taylor expansion at t = 0 and the values of the form factor at several points in the spacelike region. We do not use specific parametrizations, and the results are fully independent of the unknown phase in the inelastic region. Using for instance, from recent determinations, < r(pi)(2)> = (0.435 +/- 0.005) fm(2) and F(-1.6 GeV2) = 0.243(-0.014)(+0.022), we obtain the allowed ranges 3.75 GeV-4 less than or similar to c less than or similar to 3.98 GeV-4 and 9.91 GeV-6 less than or similar to d less than or similar to 10.46 GeV-6 for the curvature and the next Taylor coefficient, with a strong correlation between them. We also predict a large region in the complex plane where the form factor cannot have zeros.

Turbulent heat flux variation over the Monsoon-Trough region during MONTBLEX-90

Tower data collected during the Monsoon-Trough Boundary Layer Experiment (MONTBLEX-90) have been analysed to understand the observed structure of the surface layer over an arid region (Jodhpur) and a moist region (Kharagpur) during active and weak phases of the 1990 southwest monsoon. Turbulent heat and momentum fluxes are estimated by the eddy correlation method using sonic data. The turbulent momentum flux at both Jodhpur and Kharagpur was larger when the winds were stronger, reaching a maximum of the order of 0.5 N m(-2) on 5 and 6 August when a low pressure system was located over the region. The heat flux at Jodhpur is high during weak monsoon days, the maximum being 450 W m(-2), whereas during active days the flux never exceeds 200 W m(-2). At Kharagpur, the flux does not vary significantly between active and weak monsoon days, the maximum in either phase being 160 W m(-2) At Jodhpur, there is significant contrast in the near-surface air temperature, being higher during weak monsoon days as compared to active days. Cloud cover did not vary significantly in both the regions. The turbulent heat flux variation at both the sites appears to be correlated mainly with soil mixture, and less sensitive to cloud cover.

Isolation of an 800 kb contiguous DNA fragment encompassing a 3.5-cM region of chromosome 1 in Arabidopsis using YAC clones

The apetalal mutation of Arabidopsis affects floral meristem identity and the development of sepal and petal primordia of the flower. We mapped the available RFLP markers on chromosome 1 that are in the general vicinity of apetalal on a fine structure map and then chose the closest RFLP as a starting point for contiguous DNA (contig) generation. We report here a contig of about 800 kilobases (kb) that spans a 3.5 cM region of chromosome 1. We used genomic libraries of Arabidopsis prepared in yeast artificial chromosome (YAC) vectors and the detailed characterization of 19 YACs is reported. RFLPs displayed by the end fragments from the walk were mapped to align and correlate the genetic and physical maps for this region of chromosome 1. In this segment of the genome, 1 cM corresponds to a little over 200 kb of physical distance.

Probing the role of the fully conserved Cys126 in triosephosphate isomerase by site-specific mutagenesis - distal effects on dimer stability

Cys126 is a completely conserved residue in triosephosphate isomerase that is proximal to the active site but has been ascribed no specific role in catalysis. A previous study of the C126S and C126A mutants of yeast TIM reported substantial catalytic activity for the mutant enzymes, leading to the suggestion that this residue is implicated in folding and stability [Gonzalez-Mondragon E et al. (2004) Biochemistry43, 3255–3263]. We re-examined the role of Cys126 with the Plasmodium falciparum enzyme as a model. Five mutants, C126S, C126A, C126V, C126M, and C126T, were characterized. Crystal structures of the 3-phosphoglycolate-bound C126S mutant and the unliganded forms of the C126S and C126A mutants were determined at a resolution of 1.7–2.1 Å. Kinetic studies revealed an approximately five-fold drop in kcat for the C126S and C126A mutants, whereas an approximately 10-fold drop was observed for the other three mutants. At ambient temperature, the wild-type enzyme and all five mutants showed no concentration dependence of activity. At higher temperatures (> 40 °C), the mutants showed a significant concentration dependence, with a dramatic loss in activity below 15 μm. The mutants also had diminished thermal stability at low concentration, as monitored by far-UV CD. These results suggest that Cys126 contributes to the stability of the dimer interface through a network of interactions involving His95, Glu97, and Arg98, which form direct contacts across the dimer interface.

Probing the role of the fully conserved Cys126 in triosephosphate isomerase by site-specific mutagenesis - distal effects on dimer stability

DatabaseStructural data are available in the Protein Data Bank under the accession numbers 3PVF, 3PY2, and 3PWA. Structured digital abstract Tim binds to Tim by x-ray crystallography (View interaction).