Chemotypes of Achillea millefolium transferred from 14 different locations in Lithuania to the controlled environment

The composition of essential oils hydrodistilled from 19 samples of inflorescences and leaves of Achillea millefolium L. plants, which were transferred from 14 natural habitats in Lithuania to the field collection, is reported. Total content of oil was 0.15–0.55% in inflorescences and 0.06–0.19% (v/w) in leaves. In total 117 compounds were identified positively or tentatively. Data obtained clearly indicate the presence of a remarkable chemical polymorphism within the population of A. millefolium in Lithuania. The content of the major constituents in the oils from inflorescences varied in the following ranges: β-pinene, 0.33–62.29%; β-myrcene, 0.05–69.76%; α-phelandrene, 0.13–29.96%; 1,8-cineole, 2.30–21.57%; and chamazulene, 0.08–30.70%. According to the major components the essential oils' six chemotypes of A. millefolium were defined.     

Optimization of electrostatic interactions in protein-protein complexes

In this article, we present a statistical analysis of the electrostatic properties of 298 protein-protein complexes and 356 domain-domain structures extracted from the previously developed database of protein complexes (ProtCom, http://www.ces.clemson.edu/compbio/protcom). For each structure in the dataset we calculated the total electrostatic energy of the binding and its two components, Coulombic and reaction field energy. It was found that in a vast majority of the cases (>90%), the total electrostatic component of the binding energy was unfavorable. At the same time, the Coulombic component of the binding energy was found to favor the complex formation while the reaction field component of the binding energy opposed the binding. It was also demonstrated that the components in a wild-type (WT) structure are optimized/anti-optimized with respect to the corresponding distributions, arising from random shuffling of the charged side chains. The degree of this optimization was assessed through the Z-score of WT energy in respect to the random distribution. It was found that the Z-scores of Coulombic interactions peak at a considerably negative value for all 654 cases considered while the Z-score of the reaction field energy varied among different types of complexes. All these findings indicate that the Coulombic interactions within WT protein-protein complexes are optimized to favor the complex formation while the total electrostatic energy predominantly opposes the binding. This observation was used to discriminate WT structures among sets of structural decoys and showed that the electrostatic component of the binding energy is not a good discriminator of the WT; while, Coulombic or reaction field energies perform better depending upon the decoy set used.     

Structural quality of unrefined models in protein docking

Structural characterization of protein‐protein interactions is essential for understanding life processes at the molecular level. However, only a fraction of protein interactions have experimentally resolved structures. Thus, reliable computational methods for structural modeling of protein interactions (protein docking) are important for generating such structures and understanding the principles of protein recognition. Template‐based docking techniques that utilize structural similarity between target protein‐protein interaction and cocrystallized protein‐protein complexes (templates) are gaining popularity due to generally higher reliability than that of the template‐free docking. However, the template‐based approach lacks explicit penalties for intermolecular penetration, as opposed to the typical free docking where such penalty is inherent due to the shape complementarity paradigm. Thus, template‐based docking models are commonly assumed to require special treatment to remove large structural penetrations. In this study, we compared clashes in the template‐based and free docking of the same proteins, with crystallographically determined and modeled structures. The results show that for the less accurate protein models, free docking produces fewer clashes than the template‐based approach. However, contrary to the common expectation, in acceptable and better quality docking models of unbound crystallographically determined proteins, the clashes in the template‐based docking are comparable to those in the free docking, due to the overall higher quality of the template‐based docking predictions. This suggests that the free docking refinement protocols can in principle be applied to the template‐based docking predictions as well. Proteins 2016; 85:39–45. © 2016 Wiley Periodicals, Inc.     

Charge sequence coding in statistical modeling of unfolded proteins

Unfolded proteins recently attracted attention due to accumulation of experimental evidences for their significant role in different life processes. Modeling of electrostatic interactions (EI) in unfolded state of proteins is becoming increasingly important as well. In this paper, we stress on the importance of how the sequence of charged residues of a given protein is incorporated into the models for calculation of EI in the unfolded state. On the basis of the distributions of distances between titratable sites of charged residues calculated for polypeptide chains of various compositions, it was found that the distance distribution for a pair of residues, located close to each other along the sequence of a protein, depends on what residues constitute the pair in question. It was concluded that the consideration of these residue-specific distributions is essential for a statistical model to be accurate from the physical point of view. It was suggested that use of distance intervals in the spherical model of unfolded proteins accounts better for the charge sequence than the set of single distance values. This was illustrated by comparison of the pK values of the titratable groups of the unfolded N-terminal SH3 domain of the Drosophila protein drk to the available experimental data.     

How to choose templates for modeling of protein complexes: Insights from benchmarking template-based docking

Comparative docking is based on experimentally determined structures of protein-protein complexes (templates), following the paradigm that proteins with similar sequences and/or structures form similar complexes. Modeling utilizing structure similarity of target monomers to template complexes significantly expands structural coverage of the interactome. Template-based docking by structure alignment can be performed for the entire structures or by aligning targets to the bound interfaces of the experimentally determined complexes. Systematic benchmarking of docking protocols based on full and interface structure alignment showed that both protocols perform similarly, with top 1 docking success rate 26%. However, in terms of the models' quality, the interface-based docking performed marginally better. The interface-based docking is preferable when one would suspect a significant conformational change in the full protein structure upon binding, for example, a rearrangement of the domains in multidomain proteins. Importantly, if the same structure is selected as the top template by both full and interface alignment, the docking success rate increases 2-fold for both top 1 and top 10 predictions. Matching structural annotations of the target and template proteins for template detection, as a computationally less expensive alternative to structural alignment, did not improve the docking performance. Sophisticated remote sequence homology detection added templates to the pool of those identified by structure-based alignment, suggesting that for practical docking, the combination of the structure alignment protocols and the remote sequence homology detection may be useful in order to avoid potential flaws in generation of the structural templates library.     

2D finite element ecological model for the Curonian lagoon

The results of application of 2D finite element model SHYFEM to the Curonian lagoon (Baltic Sea) are considered. SHYFEM consist of a physical processes module and an eutrophication module EUTRO adapted for the SHYFEM code from well known modelling system WASP. The SHYFEM/EUTRO model calibration results were compared with the performance of various biogeochemical models analysed in other studies (153 studies published from 1990 to 2002). The performance of all corresponding state variables—dissolved oxygen, NO₃, NH₄, PO₄, phyto- and zooplankton—was slightly lower than median model performance which could be considered satisfactory given the initial state of model formulation and calibration. Model underestimates phytoplankton autumn blooms, especially for the southern part of the lagoon, where fine sediments dominate and water residence time is high. It can be concluded that, in order to increase model performance, the eutrophication module should be improved to account for the dominance of different phytoplankton groups as well as for the exchanges between the sediments and the water column. The amount and quality of the data available for the model setup and calibration are unsatisfactory and should be improved for the development of the next enhanced model version. Guest editors: A. Razinkovas, Z. R. Gasiūnaitė, J. M. Zaldivar & P. Viaroli European Lagoons and their Watersheds: Function and Biodiversity     

Decline of Phosphotransfer and Substrate Supply Metabolic Circuits Hinders ATP Cycling in Aging Myocardium

Integration of mitochondria with cytosolic ATP-consuming/ATP-sensing and substrate supply processes is critical for muscle bioenergetics and electrical activity. Whether age-dependent muscle weakness and increased electrical instability depends on perturbations in cellular energetic circuits is unknown. To define energetic remodeling of aged atrial myocardium we tracked dynamics of ATP synthesis-utilization, substrate supply, and phosphotransfer circuits through adenylate kinase (AK), creatine kinase (CK), and glycolytic/glycogenolytic pathways using ¹⁸O stable isotope-based phosphometabolomic technology. Samples of intact atrial myocardium from adult and aged rats were subjected to ¹⁸O-labeling procedure at resting basal state, and analyzed using the ¹⁸O-assisted HPLC-GC/MS technique. Characteristics for aging atria were lower inorganic phosphate Pi[¹⁸O], γ-ATP[¹⁸O], β-ADP[¹⁸O], and creatine phosphate CrP[¹⁸O] ¹⁸O-labeling rates indicating diminished ATP utilization-synthesis and AK and CK phosphotransfer fluxes. Shift in dynamics of glycolytic phosphotransfer was reflected in the diminished G6P[¹⁸O] turnover with relatively constant glycogenolytic flux or G1P[¹⁸O] ¹⁸O-labeling. Labeling of G3P[¹⁸O], an indicator of G3P-shuttle activity and substrate supply to mitochondria, was depressed in aged myocardium. Aged atrial myocardium displayed reduced incorporation of ¹⁸O into second (¹⁸O₂), third (¹⁸O₃), and fourth (¹⁸O₄) positions of Pi[¹⁸O] and a lower Pi[¹⁸O]/γ-ATP[¹⁸ O]-labeling ratio, indicating delayed energetic communication and ATP cycling between mitochondria and cellular ATPases. Adrenergic stress alleviated diminished CK flux, AK catalyzed β-ATP turnover and energetic communication in aging atria. Thus, ¹⁸O-assisted phosphometabolomics uncovered simultaneous phosphotransfer through AK, CK, and glycolytic pathways and G3P substrate shuttle deficits hindering energetic communication and ATP cycling, which may underlie energetic vulnerability of aging atrial myocardium.     

Structural templates for comparative protein docking

Structural characterization of protein‐protein interactions is important for understanding life processes. Because of the inherent limitations of experimental techniques, such characterization requires computational approaches. Along with the traditional protein‐protein docking (free search for a match between two proteins), comparative (template‐based) modeling of protein‐protein complexes has been gaining popularity. Its development puts an emphasis on full and partial structural similarity between the target protein monomers and the protein‐protein complexes previously determined by experimental techniques (templates). The template‐based docking relies on the quality and diversity of the template set. We present a carefully curated, nonredundant library of templates containing 4950 full structures of binary complexes and 5936 protein‐protein interfaces extracted from the full structures at 12 Å distance cut‐off. Redundancy in the libraries was removed by clustering the PDB structures based on structural similarity. The value of the clustering threshold was determined from the analysis of the clusters and the docking performance on a benchmark set. High structural quality of the interfaces in the template and validation sets was achieved by automated procedures and manual curation. The library is included in the Dockground resource for molecular recognition studies at http://dockground.bioinformatics.ku.edu . Proteins 2015; 83:1563–1570. © 2014 Wiley Periodicals, Inc.     

Protein Docking by the Interface Structure Similarity: How Much Structure Is Needed?

The increasing availability of co-crystallized protein-protein complexes provides an opportunity to use template-based modeling for protein-protein docking. Structure alignment techniques are useful in detection of remote target-template similarities. The size of the structure involved in the alignment is important for the success in modeling. This paper describes a systematic large-scale study to find the optimal definition/size of the interfaces for the structure alignment-based docking applications. The results showed that structural areas corresponding to the cutoff values <12 Å across the interface inadequately represent structural details of the interfaces. With the increase of the cutoff beyond 12 Å, the success rate for the benchmark set of 99 protein complexes, did not increase significantly for higher accuracy models, and decreased for lower-accuracy models. The 12 Å cutoff was optimal in our interface alignment-based docking, and a likely best choice for the large-scale (e.g., on the scale of the entire genome) applications to protein interaction networks. The results provide guidelines for the docking approaches, including high-throughput applications to modeled structures.     

GWIDD: a comprehensive resource for genome-wide structural modeling of protein-protein interactions

We live in an age of access to more information than ever before. This can be a double-edged sword. Increased access to information allows for more informed and empowered researchers, while information overload becomes an increasingly serious risk. Thus, there is a need for intelligent information retrieval systems that can summarize relevant and reliable textual sources to satisfy a user's query. Question answering is a specialized type of information retrieval with the aim of returning precise short answers to queries posed as natural language questions. We present a review and comparison of three biomedical question answering systems: askHERMES (http://www.askhermes.org/), EAGLi (http://eagl.unige.ch/EAGLi/), and HONQA (http://services.hon.ch/cgi-bin/QA10/qa.pl).