Prediction of Optimal Folding Routes of Proteins That Satisfy the Principle of Lowest Entropy Loss: Dynamic Contact Maps and Optimal Control

An optimization model is introduced in which proteins try to evade high energy regions of the folding landscape, and prefer low entropy loss routes during folding. We make use of the framework of optimal control whose convenient solution provides practical and useful insight into the sequence of events during folding. We assume that the native state is available. As the protein folds, it makes different set of contacts at different folding steps. The dynamic contact map is constructed from these contacts. The topology of the dynamic contact map changes during the course of folding and this information is utilized in the dynamic optimization model. The solution is obtained using the optimal control theory. We show that the optimal solution can be cast into the form of a Gaussian Network that governs the optimal folding dynamics. Simulation results on three examples (CI2, Sso7d and Villin) show that folding starts by the formation of local clusters. Non-local clusters generally require the formation of several local clusters. Non-local clusters form cooperatively and not sequentially. We also observe that the optimal controller prefers “zipping” or small loop closure steps during folding. The folding routes predicted by the proposed method bear strong resemblance to the results in the literature

Exploring the Levinthal limit in protein folding

According to the thermodynamic hypothesis, the native state of proteins is uniquely defined by their amino acid sequence. On the other hand, according to Levinthal, the native state is just a local minimum of the free energy and a given amino acid sequence, in the same thermodynamic conditions, can assume many, very different structures that are as thermodynamically stable as the native state. This is the Levinthal limit explored in this work. Using computer simulations, we compare the interactions that stabilize the native state of four different proteins with those that stabilize three non-native states of each protein and find that the nature of the interactions is very similar for all such 16 conformers. Furthermore, an enhancement of the degree of fluctuation of the non-native conformers can be explained by an insufficient relaxation to their local free energy minimum. These results favor Levinthal's hypothesis that protein folding is a kinetic non-equilibrium process.FCT - Foundation for Science and Technology, Portugal [UID/Multi/04326/2013]; Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP); Conselho Nacional de Desenvolvimento Cientia co e Tecnologico (CNPq

A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry

We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical derivation of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications. Moreover, much material that is found scattered in the literature has been put together here to facilitate comprehension and to serve as a handy reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and subeqn package

Prevalence of bullying and aggressive behavior and their relationship to mental health problems among 12- to 15-year-old Norwegian adolescents

The aim of this study was to examine the relationships between being bullied and aggressive behavior and self-reported mental health problems among young adolescents. A representative population sample of 2,464 young Norwegian adolescents (50.8% girls) aged 12–15 years was assessed. Being bullied was measured using three items concerning teasing, exclusion, and physical assault. Self-esteem was assessed by Harter’s self-perception profile for adolescents. Emotional and behavioral problems were measured by the Moods and Feelings Questionnaire (MFQ) and the youth self-report (YSR). Aggressive behavior was measured by four items from the YSR. One-tenth of the adolescents reported being bullied, and 5% reported having been aggressive toward others during the past 6 months. More of the students being bullied and students being aggressive toward others reported parental divorce, and they showed higher scores on all YSR subscales and on the MFQ questions, and lower scores on the global self-worth subscale (Harter) than students not being bullied or aggressive. A few differences emerged between the two groups being bullied or being aggressive toward others: those who were aggressive showed higher total YSR scores, higher aggression and delinquency scores, and lower social problems scores, and reported higher scores on the social acceptance subscale (Harter) than bullied students. However, because social problems were demonstrated in both the involved groups, interventions designed to improve social competence and interaction skills should be integrated in antibullying programs

Combining Optimal Control Theory and Molecular Dynamics for Protein Folding

A new method to develop low-energy folding routes for proteins is presented. The novel aspect of the proposed approach is the synergistic use of optimal control theory with Molecular Dynamics (MD). In the first step of the method, optimal control theory is employed to compute the force field and the optimal folding trajectory for the atoms of a Coarse-Grained (CG) protein model. The solution of this CG optimization provides an harmonic approximation of the true potential energy surface around the native state. In the next step CG optimization guides the MD simulation by specifying the optimal target positions for the atoms. In turn, MD simulation provides an all-atom conformation whose positions match closely the reference target positions determined by CG optimization. This is accomplished by Targeted Molecular Dynamics (TMD) which uses a bias potential or harmonic restraint in addition to the usual MD potential. Folding is a dynamical process and as such residues make different contacts during the course of folding. Therefore CG optimization has to be reinitialized and repeated over time to accomodate these important changes. At each sampled folding time, the active contacts among the residues are recalculated based on the all-atom conformation obtained from MD. Using the new set of contacts, the CG potential is updated and the CG optimal trajectory for the atoms is recomputed. This is followed by MD. Implementation of this repetitive CG optimization - MD simulation cycle generates the folding trajectory. Simulations on a model protein Villin demonstrate the utility of the method. Since the method is founded on the general tools of optimal control theory and MD without any restrictions, it is widely applicable to other systems. It can be easily implemented with available MD software packages

An Estimate of the Numbers and Density of Low-Energy Structures (or Decoys) in the Conformational Landscape of Proteins

The conformational energy landscape of a protein, as calculated by known potential energy functions, has several minima, and one of these corresponds to its native structure. It is however difficult to comprehensively estimate the actual numbers of low energy structures (or decoys), the relationships between them, and how the numbers scale with the size of the protein.We have developed an algorithm to rapidly and efficiently identify the low energy conformers of oligo peptides by using mutually orthogonal Latin squares to sample the potential energy hyper surface. Using this algorithm, and the ECEPP/3 potential function, we have made an exhaustive enumeration of the low-energy structures of peptides of different lengths, and have extrapolated these results to larger polypeptides.We show that the number of native-like structures for a polypeptide is, in general, an exponential function of its sequence length. The density of these structures in conformational space remains more or less constant and all the increase appears to come from an expansion in the volume of the space. These results are consistent with earlier reports that were based on other models and techniques

IFN-λ3, not IFN-λ4, likely mediates IFNL3-IFNL4 haplotype-dependent hepatic inflammation and fibrosis

Genetic variation in the IFNL3-IFNL4 (interferon-λ3-interferon-λ4) region is associated with hepatic inflammation and fibrosis. Whether IFN-λ3 or IFN-λ4 protein drives this association is not known. We demonstrate that hepatic inflammation, fibrosis stage, fibrosis progression rate, hepatic infiltration of immune cells, IFN-λ3 expression, and serum sCD163 levels (a marker of activated macrophages) are greater in individuals with the IFNL3-IFNL4 risk haplotype that does not produce IFN-λ4, but produces IFN-λ3. No difference in these features was observed according to genotype at rs117648444, which encodes a substitution at position 70 of the IFN-λ4 protein and reduces IFN-λ4 activity, or between patients encoding functionally defective IFN-λ4 (IFN-λ4-Ser70) and those encoding fully active IFN-λ4-Pro70. The two proposed functional variants (rs368234815 and rs4803217) were not superior to the discovery SNP rs12979860 with respect to liver inflammation or fibrosis phenotype. IFN-λ3 rather than IFN-λ4 likely mediates IFNL3-IFNL4 haplotype-dependent hepatic inflammation and fibrosis

How university’s activities support the development of students’ entrepreneurial abilities: case of Slovenia and Croatia

The paper reports how the offered university activities support the development of students’ entrepreneurship abilities. Data were collected from 306 students from Slovenian and 609 students from Croatian universities. The study reduces the gap between theoretical researches about the academic entrepreneurship education and individual empirical studies about the student’s estimation of the offered academic activities for development of their entrepreneurial abilities. The empirical research revealed differences in Slovenian and Croatian students’ perception about (a) needed academic activities and (b) significance of the offered university activities, for the development of their entrepreneurial abilities. Additionally, the results reveal that the impact of students’ gender and study level on their perception about the importance of the offered academic activities is not significant for most of the considered activities. The main practical implication is focused on further improvement of universities’ entrepreneurship education programs through selection and utilization of activities which can fill in the recognized gaps between the students’ needed and the offered academic activities for the development of students’ entrepreneurial abilities

A review of estimation of distribution algorithms in bioinformatics

Evolutionary search algorithms have become an essential asset in the algorithmic toolbox for solving high-dimensional optimization problems in across a broad range of bioinformatics problems. Genetic algorithms, the most well-known and representative evolutionary search technique, have been the subject of the major part of such applications. Estimation of distribution algorithms (EDAs) offer a novel evolutionary paradigm that constitutes a natural and attractive alternative to genetic algorithms. They make use of a probabilistic model, learnt from the promising solutions, to guide the search process. In this paper, we set out a basic taxonomy of EDA techniques, underlining the nature and complexity of the probabilistic model of each EDA variant. We review a set of innovative works that make use of EDA techniques to solve challenging bioinformatics problems, emphasizing the EDA paradigm's potential for further research in this domain

Examination of the efficacy of acute L-alanyl-L-glutamine ingestion during hydration stress in endurance exercise

<p>Abstract</p> <p>Background</p> <p>The effect of acute L-alanyl-L-glutamine (AG; Sustamine™) ingestion on performance changes and markers of fluid regulation, immune, inflammatory, oxidative stress, and recovery was examined in response to exhaustive endurance exercise, during and in the absence of dehydration.</p> <p>Methods</p> <p>Ten physically active males (20.8 ± 0.6 y; 176.8 ± 7.2 cm; 77.4 ± 10.5 kg; 12.3 ± 4.6% body fat) volunteered to participate in this study. During the first visit (T1) subjects reported to the laboratory in a euhydrated state to provide a baseline (BL) blood draw and perform a maximal exercise test. In the four subsequent randomly ordered trials, subjects dehydrated to -2.5% of their baseline body mass. For T2, subjects achieved their goal weight and were not rehydrated. During T3 - T5, subjects reached their goal weight and then rehydrated to 1.5% of their baseline body mass by drinking either water (T3) or two different doses (T4 and T5) of the AG supplement (0.05 g·kg^-1and 0.2 g·kg^-1, respectively). Subjects then exercised at a workload that elicited 75% of their VO₂max on a cycle ergometer. During T2 - T5 blood draws occurred once goal body mass was achieved (DHY), immediately prior to the exercise stress (RHY), and immediately following the exercise protocol (IP). Resting 24 hour (24P) blood samples were also obtained. Blood samples were analyzed for glutamine, potassium, sodium, aldosterone, arginine vasopressin (AVP), C-reactive protein (CRP), interleukin-6 (IL-6), malondialdehyde (MDA), testosterone, cortisol, ACTH, growth hormone and creatine kinase. Statistical evaluation of performance, hormonal and biochemical changes was accomplished using a repeated measures analysis of variance.</p> <p>Results</p> <p>Glutamine concentrations for T5 were significantly higher at RHY and IP than T2 - T4. When examining performance changes (difference between T2 - T5 and T1), significantly greater times to exhaustion occurred during T4 (130.2 ± 340.2 sec) and T5 (157.4 ± 263.1 sec) compared to T2 (455.6 ± 245.0 sec). Plasma sodium concentrations were greater (p < 0.05) at RHY and IP for T2 than all other trials. Aldosterone concentrations at RHY and IP were significantly lower than that at BL and DHY. AVP was significantly elevated at DHY, RHY and IP compared to BL measures. No significant differences were observed between trials in CRP, IL-6, MDA, or in any of the other hormonal or biochemical measures.</p> <p>Conclusion</p> <p>Results demonstrate that AG supplementation provided a significant ergogenic benefit by increasing time to exhaustion during a mild hydration stress. This ergogenic effect was likely mediated by an enhanced fluid and electrolyte uptake.</p