Descriptive statistics of disallowed regions and various protein secondary structures in the context of studying twisted β-hairpins

A detailed analysis of polypeptide-chain backbone conformations was carried out for polypeptide-chain segments adjacent to β-turn regions, including the sites of disallowed conformations. A cross comparison of conformations was performed for disallowed regions of the Ramachandran plot and main types of β-turns and adjacent secondary structures. Based on the results, disallowed region 2 (II, II') in the Ramachandran plot was shown to coincide mainly with β-hairpins and, more exactly, twisted β-hairpins. The frequency of residues with angles ?_i, ψ_i that fall in region 2 (II, II') in the latter is 140 times higher than in common β-hairpins.     

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.     

On residues in the disallowed region of the Ramachandran map.

An analysis of the occurrence of nonglycyl residues in conformations disallowed in the Ramachandran plot is presented. Ser, Asn, Thr, and Cys have the highest propensities to exhibit such conformations, and the branched aliphatic residues the lowest. Residues cluster in five regions and there are some trends in the types of residues and their side-chain conformations (chi(1)) occupying these. Majority of the residues are found at the edge of helices and strands and in short loops, and are involved in different types of weak, stabilizing interactions. A structural motif has been identified where a residue in disallowed conformation occurs as the first residue of a short 3(10)-helix. On the basis of the types of neighboring residues, the location in the three-dimensional structure and accessibility, there are similarities with the occurrence of cis peptide bonds in protein structures.     

Functional analysis and structure determination of alkaline protease from Aspergillus flavus

Proteases are one of the highest value commercial enzymes as they have broad applications in food, pharmaceutical, detergent, and dairy industries and serve as vital tools in determination of structure of proteins and polypeptides. Multiple application of these enzymes stimulated interest to discover them with novel properties and considerable advancement of basic research into these enzymes. A broad understanding of the active site of the enzyme and of the mechanism of its inactivation is essential for delineating its structure-function relationship. Primary structure analysis of alkaline protease showed 42% of its content to be alpha helix making it stable for three dimensional structure modeling. Homology model of alkaline protease has been constructed using the X-ray structure (3F7O) as a template and swiss model as the workspace. The model was validated by ProSA, SAVES, PROCHECK, PROSAII and RMSD. The results showed the final refined model is reliable. It has 53% amino acid sequence identity with the template, 0.24 Å as RMSD and has -7.53 as Z-score, the Ramachandran plot analysis showed that conformations for 83.4 % of amino acid residues are within the most favored regions and only 0.4% in the disallowed regions.     

Copper(II)--nucleic acid interactions--a conformational study

The effect of copper complexation upon the conformational angles in mononucleotides and dinucleotides was studied using classical potential energy expressions. The geometrical parameters were taken from crystallographic reports on the base-metal complexes. It was found that the preferred conformations of the complexes were quite different from the preferred conformations of the corresponding pure nucleotides. Notable among the disallowed conformations in case of Cu(II)-dinucleotide complexes were the regions corresponding to RRNA's, RNA's and the A, B and C forms of DNA. It is proposed that the energetics of a single strand is more important than inter-strand interactions in promoting copper-induced denaturation of DNA.     

Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety‐six of these alignment positions correspond to conserved Gly residue with (φ, ψ) values allowed for non‐glycyl residues. Reasons for this observation were investigated by in‐silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C^β atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (φ, ψ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (φ, ψ) values which are disallowed for Ala. In‐silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C^β atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.     

Vegetation of Southern Patagonia in the 1970s – Digitization of a gray-literature data set as a monitoring baseline in a changing world

Monitoring vegetation trends against objective baselines is fundamental to quantify the impacts of global change on plant biodiversity. Vegetation plot time series are a gold standard in vegetation monitoring, but such data are missing for many regions. Southern Patagonia is an example of a region strongly impacted by climate change but lacking time series data. Monitoring in this region could benefit from a comparison with vegetation survey data gathered between 1975 and 1979, as part of the multidisciplinary research program “Transecta botánica de la Patagonia austral” (hereafter Transecta). Published in 1985, it contains data on 668 vegetation plots, which were so far inaccessible to most researchers: Transecta has never been reprinted, nor fully digitized, and can only be found in specialized libraries. Here, we created a reproducible workflow, documenting how vegetation plot data from historical sources can be extracted and harmonized. The resulting open-access database we created fills a major regional gap and provides a needed baseline to assess the impacts of global change on southern Patagonia vegetation. By making these data available, we hope to inspire a new generation of vegetation scientists to resurvey the area and continue the legacy of the pioneer researchers who compiled Transecta.     

Structural study of binding of flagellin by Toll-like receptor 5

In order to predict the binding regions within the complex formed by Toll-like receptor 5 (TLR-5) and flagellin, a complementary hydropathy between the two proteins was sought. A region common to the flagellins of Salmonella enterica serovar Typhimurium, Pseudomonas aeruginosa, and Listeria monocytogenes was shown to be hydropathically complementary to the 552-to-561 fragment of TLR-5, whose sequence is EILDISRNQL. The hydrophobicity profile of this region is shared with flagellins of 377 bacterial species out of a total of 723 publicly available sequences. A conformational analysis of the predicted binding site of TLR-5, whose structure is still unknown, was carried out with a methodology already applied to similar problems. To sample the conformations available to the peptide chain, a plot of the number of conformations per unit energy interval (density of states) versus energy was built. Following a theoretical argument, conformations belonging to maxima in this plot were selected. The most stable structure obtained in this search, an alpha-helical conformation, was shown to form the electrostatic interactions Glu552-Gln89, Asp555-Arg92, and Arg558-Glu93 with the predicted binding site of the flagellin of S. enterica serovar Typhimurium, formed by the 88-to-97 chain fragment (LQRVRELAVQ), which is likewise alpha helical.     

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.     

Analysis of the code relating sequence to conformation in globular proteins. Development of a stereochemical alphabet on the basis of intra-residue information

1. The relation of primary sequence to all residue backbone conformations was explored to test out starting conformations for protein folding. 2. Information theory was used to obtain measures of information which quantitate the role of each residue in determining its own conformation; i.e. intra-residue information. 3. The information measures are plotted as a function of varphi, psi peptide-backbone angles and varphi, psi contour maps obtained for each of the 20 amino acids. These show characteristic differences between residues. 4. To find practical ways of relating sequence to varphi, psi angles, several types of stereochemical alphabet were investigated. The value of these was tested by using them to predict the varphi, psi angles of nine different proteins. 5. A difference plot was constructed to show regions of the sequence that require little or no information extra to the intra-residue information in order to predict a correct conformation. These regions are suggested to be candidates for nucleating sites in the protein.     

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. © 2013 Wiley Periodicals, Inc.     

Conformational Angles DataBase (CADB-3.0)

Transitions in amino-acid conformation angles tend to accompany various structural modifications in protein structures. Thus, to benefit the modeling of protein structures, the Conformation Angles DataBase (CADB-3.0) has been updated to visualize the conformational angles in varied regions (fully, generously, additionally and disallowed regions). In addition, options are provided to display the angles in the secondary structural elements (alpha-helix, beta-sheet and 3(10)-helix) of the Ramachandran plot. The database is being updated periodically and can be accessed over the World Wide Web at the following URL: http://cluster.physics.iisc.ernet.in/cadb/.     

Different families of double-stranded conformations of DNA as revealed by computer calculations

An algorithm has been developed that permits one to find all possible conformations of the sugar-phosphate backbone for any given disposition of DNA base pairs. For each of the conformations thus obtained, the energy of the helix was calculated by the method of atom-atom potentials. Several isolated regions in the space of the bases′ parameters (Arnott's parameters) have been found for energetically favorable helical structures. Two parameters, the distance of a base pair from the helix axis, D, and the windling angle, τ, allow one to subdivide possible conformations into the families of closely related forms. Two regions (ravines) on the (D, τ) map correspond to the know A and B families. In the B family a continuous transition has been obtained in which the double helix undergoes increasing winding, while the base pairs are moving toward the major (nonglycosidic) groove. Interrelationships between the variables, characterizing the spatial structure of the double helix, D, τ, TL and χ, when going along the bottom of the B ravine, were also obtained. Besides the Known A and B families, several new ones were found to be energetically possible. Among these the strongly underwound helices with the negative D values, as well as the forms with the C4-C5 angle in a trans position, should be mentioned. Biological roles of the different double-stranded conformations, in particular, in protein-nuclei acid interaction are discussed.     

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.     

Carbonyl-carbonyl interactions stabilize the partially allowed Ramachandran conformations of asparagine and aspartic acid

Asparagine and aspartate are known to adopt conformations in the left-handed alpha-helical region and other partially allowed regions of the Ramachandran plot more readily than any other non-glycyl amino acids. The reason for this preference has not been established. An examination of the local environments of asparagine and aspartic acid in protein structures with a resolution better than 1.5 A revealed that their side-chain carbonyls are frequently within 4 A of their own backbone carbonyl or the backbone carbonyl of the previous residue. Calculations using protein structures with a resolution better than 1.8 A reveal that this close contact occurs in more than 80% of cases. This carbonyl-carbonyl interaction offers an energetic sabilization for the partially allowed conformations of asparagine and aspartic acid with respect to all other non-glycyl amino acids. The non-covalent attractive interactions between the dipoles of two carbonyls has recently been calculated to have an energy comparable to that of a hydrogen bond. The preponderance of asparagine in the left-handed alpha-helical region, and in general of aspartic acid and asparagine in the partially allowed regions of the Ramachandran plot, may be a consequence of this carbonyl-carbonyl stacking interaction.     

锦屏水电站植被数量分类与排序

以锦屏梯级水电站大河湾段植被为研究对象,以河流流向为基础布置样方,将50个样方物种重要值及其环境因子进行TWINSPAN数量分类和DCA排序,结果显示:(1)由TWINSPAN数量分类将植被划分为7种类型,雅致雾水葛(Pouzolzia elegans)灌草群落为该区域的主要群落.(2)以TWINSPAN数量分类为基础的DCA二维排序图将样方分为6个生态类型,DCA 二维排序图横轴从左至右海拔高度逐渐降低;排序纵轴由下至上坡度逐渐升高.(3)以TWINSPAN数量分类在DCA排序图中有明确的分布范围和界限,DCA排序与TWINSPAN数量分类结合使用,能很好反映群落与环境因子间的生态关系.     

Efficient in silico exploration of RNA interhelical conformations using Euler angles and WExplore

HIV-1 TAR RNA is a two-helix bulge motif that plays a critical role in HIV viral replication and is an important drug target. However, efforts at designing TAR inhibitors have been challenged by its high degree of structural flexibility, which includes slow large-amplitude reorientations of its helices with respect to one another. Here, we use the recently introduced algorithm WExplore in combination with Euler angles to achieve unprecedented sampling of the TAR conformational ensemble. Our ensemble achieves similar agreement with experimental NMR data when compared with previous TAR computational studies, and is generated at a fraction of the computational cost. It clearly emerges from configuration space network analysis that the intermittent formation of the A22-U40 base pair acts as a reversible switch that enables sampling of interhelical conformations that would otherwise be topologically disallowed. We find that most previously determined ligand-bound structures are found in similar location in the network, and we use a sample-and-select approach to guide the construction of a set of novel conformations which can serve as the basis for future drug development efforts. Collectively, our findings demonstrate the utility of WExplore in combination with suitable order parameters as a method for exploring RNA conformational space.     

Assessment of grassland use intensity by remote sensing to support conservation schemes

Grassland is a land cover in the area of conflict between agriculture and conservation, where intensification of land use is a major threat to grassland biodiversity. Grassland use intensity is a key factor for the conservation value of grassland, and detailed spatial data on grassland use intensity is needed to improve strategies for biodiversity conservation. A new remote sensing-based approach using multi-temporal high resolution RapidEye satellite data was developed in the present study that makes a large-scale assessment of grassland use intensity possible. RapidEye is a constellation of five satellites with 6.5 m spatial resolution, which allows frequent and timely image acquisition targeted at specific growing seasons. Semi-natural grassland, extensively used grassland, intensively used grassland and tilled grassland could be reliably differentiated at the management plot level in a study area in southern Germany. Various combinations of images from different observation dates have been tested as classification input and their overall classification accuracies were validated by field data. Best results were achieved using a combination of five multi-temporal scenes with an overall accuracy of 85.7%. A three-scene combination resulted in an overall accuracy of 82.2%. The analysis showed that seasonal aspects are very important when selecting adequate observation dates. Grassland use intensity was also assessed on peatlands using a peat soil map, since land use intensity significantly affects greenhouse gas emissions from peatlands. The results demonstrate the potential of targeted multi-spectral, high spatial resolution remote sensing for the large-scale monitoring of dynamic habitats, which is of vital importance to support various environmental conservation schemes through improved monitoring and reporting capabilities.     

Introduction of a new scheme for classifying β-turns in protein structures

Protein β-turn classification remains an area of ongoing development in structural biology research. While the commonly used nomenclature defining type I, type II and type IV β-turns was introduced in the 1970s and 1980s, refinements of β-turn type definitions have been introduced as recently as 2019 by Dunbrack, Jr and co-workers who expanded the number of β-turn types to 18 (Shapovalov et al, PLOS Computat. Biol., 15, e1006844, 2019). Based on their analysis of 13 030 turns from 1074 ultrahigh resolution (≤1.2 Å) protein structures, they used a new clustering algorithm to expand the definitions used to classify protein β-turns and introduced a new nomenclature system. We recently encountered a specific problem when classifying β-turns in crystal structures of pentapeptide repeat proteins (PRPs) determined in our lab that are largely composed of β-turns that often lie close to, but just outside of, canonical β-turn regions. To address this problem, we devised a new scheme that merges the Klyne-Prelog stereochemistry nomenclature and definitions with the Ramachandran plot. The resulting Klyne-Prelog-modified Ramachandran plot scheme defines 1296 distinct potential β-turn classifications that cover all possible protein β-turn space with a nomenclature that indicates the stereochemistry of i + 1 and i + 2 backbone dihedral angles. The utility of the new classification scheme was illustrated by re-classification of the β-turns in all known protein structures in the PRP superfamily and further assessed using a database of 16 657 high-resolution protein structures (≤1.5 Å) from which 522 776 β-turns were identified and classified.