Rings and ribbons in protein structures: Characterization using helical parameters and Ramachandran plots for repeating dipeptides

Intrinsic proclivity of left-handed conformation in large Nest motif peptides inferred from molecular dynamics

The 'Nest' motif plays a functional role in protein owing to its ligand binding potential aided by geometric concavity. The presence of less favored left-handed conformation (L-state) in its structure makes this concavity possible and in shaping the native chemical environment amenable to stable binding interactions. To understand the persistent appearance of L-state torsion in the Nest motif, we analyzed 0.5μs Molecular Dynamics (MD) simulation trajectories of 35 six-residue peptides (out of a total of 50 large Nest sequences of ≥6 residues) identified in our previous study. Analysis of the MD trajectories of the individual peptides reveals initial L-state in 60% of the peptides persists for >40% of the trajectory. Further, Nests with different sequences appear to adopt a specific conformational state driven by the neighboring L-state residues. The sequences also possess short secondary structures and amino acid repeats, suggesting evolutionary conservation and the specific role of amino acids in locally predisposing the torsion angle to the L-state. These findings help us to understand how L-state conformation is an essential prerequisite in stabilizing the Nest motif and shed light on the sequence-structure-function paradigm in the rational design of peptides and peptidomimetics for therapeutics.Communicated by Ramaswamy H. Sarma.

The conserved crown bridge loop at the catalytic centre of enzymes of the haloacid dehalogenase superfamily

The crown bridge loop is hexapeptide motif in which the backbone carbonyl group at position 1 is hydrogen bonded to the backbone imino groups of positions 4 and 6. Residues at positions 1 and 4-6 are held in a tight substructure, but different orientations of the plane of the peptide bond between positions 2 and 3 result in two conformers: the 2,3-αRαR crown bridge loop - found in approximately 7% of proteins - and the 2,3-βRαL crown bridge loop, found in approximately 1-2% of proteins. We constructed a relational database in which we identified 60 instances of the 2,3-βRαL conformer, and find that about half occur in enzymes of the haloacid dehalogenase (HAD) superfamily, where they are located next to the catalytic aspartate residue. Analysis of additional enzymes of the HAD superfamily in the extensive SCOPe dataset showed this crown bridge loop to be a conserved feature. Examination of available structures showed that the 2,3-βRαL conformation - but not the 2,3-αRαR conformation - allows the backbone carbonyl group at position 2 to interact with the essential Mg2+ ion. The possibility of interconversion between the 2,3-βRαL and 2,3-αRαR conformations during catalysis is discussed.

The role of the half-turn in determining structures of Alzheimer's Aβ wild-type and mutants

Half-turns are shown to be the main determinants of many experimental Alzheimer's Aβ fibril structures. Fibril structures contain three half-turn types, βα_Rβ, βα_Lβ and βεβ which each result in a ∼90° bend in a β-strand. It is shown that only these half-turns enable cross-β stacking and thus the right-angle fold seen in fibrils is an intrinsic feature of cross-β. Encoding a strand as a conformational sequence in β, α_R, α_L and ε(β_L), pairwise combination rules for consecutive half-turns are used to decode this sequence to give the backbone path. This reveals how structures would be dramatically affected by a deletion. Using a wild-type Aβ(42) fibril structure and the pairwise combination rules, the Osaka deletion is predicted to result in exposure of surfaces that are mutually shielding from the solvent. Molecular dynamics simulations on an 11-mer β-sheet of Alzheimer's Aβ(40) of the Dutch (E22Q), Iowa (D23N), Arctic (E22G), and Osaka (E22Δ) mutants, show the crucial role glycine plays in the positioning of βα_Rβ half-turns. Their "in-phase" positions along the sequence in the wild-type, Dutch mutant and Iowa mutant means that the half-folds all fold to the same side creating the same closed structure. Their out-of-phase positions in Arctic and Osaka mutants creates a flatter structure in the former and an S-shape structure in the latter which, as predicted, exposes surfaces on the inside in the closed wild-type to the outside. This is consistent with the gain of interaction model and indicates how domain swapping might explain the Osaka mutant's unique properties.

Determination of amino acids that favour the αL region using Ramachandran propensity plots. Implications for α-sheet as the possible amyloid intermediate

In amyloid diseases an insoluble amyloid fibril forms via a soluble oligomeric intermediate. It is this intermediate that mediates toxicity and it has been suggested, somewhat controversially, that it has the α-sheet structure. Nests and α-strands are similar peptide motifs in that alternate residues lie in the α_R and γ_L regions of the Ramachandran plot for nests, or α_R and α_L regions for α-strands. In nests a concavity is formed by the main chain NH atoms whereas in α-strands the main chain is almost straight. Using "Ramachandran propensity plots" to focus on the α_L/γ_L region, it is shown that glycine favours γ_L (82% of amino acids are glycine), but disfavours α_L (3% are glycine). Most charged and polar amino acids favour α_L with asparagine having by far the highest propensity. Thus, glycine favours nests but, contrary to common expectation, should not favour α-sheet. By contrast most charged or polar amino acids should favour α-sheet by their propensity for the α_L conformation, which is more discriminating amongst amino acids than the α_R conformation. Thus, these results suggest the composition of sequences that favour α-sheet formation and point towards effective prediction of α-sheet from sequence.

Green Rust: The Simple Organizing 'Seed' of All Life?

Korenaga and coworkers presented evidence to suggest that the Earth's mantle was dry and water filled the ocean to twice its present volume 4.3 billion years ago. Carbon dioxide was constantly exhaled during the mafic to ultramafic volcanic activity associated with magmatic plumes that produced the thick, dense, and relatively stable oceanic crust. In that setting, two distinct and major types of sub-marine hydrothermal vents were active: ~400 °C acidic springs, whose effluents bore vast quantities of iron into the ocean, and ~120 °C, highly alkaline, and reduced vents exhaling from the cooler, serpentinizing crust some distance from the heads of the plumes. When encountering the alkaline effluents, the iron from the plume head vents precipitated out, forming mounds likely surrounded by voluminous exhalative deposits similar to the banded iron formations known from the Archean. These mounds and the surrounding sediments, comprised micro or nano-crysts of the variable valence FeII/FeIII oxyhydroxide known as green rust. The precipitation of green rust, along with subsidiary iron sulfides and minor concentrations of nickel, cobalt, and molybdenum in the environment at the alkaline springs, may have established both the key bio-syntonic disequilibria and the means to properly make use of them-the elements needed to effect the essential inanimate-to-animate transitions that launched life. Specifically, in the submarine alkaline vent model for the emergence of life, it is first suggested that the redox-flexible green rust micro- and nano-crysts spontaneously precipitated to form barriers to the complete mixing of carbonic ocean and alkaline hydrothermal fluids. These barriers created and maintained steep ionic disequilibria. Second, the hydrous interlayers of green rust acted as engines that were powered by those ionic disequilibria and drove essential endergonic reactions. There, aided by sulfides and trace elements acting as catalytic promoters and electron transfer agents, nitrate could be reduced to ammonia and carbon dioxide to formate, while methane may have been oxidized to methyl and formyl groups. Acetate and higher carboxylic acids could then have been produced from these C1 molecules and aminated to amino acids, and thence oligomerized to offer peptide nests to phosphate and iron sulfides, and secreted to form primitive amyloid-bounded structures, leading conceivably to protocells.

Conformations of peptoids in nanosheets result from the interplay of backbone energetics and intermolecular interactions

The conformations adopted by the molecular constituents of a supramolecular assembly influence its large-scale order. At the same time, the interactions made in assemblies by molecules can influence their conformations. Here we study this interplay in extended flat nanosheets made from nonnatural sequence-specific peptoid polymers. Nanosheets exist because individual polymers can be linear and untwisted, by virtue of polymer backbone elements adopting alternating rotational states whose twists oppose and cancel. Using molecular dynamics and quantum mechanical simulations, together with experimental data, we explore the design space of flat nanostructures built from peptoids. We show that several sets of peptoid backbone conformations are consistent with their being linear, but the specific combination observed in experiment is determined by a combination of backbone energetics and the interactions made within the nanosheet. Our results provide a molecular model of the peptoid nanosheet consistent with all available experimental data and show that its structure results from a combination of intra- and intermolecular interactions.

The Ramachandran Number: An Order Parameter for Protein Geometry

Three-dimensional protein structures usually contain regions of local order, called secondary structure, such as α-helices and β-sheets. Secondary structure is characterized by the local rotational state of the protein backbone, quantified by two dihedral angles called ϕ and ψ. Particular types of secondary structure can generally be described by a single (diffuse) location on a two-dimensional plot drawn in the space of the angles ϕ and ψ, called a Ramachandran plot. By contrast, a recently-discovered nanomaterial made from peptoids, structural isomers of peptides, displays a secondary-structure motif corresponding to two regions on the Ramachandran plot [Mannige et al., Nature 526, 415 (2015)]. In order to describe such ‘higher-order’ secondary structure in a compact way we introduce here a means of describing regions on the Ramachandran plot in terms of a single Ramachandran number, R, which is a structurally meaningful combination of ϕ and ψ. We show that the potential applications of R are numerous: it can be used to describe the geometric content of protein structures, and can be used to draw diagrams that reveal, at a glance, the frequency of occurrence of regular secondary structures and disordered regions in large protein datasets. We propose that R might be used as an order parameter for protein geometry for a wide range of applications.