Sparsely populated residue conformations in protein structures: revisiting "experimental" Ramachandran maps

Disallowed spots in protein structures

Ramachandran (ϕ, ψ) steric map was introduced in 1963 to describe available conformation space for protein structures. Subsequently, residues were observed in high-energy disallowed regions of the map. To unequivocally identify the locations of disallowed conformations of residues, we got 36 noise-free protein structures (resolution ≤1 Å, Rwork/Rfree ≤ 0.10). These stringent criteria were applied to rule out data or model errors or any crystallographic disorders. No disallowed conformation was found in the dataset. Further, we also examined disallowed conformations in a larger dataset (resolution ≤1.5 Å, devoid of any model errors, or disorders). The observed locations of disallowed residues are referred as disallowed spots. These spots include short loops of 3-5 residues, and locations where residues participate in disulfide bonding or intramolecular interactions or inter-molecular interactions with neighboring water, metals or ligands. Conformational sampling revealed that short loops in between secondary structures hardly have any opportunity to relieve from conformational strain. Residues involved in interactions, which provide energetic compensation for high-energy conformational states, were relieved from strain once the causative interaction was removed. The present study aims to identify disallowed spots in the native state of proteins, wherein residues are forced to be trapped in high-energy disallowed conformations. Moreover, it was also observed that pre-Pro, Ser, Asp, trans-Pro, Val, Asn & Gly have higher tendency to occur in disallowed conformation, which could be attributed to factors such as conformational restrictions, residue propensity of secondary structures and compensating sidechain and mainchain interactions, stabilizing turn-mimics.

Xeno Amino Acids: A Look into Biochemistry as We Do Not Know It

Would another origin of life resemble Earth's biochemical use of amino acids? Here, we review current knowledge at three levels: (1) Could other classes of chemical structure serve as building blocks for biopolymer structure and catalysis? Amino acids now seem both readily available to, and a plausible chemical attractor for, life as we do not know it. Amino acids thus remain important and tractable targets for astrobiological research. (2) If amino acids are used, would we expect the same L-alpha-structural subclass used by life? Despite numerous ideas, it is not clear why life favors L-enantiomers. It seems clearer, however, why life on Earth uses the shortest possible (alpha-) amino acid backbone, and why each carries only one side chain. However, assertions that other backbones are physicochemically impossible have relaxed into arguments that they are disadvantageous. (3) Would we expect a similar set of side chains to those within the genetic code? Many plausible alternatives exist. Furthermore, evidence exists for both evolutionary advantage and physicochemical constraint as explanatory factors for those encoded by life. Overall, as focus shifts from amino acids as a chemical class to specific side chains used by post-LUCA biology, the probable role of physicochemical constraint diminishes relative to that of biological evolution. Exciting opportunities now present themselves for laboratory work and computing to explore how changing the amino acid alphabet alters the universe of protein folds. Near-term milestones include: (a) expanding evidence about amino acids as attractors within chemical evolution; (b) extending characterization of other backbones relative to biological proteins; and (c) merging computing and laboratory explorations of structures and functions unlocked by xeno peptides.

Isolated α-turns in peptides: a selected literature survey

The results of classifying into various types the 68 examples of isolated α-turns in the X-ray diffraction crystal structures of peptides documented in the literature are presented and discussed in this review article. α-Turns characterized by the trans disposition of all ω torsion angles are common for the backbone linear peptides investigated. In contrast, the cis arrangement of the N-terminal (ω_i + 1 ) torsion angle, among those generated by the three residues internal to the α-turn, is a peculiar feature of 65% of the cyclic peptides. Among linear and cyclic peptides featuring the all-trans disposition of the ω torsion angles, only one third of the α-turns display φ,ψ values not too far from those characterizing regular α-helices. In general, our findings, taken together, suggest that a significant conformational diversity is compatible with the formation of an intramolecularly H-bonded C₁₃ -member pseudocycle (α-turn) in linear and cyclic peptides.

A gama-turn mimetic for selective sensing of Cu(II) and combinatorial multiple logic gate

We have designed and synthesized a gama-turn mimetic using fenamic acid and α-aminoisobutyricacid (Aib), the conformation and optoelectronic properties of which can be changed by appropriate external stimuli. From single-crystal...

Computational Analysis of the Mechanism of Nonenzymatic Peptide Bond Cleavage at the C-Terminal Side of an Asparagine Residue

The nonenzymatic peptide bond cleavage at the C-terminal side of Asn residues is a protein post-translational modification that occurs under physiological conditions. This reaction proceeds much slower than the deamidation of the Asn side chain and causes denaturation and hypofunction of proteins. The peptide bond cleavage of Asn is detected primarily in crystallins and aquaporin 0 in the eye lens. Therefore, cleavage is thought to be involved in age-related cataracts. In this study, to clarify the mechanism underlying succinimide formation for the peptide bond cleavage of the Asn residue, we performed quantum chemical calculations on the model compound Ace-Asn-Gly-Nme (Ace = acetyl and Nme = methylamino). The density functional theory with the B3LYP/6-31+G(d,p) level of theory was used to obtain optimized geometries. The results suggested that the reaction proceeds through two steps, cyclization and C-terminal fragment release, and the required proton transfers can be mediated by H₂PO₄ ^- and HCO₃ ^- ions. The conformational change of the main chain on the N-terminal side of Asn was needed for the C-terminal fragmentation step, and a separate conformational change at the C-terminal side was required for the cyclization step. Furthermore, the calculated activation barriers of the reactions catalyzed by the H₂PO₄ ^- ion (130 kJ mol^-1) and the HCO₃ ^- ion (123 kJ mol^-1) were sufficiently low for the reactions to occur under normal physiological conditions.

Evolution as a Guide to Designing xeno Amino Acid Alphabets

Here, we summarize a line of remarkably simple, theoretical research to better understand the chemical logic by which life’s standard alphabet of 20 genetically encoded amino acids evolved. The connection to the theme of this Special Issue, “Protein Structure Analysis and Prediction with Statistical Scoring Functions”, emerges from the ways in which current bioinformatics currently lacks empirical science when it comes to xenoproteins composed largely or entirely of amino acids from beyond the standard genetic code. Our intent is to present new perspectives on existing data from two different frontiers in order to suggest fresh ways in which their findings complement one another. These frontiers are origins/astrobiology research into the emergence of the standard amino acid alphabet, and empirical xenoprotein synthesis.

Dihedral angle preferences of amino acid residues forming various non-local interactions in proteins

In theory, a polypeptide chain can adopt a vast number of conformations, each corresponding to a set of backbone rotation angles. Many of these conformations are excluded due to steric overlaps. Ramachandran and coworkers were the first to look into this problem by plotting backbone dihedral angles in a two-dimensional plot. The conformational space in the Ramachandran map is further refined by considering the energetic contributions of various non-bonded interactions. Alternatively, the conformation adopted by a polypeptide chain may also be examined by investigating interactions between the residues. Since the Ramachandran map essentially focuses on local interactions (residues closer in sequence), out of interest, we have analyzed the dihedral angle preferences of residues that make non-local interactions (residues far away in sequence and closer in space) in the folded structures of proteins. The non-local interactions have been grouped into different types such as hydrogen bond, van der Waals interactions between hydrophobic groups, ion pairs (salt bridges), and ππ-stacking interactions. The results show the propensity of amino acid residues in proteins forming local and non-local interactions. Our results point to the vital role of different types of non-local interactions and their effect on dihedral angles in forming secondary and tertiary structural elements to adopt their native fold.

Dihedral angle preferences of DNA and RNA binding amino acid residues in proteins

A protein can interact with DNA or RNA molecules to perform various cellular processes. Identifying or analyzing DNA/RNA binding site amino acid residues is important to understand molecular recognition process. It is quite possible to accurately model DNA/RNA binding amino acid residues in experimental protein-DNA/RNA complex by using the electron density map whereas, locating/modeling the binding site amino acid residues in the predicted three dimensional structures of DNA/RNA binding proteins is still a difficult task. Considering the above facts, in the present work, we have carried out a comprehensive analysis of dihedral angle preferences of DNA and RNA binding site amino acid residues by using a classical Ramachandran map. We have computed backbone dihedral angles of non-DNA/RNA binding residues and used as control dataset to make a comparative study. The dihedral angle preference of DNA and RNA binding site residues of twenty amino acid type is presented. Our analysis clearly revealed that the dihedral angles (φ, ψ) of DNA/RNA binding amino acid residues prefer to occupy (-89° to -60°, -59° to -30°) bins. The results presented in this paper will help to model/locate DNA/RNA binding amino acid residues with better accuracy.

Extension of the classical classification of β-turns

The functional properties of a protein primarily depend on its three-dimensional (3D) structure. These properties have classically been assigned, visualized and analysed on the basis of protein secondary structures. The β-turn is the third most important secondary structure after helices and β-strands. β-turns have been classified according to the values of the dihedral angles φ and ψ of the central residue. Conventionally, eight different types of β-turns have been defined, whereas those that cannot be defined are classified as type IV β-turns. This classification remains the most widely used. Nonetheless, the miscellaneous type IV β-turns represent 1/3^rd of β-turn residues. An unsupervised specific clustering approach was designed to search for recurrent new turns in the type IV category. The classical rules of β-turn type assignment were central to the approach. The four most frequently occurring clusters defined the new β-turn types. Unexpectedly, these types, designated IV₁, IV₂, IV₃ and IV₄, represent half of the type IV β-turns and occur more frequently than many of the previously established types. These types show convincing particularities, in terms of both structures and sequences that allow for the classical β-turn classification to be extended for the first time in 25 years.

Native proteins trap high-energy transit conformations

During protein folding and as part of some conformational changes that regulate protein function, the polypeptide chain must traverse high-energy barriers that separate the commonly adopted low-energy conformations. How distortions in peptide geometry allow these barrier-crossing transitions is a fundamental open question. One such important transition involves the movement of a non-glycine residue between the left side of the Ramachandran plot (that is, ϕ < 0°) and the right side (that is, ϕ > 0°). We report that high-energy conformations with ϕ ~ 0°, normally expected to occur only as fleeting transition states, are stably trapped in certain highly resolved native protein structures and that an analysis of these residues provides a detailed, experimentally derived map of the bond angle distortions taking place along the transition path. This unanticipated information lays to rest any uncertainty about whether such transitions are possible and how they occur, and in doing so lays a firm foundation for theoretical studies to better understand the transitions between basins that have been little studied but are integrally involved in protein folding and function. Also, the context of one such residue shows that even a designed highly stable protein can harbor substantial unfavorable interactions.

Single and multiple peptide γ-turns: literature survey and recent progress

Published data on peptide isolated and repetitive γ-turns are reviewed. Advancements in our laboratories on these 3D-structures are also presented.