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

OneG: A Computational Tool for Predicting Cryptic Intermediates in the Unfolding Kinetics of Proteins under Native Conditions

Understanding the relationships between conformations of proteins and their stabilities is one key to address the protein folding paradigm. The free energy change (ΔG) of unfolding reactions of proteins is measured by traditional denaturation methods and native hydrogen-deuterium (H/D) exchange methods. However, the free energy of unfolding (ΔGU) and the free energy of exchange (ΔGHX) of proteins are not in good agreement, though the experimental conditions of both methods are well matching to each other. The anomaly is due to any one or combinations of the following reasons: (i) effects of cis-trans proline isomerisation under equilibrium unfolding reactions of proteins (ii) inappropriateness in accounting the baselines of melting curves (iii) presence of cryptic intermediates, which may elude the melting curve analysis and (iv) existence of higher energy metastable states in the H/D exchange reactions of proteins. Herein, we have developed a novel computational tool, OneG, which accounts the discrepancy between ΔGU and ΔGHX of proteins by systematically accounting all the four factors mentioned above. The program is fully automated and requires four inputs: three-dimensional structures of proteins, ΔGU, ΔGU* and residue-specific ΔGHX determined under EX2-exchange conditions in the absence of denaturants. The robustness of the program has been validated using experimental data available for proteins such as cytochrome c and apocytochrome b562 and the data analyses revealed that cryptic intermediates of the proteins detected by the experimental methods and the cryptic intermediates predicted by the OneG for those proteins were in good agreement. Furthermore, using OneG, we have shown possible existence of cryptic intermediates and metastable states in the unfolding pathways of cardiotoxin III and cobrotoxin, respectively, which are homologous proteins. The unique application of the program to map the unfolding pathways of proteins under native conditions have been brought into fore and the program is publicly available at http://sblab.sastra.edu/oneg.htm

Structure and evolutionary origin of Ca2+-dependent herring type II antifreeze protein

10.1371/journal.pone.0000548PLoS ONE26

Anopheles gambiae PGRPLC-Mediated Defense against Bacteria Modulates Infections with Malaria Parasites

Recognition of peptidoglycan (PGN) is paramount for insect antibacterial defenses. In the fruit fly Drosophila melanogaster, the transmembrane PGN Recognition Protein LC (PGRP-LC) is a receptor of the Imd signaling pathway that is activated after infection with bacteria, mainly Gram-negative (Gram−). Here we demonstrate that bacterial infections of the malaria mosquito Anopheles gambiae are sensed by the orthologous PGRPLC protein which then activates a signaling pathway that involves the Rel/NF-κB transcription factor REL2. PGRPLC signaling leads to transcriptional induction of antimicrobial peptides at early stages of hemolymph infections with the Gram-positive (Gram+) bacterium Staphylococcus aureus, but a different signaling pathway might be used in infections with the Gram− bacterium Escherichia coli. The size of mosquito symbiotic bacteria populations and their dramatic proliferation after a bloodmeal, as well as intestinal bacterial infections, are also controlled by PGRPLC signaling. We show that this defense response modulates mosquito infection intensities with malaria parasites, both the rodent model parasite, Plasmodium berghei, and field isolates of the human parasite, Plasmodium falciparum. We propose that the tripartite interaction between mosquito microbial communities, PGRPLC-mediated antibacterial defense and infections with Plasmodium can be exploited in future interventions aiming to control malaria transmission. Molecular analysis and structural modeling provided mechanistic insights for the function of PGRPLC. Alternative splicing of PGRPLC transcripts produces three main isoforms, of which PGRPLC3 appears to have a key role in the resistance to bacteria and modulation of Plasmodium infections. Structural modeling indicates that PGRPLC3 is capable of binding monomeric PGN muropeptides but unable to initiate dimerization with other isoforms. A dual role of this isoform is hypothesized: it sequesters monomeric PGN dampening weak signals and locks other PGRPLC isoforms in binary immunostimulatory complexes further enhancing strong signals

VASP: A Volumetric Analysis of Surface Properties Yields Insights into Protein-Ligand Binding Specificity

Many algorithms that compare protein structures can reveal similarities that suggest related biological functions, even at great evolutionary distances. Proteins with related function often exhibit differences in binding specificity, but few algorithms identify structural variations that effect specificity. To address this problem, we describe the Volumetric Analysis of Surface Properties (VASP), a novel volumetric analysis tool for the comparison of binding sites in aligned protein structures. VASP uses solid volumes to represent protein shape and the shape of surface cavities, clefts and tunnels that are defined with other methods. Our approach, inspired by techniques from constructive solid geometry, enables the isolation of volumetrically conserved and variable regions within three dimensionally superposed volumes. We applied VASP to compute a comparative volumetric analysis of the ligand binding sites formed by members of the steroidogenic acute regulatory protein (StAR)-related lipid transfer (START) domains and the serine proteases. Within both families, VASP isolated individual amino acids that create structural differences between ligand binding cavities that are known to influence differences in binding specificity. Also, VASP isolated cavity subregions that differ between ligand binding cavities which are essential for differences in binding specificity. As such, VASP should prove a valuable tool in the study of protein-ligand binding specificity

Compensatory Interactions between Sir3p and the Nucleosomal LRS Surface Imply Their Direct Interaction

The previously identified LRS (Loss of rDNA Silencing) domain of the nucleosome is critically important for silencing at both ribosomal DNA and telomeres. To understand the function of the LRS surface in silencing, we performed an EMS mutagenesis screen to identify suppressors of the H3 A75V LRS allele. We identified dominant and recessive mutations in histones H3, H4, and dominant mutations in the BAH (Bromo Adjacent Homology) domain of SIR3. We further characterized a surface of Sir3p critical for silencing via the LRS surface. We found that all alleles of the SIR3 BAH domain were able to 1) generally suppress the loss of telomeric silencing of LRS alleles, but 2) could not suppress SIN (Swi/Snf Independent) alleles or 3) could not suppress the telomeric silencing defect of H4 tail alleles. Moreover, we noticed a complementary trend in the electrostatic changes resulting from most of the histone mutations that gain or lose silencing and the suppressor alleles isolated in SIR3, and the genes for histones H3 and H4. Mutations in H3 and H4 genes that lose silencing tend to make the LRS surface more electronegative, whereas mutations that increase silencing make it less electronegative. Conversely, suppressors of LRS alleles in either SIR3, histone H3, or H4 also tend to make their respective surfaces less electronegative. Our results provide genetic evidence for recent data suggesting that the Sir3p BAH domain directly binds the LRS domain. Based on these findings, we propose an electrostatic model for how an extensive surface on the Sir3p BAH domain may regulate docking onto the LRS surface

Insight into the Mechanisms of Adenovirus Capsid Disassembly from Studies of Defensin Neutralization

Defensins are effectors of the innate immune response with potent antibacterial activity. Their role in antiviral immunity, particularly for non-enveloped viruses, is poorly understood. We recently found that human alpha-defensins inhibit human adenovirus (HAdV) by preventing virus uncoating and release of the endosomalytic protein VI during cell entry. Consequently, AdV remains trapped in the endosomal/lysosomal pathway rather than trafficking to the nucleus. To gain insight into the mechanism of defensin-mediated neutralization, we analyzed the specificity of the AdV-defensin interaction. Sensitivity to alpha-defensin neutralization is a common feature of HAdV species A, B1, B2, C, and E, whereas species D and F are resistant. Thousands of defensin molecules bind with low micromolar affinity to a sensitive serotype, but only a low level of binding is observed to resistant serotypes. Neutralization is dependent upon a correctly folded defensin molecule, suggesting that specific molecular interactions occur with the virion. CryoEM structural studies and protein sequence analysis led to a hypothesis that neutralization determinants are located in a region spanning the fiber and penton base proteins. This model was supported by infectivity studies using virus chimeras comprised of capsid proteins from sensitive and resistant serotypes. These findings suggest a mechanism in which defensin binding to critical sites on the AdV capsid prevents vertex removal and thereby blocks subsequent steps in uncoating that are required for release of protein VI and endosomalysis during infection. In addition to informing the mechanism of defensin-mediated neutralization of a non-enveloped virus, these studies provide insight into the mechanism of AdV uncoating and suggest new strategies to disrupt this process and inhibit infection

Selective and Irreversible Inhibitors of Mosquito Acetylcholinesterases for Controlling Malaria and Other Mosquito-Borne Diseases

New insecticides are urgently needed because resistance to current insecticides allows resurgence of disease-transmitting mosquitoes while concerns for human toxicity from current compounds are growing. We previously reported the finding of a free cysteine (Cys) residue at the entrance of the active site of acetylcholinesterase (AChE) in some insects but not in mammals, birds, and fish. These insects have two AChE genes (AP and AO), and only AP-AChE carries the Cys residue. Most of these insects are disease vectors such as the African malaria mosquito (Anopheles gambiae sensu stricto) or crop pests such as aphids. Recently we reported a Cys-targeting small molecule that irreversibly inhibited all AChE activity extracted from aphids while an identical exposure caused no effect on the human AChE. Full inhibition of AChE in aphids indicates that AP-AChE contributes most of the enzymatic activity and suggests that the Cys residue might serve as a target for developing better aphicides. It is therefore worth investigating whether the Cys-targeting strategy is applicable to mosquitocides. Herein, we report that, under conditions that spare the human AChE, a methanethiosulfonate-containing molecule at 6 µM irreversibly inhibited 95% of the AChE activity extracted from An. gambiae s. str. and >80% of the activity from the yellow fever mosquito (Aedes aegypti L.) or the northern house mosquito (Culex pipiens L.) that is a vector of St. Louis encephalitis. This type of inhibition is fast (∼30 min) and due to conjugation of the inhibitor to the active-site Cys of mosquito AP-AChE, according to our observed reactivation of the methanethiosulfonate-inhibited AChE by 2-mercaptoethanol. We also note that our sulfhydryl agents partially and irreversibly inhibited the human AChE after prolonged exposure (>4 hr). This slow inhibition is due to partial enzyme denaturation by the inhibitor and/or micelles of the inhibitor, according to our studies using atomic force microscopy, circular dichroism spectroscopy, X-ray crystallography, time-resolved fluorescence spectroscopy, and liquid chromatography triple quadrupole mass spectrometry. These results support our view that the mosquito-specific Cys is a viable target for developing new mosquitocides to control disease vectors and to alleviate resistance problems with reduced toxicity toward non-target species

Dynamic Evolution of Pathogenicity Revealed by Sequencing and Comparative Genomics of 19 Pseudomonas syringae Isolates

Closely related pathogens may differ dramatically in host range, but the molecular, genetic, and evolutionary basis for these differences remains unclear. In many Gram- negative bacteria, including the phytopathogen Pseudomonas syringae, type III effectors (TTEs) are essential for pathogenicity, instrumental in structuring host range, and exhibit wide diversity between strains. To capture the dynamic nature of virulence gene repertoires across P. syringae, we screened 11 diverse strains for novel TTE families and coupled this nearly saturating screen with the sequencing and assembly of 14 phylogenetically diverse isolates from a broad collection of diseased host plants. TTE repertoires vary dramatically in size and content across all P. syringae clades; surprisingly few TTEs are conserved and present in all strains. Those that are likely provide basal requirements for pathogenicity. We demonstrate that functional divergence within one conserved locus, hopM1, leads to dramatic differences in pathogenicity, and we demonstrate that phylogenetics-informed mutagenesis can be used to identify functionally critical residues of TTEs. The dynamism of the TTE repertoire is mirrored by diversity in pathways affecting the synthesis of secreted phytotoxins, highlighting the likely role of both types of virulence factors in determination of host range. We used these 14 draft genome sequences, plus five additional genome sequences previously reported, to identify the core genome for P. syringae and we compared this core to that of two closely related non-pathogenic pseudomonad species. These data revealed the recent acquisition of a 1 Mb megaplasmid by a sub-clade of cucumber pathogens. This megaplasmid encodes a type IV secretion system and a diverse set of unknown proteins, which dramatically increases both the genomic content of these strains and the pan-genome of the species