13Check_RNA: a tool to evaluate 13C chemical shift assignments of RNA.

Motivation: Chemical shifts (CS) are an important source of structural information of macromolecules such as RNA. In addition to the scarce availability of CS for RNA, the observed values are prone to errors due to a wrong re-calibration or miss assignments. Different groups have dedicated their efforts to correct CS systematic errors on RNA. Despite this, there are not automated and freely available algorithms for evaluating the referencing of RNA 13 C CS before their deposition to the BMRB or re-reference already deposited CS with systematic errors. Results: Based on an existent method we have implemented an open source python module to correct 13 C CS (from here on 13Cexp) systematic errors of RNAs and then return the results in 3 formats including the nmrstar one. Availability and implementation: This software is available on GitHub at https://github.com/BIOS-IMASL/13Check_RNA under a MIT license. Supplementary information: Supplementary data are available at Bioinformatics online.

Finite Dimension: A Mathematical Tool to Analise Glycans.

There is a need to develop widely applicable tools to understand glycan organization, diversity and structure. We present a graph-theoretical study of a large sample of glycans in terms of finite dimension, a new metric which is an adaptation to finite sets of the classical Hausdorff "fractal" dimension. Every glycan in the sample is encoded, via finite dimension, as a point of Glycan Space, a new notion introduced in this paper. Two major outcomes were found: (a) the existence of universal bounds that restrict the universe of possible glycans and show, for instance, that the graphs of glycans are a very special type of chemical graph, and (b) how Glycan Space is related to biological domains associated to the analysed glycans. In addition, we discuss briefly how this encoding may help to improve search in glycan databases.

Electrostatic unfolding and interactions of albumin driven by pH changes: a molecular dynamics study.

A better understanding of protein aggregation is bound to translate into critical advances in several areas, including the treatment of misfolded protein disorders and the development of self-assembling biomaterials for novel commercial applications. Because of its ubiquity and clinical potential, albumin is one of the best-characterized models in protein aggregation research; but its properties in different conditions are not completely understood. Here, we carried out all-atom molecular dynamics simulations of albumin to understand how electrostatics can affect the conformation of a single albumin molecule just prior to self-assembly. We then analyzed the tertiary structure and solvent accessible surface area of albumin after electrostatically triggered partial denaturation. The data obtained from these single protein simulations allowed us to investigate the effect of electrostatic interactions between two proteins. The results of these simulations suggested that hydrophobic attractions and counterion binding may be strong enough to effectively overcome the electrostatic repulsions between the highly charged monomers. This work contributes to our general understanding of protein aggregation mechanisms, the importance of explicit consideration of free ions in protein solutions, provides critical new insights about the equilibrium conformation of albumin in its partially denatured state at low pH, and may spur significant progress in our efforts to develop biocompatible protein hydrogels driven by electrostatic partial denaturation.

Physics-based protein-structure prediction using a hierarchical protocol based on the UNRES force field: assessment in two blind tests.

Recent improvements in the protein-structure prediction method developed in our laboratory, based on the thermodynamic hypothesis, are described. The conformational space is searched extensively at the united-residue level by using our physics-based UNRES energy function and the conformational space annealing method of global optimization. The lowest-energy coarse-grained structures are then converted to an all-atom representation and energy-minimized with the ECEPP/3 force field. The procedure was assessed in two recent blind tests of protein-structure prediction. During the first blind test, we predicted large fragments of alpha and alpha+beta proteins [60-70 residues with C(alpha) rms deviation (rmsd) <6 A]. However, for alpha+beta proteins, significant topological errors occurred despite low rmsd values. In the second exercise, we predicted whole structures of five proteins (two alpha and three alpha+beta, with sizes of 53-235 residues) with remarkably good accuracy. In particular, for the genomic target TM0487 (a 102-residue alpha+beta protein from Thermotoga maritima), we predicted the complete, topologically correct structure with 7.3-A C(alpha) rmsd. So far this protein is the largest alpha+beta protein predicted based solely on the amino acid sequence and a physics-based potential-energy function and search procedure. For target T0198, a phosphate transport system regulator PhoU from T. maritima (a 235-residue mainly alpha-helical protein), we predicted the topology of the whole six-helix bundle correctly within 8 A rmsd, except the 32 C-terminal residues, most of which form a beta-hairpin. These and other examples described in this work demonstrate significant progress in physics-based protein-structure prediction.

The protein folding problem: global optimization of the force fields.

The evolutionary development of a theoretical approach to the protein folding problem, in our laboratory, is traced. The theoretical foundations and the development of a suitable empirical all-atom potential energy function and a global optimization search are examined. Whereas the all-atom approach has thus far succeeded for relatively small molecules and for alpha-helical proteins containing up to 46 residues, it has been necessary to develop a hierarchical approach to treat larger proteins. In the hierarchical approach to single- and multiple-chain proteins, global optimization is carried out for a simplified united residue (UNRES) description of a polypeptide chain to locate the region in which the global minimum lies. Conversion of the UNRES structures in this region to all-atom structures is followed by a local search in this region. The performance of this approach in successive CASP blind tests for predicting protein structure by an ab initio physics-based method is described. Finally, a recent attempt to compute a folding pathway is discussed.

Influence of lysine content and pH on the stability of alanine-based copolypeptides.

To account for the relative contributions of lysine and alanine residues to the stability of alpha-helices of copolymers of these two residues, conformational energy calculations were carried out for several hexadecapeptides at several pHs. All the calculations considered explicitly the coupling between the conformation of the molecule and the ionization equilibria as a function of pH. The total free energy function used in these calculations included terms that account for the solvation free energy and free energy of ionization. These terms were evaluated by means of a fast multigrid boundary element method. Reasonable agreement with experimental values was obtained for the helix contents and vicinal coupling constants ((3)J(HNalpha)). The helix contents were found to depend strongly on the lysine content, in agreement with recent experimental results of Williams et al. (Journal of the American Chemical Society, 1998, Vol. 120, pp. 11033-11043) In the lowest energy conformation computed for a hexadecapeptide containing 3 lysine residues at pH 6, the lysine side chains are preferentially hydrated; this decreases the hydration of the backbone CO and NH groups, thereby forcing the latter to form hydrogen bonds with each other in the helical conformation. The lowest energy conformation computed for a hexadecapeptide containing 6 lysine residues at pH 6 shows a close proximity between the NH3(+) groups of the lysine side chains, a feature that was previously observed in calculations of short alanine-based oligopeptides. The calculation on a blocked 16-mer of alanine shows a 7% helix content based on the Boltzmann averaged vicinal coupling constants computed from the dihedral angles phi, consistent with previous experimental evidence on triblock copolymers containing a central block of alanines, and with earlier theoretical calculations.

Physical reasons for the unusual alpha-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides.

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH(2) in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge-charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent epsilon of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the alpha-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

A new approach for TU complex characterization.

In this paper, we present a new TU complex detection and characterization algorithm that consists of two stages; the first is a mathematical modeling of the electrocardiographic segment after QRS complex; the second uses classic threshold comparison techniques, over the signal and its first and second derivatives, to determine the significant points of each wave. Later, both T and U waves are morphologically classified. Amongst the principal innovations of this algorithm is the inclusion of U-wave characterization and a mathematical modeling stage, that avoids many of the problems of classic techniques when there is a low signal-to-noise ratio or when wave morphology is atypical. The results of the algorithm validation with the recently appeared QT database are also shown. For T waves these results are better when compared to other existing algorithms. U-wave results cannot be contrasted with other algorithms as, to our knowledge, none are available. Examples showing the causes of principal discrepancies between our algorithm and the QT database annotations are also given, and some ways of attempting to improve and benefit from the proposed algorithm are suggested.

On the pH-conformational dependence of the unblocked SYPYD peptide.

Simulations were carried out for an unblocked pentapeptide with the sequence Ser-Tyr-Pro-Tyr-Asp (SYPYD) with explicit consideration of the coupling between the conformation of the molecule and the ionization equilibria at a given pH. The available NMR experimental data indicate a high preference for the cis isomeric turn-like form of Tyr-Pro at intermediate pH (approximately 6) and a destabilization of the cis form at both high (approximately 9) and low (approximately 3) pH. In order to identify the source of the stability of the conformation of this pentapeptide as a function of pH, Monte Carlo simulations were used to generate an ensemble of low-energy conformations at different pH values (viz. 3, 6 and 9). The total free energy function used in these calculations includes terms that account for the solvation free energy and free energy of ionization. These terms are evaluated by means of a fast multigrid boundary element (MBE) method. In good qualitative agreement with the experiments, our results indicate that the Boltzmann averaged population of the cis isomeric form of the pentapeptide has a maximum (45 %) at pH 6 and is significantly smaller (25 % and 23 %) for higher and lower pH values, respectively, following the trend of the experimental data. Also, the degree of charge for the lowest-energy conformations, as well as the contribution of electrostatic interactions to the stability of the preferred conformations, vary widely at the different pH values. Different kinds of packing of the aromatic side-chains of Tyr2 and Tyr4 against the proline ring are observed at different pH values, indicating that their contribution to the stability of the low-energy conformations is also pH-dependent. In summary, our results provide a basis for discussing the nature of the interactions that stabilize turn-like conformations of the peptide SYPYD as a function of pH.

Role of hydrophobicity and solvent-mediated charge-charge interactions in stabilizing alpha-helices.

A theoretical study to identify the conformational preferences of lysine-based oligopeptides has been carried out. The solvation free energy and free energy of ionization of the oligopeptides have been calculated by using a fast multigrid boundary element method that considers the coupling between the conformation of the molecule and the ionization equilibria explicitly, at a given pH value. It has been found experimentally that isolated alanine and lysine residues have somewhat small intrinsic helix-forming tendencies; however, results from these simulations indicate that conformations containing right-handed alpha-helical turns are energetically favorable at low values of pH for lysine-based oligopeptides. Also, unusual patterns of interactions among lysine side chains with large hydrophobic contacts and close proximity (5-6 A) between charged NH3+ groups are observed. Similar arrangements of charged groups have been seen for lysine and arginine residues in experimentally determined structures of proteins available from the Protein Data Bank. The lowest-free-energy conformation of the sequence Ac-(LYS)6-NMe from these simulations showed large pKalpha shifts for some of the NH3+ groups of the lysine residues. Such large effects are not observed in the lowest-energy conformations of oligopeptide sequences with two, three, or four lysine residues. Calculations on the sequence Ac-LYS-(ALA)4-LYS-NMe also reveal low-energy alpha-helical conformations with interactions of one of the LYS side chains with the helix backbone in an arrangement quite similar to the one described recently by (Proc. Natl. Acad. Sci. U.S.A. 93:4025-4029). The results of this study provide a sound basis with which to discuss the nature of the interactions, such as hydrophobicity, charge-charge interaction, and solvent polarization effects, that stabilize right-handed alpha-helical conformations.

Coupling between folding and ionization equilibria: effects of pH on the conformational preferences of polypeptides.

A new approach to the conformational study of polypeptides is presented. It considers explicitly the coupling between the conformation of the molecule and the ionization equilibria at a given pH value. Calculations of the solvation free energy and free energy of ionization of a 17-residue polypeptide are carried out using a fast multigrid boundary element method (MBE). The MBE method uses an adaptive tessellation of the molecular surface by boundary elements with non-regular size to solve the Poisson equation rapidly, and with a high degree of accuracy. The MBE method is integrated into the ECEPP (Empirical Conformational Energy Program for Peptides) algorithm to compute the coupling between the ionization state and the conformation of the molecule. This approach has been applied to study the conformational preference of a short polypeptide for which the available NMR and CD experimental data indicate that conformations containing a right-handed alpha-helical segment are energetically more favorable at low values of pH. The results of calculations using the present method agree quite well with experiments, in contrast to previous applications with standard techniques (using pre-assigned charges at each pH) that were not able to reproduce the experimental findings. Also, it is shown how the coupling to the conformation leads to different degrees of ionization of a given type of residue, for example glutamic acid, at different positions in the amino acid sequence, at any given pH. The results of this study provide a sound basis to discuss the origin of the stability of polypeptide conformations, and its dependence on the environmental conditions.