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

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

Early and late outcomes after minimally invasive direct coronary artery bypass vs. full sternotomy off-pump coronary artery bypass grafting

ObjectivesMinimally-invasive direct coronary artery bypass (MIDCAB) is a less-invasive alternative to full sternotomy off-pump coronary artery bypass (FS-OPCAB) revascularization of the left anterior descending artery (LAD). Some studies suggested that MIDCAB is associated with a greater risk of graft occlusion and repeat revascularization than FS-OPCAB LIMA-to-LAD grafting. Data comparing MIDCAB to FS-OPCAB with regard to long-term follow-up is scarce. We compared short- and long-term results of MIDCAB vs. FS-OPCAB revascularization over a maximum follow-up period of 10 years.Patients and methodsFrom December 2009 to June 2020, 388 elective patients were included in our retrospective study. 229 underwent MIDCAB, and 159 underwent FS-OPCAB LIMA-to-LAD grafting. Inverse probability of treatment weighting (IPTW) was used to adjust for selection bias and to estimate treatment effects on short- and long-term outcomes. IPTW-adjusted Kaplan–Meier estimates by study group were calculated for all-cause mortality, stroke, the risk of repeat revascularization and myocardial infarction up to a maximum follow-up of 10 years.ResultsMIDCAB patients had less rethoracotomies (n = 13/3.6% vs. n = 30/8.0%, p = 0.012), fewer transfusions (0.93 units ± 1.83 vs. 1.61 units ± 2.52, p < 0.001), shorter mechanical ventilation time (7.6 ± 4.7 h vs. 12.1 ± 26.4 h, p = 0.005), and needed less hemofiltration (n = 0/0% vs. n = 8/2.4%, p = 0.004). Thirty-day mortality did not differ significantly between the two groups (n = 0/0% vs. n = 3/0.8%, p = 0.25). Long-term outcomes did not differ significantly between study groups. In the FS-OPCAB group, the probability of survival at 1, 5, and 10 years was 98.4%, 87.8%, and 71.7%, respectively. In the MIDCAB group, the corresponding values were 98.4%, 87.7%, and 68.7%, respectively (RR1.24, CI0.87–1.86, p = 0.7). In the FS group, the freedom from stroke at 1, 5, and 10 years was 97.0%, 93.0%, and 93.0%, respectively. In the MIDCAB group, the corresponding values were 98.5%, 96.9%, and 94.3%, respectively (RR0.52, CI0.25–1.09, p = 0.06). Freedom from repeat revascularization at 1, 5, and 10 years in the FS-OPCAB group was 92.2%, 84.7%, and 79.5%, respectively. In the MIDCAB group, the corresponding values were 94.8%, 90.2%, and 81.7%, respectively (RR0.73, CI0.47–1.16, p = 0.22).ConclusionMIDCAB is a safe and efficacious technique and offers comparable long-term results regarding mortality, stroke, repeat revascularization, and freedom from myocardial infarction when compared to FS-OPCAB

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

Familial hypercholesterolaemia in children and adolescents from 48 countries: a cross-sectional study

Background Approximately 450 000 children are born with familial hypercholesterolaemia worldwide every year, yet only 2·1% of adults with familial hypercholesterolaemia were diagnosed before age 18 years via current diagnostic approaches, which are derived from observations in adults. We aimed to characterise children and adolescents with heterozygous familial hypercholesterolaemia (HeFH) and understand current approaches to the identification and management of familial hypercholesterolaemia to inform future public health strategies. Methods For this cross-sectional study, we assessed children and adolescents younger than 18 years with a clinical or genetic diagnosis of HeFH at the time of entry into the Familial Hypercholesterolaemia Studies Collaboration (FHSC) registry between Oct 1, 2015, and Jan 31, 2021. Data in the registry were collected from 55 regional or national registries in 48 countries. Diagnoses relying on self-reported history of familial hypercholesterolaemia and suspected secondary hypercholesterolaemia were excluded from the registry; people with untreated LDL cholesterol (LDL-C) of at least 13·0 mmol/L were excluded from this study. Data were assessed overall and by WHO region, World Bank country income status, age, diagnostic criteria, and index-case status. The main outcome of this study was to assess current identification and management of children and adolescents with familial hypercholesterolaemia. Findings Of 63 093 individuals in the FHSC registry, 11 848 (18·8%) were children or adolescents younger than 18 years with HeFH and were included in this study; 5756 (50·2%) of 11 476 included individuals were female and 5720 (49·8%) were male. Sex data were missing for 372 (3·1%) of 11 848 individuals. Median age at registry entry was 9·6 years (IQR 5·8–13·2). 10 099 (89·9%) of 11 235 included individuals had a final genetically confirmed diagnosis of familial hypercholesterolaemia and 1136 (10·1%) had a clinical diagnosis. Genetically confirmed diagnosis data or clinical diagnosis data were missing for 613 (5·2%) of 11 848 individuals. Genetic diagnosis was more common in children and adolescents from high-income countries (9427 [92·4%] of 10 202) than in children and adolescents from non-high-income countries (199 [48·0%] of 415). 3414 (31·6%) of 10 804 children or adolescents were index cases. Familial-hypercholesterolaemia-related physical signs, cardiovascular risk factors, and cardiovascular disease were uncommon, but were more common in non-high-income countries. 7557 (72·4%) of 10 428 included children or adolescents were not taking lipid-lowering medication (LLM) and had a median LDL-C of 5·00 mmol/L (IQR 4·05–6·08). Compared with genetic diagnosis, the use of unadapted clinical criteria intended for use in adults and reliant on more extreme phenotypes could result in 50–75% of children and adolescents with familial hypercholesterolaemia not being identified. Interpretation Clinical characteristics observed in adults with familial hypercholesterolaemia are uncommon in children and adolescents with familial hypercholesterolaemia, hence detection in this age group relies on measurement of LDL-C and genetic confirmation. Where genetic testing is unavailable, increased availability and use of LDL-C measurements in the first few years of life could help reduce the current gap between prevalence and detection, enabling increased use of combination LLM to reach recommended LDL-C targets early in life. Funding Pfizer, Amgen, Merck Sharp & Dohme, Sanofi–Aventis, Daiichi Sankyo, and Regeneron

Trends in template/fragment-free protein structure prediction

Predicting the structure of a protein from its amino acid sequence is a long-standing unsolved problem in computational biology. Its solution would be of both fundamental and practical importance as the gap between the number of known sequences and the number of experimentally solved structures widens rapidly. Currently, the most successful approaches are based on fragment/template reassembly. Lacking progress in template-free structure prediction calls for novel ideas and approaches. This article reviews trends in the development of physical and specific knowledge-based energy functions as well as sampling techniques for fragment-free structure prediction. Recent physical- and knowledge-based studies demonstrated that it is possible to sample and predict highly accurate protein structures without borrowing native fragments from known protein structures. These emerging approaches with fully flexible sampling have the potential to move the field forward

Familial hypercholesterolaemia in children and adolescents from 48 countries: a cross-sectional study

Background: Approximately 450 000 children are born with familial hypercholesterolaemia worldwide every year, yet only 2·1% of adults with familial hypercholesterolaemia were diagnosed before age 18 years via current diagnostic approaches, which are derived from observations in adults. We aimed to characterise children and adolescents with heterozygous familial hypercholesterolaemia (HeFH) and understand current approaches to the identification and management of familial hypercholesterolaemia to inform future public health strategies. Methods: For this cross-sectional study, we assessed children and adolescents younger than 18 years with a clinical or genetic diagnosis of HeFH at the time of entry into the Familial Hypercholesterolaemia Studies Collaboration (FHSC) registry between Oct 1, 2015, and Jan 31, 2021. Data in the registry were collected from 55 regional or national registries in 48 countries. Diagnoses relying on self-reported history of familial hypercholesterolaemia and suspected secondary hypercholesterolaemia were excluded from the registry; people with untreated LDL cholesterol (LDL-C) of at least 13·0 mmol/L were excluded from this study. Data were assessed overall and by WHO region, World Bank country income status, age, diagnostic criteria, and index-case status. The main outcome of this study was to assess current identification and management of children and adolescents with familial hypercholesterolaemia. Findings: Of 63 093 individuals in the FHSC registry, 11 848 (18·8%) were children or adolescents younger than 18 years with HeFH and were included in this study; 5756 (50·2%) of 11 476 included individuals were female and 5720 (49·8%) were male. Sex data were missing for 372 (3·1%) of 11 848 individuals. Median age at registry entry was 9·6 years (IQR 5·8-13·2). 10 099 (89·9%) of 11 235 included individuals had a final genetically confirmed diagnosis of familial hypercholesterolaemia and 1136 (10·1%) had a clinical diagnosis. Genetically confirmed diagnosis data or clinical diagnosis data were missing for 613 (5·2%) of 11 848 individuals. Genetic diagnosis was more common in children and adolescents from high-income countries (9427 [92·4%] of 10 202) than in children and adolescents from non-high-income countries (199 [48·0%] of 415). 3414 (31·6%) of 10 804 children or adolescents were index cases. Familial-hypercholesterolaemia-related physical signs, cardiovascular risk factors, and cardiovascular disease were uncommon, but were more common in non-high-income countries. 7557 (72·4%) of 10 428 included children or adolescents were not taking lipid-lowering medication (LLM) and had a median LDL-C of 5·00 mmol/L (IQR 4·05-6·08). Compared with genetic diagnosis, the use of unadapted clinical criteria intended for use in adults and reliant on more extreme phenotypes could result in 50-75% of children and adolescents with familial hypercholesterolaemia not being identified. Interpretation: Clinical characteristics observed in adults with familial hypercholesterolaemia are uncommon in children and adolescents with familial hypercholesterolaemia, hence detection in this age group relies on measurement of LDL-C and genetic confirmation. Where genetic testing is unavailable, increased availability and use of LDL-C measurements in the first few years of life could help reduce the current gap between prevalence and detection, enabling increased use of combination LLM to reach recommended LDL-C targets early in life