C. Ramakrishnan

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. Proteins 2007.

Termination of right handed helices in proteins by residues in left handed helical conformations

An analysis of 636 helical segments, ranging in length from 4 to 32 residues, from 123 independent protein crystal structures reveals that helix termination by residues in left handed (αL) helical conformations is a common occurrence. Gly and Asn residues are the most frequent αL helix terminators, with the former having a very high propensity to adopt such conformations. The αR‐αR‐αR‐αL segment at the C termini of protein helices often possesses a 6 → 1 (π‐type) hydrogen bond between the CO of residue i and the NH of residue i + 5 with residue i + 4 occurring in the αL conformation. A stereochemical analysis of 216 examples shows that in 62 cases the 6 → 1 hydrogen bond is absent. The present analysis provides a quantitative measure of the propensity of the 20 amino acids to adopt αL helix terminating conformations.

Conformational analysis and design of cross-strand disulfides in antiparallel β-sheets

Cross‐strand disulfides bridge two cysteines in a registered pair of antiparallel β‐strands. A nonredundant data set comprising 5025 polypeptides containing 2311 disulfides was used to study cross‐strand disulfides. Seventy‐six cross‐strand disulfides were found of which 75 and 1 occurred at non‐hydrogen‐bonded (NHB) and hydrogen‐bonded (HB) registered pairs, respectively. Conformational analysis and modeling studies demonstrated that disulfide formation at HB pairs necessarily requires an extremely rare and positive χ1 value for at least one of the cysteine residues. Disulfides at HB positions also have more unfavorable steric repulsion with the main chain. Thirteen pairs of disulfides were introduced in NHB and HB pairs in four model proteins: leucine binding protein (LBP), leucine, isoleucine, valine binding protein (LIVBP), maltose binding protein (MBP), and Top7. All mutants LIVBP T247C V331C showed disulfide formation either on purification, or on treatment with oxidants. Protein stability in both oxidized and reduced states of all mutants was measured. Relative to wild type, LBP and MBP mutants were destabilized with respect to chemical denaturation, although the sole exposed NHB LBP mutant showed an increase of 3.1°C in T m . All Top7 mutants were characterized for stability through guanidinium thiocyanate chemical denaturation. Both exposed and two of the three buried NHB mutants were appreciably stabilized. All four HB Top7 mutants were destabilized (ΔΔG 0 = −3.3 to −6.7 kcal/mol). The data demonstrate that introduction of cross‐strand disulfides at exposed NHB pairs is a robust method of improving protein stability. All four exposed Top7 disulfide mutants showed mild redox activity. Proteins 2010.

Orthogonal ββ motifs in proteins

A super-secondary structural motif comprising two orthogonally oriented β-strands connected by short linking segments of ≤5 residues has been identified from a data set of 65 independent protein crystal structures. Of the 42 examples from 14 proteins, a vast majority have only a single residue as the linking element. Analysis of the conformational angles at the junction reveals that the recently described type VIII β-turn occurs frequently at the connecting hinge, while the type II β-turn is also fairly common.

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins

Short range side chain‐backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution ≤ 1.5 Å). Hydrogen bonds involving residues i − 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (Cn, where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C10 (COδi …NHi + 2), C13 (COδi …NHi + 3) and C17 (NδHi …COi − 4) structures. In contrast, Gln predominantly forms C16 (COεi …NHi − 3), C12 (NεHi …COi − 2), C15 (NεHi …COi − 3) and C18 (NεHi …COi − 4) motifs, with only the C18motif being analogous to the Asn C17structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone β‐turns and α‐turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain‐backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered. Proteins 2012;

A Novel Complexity Measure for Comparative Analysis of Protein Sequences from Complete Genomes

Abstract 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 Fc 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 Å 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).

Sparsely populated residue conformations in protein structures: Revisiting “experimental” Ramachandran maps

The Ramachandran map clearly delineates the regions of accessible conformational (φ–ψ) space for amino acid residues in proteins. Experimental distributions of φ, ψ values in high‐resolution protein structures, reveal sparsely populated zones within fully allowed regions and distinct clusters in apparently disallowed regions. Conformational space has been divided into 14 distinct bins. Residues adopting these relatively rare conformations are presented and amino acid propensities for these regions are estimated. Inspection of specific examples in a completely “arid”, fully allowed region in the top left quadrant establishes that side‐chain and backbone interactions may provide the energetic compensation necessary for populating this region of φ–ψ space. Asn, Asp, and His residues showed the highest propensities in this region. The two distinct clusters in the bottom right quadrant which are formally disallowed on strict steric considerations correspond to the gamma turn (C7 axial) conformation (Bin 12) and the i + 1 position of Type II′ β turns (Bin 13). Of the 516 non‐Gly residues in Bin 13, 384 occupied the i + 1 position of Type II′ β turns. Further examination of these turn segments revealed a high propensity to occur at the N‐terminus of helices and as a tight turn in β hairpins. The β strand–helix motif with the Type II′ β turn as a connecting element was also found in as many as 57 examples. Proteins 2014; 82:1101–1112.

Structural Compromise of Disallowed Conformations in Peptide and Protein Structures

Using a data set of 454 crystal structures of peptides and 80 crystal structures of non-homologous proteins solved at ultra high resolution of 1.2 A or better we have analyzed the occurrence of disallowed Ramachandran (phi, psi) angles. Out of 1492 and 13508 non-glycyl residues in peptides and proteins respectively 12 and 76 residues in the two datasets adopt clearly disallowed combinations of Ramachandran angles. These examples include a number of conformational points which are far away from any of the allowed regions in the Ramachandran map. According to the Ramachandran map a given (phi, psi) combination is considered disallowed when two non-bonded atoms in a system of two-linked peptide units with ideal geometry are prohibitively proximal in space. However, analysis of the disallowed conformations in peptide and protein structures reveals that none of the observations of disallowed conformations in the crystal structures correspond to a short contact between non-bonded atoms. A further analysis of deviations of bond lengths and angles, from the ideal peptide geometry, at the residue positions of disallowed conformations in the crystal structures suggest that individual bond lengths and angles are all within acceptable limits. Thus, it appears that the rare tolerance of disallowed conformations is possible by gentle and acceptable deviations in a number of bond lengths and angles, from ideal geometry, over a series of bonds resulting in a net gross effect of acceptable non-bonded inter-atomic distances.

Occurence of a single helix of the collagen type in globular proteins

The occurrence of an eight-residue long segment of polypeptide chain in collagen helical conformation has been detected in bacteriochlorophyll a-protein by the application of an algorithm for identifying secondary structures in globular proteins from their alpha-carbon positions. This segment spans residues 277 to 284 of the protein and is the longest known stretch of collagen helix to be observed in globular proteins.

Disulfide conformation and design at helix N‐termini

To understand structural and thermodynamic features of disulfides within an α‐helix, a non‐redundant dataset comprising of 5025 polypeptide chains containing 2311 disulfides was examined. Thirty‐five examples were found of intrahelical disulfides involving a CXXC motif between the N‐Cap and third helical positions. GLY and PRO were the most common amino acids at positions 1 and 2, respectively. The N‐Cap residue for disulfide bonded CXXC motifs had average (ϕ,ψ) values of (−112 ± 25.2°, 106 ± 25.4°). To further explore conformational requirements for intrahelical disulfides, CYS pairs were introduced at positions N‐Cap‐3; 1,4; 7,10 in two helices of an Escherichia coli thioredoxin mutant lacking its active site disulfide (nSS Trx). In both helices, disulfides formed spontaneously during purification only at positions N‐Cap‐3. Mutant stabilities were characterized by chemical denaturation studies (in both oxidized and reduced states) and differential scanning calorimetry (oxidized state only). All oxidized as well as reduced mutants were destabilized relative to nSS Trx. All mutants were redox active, but showed decreased activity relative to wild‐type thioredoxin. Such engineered disulfides can be used to probe helix start sites in proteins of unknown structure and to introduce redox activity into proteins. Conversely, a protein with CYS residues at positions N‐Cap and 3 of an α‐helix is likely to have redox activity. Proteins 2010.