Biomedical subjects

C Ramakrishnan

Publications and source records attributed to C Ramakrishnan.

At least 19 recordsLinked to original sources

PubMed2007

Thermodynamic effects of proline introduction on protein stability.

The amino acid Pro is more rigid than other naturally occurring amino acids and, in proteins, lacks an amide hydrogen. To understand the structural and thermodynamic effects of Pro substitutions, it was introduced at 13 different positions in four different proteins, leucine-isoleucine-valine binding protein, maltose binding protein, ribose binding protein, and thioredoxin. Three of the maltose binding protein mutants were characterized by X-ray crystallography to confirm that no structural changes had occurred upon mutation. In the remaining cases, fluorescence and CD spectroscopy were used to show the absence of structural change. Stabilities of wild type and mutant proteins were characterized by chemical denaturation at neutral pH and by differential scanning calorimetry as a function of pH. The mutants did not show enhanced stability with respect to chemical denaturation at room temperature. However, 6 of the 13 single mutants showed a small but significant increase in the free energy of thermal unfolding in the range of 0.3-2.4 kcal/mol, 2 mutants showed no change, and 5 were destabilized. In five of the six cases, the stabilization was because of reduced entropy of unfolding. However, the magnitude of the reduction in entropy of unfolding was typically several fold larger than the theoretical estimate of -4 cal K(-1) mol(-1) derived from the relative areas in the Ramachandran map accessible to Pro and Ala residues, respectively. Two double mutants were constructed. In both cases, the effects of the single mutations on the free energy of thermal unfolding were nonadditive.

Amino Acid Substitution↗

PubMed2005

Protein stabilization by introduction of cross-strand disulfides.

Disulfides cross-link residues in a protein that are separated in primary sequence and stabilize the protein through entropic destabilization of the unfolded state. While the removal of naturally occurring disulfides leads to protein destabilization, introduction of engineered disulfides does not always lead to significant stabilization of a protein. We have analyzed naturally occurring disulfides that span adjacent antiparallel strands of beta sheets (cross-strand disulfides). Cross-strand disulfides have recently been implicated as redox-based conformational switches in proteins such as gp120 and CD4. The propensity of these disulfides to act as conformational switches was postulated on the basis of the hypothesis that this class of disulfide is conformationally strained. In the present analysis, there was no evidence to suggest that cross-strand disulfides are more strained compared to other disulfides as assessed by their torsional energy. It was also observed that these disulfides occur solely at non-hydrogen-bonded (NHB) registered pairs of adjacent antiparallel strands and not at hydrogen-bonded (HB) positions as suggested previously. One of the half-cystines involved in cross-strand disulfide formation often occurs at an edge strand. Experimental confirmation of the stabilizing effects of such disulfides was carried out in Escherichia coli thioredoxin. Four pairs of cross-strand cysteines were introduced, two at HB and two at NHB pairs. Disulfides were formed in all four cases. However, as predicted from our analysis, disulfides at NHB positions resulted in an increase in melting temperature of 7-10 degrees C, while at HB positions there was a corresponding decrease of -7 degrees C. The reduced state of all proteins had similar stability.

Animals↗

PubMed2004

Conformational analysis of invariant peptide sequences in bacterial genomes.

The functional significance of evolutionarily conserved motifs/patterns of short regions in proteins is well documented. Although a large number of sequences are conserved, only a small fraction of these are invariant across several organisms. Here, we have examined the structural features of the functionally important peptide sequences, which have been found invariant across diverse bacterial genera. Ramachandran angles (phi,psi) have been used to analyze the conformation, folding patterns and geometrical location (buried/exposed) of these invariant peptides in different crystal structures harboring these sequences. The analysis indicates that the peptides preferred a single conformation in different protein structures, with the exception of only a few longer peptides that exhibited some conformational variability. In addition, it is noticed that the variability of conformation occurs mainly due to flipping of peptide units about the virtual C(alpha)...C(alpha) bond. However, for a given invariant peptide, the folding patterns are found to be similar in almost all the cases. Over and above, such peptides are found to be buried in the protein core. Thus, we can safely conclude that these invariant peptides are structurally important for the proteins, since they acquire unique structures across different proteins and can act as structural determinants (SD) of the proteins. The location of these SD peptides on the protein chain indicated that most of them are clustered towards the N-terminal and middle region of the protein with the C-terminal region exhibiting low preference. Another feature that emerges out of this study is that some of these SD peptides can also play the roles of "fold boundaries" or "hinge nucleus" in the protein structure. The study indicates that these SD peptides may act as chain-reversal signatures, guiding the proteins to adopt appropriate folds. In some cases the invariant signature peptides may also act as folding nuclei (FN) of the proteins.

Amino Acid Sequence↗

PubMed2003

A novel complexity measure for comparative analysis of protein sequences from complete genomes.

Analysis of sequence complexities of proteins is an important step in the characterization and classification of new genomes. A new measure has been proposed to compute sequence complexity in protein sequences based on linguistic complexity. The algorithm requires a single parameter, is computationally simple and provides a framework for comparative genomic analysis. Protein sequences were classified into groups of high or low complexity based on a quantitative measure termed F(c), which is proportional to the fraction of low complexity sequence present in the protein. The algorithm was tested on sequences of 196 non-homologous proteins whose crystal structures are available at </=2.0 A resolution. Protein sequences of high complexity had 'globular' structures (95% agreement), whereas those of low complexity had non-globular structures (80% agreement). Application of this measure to proteins of unknown structure/function from different genomes revealed that the sequences of high complexity constitute the majority in all genomes (about 90% in Archaea, about 93% in Eubacteria, 89% in Saccharomyces cerevisiae and 90% in Caenorhabditis elegans). Aeropyrum pernix among Archaeae and Deinococcus radiodurans among Eubacteria have the lowest fraction of high complexity proteins (75% and 80% respectively). Further, it was observed that a few bacterial pathogens (Mycobacterium tuberculosis, Pseudomonas aeruginosa) have high fraction of low complexity proteins. The program ScanCom is available from the authors as a PERL script (UNIX system).

Algorithms↗

PubMed2003

MODIP revisited: re-evaluation and refinement of an automated procedure for modeling of disulfide bonds in proteins.

There have been several attempts to stabilize proteins through the introduction of engineered disulfide bonds. For reasons that are currently unclear, these have met with mixed success. Hence identification of locations where introduction of a disulfide cross-link will lead to protein stabilization is still a challenging task. A computational procedure, MODIP, was introduced more than a decade ago to select sites in protein structures that have the correct geometry for disulfide formation when replaced by Cys. In this study, we re-evaluated the stereochemical criteria used by MODIP for the selection and gradation of sites for modeling disulfides. We introduced steric criteria to check for energetically unfavorable non-bonded contacts with the modeled disulfide, since these can considerably offset the stabilizing effect of the cross-link. The performance of the refined procedure was checked for its ability to correctly predict naturally occurring disulfide bonds in proteins. A set of proteins in which disulfide bonds were introduced experimentally were analyzed with respect to MODIP predictions, stability and other parameters such as accessibility, residue depth, B-factors of the mutated sites, change in volume upon mutation and loop length enclosed by the disulfide. The analysis suggests that in addition to proper stereochemistry, stabilizing disulfides occur in regions of low depth, relatively high mobility, have a loop length greater than 25 and where the disulfide typically occupies a volume less than or equal to that of the original residues.

Computational Biology↗

PubMed2003

Stranded in isolation: structural role of isolated extended strands in proteins.

Reasons for the formation of extended-strands (E-strands) in proteins are often associated with the formation of beta-sheets. However E-strands, not part of beta-sheets, commonly occur in proteins. This raises questions about the structural role and stability of such isolated E-strands. Using a dataset of 250 largely non-homologous and high-resolution (<2 A) crystal structures of proteins, we have identified 518 isolated E-strands from 187 proteins. The two most distinguishing features of isolated E-strands from beta-strands in beta-sheets are the high preponderance of prolyl residues occuring in isolated E-strands and their high exposure to the surroundings. Removal of regions with polyproline conformation from the dataset did not significantly reduce the propensity of prolyl residues to occur in isolated E-strands. Isolated E-strands are often characterized by their main-chain amide and carbonyl groups involved in hydrogen bonding with polar side chains or water. They are often flanked by irregular loop structures and are less well conserved, than beta-sheet forming beta-strands, among homologous protein structures. It is suggested that isolated beta-strands have many characteristics of loop segments but with repetitive (phi,psi) values falling within the beta-region of the Ramachandran map.

Crystallography, X-Ray↗

PubMed2002

Influence of solvent molecules on the stereochemical code of glycyl residues in proteins.

The Ramachandran steric map and energy diagrams of the glycyl residue are symmetric. A plot of (phi,psi) angles of glycyl residues in 250 nonhomologous and high-resolution protein structures is also largely symmetric. However, there is a clear aberration in the symmetry. Although there is a cluster of points corresponding to the right-handed alpha-helical region, the "equivalent" cluster is clearly shifted to in and around the (phi,psi) values of (90 degrees, 0 degrees ) instead of being centered at the left-handed alpha-helical region of (60 degrees, 40 degrees ). This lack of symmetry exists even in the (phi,psi) distribution of residues from non-alpha-helical regions in proteins. Here we provide an explanation for this observation. An analysis of glycyl conformations in small peptide structures and in "coil" proteins, which are largely devoid of helical and sheet regions, shows that glycyl residues prefer to adopt conformations around (+/-90 degrees, 0 degrees ) instead of right- and left-handed alpha-helical regions. By using theoretical calculations, such conformations are shown to have highest solvent accessibility in a system of two-linked peptide units with glycyl residue at the central C(alpha) atom. This finding is consistent with the observations from 250 nonhomologous protein structures where glycyl residues with conformations close to (+/-90 degrees, 0 degrees ) are seen to have high solvent accessibility. Analysis of a subset of nonhomologous structures with very high resolution (1.5 A or better) shows that water molecules are indeed present at distances suitable for hydrogen bond interaction with glycyl residues possessing conformations close to (+/-90 degrees, 0 degrees ). It is suggested that water molecules play a key role in determining and stabilizing these conformations of glycyl residues and explain the aberration in the symmetry of glycyl conformations in proteins.

Animals↗

PubMed2002

A conformational analysis of Walker motif A [GXXXXGKT (S)] in nucleotide-binding and other proteins.

The sequence GXXXXGKT/S, popularly known as Walker motif A, is widely believed to be the site for binding nucleotides in many proteins. Examination of the crystal structures in the Protein Data Bank showed that about half of the examples having these sequences do not bind or use nucleotides. Data analyses showed 92 different Walker sequences of the variable quartet (XXXX). Ramachandran angles in this segment revealed conformational similarity in the group of 45 proteins, known to bind or utilize nucleotides. The conformations of this segment in other proteins differ widely and it is not known whether they play any role in their functions. A flip of a peptide unit at different locations, with little change in the backbone conformation was noted in nine pairs of these proteins having same Walker sequence. An examination of the immediate neighborhood of the Walker sequence indicates that this region is preceded by a beta-strand and followed by an alpha-helix, resulting in the motif beta-W-alpha, an invariant feature amongst nucleotide-binding proteins.

Amino Acid Motifs↗

PubMed2001

In memoriam: Professor G.N. Ramachandran (1922-2001).

Few scientists contribute an idea of such clarity and power that it appears in all the discipline's textbooks and bears the author's name. For the contribution to be relevant and universally employed almost forty years after it first appeared is even less common. Structural biology lost the author of such an idea with the death of G.N. Ramachandran, whose picture appears on the cover of this issue of Protein Science. His seminal contribution is described in remembrances of Professor Ramachandran's life and career by colleague and co-author of the 1963 paper, C. Ramakrishnan. A perspective by George D. Rose follows, which articulates the enduring impact of that work.

Biophysics↗

PubMed2000

Deterministic features of side-chain main-chain hydrogen bonds in globular protein structures.

A total of 19 835 polar residues from a data set of 250 non-homologous and highly resolved protein crystal structures were used to identify side-chain main-chain (SC-MC) hydrogen bonds. The ratio of the number of SC-MC hydrogen bonds to the total number of polar residues is close to 1:2, indicating the ubiquitous nature of such hydrogen bonds. Close to 56% of the SC-MC hydrogen bonds are local involving side-chain acceptor/donor ('i') and a main-chain donor/acceptor within the window i-5 to i+5. These short-range hydrogen bonds form well defined conformational motifs characterized by specific combinations of backbone and side-chain torsion angles. (a) The Ser/Thr residues show the greatest preference in forming intra-helical hydrogen bonds between the atoms O(gamma)(i) and O(i-4). More than half the examples of such hydrogen bonds are found at the middle of alpha-helices rather than at their ends. The most favoured motif of these examples is alpha(R)alpha(R)alpha(R)alpha(R)(g(-)). (b) These residues also show great preference to form hydrogen bonds between O(gamma)(i) and O(i-3), which are closely related to the previous type and though intra-helical, these hydrogen bonds are more often found at the C-termini of helices than at the middle. The motif represented by alpha(R)alpha(R)alpha(R)alpha(R)(g(+)) is most preferred in these cases. (c) The Ser, Thr and Glu are the most frequently found residues participating in intra-residue hydrogen bonds (between the side-chain and main-chain of the same residue) which are characterized by specific motifs of the form beta(g(+)) for Ser/Thr residues and alpha(R)(g(-)g(+)t) for Glu/Gln. (d) The side-chain acceptor atoms of Asn/Asp and Ser/Thr residues show high preference to form hydrogen bonds with acceptors two residues ahead in the chain, which are characterized by the motifs beta (tt')alphaR and beta(t)alpha(R), respectively. These hydrogen bonded segments, referred to as Asx turns, are known to provide stability to type I and type I' beta-turns. (e) Ser/Thr residues often form a combination of SC-MC hydrogen bonds, with the side-chain donor hydrogen bonded to the carbonyl oxygen of its own peptide backbone and the side-chain acceptor hydrogen bonded to an amide hydrogen three residues ahead in the sequence. Such motifs are quite often seen at the beginning of alpha-helices, which are characterized by the beta(g(+))alpha(R)alpha(R) motif. A remarkable majority of all these hydrogen bonds are buried from the protein surface, away from the surrounding solvent. This strongly indicates the possibility of side-chains playing the role of the backbone, in the protein interiors, to satisfy the potential hydrogen bonding sites and maintaining the network of hydrogen bonds which is crucial to the structure of the protein.

Amino Acids↗

PubMed1999

Prediction of the maximal stability temperature of monomeric globular proteins solely from amino acid sequence.

Globular protein thermostability is characterized the cold denaturation, maximal stability (Tms) and heat denaturation temperatures. For mesophilic globular proteins, Tms typically ranges from -25 degrees C to +35 degrees C. We show that the indirect estimate of Tms from calorimetry and the direct estimate from chemical denaturation performed in a range of temperatures are in close agreement. The heat capacity change of unfolding per mol residue (delta Cp) alone is shown to accurately predict Tms. Delta Cp and hence Tms can be predicted solely from the protein sequence. The average difference in free energy of unfolding at the observed and predicted values of Tms is 1.0 kcal mol(-1), which is small compared to typical values of the total free energy of unfolding.

Algorithms↗

PubMed1999

Secondary structures without backbone: an analysis of backbone mimicry by polar side chains in protein structures.

Backbone mimicry by the formation of closed-loop C7, C10 and C13 (mimics of gamma-, beta- and alpha-turns) conformations through side chain-main chain hydrogen bonds by polar groups is a frequent observation in protein structures. A data set of 250 non-homologous and high-resolution protein crystal structures was used to analyze these conformations for their characteristic features. Seven out of the nine polar residues (Ser, Thr, Asn, Asp, Gln, Glu and His) have hydrogen bonding groups in their side chains which can participate in such mimicry and as many as 15% of all these polar residues engage in such conformations. The distributions of dihedral angles of these mimics indicate that only certain combinations of the dihedral angles involved aid the formation of these mimics. The observed examples were categorized into various classes based on these combinations, resulting in well defined motifs. Asn and Asp residues show a very high capability to perform such backbone secondary structural mimicry. The most highly mimicked backbone structure is of the C10 conformation by the Asx residues. The mimics formed by His, Ser, Thr and Glx residues are also discussed. The role of such conformations in initiating the formation of regular secondary structures during the course of protein folding seems significant.

Databases, Factual↗

PubMed1999

Knowledge-based modeling of the serine protease triad into non-proteases.

The Asp-His-Ser triad of serine proteases has been regarded, in the present study, as an independent catalytic motif, because in nature it has been incorporated at the active sites of enzymes as diverse as the serine proteases and the lipases. Incorporating this motif into non-protease scaffolds, by rational design and mutagenesis, might lead to the generation of novel catalysts. As an aid to such experiments, a knowledge-based computer modeling procedure has been developed to model the protease Asp-His-Ser triad into non-proteases. Catalytic triads from a set of trypsin family proteases have been analyzed and criteria that characterize the geometry of the triads have been obtained. Using these criteria, the modeling procedure first identifies sites in non-proteases that are suitable for modeling the protease triad. H-bonded Asp-His-Ser triads, that mimic the protease catalytic triad in geometry, are then modeled in at these sites, provided it is stereochemically possible to do so. Thus non-protease sites at which H-bonded Asp-His-Ser triads are successfully modeled in may be considered for mutagenesis experiments that aim at introducing the protease triad into non-proteases. The triad modeling procedure has been used to identify sites for introducing the protease triad in three binding proteins and an immunoglobulin. A scoring function, depending on inter-residue distances, solvent accessibility and the substitution potential of amino acid residues at the modeling sites in the host proteins, has been used to assess the quality of the model triads.

Animals↗

PubMed1998

Conformational interconversions in peptide beta-turns: analysis of turns in proteins and computational estimates of barriers.

The two most important beta-turn features in peptides and proteins are the type I and type II turns, which differ mainly in the orientation of the central peptide unit. Facile conformational interconversion is possible, in principle, by a flip of the central peptide unit. Homologous crystal structures afford an opportunity to structurally characterize both possible conformational states, thus allowing identification of sites that are potentially stereochemically mobile. A representative data set of 250 high-resolution (</=2.0 A), non-homologous protein crystal structures and corresponding variant and homologous entries, obtained from the Brookhaven Protein Data Bank, was examined to identify turns that are assigned different conformational types (type I/type II) in related structures. A total of 55 examples of beta-turns were identified as possible candidates for a stereochemically mobile site. Of the 55 examples, 45 could be classified as a potential site for interconversion between type I and type II beta-turns, while ten correspond to flips from type I' to type II' structures. As a further check, the temperature factors of the central peptide unit carbonyl oxygen atom of the 55 examples were examined. The analysis reveals that the turn assignments are indeed reliable. Examination of the secondary structures at the flanking positions of the flippable beta-turns reveals that seven examples occur in the loop region of beta-hairpins, indicating that the formation of ordered secondary structures on either side of the beta-turn does not preclude local conformational variations. In these beta-turns, Pro (11 examples), Lys (nine examples) and Ser (seven examples) were most often found at the i+1 position. Glycine was found to occur overwhelmingly at position i+2 (28 examples), while Ser (seven examples) and Asn (six examples) were amongst the most frequent residues. Activation energy barriers for the interconversion between type I and type II beta-turns were computed using the peptide models Ac-Pro-Aib-NHMe and Ac-Pro-Gly-NHMe within the framework of the AM1 semi-empirical molecular orbital procedure. In order to have a uniform basis for comparison and to eliminate the distracting influence of the deviation of backbone dihedral angles from that expected for ideal beta-turns, the dihedral angles phii+1 and psii+2 were fixed at the ideal values (phii+1=-60 degrees and psii+2=0 degrees). The other two angles (psii+1 and phii+2) were varied systematically to go from type II to type I beta-turn structures. The computational results suggest that there exists one stereospecific, concerted flip of the central peptide unit involving correlated single bond rotation that can occur with an activation barrier of the order of 3 kcal/mol. The results presented here suggest that conformational variations in beta-turns are observed in protein crystal structures and such changes may be an important dynamic feature in solution.

Models, Chemical↗

PubMed1998

Stereochemical punctuation marks in protein structures: glycine and proline containing helix stop signals.

An analysis on the nature of alpha-helix stop signals has been carried out, using a dataset of 1057 helices identified from 250 high resolution (</=2.0 A), non-homologous, protein crystal structures. The backbone dihedral angles (phi, psi) of the terminating residue (T) were found to cluster either in the left-handed helical region (alphaL: phi=20 degrees to 125 degrees and psi=-45 degrees to 90 degrees ; 469 helices (44%)) or in the extended region (E: phi=-180 degrees to -30 degrees and psi=60 degrees to 180 degrees and -180 degrees to -150 degrees ; 459 helices (43%)) of the Ramachandran map. These two broad categories of helix stop signals, alphaL and E-terminated helices, were further examined for sequence preferences. Gly residues were found to have an overwhelming preference to occur as the "alphaL-terminator (T)" resulting in the classical Schellman motif, with a strong preference for hydrophobic residues at position T-4 and T+1. In the case of E-terminated helices His, Asn, Leu and Phe were found to occur with high propensity at position T. Quite remarkably Pro residues, with single exception, were absent at position T, but had the highest propensity at position T+1. Examination of the frequencies of hydrophobic (h) and polar (p) residues at positions flanking Gly/Pro permitted delineation of exclusive patterns and predictive rules for Gly-terminated helices and Pro-terminated helices. The analysis reveals that Pro residues flanked by polar amino acids have a very strong tendency to terminate helices. Examination of a segment ranging from T-4 to T+3 appeared to be necessary to determine whether helix termination or continuation occur at Gly residues. The two types of helix termination (alphaL, E) signals also differed dramatically in their solvent accessibility. Gly and Pro residues at helix termini appeared to be strongly conserved in homologous sequences.

Forecasting↗

PubMed1998

Energy minimization studies on alpha-turns.

Using a grid search technique, the entire conformational space of a system of four linked peptide units (tetrapeptide) was scanned to pick out geometrically possible 5-->1 type hydrogen-bonded conformations defined as an alpha-turn. The energy minimization of these conformations led to 23 distinct minimum energy conformations (MECs) falling in 13 different classes. The presence of beta and gamma turn type hydrogen bonds along with 5-->1 type hydrogen bond gave conformational variability in a given class. The occurrence of bifurcated hydrogen bonding network was a characteristic feature of most of the MECs. In many prototype MECs non-glycyl residues such as Ala and Pro could be accommodated. Comparison of MECs with the alpha-turn examples that are observed in proteins showed that the conformationally worked out MECs occurred in isolation in proteins, with the alpha-helical alpha-turn being distinctly the most predominant.

Energy Metabolism↗

PubMed1998

Molecular chaperones and subcellular trafficking of steroid receptors.

Unliganded steroid receptors exist as heteromeric complexes comprised of heat shock and immunophilin proteins that associate either directly or indirectly with receptor carboxyl-terminal ligand-binding domains. Molecular chaperons, and other proteins associated with steroid receptors, play an important role in the maturation of receptors to a hormone-binding competent state. Steroid receptor-associated 90 and 70 kDa heat shock proteins, hsp90 and hsp70, respectively, have well established roles in protein folding in addition to participating in numerous subcellular trafficking pathways. In this review, we discuss the possible roles that molecular chaperons, such as hsp90, hsp70 and DnaJ proteins, have in steroid receptor trafficking within two distinct subcellular compartments, i.e. the cytoplasm and nucleus.

Amino Acid Sequence↗

PubMed1997

Steroid hormone responsiveness of a family of closely related mouse proviral elements.

Regulation of the mouse sex-limited protein (Slp) gene in unusual in that hormone response is conferred by the 5' LTR of an upstream inserted provirus, dubbed the imposon (imp1). In a search for additional genes whose regulation has been affected by retrotransposition events, we isolated two partial proviral elements by stringent screening of a mouse genomic library. One clone (imp2) contained a portion of the envelope gene and a 3' LTR that was nearly identical to the 3' LTR of imp1; this similarity extended to insertion into a B1 repetitive element. The second proviral clone (imp3) contained a 5' LTR and associated coding sequences, but lacked its 3' LTR; the LTR of imp3 differed by 12% from the imp1 sequence. To assess potential hormone response, proviral enhancer regions cloned into reporter vectors were tested in transfection. The imp2 enhancer was similar in behavior to imp1, conferring both androgen and glucocorticoid induction in one fragment context and an androgen-specific response in another. In contrast, the imp3 enhancer allowed high expression in the absence of hormone and was less responsive to steroids in general and androgen in particular. These three proviral elements define a small family of steroid responsive proviruses in the mouse genome, and at least one member has had a lasting impact on an endogenous gene's regulation.

Animals↗