New genes of Xanthomonas citri subsp. citri involved in pathogenesis and adaptation revealed by a transposon-based mutant library

Background  

Citrus canker is a disease caused by the phytopathogens Xanthomonas citri subsp. citri, Xanthomonas fuscans subsp. aurantifolli and Xanthomonas alfalfae subsp. citrumelonis. The first of the three species, which causes citrus bacterial canker type A, is the most widely spread and severe, attacking all citrus species. In Brazil, this species is the most important, being found in practically all areas where citrus canker has been detected. Like most phytobacterioses, there is no efficient way to control citrus canker. Considering the importance of the disease worldwide, investigation is needed to accurately detect which genes are related to the pathogen-host adaptation process and which are associated with pathogenesis.     

Fast protein folding on downhill energy landscape

Proteins fold in a time range of microseconds to minutes despite the large amount of possible conformers. Molecular dynamics simulations of a three-stranded antiparallel beta-sheet peptide (for a total of 12.6 microsec and 72 folding events) show that at the melting temperature the unfolded state ensemble contains many more conformers than those sampled during a folding event.     

The role of aromaticity, exposed surface, and dipole moment in determining protein aggregation rates

The mechanisms by which peptides and proteins form ordered aggregates are not well understood. Here we focus on the physicochemical properties of amino acids that favor ordered aggregation and suggest a parameter-free model that is able to predict the change of aggregation rates over a large set of natural sequences. Furthermore, the results of the parameter-free model correlate well with the aggregation propensities of a set of peptides designed by computer simulations.     

Computer simulations of protein folding by targeted molecular dynamics

We have performed 128 folding and 45 unfolding molecular dynamics runs of chymotrypsin inhibitor 2 (CI2) with an implicit solvation model for a total simulation time of 0.4 microseconds. Folding requires that the three-dimensional structure of the native state is known. It was simulated at 300 K by supplementing the force field with a harmonic restraint which acts on the root-mean-square deviation and allows to decrease the distance to the target conformation. High temperature and/or the harmonic restraint were used to induce unfolding. Of the 62 folding simulations started from random conformations, 31 reached the native structure, while the success rate was 83% for the 66 trajectories which began from conformations unfolded by high-temperature dynamics. A funnel-like energy landscape is observed for unfolding at 475 K, while the unfolding runs at 300 K and 375 K as well as most of the folding trajectories have an almost flat energy landscape for conformations with less than about 50% of native contacts formed. The sequence of events, i.e., secondary and tertiary structure formation, is similar in all folding and unfolding simulations, despite the diversity of the pathways. Previous unfolding simulations of CI2 performed with different force fields showed a similar sequence of events. These results suggest that the topology of the native state plays an important role in the folding process.     

Hydrophobicity at the surface of proteins

A new method is presented to quantitatively estimate and graphically display the propensity of nonpolar groups to bind at the surface of proteins. It is based on the calculation of the binding energy, i.e., van der Waals interaction plus protein electrostatic desolvation, of a nonpolar probe sphere rolled over the protein surface, and on the color coding of this quantity on a smooth molecular surface (hydrophobicity map). The method is validated on ten protein-ligand complexes and is shown to distinguish precisely where polar and nonpolar groups preferentially bind. Comparisons with existing approaches, like the display of the electrostatic potential or the curvature, illustrate the advantages and the better predictive power of the present method. Hydrophobicity maps will play an important role in the characterization of binding sites for the large number of proteins emerging from the genome projects and structure modeling approaches.     

Computational analysis of the S. cerevisiae proteome reveals the function and cellular localization of the least and most amyloidogenic proteins

Protein sequences have evolved to optimize biological function that usually requires a well-defined three-dimensional structure and a monomeric (or oligomeric) state. These two requirements may be in conflict as the propensity for beta-sheet structure, which is one of the two most common regular conformations of the polypeptide chain in folded proteins, favors also the formation of ordered aggregates of multiple copies of the same protein (fibril, i.e., polymeric state). Such beta-aggregation is typical of amyloid diseases that include Alzheimer's, Parkinson's, and type II diabetes as well as the spongiform encephalopathies. Here, an analytical model previously developed for evaluating the amyloidogenic potential of polypeptides is applied to the proteome of the budding yeast (Saccharomyces cerevisiae). The model is based on the physicochemical properties that are relevant for beta-aggregation and requires only the protein sequence as input. It is shown that beta-aggregation prone proteins in yeast are accrued in molecular transport, protein biosynthesis, and cell wall organization processes while they are underrepresented in ribosome biogenesis, RNA metabolism, and vitamin metabolism. Furthermore, beta-aggregation prone proteins are much more abundant in the cell wall, endoplasmic reticulum, and plasma membrane than in the nucleolus, ribosome, and nucleus. Thus, this study indicates that evolution has not only prevented the selection of amyloidogenic sequences in cellular compartments characterized by a high concentration of unfolded proteins but also tried to exploit the beta-aggregated state for certain functions (e.g. molecular transport) and in well-confined cellular environments or organelles to protect the rest of the cell from toxic (pre-)fibrillar species.     

Disulfide bridge formation influences ligand recognition by the ATAD2 bromodomain

The ATPase family, AAA domain-containing protein 2 (ATAD2) has a C-terminal bromodomain, which functions as a chromatin reader domain recognizing acetylated lysine on the histone tails within the nucleosome. ATAD2 is overexpressed in many cancers and its expression is correlated with poor patient outcomes, making it an attractive therapeutic target and potential biomarker. We solved the crystal structure of the ATAD2 bromodomain and found that it contains a disulfide bridge near the base of the acetyllysine binding pocket (Cys1057-Cys1079). Site-directed mutagenesis revealed that removal of a free C-terminal cysteine (C1101) residue greatly improved the solubility of the ATAD2 bromodomain in vitro. Isothermal titration calorimetry experiments in combination with the Ellman's assay demonstrated that formation of an intramolecular disulfide bridge negatively impacts the ligand binding affinities and alters the thermodynamic parameters of the ATAD2 bromodomain interaction with a histone H4K5ac peptide as well as a small molecule bromodomain ligand. Molecular dynamics simulations indicate that the formation of the disulfide bridge in the ATAD2 bromodomain does not alter the structure of the folded state or flexibility of the acetyllysine binding pocket. However, consideration of this unique structural feature should be taken into account when examining ligand-binding affinity, or in the design of new bromodomain inhibitor compounds that interact with this acetyllysine reader module.