Potentials of Mean Force for Protein Structure Prediction Vindicated, Formalized and Generalized

Understanding protein structure is of crucial importance in science, medicine and biotechnology. For about two decades, knowledge based potentials based on pairwise distances -- so-called "potentials of mean force" (PMFs) -- have been center stage in the prediction and design of protein structure and the simulation of protein folding. However, the validity, scope and limitations of these potentials are still vigorously debated and disputed, and the optimal choice of the reference state -- a necessary component of these potentials -- is an unsolved problem. PMFs are loosely justified by analogy to the reversible work theorem in statistical physics, or by a statistical argument based on a likelihood function. Both justifications are insightful but leave many questions unanswered. Here, we show for the first time that PMFs can be seen as approximations to quantities that do have a rigorous probabilistic justification: they naturally arise when probability distributions over different features of proteins need to be combined. We call these quantities reference ratio distributions deriving from the application of the reference ratio method. This new view is not only of theoretical relevance, but leads to many insights that are of direct practical use: the reference state is uniquely defined and does not require external physical insights; the approach can be generalized beyond pairwise distances to arbitrary features of protein structure; and it becomes clear for which purposes the use of these quantities is justified. We illustrate these insights with two applications, involving the radius of gyration and hydrogen bonding. In the latter case, we also show how the reference ratio method can be iteratively applied to sculpt an energy funnel. Our results considerably increase the understanding and scope of energy functions derived from known biomolecular structures

Knowledge-based energy functions for computational studies of proteins

This chapter discusses theoretical framework and methods for developing knowledge-based potential functions essential for protein structure prediction, protein-protein interaction, and protein sequence design. We discuss in some details about the Miyazawa-Jernigan contact statistical potential, distance-dependent statistical potentials, as well as geometric statistical potentials. We also describe a geometric model for developing both linear and non-linear potential functions by optimization. Applications of knowledge-based potential functions in protein-decoy discrimination, in protein-protein interactions, and in protein design are then described. Several issues of knowledge-based potential functions are finally discussed.Comment: 57 pages, 6 figures. To be published in a book by Springe

Early prevention of diabetes microvascular complications in people with hyperglycaemia in Europe. ePREDICE randomized trial. Study protocol, recruitment and selected baseline data

Objectives To assess the effects of early management of hyperglycaemia with antidiabetic drugs plus lifestyle intervention compared with lifestyle alone, on microvascular function in adults with pre-diabetes. Methods Trial design: International, multicenter, randomised, partially double-blind, placebo-controlled, clinical trial. Participants Males and females aged 45-74 years with IFG, IGT or IFG+IGT, recruited from primary care centres in Australia, Austria, Bulgaria, Greece, Kuwait, Poland, Serbia, Spain and Turkey. Intervention Participants were randomized to placebo; metformin 1.700 mg/day; linagliptin 5 mg/day or fixed-dose combination of linagliptin/metformin. All patients were enrolled in a lifestyle intervention program (diet and physical activity). Drug intervention will last 2 years. Primary Outcome: Composite end-point of diabetic retinopathy estimated by the Early Treatment Diabetic Retinopathy Study Score, urinary albumin to creatinine ratio, and skin conductance in feet estimated by the sudomotor index. Secondary outcomes in a subsample include insulin sensitivity, beta-cell function, biomarkers of inflammation and fatty liver disease, quality of life, cognitive function, depressive symptoms and endothelial function. Results One thousand three hundred ninety one individuals with hyperglycaemia were assessed for eligibility, 424 excluded after screening, 967 allocated to placebo, metformin, linagliptin or to fixed-dose combination of metformin + linagliptin. A total of 809 people (91.1%) accepted and initiated the assigned treatment. Study sample after randomization was well balanced among the four groups. No statistical differences for the main risk factors analysed were observed between those accepting or rejecting treatment initiation. At baseline prevalence of diabetic retinopathy was 4.2%, severe neuropathy 5.3% and nephropathy 5.7%. Conclusions ePREDICE is the first -randomized clinical trial with the aim to assess effects of different interventions (lifestyle and pharmacological) on microvascular function in people with prediabetes. The trial will provide novel data on lifestyle modification combined with glucose lowering drugs for the prevention of early microvascular complications and diabetes

PoPMuSiC 2.1: a web server for the estimation of protein stability changes upon mutation and sequence optimality

ABSTRACT:Journal ArticleSCOPUS: ar.jinfo:eu-repo/semantics/publishe

Inferring stabilizing mutations from protein phylogenies : application to influenza hemagglutinin

One selection pressure shaping sequence evolution is the requirement that a protein fold with sufficient stability to perform its biological functions. We present a conceptual framework that explains how this requirement causes the probability that a particular amino acid mutation is fixed during evolution to depend on its effect on protein stability. We mathematically formalize this framework to develop a Bayesian approach for inferring the stability effects of individual mutations from homologous protein sequences of known phylogeny. This approach is able to predict published experimentally measured mutational stability effects (ΔΔG values) with an accuracy that exceeds both a state-of-the-art physicochemical modeling program and the sequence-based consensus approach. As a further test, we use our phylogenetic inference approach to predict stabilizing mutations to influenza hemagglutinin. We introduce these mutations into a temperature-sensitive influenza virus with a defect in its hemagglutinin gene and experimentally demonstrate that some of the mutations allow the virus to grow at higher temperatures. Our work therefore describes a powerful new approach for predicting stabilizing mutations that can be successfully applied even to large, complex proteins such as hemagglutinin. This approach also makes a mathematical link between phylogenetics and experimentally measurable protein properties, potentially paving the way for more accurate analyses of molecular evolution

Chemical Synthesis of Staphyloferrin B Affords Insight into the Molecular Structure, Iron Chelation, and Biological Activity of a Polycarboxylate Siderophore Deployed by the Human Pathogen

Staphyloferrin B (SB) is a citrate-based polycarboxylate siderophore produced and utilized by the human pathogen Staphylococcus aureus for acquiring iron when colonizing the vertebrate host. The first chemical synthesis of SB is reported, which enables further molecular and biological characterization and provides access to structural analogues of the siderophore. Under conditions of iron limitation, addition of synthetic SB to bacterial growth medium recovered the growth of the antibiotic resistant community isolate S. aureus USA300 JE2. Two structural analogues of SB, epiSB and SBimide, were also synthesized and employed to investigate how epimerization of the citric acid moiety or imide formation influence its function as a siderophore. Epimerization of the citric acid stereocenter perturbed the iron-binding properties and siderophore function of SB as evidenced by experimental and computational modeling studies. Although epiSB provided growth recovery to S. aureus USA300 JE2 cultured in iron-deficient medium, the effect was attenuated relative to that of SB. Moreover, SB more effectively sequestered the Fe(III) bound to human holo-transferrin, an iron source of S. aureus, than epiSB. SBimide is an imide analogous to the imide forms of other citric acid siderophores that are often observed when these molecules are isolated from natural sources. Here, SBimide is shown to be unstable, converting to native SB at physiological pH. SB is considered to be a virulence factor of S. aureus, a pathogen that poses a particular threat to public health because of the number of drug-resistant strains emerging in hospital and community settings. Iron acquisition by S. aureus is important for its ability to colonize the human host and cause disease, and new chemical insights into the structure and function of SB will inform the search for new therapeutic strategies for combating S. aureus infections.Alfred Benzon Foundation (Postdoctoral fellowship)Pacific Southwest Regional Center of ExcellenceAlfred P. Sloan Foundatio

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

Limpet Shells from the Aterian Level 8 of El Harhoura 2 Cave (Témara, Morocco): Preservation State of Crossed-Foliated Layers

International audienceThe exploitation of mollusks by the first anatomically modern humans is a central question for archaeologists. This paper focuses on level 8 (dated around * 100 ka BP) of El Har-houra 2 Cave, located along the coastline in the Rabat-Témara region (Morocco). The large quantity of Patella sp. shells found in this level highlights questions regarding their origin and preservation. This study presents an estimation of the preservation status of these shells. We focus here on the diagenetic evolution of both the microstructural patterns and organic components of crossed-foliated shell layers, in order to assess the viability of further investigations based on shell layer minor elements, isotopic or biochemical compositions. The results show that the shells seem to be well conserved, with microstructural patterns preserved down to sub-micrometric scales, and that some organic components are still present in situ. But faint taphonomic degradations affecting both mineral and organic components are nonetheless evidenced, such as the disappearance of organic envelopes surrounding crossed-foliated lamellae, combined with a partial recrystallization of the lamellae. Our results provide a solid case-study of the early stages of the diagenetic evolution of crossed-foliated shell layers. Moreover, they highlight the fact that extreme caution must be taken before using fossil shells for palaeoenvironmental or geochronological reconstructions. Without thorough investigation, the alteration patterns illustrated here would easily have gone unnoticed. However, these degradations are liable to bias any proxy based on the elemental, isotopic or biochemical composition of the shells. This study also provides significant data concerning human subsistence behavior: the presence of notches and the good preservation state of limpet shells (no dissolution/recrystallization, no bioerosion and no abrasion/fragmentation aspects) would attest that limpets were gathered alive with tools by Middle Palaeolithic (Aterian) populations in North Africa for consumption

Codon Size Reduction as the Origin of the Triplet Genetic Code

The genetic code appears to be optimized in its robustness to missense errors and frameshift errors. In addition, the genetic code is near-optimal in terms of its ability to carry information in addition to the sequences of encoded proteins. As evolution has no foresight, optimality of the modern genetic code suggests that it evolved from less optimal code variants. The length of codons in the genetic code is also optimal, as three is the minimal nucleotide combination that can encode the twenty standard amino acids. The apparent impossibility of transitions between codon sizes in a discontinuous manner during evolution has resulted in an unbending view that the genetic code was always triplet. Yet, recent experimental evidence on quadruplet decoding, as well as the discovery of organisms with ambiguous and dual decoding, suggest that the possibility of the evolution of triplet decoding from living systems with non-triplet decoding merits reconsideration and further exploration. To explore this possibility we designed a mathematical model of the evolution of primitive digital coding systems which can decode nucleotide sequences into protein sequences. These coding systems can evolve their nucleotide sequences via genetic events of Darwinian evolution, such as point-mutations. The replication rates of such coding systems depend on the accuracy of the generated protein sequences. Computer simulations based on our model show that decoding systems with codons of length greater than three spontaneously evolve into predominantly triplet decoding systems. Our findings suggest a plausible scenario for the evolution of the triplet genetic code in a continuous manner. This scenario suggests an explanation of how protein synthesis could be accomplished by means of long RNA-RNA interactions prior to the emergence of the complex decoding machinery, such as the ribosome, that is required for stabilization and discrimination of otherwise weak triplet codon-anticodon interactions

Investigating the Structural Impacts of I64T and P311S Mutations in APE1-DNA Complex: A Molecular Dynamics Approach

Elucidating the molecular dynamic behavior of Protein-DNA complex upon mutation is crucial in current genomics. Molecular dynamics approach reveals the changes on incorporation of variants that dictate the structure and function of Protein-DNA complexes. Deleterious mutations in APE1 protein modify the physicochemical property of amino acids that affect the protein stability and dynamic behavior. Further, these mutations disrupt the binding sites and prohibit the protein to form complexes with its interacting DNA.In this study, we developed a rapid and cost-effective method to analyze variants in APE1 gene that are associated with disease susceptibility and evaluated their impacts on APE1-DNA complex dynamic behavior. Initially, two different in silico approaches were used to identify deleterious variants in APE1 gene. Deleterious scores that overlap in these approaches were taken in concern and based on it, two nsSNPs with IDs rs61730854 (I64T) and rs1803120 (P311S) were taken further for structural analysis.Different parameters such as RMSD, RMSF, salt bridge, H-bonds and SASA applied in Molecular dynamic study reveals that predicted deleterious variants I64T and P311S alters the structure as well as affect the stability of APE1-DNA interacting functions. This study addresses such new methods for validating functional polymorphisms of human APE1 which is critically involved in causing deficit in repair capacity, which in turn leads to genetic instability and carcinogenesis