Permalink
Cannot retrieve contributors at this time
Fetching contributors…
Cannot retrieve contributors at this time
| <?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1d1 20130915//EN" "JATS-archivearticle1.dtd"><article article-type="research-article" dtd-version="1.1d1" xmlns:xlink="http://www.w3.org/1999/xlink"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="hwp">eLife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">03532</article-id><article-id pub-id-type="doi">10.7554/eLife.03532</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research article</subject></subj-group><subj-group subj-group-type="heading"><subject>Biophysics and structural biology</subject></subj-group></article-categories><title-group><article-title>Allosteric signalling in the outer membrane translocation domain of PapC usher</article-title></title-group><contrib-group><contrib contrib-type="author" id="author-14824"><name><surname>Farabella</surname><given-names>Irene</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="other" rid="par-1"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-14825"><name><surname>Pham</surname><given-names>Thieng</given-names></name><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-14826"><name><surname>Henderson</surname><given-names>Nadine S</given-names></name><xref ref-type="aff" rid="aff3"/><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-14827"><name><surname>Geibel</surname><given-names>Sebastian</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="pa1">†</xref><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-14828"><name><surname>Phan</surname><given-names>Gilles</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="pa2">‡</xref><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-14829"><name><surname>Thanassi</surname><given-names>David G</given-names></name><xref ref-type="aff" rid="aff3"/><xref ref-type="other" rid="par-4"/><xref ref-type="fn" rid="con6"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-14830"><name><surname>Delcour</surname><given-names>Anne H</given-names></name><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="con7"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-11280"><name><surname>Waksman</surname><given-names>Gabriel</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="con8"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-15325"><name><surname>Topf</surname><given-names>Maya</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="corresp" rid="cor1">*</xref><xref ref-type="other" rid="par-2"/><xref ref-type="other" rid="par-3"/><xref ref-type="fn" rid="con9"/><xref ref-type="fn" rid="conf1"/></contrib><aff id="aff1"><institution content-type="dept">Institute of Structural and Molecular Biology</institution>, <institution>Birkbeck College and University College London</institution>, <addr-line><named-content content-type="city">London</named-content></addr-line>, <country>United Kingdom</country></aff><aff id="aff2"><institution content-type="dept">Department of Biology and Biochemistry</institution>, <institution>University of Houston</institution>, <addr-line><named-content content-type="city">Houston</named-content></addr-line>, <country>United States</country></aff><aff id="aff3"><institution content-type="dept">Department of Molecular Genetics and Microbiology, Center for Infectious Diseases</institution>, <institution>Stony Brook University</institution>, <addr-line><named-content content-type="city">Stony Brook</named-content></addr-line>, <country>United States</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Dötsch</surname><given-names>Volker</given-names></name><role>Reviewing editor</role><aff><institution>Goethe University</institution>, <country>Germany</country></aff></contrib></contrib-group><author-notes><corresp id="cor1"><label>*</label>For correspondence: <email>m.topf@cryst.bbk.ac.uk</email></corresp><fn fn-type="present-address" id="pa1"><label>†</label><p>Institut für Molekulare Infektionsbiologie, University of Würzburg, Würzburg, Germany</p></fn><fn fn-type="present-address" id="pa2"><label>‡</label><p>Faculté des Sciences, Université Paris Descartes, Paris, France</p></fn></author-notes><pub-date date-type="pub" publication-format="electronic"><day>01</day><month>10</month><year>2014</year></pub-date><pub-date pub-type="collection"><year>2014</year></pub-date><volume>3</volume><elocation-id>e03532</elocation-id><history><date date-type="received"><day>30</day><month>05</month><year>2014</year></date><date date-type="accepted"><day>29</day><month>09</month><year>2014</year></date></history><permissions><copyright-statement>© 2014, Farabella et al</copyright-statement><copyright-year>2014</copyright-year><copyright-holder>Farabella et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife03532.pdf"/><abstract><object-id pub-id-type="doi">10.7554/eLife.03532.001</object-id><p>PapC ushers are outer-membrane proteins enabling assembly and secretion of P pili in uropathogenic <italic>E. coli</italic>. Their translocation domain is a large β-barrel occluded by a plug domain, which is displaced to allow the translocation of pilus subunits across the membrane. Previous studies suggested that this gating mechanism is controlled by a β-hairpin and an α-helix. To investigate the role of these elements in allosteric signal communication, we developed a method combining evolutionary and molecular dynamics studies of the native translocation domain and mutants lacking the β-hairpin and/or the α-helix. Analysis of a hybrid residue interaction network suggests distinct regions (residue ‘communities’) within the translocation domain (especially around β12–β14) linking these elements, thereby modulating PapC gating. Antibiotic sensitivity and electrophysiology experiments on a set of alanine-substitution mutants confirmed functional roles for four of these communities. This study illuminates the gating mechanism of PapC ushers and its importance in maintaining outer-membrane permeability.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.001">http://dx.doi.org/10.7554/eLife.03532.001</ext-link></p></abstract><abstract abstract-type="executive-summary"><object-id pub-id-type="doi">10.7554/eLife.03532.002</object-id><title>eLife digest</title><p><italic>Escherichia coli</italic> is a bacterium that commonly lives in the intestines of mammals, including humans, where it is usually harmless and can even be beneficial to its host. However, some types of <italic>E. coli</italic> produce hair-like filaments called P pili that allow the bacteria to attach to the human urinary tract and cause disease. To pass through the outer membrane of the <italic>E. coli</italic> cell, the filaments have to travel through a protein in the membrane called PapC usher.</p><p>The PapC usher protein—which is also involved in the assembly of the P pili filaments—contains a tube-like part called a β-barrel that is usually blocked by another part of the protein called the ‘plug domain’. For the P pili to pass through the β-barrel, the plug domain has to move. This movement is controlled by two parts of the PapC protein, known as the α-helix and the β-hairpin, but it is not clear how.</p><p>To address this question, Farabella et al. made computer models of the normal PapC protein and versions that lacked the α-helix and/or the β-hairpin. Looking at these structural models and analyzing the evolution of PapC proteins helped to predict that certain regions of the β-barrel may be involved in controlling the movement of the plug domain, and this was then confirmed experimentally. Farabella et al. propose that these regions—together with the α-helix and β-hairpin—control the opening and closing of the β-barrel.</p><p>Further work is needed to investigate how other parts of the PapC protein are involved in P pili formation. These new insights could prove useful in the development of alternative treatments to fight bacterial infection.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.002">http://dx.doi.org/10.7554/eLife.03532.002</ext-link></p></abstract><kwd-group kwd-group-type="author-keywords"><title>Author keywords</title><kwd>dynamics</kwd><kwd>outer membrane protein</kwd><kwd>evolution</kwd><kwd>structure</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd><italic>E. coli</italic></kwd></kwd-group><funding-group><award-group id="par-1"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100004440</institution-id><institution>Wellcome Trust</institution></institution-wrap></funding-source><award-id>PhD training programme 086711/z/08/z</award-id><principal-award-recipient><name><surname>Farabella</surname><given-names>Irene</given-names></name></principal-award-recipient></award-group><award-group id="par-2"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/501100000265</institution-id><institution>Medical Research Council</institution></institution-wrap></funding-source><award-id>G0600084</award-id><principal-award-recipient><name><surname>Topf</surname><given-names>Maya</given-names></name></principal-award-recipient></award-group><award-group id="par-3"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/501100000268</institution-id><institution>Biotechnology and Biological Sciences Research Council</institution></institution-wrap></funding-source><award-id>BB/K01692X/1</award-id><principal-award-recipient><name><surname>Topf</surname><given-names>Maya</given-names></name></principal-award-recipient></award-group><award-group id="par-4"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000002</institution-id><institution>National Institutes of Health</institution></institution-wrap></funding-source><award-id>GM062987</award-id><principal-award-recipient><name><surname>Thanassi</surname><given-names>David G</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta><meta-name>elife-xml-version</meta-name><meta-value>2</meta-value></custom-meta><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>Input from computational models has enabled the detection of allosteric communication that modulates the gating mechanism of a bacterial outer-membrane protein.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>Gram-negative pathogens commonly express a vast variety of complex surface organelles that are involved in different cellular processes. One of these organelles, known as pili (or fimbriae), forms a class of virulence factors involved in host cell adhesion and recognition, invasion, cell mobility, and biofilm formation. P pili from uropathogenic <italic>Escherichia coli</italic> are specifically required for the colonization of the human kidney epithelium, a critical event in the kidney infection process (pyelonephritis) (<xref ref-type="bibr" rid="bib51">Roberts et al., 1994</xref>). P pili are assembled on the bacterial outer membrane (OM) via the chaperone/usher (CU) pathway (<xref ref-type="bibr" rid="bib59">Thanassi et al., 1998</xref>), which is often used as a model system to elucidate the mechanism of pilus biogenesis (<xref ref-type="bibr" rid="bib64">Waksman and Hultgren, 2009</xref>).</p><p>The biogenesis of pili via the CU pathway is a highly ordered process that comprises sequential steps. The chaperone protein (PapD) brings the pilins to the bacterial OM where they are assembled into a pilus at a transmembrane pore protein known as the usher (PapC). The usher (∼800 residues) is composed of five domains (<xref ref-type="fig" rid="fig1">Figure 1A</xref>): a periplasmic N-terminal domain (NTD), an OM central translocation domain (TD) that comprises a translocation pore domain (TP), interrupted by a conserved Ig-like plug domain (PD), and two domains at the periplasmic C-terminal end (CTD1 and CTD2) (<xref ref-type="bibr" rid="bib60">Thanassi et al., 2002</xref>; <xref ref-type="bibr" rid="bib45">Ng et al., 2004</xref>; <xref ref-type="bibr" rid="bib9">Capitani et al., 2006</xref>; <xref ref-type="bibr" rid="bib49">Phan et al., 2011</xref>; <xref ref-type="bibr" rid="bib18">Geibel et al., 2013</xref>). The structure of the apo TD (<xref ref-type="fig" rid="fig1">Figure 1B,C</xref>) consists of a 24-stranded kidney-shaped β-barrel where the PD is inserted into the loop connecting two β-strands (β6–β7), occluding the luminal volume of the pore (<xref ref-type="bibr" rid="bib50">Remaut et al., 2008</xref>; <xref ref-type="bibr" rid="bib27">Huang et al., 2009</xref>). In the activated form of another archetypal member of the usher family, FimD, the PD is located outside the pore lumen in the periplasm, next to the NTD (<xref ref-type="bibr" rid="bib49">Phan et al., 2011</xref>; <xref ref-type="bibr" rid="bib18">Geibel et al., 2013</xref>). In addition to the PD, there are two secondary structure elements that uniquely characterize the large β-barrel structures of the usher TD (<xref ref-type="fig" rid="fig1">Figure 1B</xref>). The first element is a <italic>β-hairpin</italic> that creates a large gap in the side of the β-barrel, a feature unprecedented in previously known OM β-barrel structures (<xref ref-type="bibr" rid="bib50">Remaut et al., 2008</xref>). This element (located between strands β5 and β6 of the barrel, <xref ref-type="fig" rid="fig1">Figure 1C</xref>) folds into the barrel lumen and constrains the PD laterally inside the barrel pore. Mutants lacking the β-hairpin show an increased pore permeability suggesting that the β-hairpin has a role in maintaining the PD in a closed conformation (<xref ref-type="bibr" rid="bib63">Volkan et al., 2013</xref>). The second element is an <italic>α-helix</italic> (located on the loop between β13 and β14, <xref ref-type="fig" rid="fig1">Figure 1B</xref>), which caps the β-hairpin from the extracellular side. Mutants lacking the α-helix, or in which the interface between the helix and the PD is disrupted, present a remarkable increase in pore permeability, comparable with that of the mutant lacking the PD, suggesting a role for the helix in maintaining the PD in a closed state (<xref ref-type="bibr" rid="bib36">Mapingire et al., 2009</xref>; <xref ref-type="bibr" rid="bib63">Volkan et al., 2013</xref>).<fig-group><fig id="fig1" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.003</object-id><label>Figure 1.</label><caption><title>PapC usher organization and detail of its translocation domain.</title><p>(<bold>A</bold>) A diagram of the domain organization of PapC usher. NTD (dark-blue) represents the N-terminal domain, CTD1 (light-violet) and CTD2 (dark-violet) represent the C-terminal domains; TD represents the translocation domain, comprising the TP (translocation pore, light-blue) and the PD (plug domain, magenta). (<bold>B</bold> and <bold>C</bold>): Ribbon representation of the starting model of the native translocation domain (TD) of PapC with the labels ‘N’ and ‘C’ indicating the N and C termini of the translocation channel. The β-barrel, PD (including the P-linkers), β-hairpin, and α-helix (including the H-linkers) are coloured blue, magenta, orange, and yellow, respectively. The outer membrane position is represented schematically with the labels ‘E’, ’M’, and ‘P’ indicating the extracellular side, the membrane, and the periplasmic side, respectively. Side view of the TD (<bold>B</bold>) is shown with the α-helix, β-hairpin, H-linker1, H-linker2, P-linker1, P-linker2, and PD, labelled. Extracellular top view of the TD (<bold>C</bold>) is shown with the barrel β strands labelled β1 through β24 and with the PD strands labelled βA through βF. The figures were created with Chimera (<xref ref-type="bibr" rid="bib48">Pettersen et al., 2004</xref>).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.003">http://dx.doi.org/10.7554/eLife.03532.003</ext-link></p></caption><graphic xlink:href="elife03532f001"/></fig><fig id="fig1s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.03532.004</object-id><label>Figure 1—figure supplement 1.</label><caption><title>MD simulations of the native PapC TD and its mutants.</title><p>(<bold>A</bold>) Cutaway view across the membrane plane of the native PapC TD starting model in a POPE/POPG lipid bilayer (sim1, t = 0). Molecular surface of PapC TD is coloured as in <xref ref-type="fig" rid="fig1">Figure 1</xref>, the lipids are shown in grey with the lipid head group coloured by element, the water is coloured by element, and the ions (the Na<sup>+</sup> in blue and the CL<sup>−</sup> in yellow) are represented as sphere. The Cα-RMSD values for each system from the starting structure (t = 0) for the native TD (<bold>B</bold>), the hairpin mutant (<bold>C</bold>), helix mutant (<bold>D</bold>), and helix-hairpin mutant (<bold>E</bold>) are plotted as a function of time.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.004">http://dx.doi.org/10.7554/eLife.03532.004</ext-link></p></caption><graphic xlink:href="elife03532fs001"/></fig></fig-group></p><p>The mutant lacking both the β-hairpin and the α-helix is defective for pilus biogenesis (<xref ref-type="bibr" rid="bib36">Mapingire et al., 2009</xref>). It has been observed in other OMP β-barrels that such secondary structure elements (e.g., an α-helix that protrudes inside the barrel or packs against the transmembrane strands) can use complex allosteric mechanisms to mediate their function (<xref ref-type="bibr" rid="bib43">Naveed et al., 2009</xref>). These are often combinations of large conformational changes (‘global motions’) dictated by the overall architecture (including movement of secondary structure elements) and smaller changes (‘local motions’, such as the motion of recognition loops and side-chain fluctuations) (<xref ref-type="bibr" rid="bib34">Liu and Bahar, 2012</xref>). Additionally, it has been shown that important residues in terms of evolution (highly-coevolved or conserved) could have a pivotal role in mediating such allosteric communications (<xref ref-type="bibr" rid="bib56">Suel et al., 2003</xref>; <xref ref-type="bibr" rid="bib57">Tang et al., 2007</xref>).</p><p>In this study, to understand the allosteric mechanism leading to the plug displacement in PapC and the involvement of the α-helix and β-hairpin, we used a hybrid computational approach and verified our results experimentally. By combining sequence conservation analysis, mutual information-based coevolution analysis, and all-atom molecular dynamics (AA-MD), we modelled the interaction network within the native PapC TD as well as within different mutants lacking the α-helix, β-hairpin, and both. This unique computational approach allowed us to identify residues that are likely to be involved in the transmission of the allosteric signal between the α-helix, β-hairpin elements and the plug. These residues were investigated by site-directed mutagenesis, functional studies, and planar lipid bilayer electrophysiology. The results confirmed the involvement of 4 of the 5 distinct communities of residues in modulating the usher's channel activity and gating, suggesting that they all participate in the allosteric mechanism controlling plug displacement.</p></sec><sec id="s2" sec-type="results"><title>Results</title><p>To investigate if the β-hairpin or α-helix (or both) of the TD (residues 146–637 in the full length PapC) have a role in the allosteric communication leading to the displacement of the PD (residues 264–324), we performed four independent MD simulations, corresponding to the PapC TD model (sim1, <xref ref-type="table" rid="tbl1">Table 1</xref>) and three mutants embedded in a mixed lipid bilayer (<xref ref-type="table" rid="tbl1">Table 1</xref>): (i) where the region corresponding to the hairpin between β5 and β6 (residues 233–240) is deleted (sim2); (ii) where the α-helix between β13 and β14 (residues 447–460) is removed; and where both the regions were removed (sim4). The last 50 ns of simulation were considered for analysis, where the averaged root-mean-square deviation of C<sub>α</sub> atoms (C<sub>α</sub>-RMSD) from the averaged structures stabilized around 2.00 ± 0.09 Å, 1.80 ± 0.09 Å, 1.86 ± 0.11 Å, and 2.03 ± 0.10 Å, for the <italic>native</italic> (sim1), <italic>hairpin mutant</italic> (sim2), <italic>helix mutant</italic> (sim3), and <italic>helix-hairpin mutant</italic> (sim4), respectively (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1</xref>). This timescale, although limited for a full exploration of the structural changes induced by the mutations, was informative in revealing how local structural perturbations may affect allosteric changes leading to the plug displacement in PapC TD.<table-wrap id="tbl1" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.005</object-id><label>Table 1.</label><caption><p>Summary of the simulations</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.005">http://dx.doi.org/10.7554/eLife.03532.005</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Simulation</th><th>Model systems</th><th>Length (ns)</th></tr></thead><tbody><tr><td>Sim1</td><td>Native PapC TD</td><td>72</td></tr><tr><td>Sim2</td><td>Hairpin mutant</td><td>70</td></tr><tr><td>Sim3</td><td>Helix mutant</td><td>70</td></tr><tr><td>Sim4</td><td>Helix-hairpin mutant</td><td>70</td></tr></tbody></table><table-wrap-foot><fn><p>Descriptions of the items are: Simulation, the name of the simulation; Model systems, PapC TD model systems simulated; and Length, the length of the simulation.</p></fn></table-wrap-foot></table-wrap></p><sec id="s2-1"><title>Non-covalent interaction network in the native PapC translocation domain and its perturbation in the absence of the β-hairpin, α-helix, or both</title><p>The changes in the non-covalent interactions (hydrogen bonds and salt bridges) between all residue pairs were analysed within the native TD by calculating their <italic>non-covalent interaction score</italic> (NCI score) (see ‘Materials and methods’). A non-covalent residue–residue interaction network (RIN) comprising 492 nodes (residues) and 1350 edges (interactions) was then constructed as a weighted undirected graph for the native TD (<xref ref-type="fig" rid="fig2">Figure 2</xref>) and the three mutant systems (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A–C</xref>), with the weight for each edge given by the corresponding NCI score (<xref ref-type="table" rid="tbl2">Table 2</xref>). All four RINs have properties typical of small-world networks (<xref ref-type="bibr" rid="bib5">Atilgan et al., 2004</xref>; <xref ref-type="bibr" rid="bib21">Haiyan and Jihua, 2009</xref>; <xref ref-type="bibr" rid="bib58">Taylor, 2013</xref>), with significant higher clustering coefficient compared to a corresponding random network and a higher mean short path length (<xref ref-type="table" rid="tbl2">Table 2</xref>). Within the constructed non-covalent native RIN, we identified 246 weak-to-strong interactions (connecting 362 nodes) with an NCI score of at least 0.3. Among these, 231 nodes connected by 133 edges showed an NCI score greater than 0.6 (i.e., strong interaction) of which 78 involve residues that are part of the barrel strands (58.6%).<fig-group><fig id="fig2" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.006</object-id><label>Figure 2.</label><caption><title>The native TD and its non-covalent interaction network (non-covalent RIN).</title><p>(<bold>A</bold>) Ribbon representation of the starting model of the native translocation domain (TD) of PapC with the labels ‘N’ and ‘C’ indicating the N and C termini of the translocation channel. The β-barrel, PD, P-linker1, P-linker2, β-hairpin, and α-helix (including the H-linkers) are coloured grey, magenta, light purple, dark purple, orange, and yellow, respectively. The α-helix, β-hairpin, P-linker1, P-linker2, and PD are labelled. (<bold>B</bold>) Non-covalent RIN representation of the native translocation domain (TD) of PapC visualized with Cytoscape 2.8.2 (<xref ref-type="bibr" rid="bib55">Smoot et al., 2011</xref>) based on RINalyzer plug-in analysis (<xref ref-type="bibr" rid="bib15">Doncheva et al., 2011</xref>) (see <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref> for the RINs of the TD mutants). The nodes (representing residues) are coloured by structural element as in (<bold>A</bold>) Edges (connecting two residues) are shown in blue, the edge width is proportional to its NCI score from lower to higher values.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.006">http://dx.doi.org/10.7554/eLife.03532.006</ext-link></p></caption><graphic xlink:href="elife03532f002"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.03532.007</object-id><label>Figure 2—figure supplement 1.</label><caption><title>Non-covalent interaction network (non-covalent RIN) for PapC TD mutants.</title><p>Non-covalent RIN representation of the hairpin mutant (<bold>A</bold>), helix mutant (<bold>B</bold>), and helix-hairpin mutant (<bold>C</bold>) translocation domain (TD) of PapC visualized as in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.007">http://dx.doi.org/10.7554/eLife.03532.007</ext-link></p></caption><graphic xlink:href="elife03532fs002"/></fig><fig id="fig2s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.03532.008</object-id><label>Figure 2—figure supplement 2.</label><caption><title>Combined RIN of the difference in non-covalent interaction score (ΔNCI score).</title><p>RINs showing the difference in non-covalent interaction score (ΔNCI score) between the native TD system and the hairpin mutant (<bold>A</bold>), helix mutant (<bold>B</bold>), and helix-hairpin mutant (<bold>C</bold>). The combined RIN (<bold>D</bold>) was created by merging the three RINs. The nodes (representing residues) are coloured as in <xref ref-type="fig" rid="fig2">Figure 2A</xref>. Edges (connecting two residues) are shown in blue, with edge width proportional to its corresponding ΔNCI score (from lower to higher values).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.008">http://dx.doi.org/10.7554/eLife.03532.008</ext-link></p></caption><graphic xlink:href="elife03532fs003"/></fig></fig-group><table-wrap id="tbl2" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.009</object-id><label>Table 2.</label><caption><p>Summary of the residue–residue interaction networks (RINs) parameter</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.009">http://dx.doi.org/10.7554/eLife.03532.009</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>RIN</th><th>Full RIN</th><th>C</th><th>Cr</th><th>C/Cr</th><th>L</th><th>Lr</th><th>L/Lr</th></tr></thead><tbody><tr><td>Native PapC TD</td><td>1350 (492)</td><td>0.384</td><td>0.012</td><td>32.00</td><td>6.67</td><td>3.78</td><td>1.76</td></tr><tr><td>Hairpin mutant</td><td>1196 (485)</td><td>0.368</td><td>0.011</td><td>33.45</td><td>7.20</td><td>3.90</td><td>1.84</td></tr><tr><td>Helix mutant</td><td>1225 (476)</td><td>0.362</td><td>0.011</td><td>32.90</td><td>6.67</td><td>3.90</td><td>1.71</td></tr><tr><td>Helix-hairpin mutant</td><td>854 (466)</td><td>0.262</td><td>0.008</td><td>32.75</td><td>8.10</td><td>4.70</td><td>1.72</td></tr></tbody></table><table-wrap-foot><fn><p>Descriptions of the items are: RIN, residue–residue interaction networks of the different model systems; Full RIN, number of edges in the RIN, in parenthesis the number of node; C, average clustering coefficient; Cr, average clustering coefficient for the random networks with the same size; C/Cr, average clustering coefficient ratio (as used in <xref ref-type="bibr" rid="bib5">Atilgan et al., 2004</xref>); L, average shortest path length; Lr, average shortest path length for the random networks with the same size; L/Lr, average shortest path length ratio (as used in <xref ref-type="bibr" rid="bib5">Atilgan et al., 2004</xref>).</p></fn></table-wrap-foot></table-wrap></p><p>Comparative analysis between the RINs of native and mutants systems revealed slight changes, suggesting a rearrangement in the interaction network. To better understand the mutation-induced changes in network components, we calculated the difference in non-covalent interaction score (ΔNCI score) between the native TD system and each of the mutant systems (the weakened interactions are shown in <xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2A–C</xref>). This information was then added as a weighted undirected edge to the pre-existing native non-covalent RIN (the ΔNCI edges are shown in <xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2D</xref>). Interestingly, 24% of the strong interactions in the native RIN were weakened relative to the RIN of the mutant lacking the β-hairpin, 22.6% relative to the mutant lacking the α-helix, and 23.3% relative to the mutant lacking both, suggesting that interactions between nodes that are not part of the deleted secondary structure elements were consistently weakened in the absence of these elements.</p></sec><sec id="s2-2"><title>Evolutionary analysis of PapC TD</title><p>We first extracted evolutionary information from a multiple sequence alignment of the PapC TD family. The patterns of conservation in the TD using Consurf (<xref ref-type="bibr" rid="bib3">Ashkenazy et al., 2010</xref>) analysis suggested that the highly conserved residues (score 9) tend to be clustered in two specific regions of the usher (<xref ref-type="fig" rid="fig3">Figure 3A</xref>). The first cluster mapped onto the PD and the P-linkers (P-linker1 residues 248–263; P-linker2 residues 325–335) connecting it to the TP. The second cluster (which included the majority of the highly-conserved residues) mapped onto one side of the TP (strand β1–14 and β24). It includes residues: (i) near the periplasmic side of the β-barrel within β1–4 strands and β24 strand; (ii) on the extracellular side of the barrel (within β5–10); (iii) in the β-hairpin region (β-hairpin and β7–9); and (iv) in the area of β10–14 capped by the α-helix region, which comprises the α-helix and its linkers—H-linker1 (residues 445–450) and H-linker2 (461–468, respectively). Surface representation of the TD reveals a continuous patch of conserved residues facing the lipid bilayer, including β13, the extracellular half of β14 and the periplasmic half of β12 (<xref ref-type="fig" rid="fig3">Figure 3B</xref>). Intriguingly, this patch (‘β13 conserved patch’) reaches the full height of the pore from the α-helix region to a functionally important loop located between β12 and β13 strands (<xref ref-type="bibr" rid="bib16">Farabella, 2013</xref>; <xref ref-type="bibr" rid="bib63">Volkan et al., 2013</xref>).<fig id="fig3" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.010</object-id><label>Figure 3.</label><caption><title>Evolutionary analysis of PapC TD.</title><p>(<bold>A</bold>–<bold>B</bold>) Sequence conservation calculated with Consurf (<xref ref-type="bibr" rid="bib3">Ashkenazy et al., 2010</xref>) and mapped onto the initial model of the native PapC TD (sim1, t = 0). Amino acid conservation scores are classified into nine levels. The colour scale for residue conservation goes from cyan (non-conserved: grade 1) to maroon (highly conserved: grade 9), unreliable positions are coloured light yellow. (<bold>A</bold>) Ribbon representation of the model with the highly conserved residues (grade9) shown as spheres and key elements labelled. (<bold>B</bold>) Molecular surface of the model with β12–β14 labelled. (<bold>C</bold>–<bold>D</bold>) Sequence co-evolution calculated with PyCogent (<xref ref-type="bibr" rid="bib32">Knight et al., 2007</xref>; <xref ref-type="bibr" rid="bib10">Caporaso et al., 2008</xref>). (<bold>C</bold>) The co-evolving residues are mapped onto the initial model of the native PapC TD (sim1, t = 0). (<bold>D</bold>) The co-evolution network as visualized with Cytoscape 2.8.2 Cytoscape 2.8.2 (<xref ref-type="bibr" rid="bib55">Smoot et al., 2011</xref>) based on RINalyzer plug-in analysis (<xref ref-type="bibr" rid="bib15">Doncheva et al., 2011</xref>). Edges (connecting two co-evolved residues) are shown in blue, and nodes (representing coevolved residues) are coloured by structural element. The PD, P-linker1, P-linker2, β-hairpin, and α-helix are indicated schematically and coloured as in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The node size is proportional to its degree of connectivity.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.010">http://dx.doi.org/10.7554/eLife.03532.010</ext-link></p></caption><graphic xlink:href="elife03532f003"/></fig></p><p>In addition to investigating conservation, we performed an analysis to identify the coevolutionary relationships between residues in the structure. Using normalized mutual information (NMI) analysis (<xref ref-type="bibr" rid="bib38">Martin et al., 2005</xref>) with a Z-score cut-off = 4 (see ‘Materials and methods’) to detect the intra-molecular coevolved residues within PapC TD, a coevolutionary RIN containing 100 coevolved residues (nodes) and 357 connections (edges) was derived (<xref ref-type="fig" rid="fig3">Figure 3D</xref>). Mapping the network onto the PapC TD structure showed that many of the residues involved are also connected spatially and are clustered in the same regions where the highly-conserved residues were found (P-linkers, the PD, and the barrel wall capped by the α-helix, in close proximity to the β-hairpin) (<xref ref-type="fig" rid="fig3">Figure 3C,D</xref>). The obtained coevolutionary RIN showed a significant clustering coefficient compared to a corresponding random network (of 0.493 vs 0.187, respectively) and a comparable mean short path length (3.15 vs 2.57, respectively) (<xref ref-type="bibr" rid="bib11">Daily et al., 2008</xref>).</p></sec><sec id="s2-3"><title>Identifying allosteric ‘hot spots’ from a hybrid residue interaction network</title><p>We constructed one hybrid RIN in which the attributes for the nodes and edges are defined by the properties described above (non-covalent networks and evolutionary analysis, see ‘Materials and methods’). Starting from the secondary structure elements (that uniquely characterise the barrel–the α-helix and β-hairpin) in this hybrid RIN, we used a multi-step procedure to reconstruct a pathway of communication between them (<xref ref-type="fig" rid="fig4">Figure 4</xref>).<fig-group><fig id="fig4" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.011</object-id><label>Figure 4.</label><caption><title>Detection of allosteric hot spots.</title><p>A flowchart representing the multistep procedure used to identify allosteric hot spots. First, a sub-network of the protein hybrid RIN was generated starting from the α-helix and β-hairpin. Then, filters based on the evolutionary information and on the interactions analysis were applied (see <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1</xref>) resulting in a sub-network of ‘hot spot’ residues.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.011">http://dx.doi.org/10.7554/eLife.03532.011</ext-link></p></caption><graphic xlink:href="elife03532f004"/></fig><fig id="fig4s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.03532.012</object-id><label>Figure 4—figure supplement 1.</label><caption><title>Contribution of each filter in the detection of allosteric hot spots.</title><p>Venn diagrams illustrating the contribution of each filter in the node set detection. (<bold>A</bold>) The dynamic filter resulting from the intersection of the ΔNCI score between the native TD and each of the mutants. (<bold>B</bold>) The relative combination of the dynamic filter set and the set from the evolutionary filter to the final hot spots sub-network.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.012">http://dx.doi.org/10.7554/eLife.03532.012</ext-link></p></caption><graphic xlink:href="elife03532fs004"/></fig></fig-group></p><p>This initial large sub-network is formed by 208 nodes (residues) connected by 456 NCI edges (in the native RIN). Applying the dynamic filter (independently) on the edges, based on the difference in non-covalent interaction score between the native TD and each of the mutants (ΔNCI > 0), revealed that in each case a large part of network has weakened interactions (hairpin mutant: 200 nodes, 438 NCI edges; helix mutant: 199 nodes, 437 NCI edges; helix-hairpin: 202 nodes, 443 NCI edges) (<xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1A</xref>). The application of the evolutionary filter revealed that only a small part of the sub-network is made of evolutionary important residues (75 nodes connected by 104 native NCI edges). Combining the filters (<xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1B</xref>) resulted in 69 nodes connected by 100 NCI edges (thus representing interacting residues in the native PapC network). The residues of this sub-network (14% of all residues in the TD) were considered ‘hot spots’ in the communication pathway of PapC TD (‘hot spot’ sub-network). Mapping them onto the structure revealed that they are located close together in a continuous area within PapC TD.</p><p>We analysed the community structure of the hot spot sub-network using the edge-betweenness clustering algorithm (<xref ref-type="bibr" rid="bib19">Girvan and Newman, 2002</xref>; <xref ref-type="bibr" rid="bib41">Morris et al., 2010</xref>). This analysis shows that the sub-network has a modular structure, with a modularity index of 0.73 (maximum value of the modularity index is 1), which is typical of 3D-structure based RIN (<xref ref-type="bibr" rid="bib44">Newman and Girvan, 2004</xref>; <xref ref-type="bibr" rid="bib53">Sethi et al., 2009</xref>). Here, a total of 11 communities containing two or more residues were identified, from which only five communities are composed by more than five residues. For further analysis, we chose to consider only these five largest communities, which are located: between β<italic>7–</italic>9 and the P-linkers (C1); between the β-hairpin and the conserved region at the base of the α-helix (β12–14) (C2); between β12<italic>–</italic>β13_loop and the P-linker1 (C3); between the β-hairpin, P-linker2 and the PD (βE–F) (C4); and on the tip of the PD (βE–F loop and βA–B loop) (C5) (<xref ref-type="table" rid="tbl3">Table 3</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref>).<table-wrap id="tbl3" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.013</object-id><label>Table 3.</label><caption><p>Communities in the hot spot sub-network</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.013">http://dx.doi.org/10.7554/eLife.03532.013</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Community</th><th>Residues</th></tr></thead><tbody><tr><td>C1</td><td>E247, D249, Y329, L330, T331, G334, Q335, R337, K339, E361, S363, W364, G365, L366, S371, L372</td></tr><tr><td>C2</td><td>R237, D402, S420, Y441, R442, F443, S444, K468, E469, M470, E475, W496</td></tr><tr><td>C3</td><td>Y260, Y425, S426, K427, T437, F438, A439</td></tr><tr><td>C4</td><td>S233, R303, G304, L306, V308, F320, T324, A325, V327</td></tr><tr><td>C5</td><td>E269, E312, N314, G315, R316, K318</td></tr></tbody></table><table-wrap-foot><fn><p>Descriptions of the items are: Community, the name of the community; Residues, residues that are part of the community.</p></fn></table-wrap-foot></table-wrap><fig id="fig5" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.014</object-id><label>Figure 5.</label><caption><title>PapC TD communities.</title><p>The communities of the hot spot sub-network are shown as surface by colours and indicated schematically (C1 to C5). The inset shows a close up of the identified core residues located in β<italic>7,</italic> β8, the P-linkers, the β-hairpin, the conserved region at the base of the α-helix, in the junction between β12–β13_loop. The core residues are labelled in bold and numbered according to the X-ray structure of the apo PapC TD (PDB id: 2vqi).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.014">http://dx.doi.org/10.7554/eLife.03532.014</ext-link></p></caption><graphic xlink:href="elife03532f005"/></fig></p><p>We selected a number of key residues from the communities (core hot spot residues), which link different elements within each community, for further experimental investigations (<xref ref-type="fig" rid="fig5">Figure 5</xref>). These were found in communities C1–C4: in C1, residues linking the P-linkers and the barrel wall that possibly help in maintaining the P-linkers in a closed configuration (P-linker1:D249, P-linker2:Y329, P-linker2:T331, β7:R337, and β8:S363); in C2, residues that bridge the base of the α-helix (the extracellular end of the β13 conserved patch) and β-hairpin (β<italic>-</italic>hairpin:R237, β12:S420, β13:R442, β13:S444); in C3, residues on the interface between P-linker1 and β12–β13_loop (P-linker1:Y260, β12–β13_loop:K427, β13:T437, β13:F438); and in C4, residues that are part of the interface with P-linker1 and the periplasmic end of the β13 conserved patch (P-linker2:A325, P-linker2:V327).</p></sec><sec id="s2-4"><title>Experimental analysis of residues in the hot spot sub-network</title><p>To test experimentally if the key hot spot residues identified above (linking elements within each community) contribute to allosteric signalling within PapC, we constructed a set of single alanine substitution mutations (<xref ref-type="table" rid="tbl4">Table 4</xref>). Each of the mutants was present at a similar level in the OM compared to the wild-type PapC usher, and the mutations did not affect the ability of the usher to form a stable β-barrel in the OM (data not shown). The functionality of the PapC substitution mutants was assessed by ability to assemble P pili on the bacterial surface. P pili bind to receptors on human red blood cells, and assembly of functional P pili was determined using a hemagglutination assay (HA). Seven (D249A, T331A, R442A, S444A, Y260A, K427A, T437A) of the 14 tested mutants exhibited greater than twofold defects in agglutination titers compared with wild-type PapC, with 4 of the mutants (D249A, R442A, Y260A, and K427A) exhibiting no agglutination activity (HA titer = 0) (<xref ref-type="table" rid="tbl4">Table 4</xref>). The defective mutants were in key residues from communities C1, C2, and C3, confirming roles for these communities in proper usher function.<table-wrap id="tbl4" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.015</object-id><label>Table 4.</label><caption><p>Analysis of PapC substitution mutants</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.015">http://dx.doi.org/10.7554/eLife.03532.015</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>PapC</th><th>Community</th><th>HA titer</th><th colspan="3">Antibiotic sensitivity</th></tr><tr><th/><th/><th/><th>SDS</th><th>Erythromycin</th><th>Vancomycin</th></tr></thead><tbody><tr><td>WT</td><td/><td>64</td><td>15</td><td>6</td><td>6</td></tr><tr><td>D249A</td><td>C1</td><td>0</td><td>15</td><td>6</td><td>6</td></tr><tr><td>Y329A</td><td>C1</td><td>32</td><td>15</td><td>6</td><td>6</td></tr><tr><td>T331A</td><td>C1</td><td>24</td><td>15</td><td>15</td><td>10</td></tr><tr><td>R337A</td><td>C1</td><td>64</td><td>15</td><td>6</td><td>6</td></tr><tr><td>S363A</td><td>C1</td><td>64</td><td>15</td><td>6</td><td>6</td></tr><tr><td>R237A</td><td>C2</td><td>64</td><td>16</td><td>6</td><td>15</td></tr><tr><td>S420A</td><td>C2</td><td>32</td><td>15</td><td>6</td><td>6</td></tr><tr><td>R442A</td><td>C2</td><td>0</td><td>15</td><td>6</td><td>6</td></tr><tr><td>S444A</td><td>C2</td><td>24</td><td>14</td><td>6</td><td>6</td></tr><tr><td>Y260A</td><td>C3</td><td>0</td><td>15</td><td>14</td><td>6</td></tr><tr><td>K427A</td><td>C3</td><td>0</td><td>14</td><td>6</td><td>6</td></tr><tr><td>T437A</td><td>C3</td><td>24</td><td>14</td><td>6</td><td>6</td></tr><tr><td>F438A</td><td>C3</td><td>32</td><td>14</td><td>12</td><td>6</td></tr><tr><td>V327A</td><td>C4</td><td>64</td><td>20</td><td>14</td><td>16</td></tr></tbody></table><table-wrap-foot><fn><p>Descriptions of the items are: PapC, the PapC construct tested; Community, the name of the community to which the mutated residue belongs; HA (hemagglutination assay) titer, the maximum fold dilution of bacteria able to agglutinate human red blood cells; Antibiotic sensitivity, the diameter of zone of inhibition (mm) around filter disc impregnated with SDS (750 µg), erythromycin (15 µg), or vancomycin (20 µg). The antibiotic sensitivity measurement includes the filter disc (6 mm diameter).</p></fn></table-wrap-foot></table-wrap></p><p>We next used an antibiotic sensitivity assay to screen the PapC substitution mutants for effects on channel activity of the usher. The OM of Gram-negative bacteria has low permeability to detergents such as SDS and to antibiotics such as erythromycin and vancomycin, providing resistance to these molecules. In its resting state, the usher TP is gated closed by the PD, preserving integrity of the OM. Mutations that disrupt channel gating by the PD will result in the opening of the large TD channel, leading to increased sensitivity of the bacteria to antibiotics. Bacteria expressing five of the PapC substitution mutants (T331A, R237A, Y260A, F438A, and V327A) exhibited increased sensitivity to one or more of the tested molecules (<xref ref-type="table" rid="tbl4">Table 4</xref>). Y329A, R337A, S363A, and S420A did not appear to perturb the allosteric signalling within PapC, showing the same antibiotic sensitivity phenotype and ability to assemble pili of the native PapC (<xref ref-type="table" rid="tbl4">Table 4</xref>). However, the hemagglutination assay and antibiotic sensitivity assay are screening tools, and as such, they lack in sensitivity to pick up smaller changes in the ability to assemble pili or channel activity of the usher.</p></sec><sec id="s2-5"><title>Electrophysiological analysis of selected mutants</title><p>In total, 10 mutated PapC TDs were found to be affected either in their ability to trigger hemagglutination or in their permeability to SDS or antibiotics (T331A, D249A in C1; R237A, R442A, S444A in C2; Y260A, K427A, T437A, F438A in C3; and V327A in C4). We attempted to purify those mutants in view of examining their channel activity using planar lipid bilayer electrophysiology (which is a more sensitive assay). Unfortunately, only seven of these 10 mutants yielded protein stable enough (as wild-type) in detergent solutions to carry out the planned experiments (T331A in C1; R237A, S444A in C2; K427A, T437A, F438A in C3; and V327A in C4). During the OM extraction procedure, D249A, R442A, and Y260A were not stable enough due to the loss of the membrane bilayer environment and inability to maintain their native conformation in detergents. Insertion of PapC purified proteins (see ‘Materials and methods’) was promoted in planar lipid bilayers by clamping the membrane potential to −90 mV. As soon as channel activity was observed, the potential was briefly returned to zero and the chamber stirring stopped to minimize further insertions. 10-min long recordings of channel activity at + and −90 mV, and at + and −50 mV were performed.</p><p>The typical electrophysiological signature of the wild-type PapC usher is characterized by prolonged dwell times at a low current level, representing the closed state of the usher, and brief transitions of various current amplitudes. These transitions represent short-lived openings of various conductance, ranging from 50 to 600 pS (‘transient-mixed’ behaviour, TM) (<xref ref-type="fig" rid="fig6">Figure 6A</xref>). Although it is not possible to know exactly how many individual pores were inserted into the bilayer, the observed fluctuations of various sizes are taken to represent various conformational states of a single pore. As documented previously, the openings of the ‘transient-mixed’ behaviour appear rather small and may be due to the jiggling of the plug within the TP and/or the thermally induced mobility of various domains of the protein, such as the NTD and CTDs or loops (<xref ref-type="bibr" rid="bib36">Mapingire et al., 2009</xref>). Occasionally, and more so at higher membrane potential, very large and sustained openings (‘large-open’ behaviour, LO) are observed in wild-type PapC usher (<xref ref-type="fig" rid="fig6">Figure 6B</xref>). These openings have a conductance of ∼3–4 nS, which is similar to the monomeric conductance of the mutant lacking the PD (<xref ref-type="bibr" rid="bib36">Mapingire et al., 2009</xref>) and are interpreted as representing a full displacement of the PD from a single monomer. Prolonged opening of intermediate conductance (0.5–1 nS) can also be observed and may represent partial PD displacement.<fig id="fig6" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.016</object-id><label>Figure 6.</label><caption><title>Kinetic signatures of channel activity in wildtype and mutant PapC ushers and frequency of PD displacement.</title><p>Fifty-second segments of recordings obtained in planar lipid bilayers were selected to illustrate the behaviour of the different proteins. (<bold>A</bold>) Recording from the wild-type PapC usher showing the characteristic ‘transient-mixed’ behaviour. (<bold>B</bold>) Recording from the wild-type PapC usher showing an example of spontaneous large openings due to plug displacement. Note the large amount of current fluctuations during the openings, and the ‘transient-mixed’ behaviour in between such events. Examples of similar large openings (<bold>C</bold> and <bold>D</bold>) are shown for the V327A and T331A mutants, respectively. (<bold>E</bold>) A recording from the K427A shows that the channel barely displays any activity at this voltage. The voltage was +90 mV for all panels. The current level for the closed channels is marked as ‘C’, and openings are seen as upward deflections of the traces; current levels corresponding to fully open monomeric or dimeric forms are denoted by ‘M’ and ‘D’, respectively. Note that the traces are plotted as conductance, rather than current, vs time and the scale bars are given in nS. (<bold>F</bold>) The percent of sweeps displaying ‘large-open behaviour’ (LO) indicative of PD displacement is shown for WT and each mutant at the indicated voltages. The number of individual bilayers investigated in each case is given above the bars.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.016">http://dx.doi.org/10.7554/eLife.03532.016</ext-link></p></caption><graphic xlink:href="elife03532f006"/></fig></p><p>Because the electrophysiological behaviour of wild-type PapC usher is quite variable, and in attempt to quantify the propensity at spontaneous PD displacement, we have counted the number of 10-min long recordings (sweeps) that show ‘large-open’ behaviour, and we report the percent of such sweeps in various conditions. The frequency of observing these large openings in wild-type PapC usher is ∼20% at ±50 mV, but increases to ∼60% at ±90 mV. The application of a larger transmembrane voltage is likely to disrupt the interactions between key residues involved in keeping the PD in place, leading to a more frequent spontaneous displacement of the latter.</p><p>3 of the 7 analysed mutants, V327A, T331A, and F438A, showed an increased propensity at displaying large openings, relative to the wild-type PapC usher, as illustrated for V327A and T331A (<xref ref-type="fig" rid="fig6">Figure 6C,D</xref>). This was particularly true at ±90 mV where the percent of sweeps with large openings reaches values of 75–90% (<xref ref-type="fig" rid="fig6">Figure 6F</xref>). The T331A mutant was consistently more prone to open than WT and any other mutants, which led to the occasional simultaneous opening of several monomers (<xref ref-type="fig" rid="fig6">Figure 6D</xref>). Two of the seven mutants R237A and T437A still opened occasionally to the 3 nS level, but the frequency of sweeps with such events was slightly diminished relative to the wild-type PapC usher at ±90 mV (<xref ref-type="fig" rid="fig6">Figure 6F</xref>) suggesting that these mutants are likely to be insignificantly different from the wild-type PapC usher.</p><p>The K427A and S444A mutants showed a decreased frequency (or complete absence) of large openings. The K427A mutant almost never showed ‘large-open’ behaviour at ±50 and ±90 mV, indicating an extremely closed channel (<xref ref-type="fig" rid="fig6">Figure 6E</xref>, <xref ref-type="fig" rid="fig6">Figure 6E,F</xref>). The S444A mutant was even less prone to open, with 0% occurrence of PD displacement at ±50 mV in the 11 bilayers that we investigated (<xref ref-type="fig" rid="fig6">Figure 6F</xref>). However, increased activity with fast flickers and occasional more prolonged openings could be seen for both mutants if the membrane potential was switched to voltages in the ±100–150 mV range, indicating that the channels are present in the bilayer, but require higher voltages for activation.</p></sec></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>PapC usher catalyses the translocation across the outer membrane of P pili, and its gating mechanism is important for bacterial homeostasis and for catalysis of pilus assembly. The TD of PapC is formed by the largest β-barrel pore known to be formed by a single chain. The PD occludes the pore in an inactive state and maintains the permeability of the channel. As previously documented, the native PapC channel is highly dynamic and is characterized by spontaneous short-lived openings of various conductance levels (<xref ref-type="bibr" rid="bib36">Mapingire et al., 2009</xref>). Two distinct structural elements, the β-hairpin and the α-helix, play an important role in maintaining the PD in a closed conformation. In the absence of both elements the usher is defective for pilus biogenesis. Our analysis of the non-covalent interaction RIN in the native PapC TD shows that the interactions found between the TP and the PD are mostly weak, possibly to allow an easy release of the PD. This finding supports the idea that the highly dynamic behaviour of the native PapC channel is originated from the ‘jiggling’ of the PD within the TP (<xref ref-type="bibr" rid="bib36">Mapingire et al., 2009</xref>). On the other hand, we find that the interactions between the TP and the P-linkers, between the TP and the β-hairpin, and between the TP and the α-helix region are mainly stable, supporting their role in maintaining the PD in a closed conformation.</p><p>Analysis of the mutation-induced perturbation of the non-covalent interaction RIN of the native PapC TD in the absence of the β-hairpin, α-helix, or both, shows that the interactions between nodes that are not part of the deleted secondary structure elements are consistently weakened (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref>). This feature suggests that the two elements are not independent and they are part of a complex allosteric process regulating the PD gating mechanism. It has been proposed that only a few residues play essential roles during allosteric processes and that perturbing the interactions between these residues can facilitate the population shift of the conformational ensembles (<xref ref-type="bibr" rid="bib14">Del Sol et al., 2009</xref>; <xref ref-type="bibr" rid="bib61">Tsai et al., 2009</xref>). Additionally, it has been shown that residues with a pivotal role in mediating such allosteric communications are also important in terms of evolution (both highly coevolved residues and conserved residues) (<xref ref-type="bibr" rid="bib56">Suel et al., 2003</xref>; <xref ref-type="bibr" rid="bib17">Ferguson et al., 2007</xref>; <xref ref-type="bibr" rid="bib57">Tang et al., 2007</xref>). Remarkably, we show here that PapC is characterised by an uneven distribution of the evolutionary important residues, clustered in the P-linkers, the PD, and the barrel wall capped by the α-helix, in close proximity to the β-hairpin. Another interesting finding is the presence of the ‘β13 conserved patch’ that reaches the full height of the pore from the α-helix region to the loop located between β12 and β13 strands (β12–β13_loop), which has been recently identified as important (<xref ref-type="bibr" rid="bib16">Farabella, 2013</xref>; <xref ref-type="bibr" rid="bib63">Volkan et al., 2013</xref>).</p><p>To detect the allosteric network, we implement a new method that integrates dynamic and evolutionary information in a hybrid RIN and then apply network analysis. This approach allows us to explore a large part of the protein, resulting in the detection of only 14% of all residues in the TD as potential candidates. It has been shown that detecting ‘residue communities’ in protein structure networks leads to the identification of key residues that are often part of a signal transduction pathway (<xref ref-type="bibr" rid="bib7">Bode et al., 2007</xref>; <xref ref-type="bibr" rid="bib13">Del Sol et al., 2007</xref>). The interconnection within and between the communities is pivotal for the flow of allosteric signalling. Residues in the same community are densely interconnected and have multiple routes to communicate with one another. However, the interconnections between communities involve only a few edges, which form the bottleneck for the flow of the signal in the network (<xref ref-type="bibr" rid="bib7">Bode et al., 2007</xref>; <xref ref-type="bibr" rid="bib13">Del Sol et al., 2007</xref>; <xref ref-type="bibr" rid="bib53">Sethi et al., 2009</xref>). Here, all the identified communities (C1–C5, <xref ref-type="fig" rid="fig5">Figure 5</xref>) comprise residues from multiple elements (e.g., the β-hairpin, P-linker1, P-linker2, and distinct part of the TP) of PapC TD, except C5 that is composed only by PD residues. Mutation of key residues linking elements within each of the communities C1–C4 shows an altered antibiotic sensitivity phenotype, confirming a role of these communities in the pore-gating mechanism. Additionally, communities C1–C3 are required for proper usher function (as estimated by the hemagglutination assay), suggesting dependency of the pore gating function and pilus assembly (based on some of the mutations). For example, mutations in β12–β13_loop:K427 in C3 and β13:S444 in C2, which lead to a drastic decrease in the frequency of plug displacement (stabilizing the closed state, based on channel activity analysis), also show defective pilus assembly (possibly due to their deficiency in relocating the plug in a functional conformation).</p><p>Additionally, our study identifies two residues with opposite effects on plug displacement in the same community: C3: β12–β13_loop:K427, which is located on β12–β13_loop, and β13:F438, which is located at the periplasmic end of the β13 conserved patch (linking β12–β13_loop with the α-helix region). This observation indicates that β12–β13_loop has a pivotal role in the regulation of plug displacement by acting as a ‘latch’, as proposed previously (<xref ref-type="bibr" rid="bib16">Farabella, 2013</xref>) and was also shown by mutagenesis studies (<xref ref-type="bibr" rid="bib63">Volkan et al., 2013</xref>). Intriguingly, at the extracellular end of the β13 conserved patch there is another mutation (β13:S444) leading to a channel more reluctant to open, suggesting a regulatory role for the β13 conserved patch (<xref ref-type="fig" rid="fig7">Figure 7</xref>) in modulating the latch.<fig id="fig7" position="float"><object-id pub-id-type="doi">10.7554/eLife.03532.017</object-id><label>Figure 7.</label><caption><title>Residues involved in the allosteric signalling to control PapC gating.</title><p>A schematic model summarizing the location of the detected hot spots involved in the gating mechanism. The β12–β13_loop (the ‘latch’) and the β13 conserved patch are coloured in dark grey and light grey respectively. The PD, P-linker1, P-linker2, β-hairpin, and α-helix are coloured as in <xref ref-type="fig" rid="fig2">Figure 2</xref>. Hot spot residues are colour-coded based on their communities (C1–C4, as in <xref ref-type="fig" rid="fig5">Figure 5</xref>). The different symbols indicate the mutant's electrophysiological behaviour (‘X’ where no data were available). Mutants that show a pilus assembly defect or an increased antibiotic sensitivity, or both, are represented by a triangle, pentagon and rectangle, respectively.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.017">http://dx.doi.org/10.7554/eLife.03532.017</ext-link></p></caption><graphic xlink:href="elife03532f007"/></fig></p><p>Interestingly, our community residues in some cases are found to have very different patterns of interactions in an alternative conformation. For example, the functionally important latch (β12–β13_loop) is shown to be in a different conformation in the open state (based on FimD:FimC:FimH structure [<xref ref-type="bibr" rid="bib49">Phan et al., 2011</xref>]) as well as other residues in C3 and C4 communities. As a result, some residues that are identified by our method, such as V327, could in principle be selected also by visual inspection for further experimental investigation. V327 (located on P-linker2 in C4 at interface with P-linker1 and the barrel wall) is found to have a different interaction pattern in the open state of the usher and is also shown to promote plug displacement. However, more importantly, the method is able to predict residues in regions that would not attract our attention at all but may contain important information in respect to the allosteric pathway. (This aspect can become even more significant in the absence of an alternative conformation.) For example, we predict two such residues to be functionally important—T331 and S444. T331 is located on P-linker2 at the interface between P-linker2 and the barrel wall, which stays intact in the open conformation of the usher; S444 is located on β13 at the base of the α-helix, that is, on the barrel itself, and it lacks any direct contact with the plug domain. Surprisingly, and in support to our prediction, both residues are found to have an effect on plug displacement. Thus, the main strength of our method is that the knowledge of the structure of the protein (in one conformation only) and enough sequence information (to extract evolutionary information) are sufficient for the detection of functionally important residues that are pivotal for transferring the regulatory information within the protein.</p><sec id="s3-1"><title>Conclusion</title><p>In this study, we provide a first deep insight into the allosteric regulation of the gating mechanism of the usher family. Using PapC TD as model system, we developed an integrative approach combining computational modelling, sequence conservation analysis, mutual information-based coevolution analysis and information from AA-MD simulations, to study the potential involvement of particular secondary structure elements (the α-helix and β-hairpin) in the allosteric communication. The construction of a hybrid interaction network and the use of network analysis allowed us to identify communities of residues within the TD that potentially mediate this process. Antibiotic sensitivity and electrophysiology experiments on a set of alanine-substitution mutants confirmed that residues located in the P-linkers, the β-hairpin, and β13 conserved patch (part of four communities) alter channel gating and that residues located in P-linker2, β12–β13_loop, and β13 conserved patch (both periplasmic end and extracellular end) are sensitive to plug displacement. Therefore, we suggest that the β13 conserved patch acts as a regulator of β12–β13_loop (the latch), mediating channel opening. Furthermore, our study shows how the integration of different computational approaches based on evolution, structure, and dynamics of proteins, into a hybrid network can unveil communication pathways within proteins. Such an integrative approach can guide the experimental investigation by pinpointing key candidates involved in the transmission of the allosteric signal.</p></sec></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><sec id="s4-1"><title>Systems modelling</title><p>We built four model systems based on the X-ray structure of the TD (residues 1–492) of PapC at 3.2 Å resolution (PDB ID: 2vqi [<xref ref-type="bibr" rid="bib50">Remaut et al., 2008</xref>]). The starting model for the simulation of the native TD (sim1, native) was generated by adding the missing loops to the X-ray structure using the <italic>dope_loopmodel</italic> method (<xref ref-type="bibr" rid="bib54">Shen and Sali, 2006</xref>) in MODELLER-9v7 (<xref ref-type="bibr" rid="bib52">Sali and Blundell, 1993</xref>). Additionally, three mutant model systems were constructed based on the native model: a mutant lacking the β-hairpin (Δ233–240) (sim2, hairpin mutant), a mutant lacking the α-helix (Δ447–460) (sim3, helix mutant), and a mutant lacking both (Δ233–240 and Δ447–460) (sim4, helix-hairpin mutant).</p><p>Each of the systems was oriented with respect to the membrane normal (the Z axis by definition) using the database (<xref ref-type="bibr" rid="bib35">Lomize et al., 2006</xref>). For the native model (sim1), a mixed lipid bilayer (POPE/POPG 3:1) was generated around the protein using the replacement method (<xref ref-type="bibr" rid="bib28">Jo et al., 2007</xref>). To obtain a mixed lipid bilayer that reproduces an estimated surface area per lipid of 61.5 ± 0.2 Å<sup>2</sup> (<xref ref-type="bibr" rid="bib42">Murzyn et al., 2005</xref>), we used the InflateGRO method (<xref ref-type="bibr" rid="bib31">Kandt et al., 2007</xref>).</p></sec><sec id="s4-2"><title>Systems set up for MD simulations</title><p>All MD simulations were performed using Gromacs 4.0.5 (<xref ref-type="bibr" rid="bib62">Van Der Spoel et al., 2005</xref>). TIP3P parameters were used for water molecules (<xref ref-type="bibr" rid="bib29">Jorgensen et al., 1983</xref>), the OPLSA-AA force-field (<xref ref-type="bibr" rid="bib30">Kaminski et al., 2001</xref>) was applied to the protein and ions, and the Berger force-field (<xref ref-type="bibr" rid="bib6">Berger et al., 1997</xref>) to the lipids. All four systems were solvated in water and ions were added to neutralize the total charge (0.15 M NaCl), resulting in more than 75,000 atoms in total. Next, each system was energy-minimized using a steepest descent algorithm in the presence of different position restraints on the protein and the lipid bilayer head-groups, which were gradually removed. In the simulations of the mutant systems (sim2 to sim4), the protein models were embedded in the pre-equilibrated membrane obtained after 15 ns of unrestrained equilibration of the native TD (sim1).</p></sec><sec id="s4-3"><title>Equilibration procedure and production run</title><p>The assembled systems were equilibrated in a multistage process using periodic boundary conditions and a 2 fs time step. Short-range interactions were used with a cut-off of 1 nm. The PME algorithm (<xref ref-type="bibr" rid="bib12">Darden et al., 1993</xref>) was used for long-range electrostatic interactions. All bonds were constrained using the LINCS algorithm (<xref ref-type="bibr" rid="bib25">Hess et al., 1997</xref>; <xref ref-type="bibr" rid="bib24">Hess, 2008</xref>). The first equilibration step was performed in the NVT ensemble, using a restraining force of 1000 kJ/(mol nm<sup>2</sup>) for 0.1 ns on the protein and lipids. The V-rescale thermostat (<xref ref-type="bibr" rid="bib8">Bussi et al., 2007</xref>) was employed to couple the temperature of the system to 310 K with a time constant of t<sub>T</sub> = 0.1 ps. All the following equilibration steps were performed in the NPT ensemble. During the next three steps, different parts of the system were restrained using a force constant of 1000 kJ/(mol nm<sup>2</sup>): the protein and lipids, the protein atoms only, and the protein backbone atoms.</p><p>The resulting model of each system was then simulated without restraints. Constant temperature of 310 K was maintained using the Nose–Hoover thermostat (<xref ref-type="bibr" rid="bib26">Hoover, 1985</xref>; <xref ref-type="bibr" rid="bib46">Nose, 1984</xref>) with a time constant of t<sub>T</sub> = 0.1 ps. Using semi-isotropic coupling with a Parrinello–Rahman barostat (<xref ref-type="bibr" rid="bib47">Parrinello, 1981</xref>), a constant pressure of 1 bar was applied with a coupling constant (t<sub>P</sub>) of 1 ps and a compressibility 4.5 e<sup>−5</sup> bar<sup>−1</sup>. Each unrestrained simulation was performed for ∼70–72 ns. The last 50 ns of simulation was used for analysis.</p></sec><sec id="s4-4"><title>Non-covalent residue interaction network</title><p>Hydrogen bonds were defined using a cut-off of 30° for the acceptor–donor–hydrogen angle and a cut-off of 3.5 Å for the hydrogen-acceptor distance. The definition of salt-bridges was based on a 4 Å distance cut-off between any oxygen atoms of acidic residues and nitrogen atoms of basic residues. The non-covalent interaction score (NCI score) of the identified bonds was defined as the percentage of simulation time during which a bond occurs between two amino acids normalised by the number of bonds. Using the normalized score, a non-covalent residue interaction network (RIN) was built for each of the simulated systems as a weighted undirected graph, in which each node represents a residue and each edge is weighted by the normalized score. The difference in non-covalent interaction score (ΔNCI score) between the native TD system and each of the mutant systems was then calculated and added as a weighted undirected edge to the pre-existing non-covalent RIN.</p></sec><sec id="s4-5"><title>Sequence conservation analysis</title><p>The sequence corresponding to the X-ray structures of PapC TD (PDB id: 2vqi; Uniprot id: P07110) was used as input to psiBlast resulting in a set of unique related sequences from the non-redundant NCBI data set (<xref ref-type="bibr" rid="bib1">Altschul et al., 1997</xref>). The E-value threshold was set as 10<sup>−3</sup> and sequences with id >90% and <30% sequence identity were excluded. A structure-based multiple sequence alignment was calculated using Expresso (3DCoffee) (<xref ref-type="bibr" rid="bib2">Armougom et al., 2006</xref>). Finally, an evolutionary conservation score was calculated for each residue an empirical Baysian inference method (<xref ref-type="bibr" rid="bib39">Mayrose et al., 2004</xref>) as implemented in the ConSurf web server (<xref ref-type="bibr" rid="bib3">Ashkenazy et al., 2010</xref>).</p></sec><sec id="s4-6"><title>Sequence coevolution analysis</title><p>To estimate the coevolution within the residues in the usher TD, we used the normalized mutual information (NMI) (<xref ref-type="bibr" rid="bib38">Martin et al., 2005</xref>) over all position pairs in the multiple sequence alignment obtained as described above. NMI calculations were performed using PyCogent (<xref ref-type="bibr" rid="bib32">Knight et al., 2007</xref>; <xref ref-type="bibr" rid="bib10">Caporaso et al., 2008</xref>) and a Z-score was calculated for each residue pair based on the standard deviation from the mean NMI values. Only residue pairs that had Z-score > 4 were identified as coevolved pairs (<xref ref-type="bibr" rid="bib20">Gloor et al., 2005</xref>; <xref ref-type="bibr" rid="bib38">Martin et al., 2005</xref>).</p><p>Next, a coevolutionary RIN was built as a weighted undirected graph where each node represents a residue (as in the non-covalent RIN) and an edge connecting two nodes is the NMI score. The network was visualized and analszed in Cytoscape 2.8.2 (<xref ref-type="bibr" rid="bib55">Smoot et al., 2011</xref>) using NetworkAnalyzer plug-in (<xref ref-type="bibr" rid="bib4">Assenov et al., 2008</xref>) for calculating degrees of connectivity and RINalyzer plug-in (<xref ref-type="bibr" rid="bib15">Doncheva et al., 2011</xref>) for mapping the network on the PapC structure.</p></sec><sec id="s4-7"><title>Hybrid RIN</title><p>To store the entire information, we combined the coevolutionary RIN, the non-covalent RIN, and the conservation analysis into one network. In this hybrid network each node represents each PapC TD residue and is associated with the corresponding conservation score. Two nodes can have multiple edges, each weighted according to the information it carries (NMI score, NCI score, or ΔNCI).</p></sec><sec id="s4-8"><title>Allosteric ‘hot spots’ sub-network reconstruction</title><p>Using the hybrid network and starting from the α-helix region (residues 445–468) and β-hairpin residues (residues 230–240), we generated a sub-network of first neighbours residues based only on the NCI score higher than 0.3 in the native RIN. This sub-network was expanded by again adding only neighbouring residues connected by NCI score higher than 0.3. The procedure was repeated until no new residues could be added to the sub-network. A first set was generated by filtering the sub-network based on the evolutionary information. The filtering was done by selecting nodes with a conservation score of 9 (i.e., highly conserved) or nodes that are connected by an NMI edge (i.e., coevolved significantly). A second set was generated by filtering out the nodes that are not connected by a weakened interaction (ΔNCI > 0) in each of the mutant systems. Intersecting the identified sets resulted in a ‘hot spots’ sub-network. The ‘hot spots’ sub-network was then decomposed into communities using the edge-betweenness clustering algorithm (GLay) as implemented in clusterMaker (<xref ref-type="bibr" rid="bib19">Girvan and Newman, 2002</xref>; <xref ref-type="bibr" rid="bib41">Morris et al., 2010</xref>). Cytoscape 2.8.2 (<xref ref-type="bibr" rid="bib55">Smoot et al., 2011</xref>), RINalyzer plug-in (<xref ref-type="bibr" rid="bib15">Doncheva et al., 2011</xref>), and Chimera (<xref ref-type="bibr" rid="bib48">Pettersen et al., 2004</xref>) were used for mapping the network on the PapC structure.</p></sec><sec id="s4-9"><title>PapC substitution mutants</title><p>The PapC alanine substitution mutants were derived from plasmid pDG2 using the QuikChange Site-Directed Mutagenesis Kit (Stratagene/Agilent Technologies, Santa Clara, CA) and the primers listed in <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>. Plasmid pDG2 encodes wild-type <italic>papC</italic> with a C-terminal, thrombin-cleavable His-tag (<xref ref-type="bibr" rid="bib33">Li et al., 2004</xref>). All mutants were sequenced to verify the intended mutation.</p></sec><sec id="s4-10"><title>Expression and folding of the PapC substitution mutants in the outer membrane</title><p>Each of the PapC mutants was compared with wild-type PapC for expression levels and ability to fold into a stable β-barrel in the OM. OM isolation, analysis of usher protein levels, and the heat-modifiable mobility assay for β-barrel stability were done as previously described (<xref ref-type="bibr" rid="bib22">Henderson et al., 2011</xref>).</p></sec><sec id="s4-11"><title>Hemagglutination assay</title><p>HA assays were performed to test the ability of each of the PapC substitution mutants to assemble functional P pili on the bacterial surface. HA assays were performed by serial dilution in microtiter plates as previously described (<xref ref-type="bibr" rid="bib22">Henderson et al., 2011</xref>). HA titers were determined visually as the highest fold dilution of bacteria still able to agglutinate human red blood cells. Each assay was performed in triplicate; each mutant was analszed twice and the values were averaged.</p></sec><sec id="s4-12"><title>Top soft agar assay for antibiotic sensitivity</title><p>Bacteria were grown in LB medium supplemented with 100 µg/ml ampicillin (Amp) to an OD<sub>600</sub> of 0.6 and then induced for PapC expression with 0.1% arabinose for 1 hr. Aliquots of 0.1 ml bacteria were added to 3 ml melted soft top agar (0.75% LB agar) cooled to 45°C and supplemented with 100 µg/ml Amp and 0.1% arabinose. The bacteria and melted agar were mixed well and poured on top of 1.5% solid LB agar plates containing 100 µg/ml Amp and 0.1% arabinose. Once the top agar solidified, sterile 6-mm filter discs were placed on top and 10 μl of the following antibiotics were added: 75 mg/ml SDS, 2 mg/ml vancomycin, or 1.5 mg/ml erythromycin. The diameter of the growth inhibition zone around the antibiotics, including the filter disc, was measured after overnight growth at 37°C. Each PapC mutant strain was tested twice and the values were averaged.</p></sec><sec id="s4-13"><title>Electrophysiological analysis of selected PapC substitution mutants</title><p>PapC mutants (R237A, V327A, T331A, K427A, T437A, F438A, and S444A) were purified according to published protocols (<xref ref-type="bibr" rid="bib23">Henderson and Thanassi, 2013</xref>) and investigated by planar lipid bilayer electrophysiology. Planar bilayers were made from a preparation of L-α-phosphatidylcholine Type II-S from Sigma (also known as asolectin) according to the Montal and Mueller technique (<xref ref-type="bibr" rid="bib40">Montal and Mueller, 1972</xref>) following a published protocol (<xref ref-type="bibr" rid="bib37">Mapingire et al., 2013</xref>). Protein aliquots were diluted 1:1 in buffer T (1 M KCl, 5 mM Hepes, pH 7.2) containing either 1 or 2% <italic>N</italic>-Octyl-oligo-oxyethylene (octyl-POE, Axxora). Eight micrograms of protein from the diluted sample was added to the <italic>cis</italic> side of a planar lipid bilayer chamber containing ∼1.5 ml of buffer T. Gentle stirring was applied to promote spontaneous insertions of the protein into the bilayer. Channel activity was monitored by measuring current under voltage clamp conditions using an Axopatch 1D amplifier with a CV4B headstage or an Axopatch 200B amplifier (Axon Instruments). The current was digitized (ITC-18; Instrutech) and stored on a PC computer using the Acquire software (Bruxton). 10-min long traces were sampled at 1.25 ms intervals and filtered at 500 Hz. Both chambers contained buffer T and Ag/AgCl electrodes with pellet. The <italic>trans</italic> side of the bilayer was set as ground. Insertions were typically performed at −90 mV. Data display and analysis were done with pCLAMP software (Axon Instruments).</p></sec></sec></body><back><ack id="ack"><title>Acknowledgements</title><p>The authors are grateful to Dr David Houldershaw for computing support, and Drs Tsjerk Wassenaar and Paula Petrone for useful discussions. The authors acknowledge the use of the UCL Legion High Performance Computing facility, and associated services, in the completion of this work. The work was supported in part with MRC grant 018434 (to GW), a Wellcome Trust PhD Training Programme (086711/z/08/z) (to IF), MRC (G0600084) and BBSRC (BB/K01692X/1) (to MT) as well as National Institutes of Health grant GM062987 (to DGT).</p></ack><sec sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="conflict" id="conf1"><p>The authors declare that no competing interests exist.</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>IF, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con2"><p>TP, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con3"><p>NSH, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con4"><p>SG, Acquisition of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con5"><p>GP, Acquisition of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con6"><p>DGT, Acquisition of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con7"><p>AHD, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con8"><p>GW, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con9"><p>MT, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn></fn-group></sec><sec sec-type="supplementary-material"><title>Additional files</title><supplementary-material id="SD1-data"><object-id pub-id-type="doi">10.7554/eLife.03532.018</object-id><label>Supplementary file 1.</label><caption><p>Primers used in this study to generate PapC substitution mutations.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.03532.018">http://dx.doi.org/10.7554/eLife.03532.018</ext-link></p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife03532s001.xlsx"/></supplementary-material><sec sec-type="datasets"><title>Major dataset</title><p>The following previously published dataset was used:</p><p><related-object content-type="existing-dataset" source-id="http://www.pdb.org/pdb/explore/explore.do?structureId=2vqi" source-id-type="uri" id="dataro1"><collab collab-type="author">Remaut H</collab>, <collab collab-type="author">Tang C</collab>, <collab collab-type="author">Henderson NS</collab>, <collab collab-type="author">Pinkner JS</collab>, <collab collab-type="author">Wang T</collab>, <collab collab-type="author">Hultgren SJ</collab>, <collab collab-type="author">Thanassi DG</collab>, <collab collab-type="author">Waksman G</collab>, <collab collab-type="author">Li H</collab>, <year>2008</year><x>, </x><source>Structure of the p pilus usher (PapC) translocation pore</source><x>, </x><ext-link ext-link-type="uri" xlink:href="http://www.pdb.org/pdb/explore/explore.do?structureId=2vqi">http://www.pdb.org/pdb/explore/explore.do?structureId=2vqi</ext-link><x>, </x><comment>Publicly available at RCSB Protein Data Bank.</comment></related-object></p></sec></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Altschul</surname><given-names>SF</given-names></name><name><surname>Madden</surname><given-names>TL</given-names></name><name><surname>Schaffer</surname><given-names>AA</given-names></name><name><surname>Zhang</surname><given-names>J</given-names></name><name><surname>Zhang</surname><given-names>Z</given-names></name><name><surname>Miller</surname><given-names>W</given-names></name><name><surname>Lipman</surname><given-names>DJ</given-names></name></person-group><year>1997</year><article-title>Gapped BLAST and PSI-BLAST: a new generation of protein database search programs</article-title><source>Nucleic Acids Research</source><volume>25</volume><fpage>3389</fpage><lpage>3402</lpage><pub-id pub-id-type="doi">10.1093/nar/25.17.3389</pub-id></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Armougom</surname><given-names>F</given-names></name><name><surname>Moretti</surname><given-names>S</given-names></name><name><surname>Poirot</surname><given-names>O</given-names></name><name><surname>Audic</surname><given-names>S</given-names></name><name><surname>Dumas</surname><given-names>P</given-names></name><name><surname>Schaeli</surname><given-names>B</given-names></name><name><surname>Keduas</surname><given-names>V</given-names></name><name><surname>Notredame</surname><given-names>C</given-names></name></person-group><year>2006</year><article-title>Expresso: automatic incorporation of structural information in multiple sequence alignments using 3D-Coffee</article-title><source>Nucleic Acids Research</source><volume>34</volume><fpage>W604</fpage><lpage>W608</lpage><pub-id pub-id-type="doi">10.1093/nar/gkl092</pub-id></element-citation></ref><ref id="bib3"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ashkenazy</surname><given-names>H</given-names></name><name><surname>Erez</surname><given-names>E</given-names></name><name><surname>Martz</surname><given-names>E</given-names></name><name><surname>Pupko</surname><given-names>T</given-names></name><name><surname>Ben-Tal</surname><given-names>N</given-names></name></person-group><year>2010</year><article-title>ConSurf 2010: calculating evolutionary conservation in sequence and structure of proteins and nucleic acids</article-title><source>Nucleic Acids Research</source><volume>38</volume><fpage>W529</fpage><lpage>W533</lpage><pub-id pub-id-type="doi">10.1093/nar/gkq399</pub-id></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Assenov</surname><given-names>Y</given-names></name><name><surname>Ramirez</surname><given-names>F</given-names></name><name><surname>Schelhorn</surname><given-names>SE</given-names></name><name><surname>Lengauer</surname><given-names>T</given-names></name><name><surname>Albrecht</surname><given-names>M</given-names></name></person-group><year>2008</year><article-title>Computing topological parameters of biological networks</article-title><source>Bioinformatics</source><volume>24</volume><fpage>282</fpage><lpage>284</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btm554</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Atilgan</surname><given-names>AR</given-names></name><name><surname>Akan</surname><given-names>P</given-names></name><name><surname>Baysal</surname><given-names>C</given-names></name></person-group><year>2004</year><article-title>Small-world communication of residues and significance for protein dynamics</article-title><source>Biophysical Journal</source><volume>86</volume><fpage>85</fpage><lpage>91</lpage><pub-id pub-id-type="doi">10.1016/S0006-3495(04)74086-2</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berger</surname><given-names>O</given-names></name><name><surname>Edholm</surname><given-names>O</given-names></name><name><surname>Jahnig</surname><given-names>F</given-names></name></person-group><year>1997</year><article-title>Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature</article-title><source>Biophysical Journal</source><volume>72</volume><fpage>2002</fpage><lpage>2013</lpage><pub-id pub-id-type="doi">10.1016/S0006-3495(97)78845-3</pub-id></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bode</surname><given-names>C</given-names></name><name><surname>Kovacs</surname><given-names>IA</given-names></name><name><surname>Szalay</surname><given-names>MS</given-names></name><name><surname>Palotai</surname><given-names>R</given-names></name><name><surname>Korcsmaros</surname><given-names>T</given-names></name><name><surname>Csermely</surname><given-names>P</given-names></name></person-group><year>2007</year><article-title>Network analysis of protein dynamics</article-title><source>FEBS Letters</source><volume>581</volume><fpage>2776</fpage><lpage>2782</lpage><pub-id pub-id-type="doi">10.1016/j.febslet.2007.05.021</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bussi</surname><given-names>G</given-names></name><name><surname>Donadio</surname><given-names>D</given-names></name><name><surname>Parrinello</surname><given-names>M</given-names></name></person-group><year>2007</year><article-title>Canonical sampling through velocity rescaling</article-title><source>The Journal of Chemical Physics</source><volume>126</volume><fpage>014101</fpage><pub-id pub-id-type="doi">10.1063/1.2408420</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Capitani</surname><given-names>G</given-names></name><name><surname>Eidam</surname><given-names>O</given-names></name><name><surname>Glockshuber</surname><given-names>R</given-names></name><name><surname>Grutter</surname><given-names>MG</given-names></name></person-group><year>2006</year><article-title>Structural and functional insights into the assembly of type 1 pili from <italic>Escherichia coli</italic></article-title><source>Microbes and Infection</source><volume>8</volume><fpage>2284</fpage><lpage>2290</lpage><pub-id pub-id-type="doi">10.1016/j.micinf.2006.03.013</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Caporaso</surname><given-names>JG</given-names></name><name><surname>Smit</surname><given-names>S</given-names></name><name><surname>Easton</surname><given-names>BC</given-names></name><name><surname>Hunter</surname><given-names>L</given-names></name><name><surname>Huttley</surname><given-names>GA</given-names></name><name><surname>Knight</surname><given-names>R</given-names></name></person-group><year>2008</year><article-title>Detecting coevolution without phylogenetic trees? Tree-ignorant metrics of coevolution perform as well as tree-aware metrics</article-title><source>BMC Evolutionary Biology</source><volume>8</volume><fpage>327</fpage><pub-id pub-id-type="doi">10.1186/1471-2148-8-327</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Daily</surname><given-names>MD</given-names></name><name><surname>Upadhyaya</surname><given-names>TJ</given-names></name><name><surname>Gray</surname><given-names>JJ</given-names></name></person-group><year>2008</year><article-title>Contact rearrangements form coupled networks from local motions in allosteric proteins</article-title><source>Proteins</source><volume>71</volume><fpage>455</fpage><lpage>466</lpage><pub-id pub-id-type="doi">10.1002/prot.21800</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Darden</surname><given-names>T</given-names></name><name><surname>York</surname><given-names>D</given-names></name><name><surname>Pedersen</surname><given-names>L</given-names></name></person-group><year>1993</year><article-title>Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems</article-title><source>The Journal of Chemical Physics</source><volume>98</volume><fpage>10089</fpage><pub-id pub-id-type="doi">10.1063/1.464397</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Del Sol</surname><given-names>A</given-names></name><name><surname>Arauzo-Bravo</surname><given-names>MJ</given-names></name><name><surname>Amoros</surname><given-names>D</given-names></name><name><surname>Nussinov</surname><given-names>R</given-names></name></person-group><year>2007</year><article-title>Modular architecture of protein structures and allosteric communications: potential implications for signaling proteins and regulatory linkages</article-title><source>Genome Biology</source><volume>8</volume><fpage>R92</fpage><pub-id pub-id-type="doi">10.1186/gb-2007-8-5-r92</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Del Sol</surname><given-names>A</given-names></name><name><surname>Tsai</surname><given-names>CJ</given-names></name><name><surname>Ma</surname><given-names>B</given-names></name><name><surname>Nussinov</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>The origin of allosteric functional modulation: multiple pre-existing pathways</article-title><source>Structure</source><volume>17</volume><fpage>1042</fpage><lpage>1050</lpage><pub-id pub-id-type="doi">10.1016/j.str.2009.06.008</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Doncheva</surname><given-names>NT</given-names></name><name><surname>Klein</surname><given-names>K</given-names></name><name><surname>Domingues</surname><given-names>FS</given-names></name><name><surname>Albrecht</surname><given-names>M</given-names></name></person-group><year>2011</year><article-title>Analyzing and visualizing residue networks of protein structures</article-title><source>Trends in Biochemical Sciences</source><volume>36</volume><fpage>179</fpage><lpage>182</lpage><pub-id pub-id-type="doi">10.1016/j.tibs.2011.01.002</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Farabella</surname><given-names>I</given-names></name></person-group><year>2013</year><article-title>Computational studies of the usher outer-membrane assembly platform involved in pilus biogenesis</article-title><comment>Doctor of Philosophy, Birkbeck College, University of London</comment></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ferguson</surname><given-names>AD</given-names></name><name><surname>Amezcua</surname><given-names>CA</given-names></name><name><surname>Halabi</surname><given-names>NM</given-names></name><name><surname>Chelliah</surname><given-names>Y</given-names></name><name><surname>Rosen</surname><given-names>MK</given-names></name><name><surname>Ranganathan</surname><given-names>R</given-names></name><name><surname>Deisenhofer</surname><given-names>J</given-names></name></person-group><year>2007</year><article-title>Signal transduction pathway of TonB-dependent transporters</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>104</volume><fpage>513</fpage><lpage>518</lpage><pub-id pub-id-type="doi">10.1073/pnas.0609887104</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Geibel</surname><given-names>S</given-names></name><name><surname>Procko</surname><given-names>E</given-names></name><name><surname>Hultgren</surname><given-names>SJ</given-names></name><name><surname>Baker</surname><given-names>D</given-names></name><name><surname>Waksman</surname><given-names>G</given-names></name></person-group><year>2013</year><article-title>Structural and energetic basis of folded-protein transport by the FimD usher</article-title><source>Nature</source><volume>496</volume><fpage>243</fpage><lpage>246</lpage><pub-id pub-id-type="doi">10.1038/nature12007</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Girvan</surname><given-names>M</given-names></name><name><surname>Newman</surname><given-names>ME</given-names></name></person-group><year>2002</year><article-title>Community structure in social and biological networks</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>99</volume><fpage>7821</fpage><lpage>7826</lpage><pub-id pub-id-type="doi">10.1073/pnas.122653799</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gloor</surname><given-names>GB</given-names></name><name><surname>Martin</surname><given-names>LC</given-names></name><name><surname>Wahl</surname><given-names>LM</given-names></name><name><surname>Dunn</surname><given-names>SD</given-names></name></person-group><year>2005</year><article-title>Mutual information in protein multiple sequence alignments reveals two classes of coevolving positions</article-title><source>Biochemistry</source><volume>44</volume><fpage>7156</fpage><lpage>7165</lpage><pub-id pub-id-type="doi">10.1021/bi050293e</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Haiyan</surname><given-names>L</given-names></name><name><surname>Jihua</surname><given-names>W</given-names></name></person-group><year>2009</year><article-title>Network properties of protein structures at three different length scales</article-title><comment>Proceedings of the 2009 Second Pacific-Asia Conference on Web Mining and Web-based Application, IEEE Computer Society</comment><fpage>161</fpage><lpage>164</lpage><pub-id pub-id-type="doi">10.1109/WMWA.2009.6</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Henderson</surname><given-names>NS</given-names></name><name><surname>Ng</surname><given-names>TW</given-names></name><name><surname>Talukder</surname><given-names>I</given-names></name><name><surname>Thanassi</surname><given-names>DG</given-names></name></person-group><year>2011</year><article-title>Function of the usher N-terminus in catalysing pilus assembly</article-title><source>Molecular Microbiology</source><volume>79</volume><fpage>954</fpage><lpage>967</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2010.07505.x</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Henderson</surname><given-names>NS</given-names></name><name><surname>Thanassi</surname><given-names>DG</given-names></name></person-group><year>2013</year><article-title>Purification of the outer membrane usher protein and periplasmic chaperone-subunit complexes from the P and type 1 pilus systems</article-title><source>Methods in Molecular Biology</source><volume>966</volume><fpage>37</fpage><lpage>52</lpage><pub-id pub-id-type="doi">10.1007/978-1-62703-245-2_3</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hess</surname><given-names>B</given-names></name></person-group><year>2008</year><article-title>P-LINCS: a parallel linear constraint solver for molecular simulation</article-title><source>Journal of Chemical Theory and Computation</source><volume>4</volume><fpage>116</fpage><lpage>122</lpage><pub-id pub-id-type="doi">10.1021/ct700200b</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hess</surname><given-names>B</given-names></name><name><surname>Bekker</surname><given-names>H</given-names></name><name><surname>Berendsen</surname><given-names>HJC</given-names></name><name><surname>Fraaije</surname><given-names>JGEM</given-names></name></person-group><year>1997</year><article-title>LINCS: a linear constraint solver for molecular simulations</article-title><source>Journal of Computational Chemistry</source><volume>18</volume><fpage>1463</fpage><lpage>1472</lpage><pub-id pub-id-type="doi">10.1002/(SICI)1096-987X(199709)18:123.0.CO;2-H</pub-id></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hoover</surname><given-names>WG</given-names></name></person-group><year>1985</year><article-title>Canonical dynamics: Equilibrium phase-space distributions</article-title><source>Physical Review A</source><volume>31</volume><fpage>1695</fpage><lpage>1697</lpage><pub-id pub-id-type="doi">10.1103/PhysRevA.31.1695</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname><given-names>Y</given-names></name><name><surname>Smith</surname><given-names>BS</given-names></name><name><surname>Chen</surname><given-names>LX</given-names></name><name><surname>Baxter</surname><given-names>RH</given-names></name><name><surname>Deisenhofer</surname><given-names>J</given-names></name></person-group><year>2009</year><article-title>Insights into pilus assembly and secretion from the structure and functional characterization of usher PapC</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>106</volume><fpage>7403</fpage><lpage>7407</lpage><pub-id pub-id-type="doi">10.1073/pnas.0902789106</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jo</surname><given-names>S</given-names></name><name><surname>Kim</surname><given-names>T</given-names></name><name><surname>Im</surname><given-names>W</given-names></name></person-group><year>2007</year><article-title>Automated builder and database of protein/membrane complexes for molecular dynamics simulations</article-title><source>PLOS ONE</source><volume>2</volume><fpage>e880</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0000880</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jorgensen</surname><given-names>WL</given-names></name><name><surname>Chandrasekhar</surname><given-names>J</given-names></name><name><surname>Madura</surname><given-names>JD</given-names></name><name><surname>Impey</surname><given-names>RW</given-names></name><name><surname>Klein</surname><given-names>ML</given-names></name></person-group><year>1983</year><article-title>Comparison of simple potential functions for simulating liquid water</article-title><source>The Journal of Chemical Physics</source><volume>79</volume><fpage>926</fpage><pub-id pub-id-type="doi">10.1063/1.445869</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kaminski</surname><given-names>GA</given-names></name><name><surname>Friesner</surname><given-names>RA</given-names></name><name><surname>Tirado-Rives</surname><given-names>J</given-names></name><name><surname>Jorgensen</surname><given-names>WL</given-names></name></person-group><year>2001</year><article-title>Evaluation and reparametrization of the OPLS-AA force field for proteins via comparison with accurate quantum chemical calculations on peptides</article-title><source>The Journal of Physical Chemistry B</source><volume>105</volume><fpage>6474</fpage><lpage>6487</lpage><pub-id pub-id-type="doi">10.1021/jp003919d</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kandt</surname><given-names>C</given-names></name><name><surname>Ash</surname><given-names>WL</given-names></name><name><surname>Tieleman</surname><given-names>DP</given-names></name></person-group><year>2007</year><article-title>Setting up and running molecular dynamics simulations of membrane proteins</article-title><source>Methods</source><volume>41</volume><fpage>475</fpage><lpage>488</lpage><pub-id pub-id-type="doi">10.1016/j.ymeth.2006.08.006</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Knight</surname><given-names>R</given-names></name><name><surname>Maxwell</surname><given-names>P</given-names></name><name><surname>Birmingham</surname><given-names>A</given-names></name><name><surname>Carnes</surname><given-names>J</given-names></name><name><surname>Caporaso</surname><given-names>JG</given-names></name><name><surname>Easton</surname><given-names>BC</given-names></name><name><surname>Eaton</surname><given-names>M</given-names></name><name><surname>Hamady</surname><given-names>M</given-names></name><name><surname>Lindsay</surname><given-names>H</given-names></name><name><surname>Liu</surname><given-names>Z</given-names></name><name><surname>Lozupone</surname><given-names>C</given-names></name><name><surname>Mcdonald</surname><given-names>D</given-names></name><name><surname>Robeson</surname><given-names>M</given-names></name><name><surname>Sammut</surname><given-names>R</given-names></name><name><surname>Smit</surname><given-names>S</given-names></name><name><surname>Wakefield</surname><given-names>MJ</given-names></name><name><surname>Widmann</surname><given-names>J</given-names></name><name><surname>Wikman</surname><given-names>S</given-names></name><name><surname>Wilson</surname><given-names>S</given-names></name><name><surname>Ying</surname><given-names>H</given-names></name><name><surname>Huttley</surname><given-names>GA</given-names></name></person-group><year>2007</year><article-title>PyCogent: a toolkit for making sense from sequence</article-title><source>Genome Biology</source><volume>8</volume><fpage>R171</fpage><pub-id pub-id-type="doi">10.1186/gb-2007-8-8-r171</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>H</given-names></name><name><surname>Qian</surname><given-names>L</given-names></name><name><surname>Chen</surname><given-names>Z</given-names></name><name><surname>Thahbot</surname><given-names>D</given-names></name><name><surname>Liu</surname><given-names>G</given-names></name><name><surname>Liu</surname><given-names>T</given-names></name><name><surname>Thanassi</surname><given-names>DG</given-names></name></person-group><year>2004</year><article-title>The outer membrane usher forms a twin-pore secretion complex</article-title><source>Journal of Molecular Biology</source><volume>344</volume><fpage>1397</fpage><lpage>1407</lpage><pub-id pub-id-type="doi">10.1016/j.jmb.2004.10.008</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>Y</given-names></name><name><surname>Bahar</surname><given-names>I</given-names></name></person-group><year>2012</year><article-title>Sequence evolution correlates with structural dynamics</article-title><source>Molecular Biology and Evolution</source><volume>29</volume><fpage>2253</fpage><lpage>2263</lpage><pub-id pub-id-type="doi">10.1093/molbev/mss097</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lomize</surname><given-names>MA</given-names></name><name><surname>Lomize</surname><given-names>AL</given-names></name><name><surname>Pogozheva</surname><given-names>ID</given-names></name><name><surname>Mosberg</surname><given-names>HI</given-names></name></person-group><year>2006</year><article-title>OPM: orientations of proteins in membranes database</article-title><source>Bioinformatics</source><volume>22</volume><fpage>623</fpage><lpage>625</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btk023</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mapingire</surname><given-names>OS</given-names></name><name><surname>Henderson</surname><given-names>NS</given-names></name><name><surname>Duret</surname><given-names>G</given-names></name><name><surname>Thanassi</surname><given-names>DG</given-names></name><name><surname>Delcour</surname><given-names>AH</given-names></name></person-group><year>2009</year><article-title>Modulating effects of the plug, helix, and N- and C-terminal domains on channel properties of the PapC usher</article-title><source>The Journal of Biological Chemistry</source><volume>284</volume><fpage>36324</fpage><lpage>36333</lpage><pub-id pub-id-type="doi">10.1074/jbc.M109.055798</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mapingire</surname><given-names>OS</given-names></name><name><surname>Wager</surname><given-names>B</given-names></name><name><surname>Delcour</surname><given-names>AH</given-names></name></person-group><year>2013</year><article-title>Electrophysiological characterization of bacterial pore-forming proteins in planar lipid bilayers</article-title><source>Methods in Molecular Biology</source><volume>966</volume><fpage>381</fpage><lpage>396</lpage><pub-id pub-id-type="doi">10.1007/978-1-62703-245-2_24</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martin</surname><given-names>LC</given-names></name><name><surname>Gloor</surname><given-names>GB</given-names></name><name><surname>Dunn</surname><given-names>SD</given-names></name><name><surname>Wahl</surname><given-names>LM</given-names></name></person-group><year>2005</year><article-title>Using information theory to search for co-evolving residues in proteins</article-title><source>Bioinformatics</source><volume>21</volume><fpage>4116</fpage><lpage>4124</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/bti671</pub-id></element-citation></ref><ref id="bib39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mayrose</surname><given-names>I</given-names></name><name><surname>Graur</surname><given-names>D</given-names></name><name><surname>Ben-Tal</surname><given-names>N</given-names></name><name><surname>Pupko</surname><given-names>T</given-names></name></person-group><year>2004</year><article-title>Comparison of site-specific rate-inference methods for protein sequences: empirical Bayesian methods are superior</article-title><source>Molecular Biology and Evolution</source><volume>21</volume><fpage>1781</fpage><lpage>1791</lpage><pub-id pub-id-type="doi">10.1093/molbev/msh194</pub-id></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Montal</surname><given-names>M</given-names></name><name><surname>Mueller</surname><given-names>P</given-names></name></person-group><year>1972</year><article-title>Formation of bimolecular membranes from lipid monolayers and a study of their electrical properties</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>69</volume><fpage>3561</fpage><lpage>3566</lpage><pub-id pub-id-type="doi">10.1073/pnas.69.12.3561</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Morris</surname><given-names>JH</given-names></name><name><surname>Meng</surname><given-names>EC</given-names></name><name><surname>Ferrin</surname><given-names>TE</given-names></name></person-group><year>2010</year><article-title>Computational tools for the interactive exploration of proteomic and structural data</article-title><source>Molecular & Cellular Proteomics</source><volume>9</volume><fpage>1703</fpage><lpage>1715</lpage><pub-id pub-id-type="doi">10.1074/mcp.R000007-MCP201</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Murzyn</surname><given-names>K</given-names></name><name><surname>Rog</surname><given-names>T</given-names></name><name><surname>Pasenkiewicz-Gierula</surname><given-names>M</given-names></name></person-group><year>2005</year><article-title>Phosphatidylethanolamine-phosphatidylglycerol bilayer as a model of the inner bacterial membrane</article-title><source>Biophysical Journal</source><volume>88</volume><fpage>1091</fpage><lpage>1103</lpage><pub-id pub-id-type="doi">10.1529/biophysj.104.048835</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Naveed</surname><given-names>H</given-names></name><name><surname>Jackups</surname><given-names>R</given-names><suffix>Jnr</suffix></name><name><surname>Liang</surname><given-names>J</given-names></name></person-group><year>2009</year><article-title>Predicting weakly stable regions, oligomerization state, and protein-protein interfaces in transmembrane domains of outer membrane proteins</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>106</volume><fpage>12735</fpage><lpage>12740</lpage><pub-id pub-id-type="doi">10.1073/pnas.0902169106</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Newman</surname><given-names>ME</given-names></name><name><surname>Girvan</surname><given-names>M</given-names></name></person-group><year>2004</year><article-title>Finding and evaluating community structure in networks</article-title><source>Physical Review E, Statistical, Nonlinear, and Soft Matter Physics</source><volume>69</volume><fpage>026113</fpage><pub-id pub-id-type="doi">10.1103/PhysRevE.69.026113</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ng</surname><given-names>TW</given-names></name><name><surname>Akman</surname><given-names>L</given-names></name><name><surname>Osisami</surname><given-names>M</given-names></name><name><surname>Thanassi</surname><given-names>DG</given-names></name></person-group><year>2004</year><article-title>The usher N terminus is the initial targeting site for chaperone-subunit complexes and participates in subsequent pilus biogenesis events</article-title><source>Journal of Bacteriology</source><volume>186</volume><fpage>5321</fpage><lpage>5331</lpage><pub-id pub-id-type="doi">10.1128/JB.186.16.5321-5331.2004</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nose</surname><given-names>S</given-names></name></person-group><year>1984</year><article-title>A molecular dynamics method for simulations in the canonical ensemble</article-title><source>Molecular Physics</source><volume>52</volume><fpage>255</fpage><lpage>268</lpage><pub-id pub-id-type="doi">10.1080/00268978400101201</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Parrinello</surname><given-names>M</given-names></name></person-group><year>1981</year><article-title>Polymorphic transitions in single crystals: a new molecular dynamics method</article-title><source>Journal of Applied Physics</source><volume>52</volume><fpage>7182</fpage><pub-id pub-id-type="doi">10.1063/1.328693</pub-id></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pettersen</surname><given-names>EF</given-names></name><name><surname>Goddard</surname><given-names>TD</given-names></name><name><surname>Huang</surname><given-names>CC</given-names></name><name><surname>Couch</surname><given-names>GS</given-names></name><name><surname>Greenblatt</surname><given-names>DM</given-names></name><name><surname>Meng</surname><given-names>EC</given-names></name><name><surname>Ferrin</surname><given-names>TE</given-names></name></person-group><year>2004</year><article-title>UCSF Chimera–a visualization system for exploratory research and analysis</article-title><source>Journal of Computational Chemistry</source><volume>25</volume><fpage>1605</fpage><lpage>1612</lpage><pub-id pub-id-type="doi">10.1002/jcc.20084</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Phan</surname><given-names>G</given-names></name><name><surname>Remaut</surname><given-names>H</given-names></name><name><surname>Wang</surname><given-names>T</given-names></name><name><surname>Allen</surname><given-names>WJ</given-names></name><name><surname>Pirker</surname><given-names>KF</given-names></name><name><surname>Lebedev</surname><given-names>A</given-names></name><name><surname>Henderson</surname><given-names>NS</given-names></name><name><surname>Geibel</surname><given-names>S</given-names></name><name><surname>Volkan</surname><given-names>E</given-names></name><name><surname>Yan</surname><given-names>J</given-names></name><name><surname>Kunze</surname><given-names>MB</given-names></name><name><surname>Pinkner</surname><given-names>JS</given-names></name><name><surname>Ford</surname><given-names>B</given-names></name><name><surname>Kay</surname><given-names>CW</given-names></name><name><surname>Li</surname><given-names>H</given-names></name><name><surname>Hultgren</surname><given-names>SJ</given-names></name><name><surname>Thanassi</surname><given-names>DG</given-names></name><name><surname>Waksman</surname><given-names>G</given-names></name></person-group><year>2011</year><article-title>Crystal structure of the FimD usher bound to its cognate FimC-FimH substrate</article-title><source>Nature</source><volume>474</volume><fpage>49</fpage><lpage>53</lpage><pub-id pub-id-type="doi">10.1038/nature10109</pub-id></element-citation></ref><ref id="bib50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Remaut</surname><given-names>H</given-names></name><name><surname>Tang</surname><given-names>C</given-names></name><name><surname>Henderson</surname><given-names>NS</given-names></name><name><surname>Pinkner</surname><given-names>JS</given-names></name><name><surname>Wang</surname><given-names>T</given-names></name><name><surname>Hultgren</surname><given-names>SJ</given-names></name><name><surname>Thanassi</surname><given-names>DG</given-names></name><name><surname>Waksman</surname><given-names>G</given-names></name><name><surname>Li</surname><given-names>H</given-names></name></person-group><year>2008</year><article-title>Fiber formation across the bacterial outer membrane by the chaperone/usher pathway</article-title><source>Cell</source><volume>133</volume><fpage>640</fpage><lpage>652</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2008.03.033</pub-id></element-citation></ref><ref id="bib51"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Roberts</surname><given-names>JA</given-names></name><name><surname>Marklund</surname><given-names>BI</given-names></name><name><surname>Ilver</surname><given-names>D</given-names></name><name><surname>Haslam</surname><given-names>D</given-names></name><name><surname>Kaack</surname><given-names>MB</given-names></name><name><surname>Baskin</surname><given-names>G</given-names></name><name><surname>Louis</surname><given-names>M</given-names></name><name><surname>Mollby</surname><given-names>R</given-names></name><name><surname>Winberg</surname><given-names>J</given-names></name><name><surname>Normark</surname><given-names>S</given-names></name></person-group><year>1994</year><article-title>The Gal(alpha 1-4)Gal-specific tip adhesin of <italic>Escherichia coli</italic> P-fimbriae is needed for pyelonephritis to occur in the normal urinary tract</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>91</volume><fpage>11889</fpage><lpage>11893</lpage><pub-id pub-id-type="doi">10.1073/pnas.91.25.11889</pub-id></element-citation></ref><ref id="bib52"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sali</surname><given-names>A</given-names></name><name><surname>Blundell</surname><given-names>TL</given-names></name></person-group><year>1993</year><article-title>Comparative protein modelling by satisfaction of spatial restraints</article-title><source>Journal of Molecular Biology</source><volume>234</volume><fpage>779</fpage><lpage>815</lpage><pub-id pub-id-type="doi">10.1006/jmbi.1993.1626</pub-id></element-citation></ref><ref id="bib53"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sethi</surname><given-names>A</given-names></name><name><surname>Eargle</surname><given-names>J</given-names></name><name><surname>Black</surname><given-names>AA</given-names></name><name><surname>Luthey-Schulten</surname><given-names>Z</given-names></name></person-group><year>2009</year><article-title>Dynamical networks in tRNA:protein complexes</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>106</volume><fpage>6620</fpage><lpage>6625</lpage><pub-id pub-id-type="doi">10.1073/pnas.0810961106</pub-id></element-citation></ref><ref id="bib54"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shen</surname><given-names>MY</given-names></name><name><surname>Sali</surname><given-names>A</given-names></name></person-group><year>2006</year><article-title>Statistical potential for assessment and prediction of protein structures</article-title><source>Protein Science</source><volume>15</volume><fpage>2507</fpage><lpage>2524</lpage><pub-id pub-id-type="doi">10.1110/ps.062416606</pub-id></element-citation></ref><ref id="bib55"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smoot</surname><given-names>ME</given-names></name><name><surname>Ono</surname><given-names>K</given-names></name><name><surname>Ruscheinski</surname><given-names>J</given-names></name><name><surname>Wang</surname><given-names>PL</given-names></name><name><surname>Ideker</surname><given-names>T</given-names></name></person-group><year>2011</year><article-title>Cytoscape 2.8: new features for data integration and network visualization</article-title><source>Bioinformatics</source><volume>27</volume><fpage>431</fpage><lpage>432</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btq675</pub-id></element-citation></ref><ref id="bib56"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suel</surname><given-names>GM</given-names></name><name><surname>Lockless</surname><given-names>SW</given-names></name><name><surname>Wall</surname><given-names>MA</given-names></name><name><surname>Ranganathan</surname><given-names>R</given-names></name></person-group><year>2003</year><article-title>Evolutionarily conserved networks of residues mediate allosteric communication in proteins</article-title><source>Nature Structural Biology</source><volume>10</volume><fpage>59</fpage><lpage>69</lpage><pub-id pub-id-type="doi">10.1038/nsb881</pub-id></element-citation></ref><ref id="bib57"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tang</surname><given-names>S</given-names></name><name><surname>Liao</surname><given-names>JC</given-names></name><name><surname>Dunn</surname><given-names>AR</given-names></name><name><surname>Altman</surname><given-names>RB</given-names></name><name><surname>Spudich</surname><given-names>JA</given-names></name><name><surname>Schmidt</surname><given-names>JP</given-names></name></person-group><year>2007</year><article-title>Predicting allosteric communication in myosin via a pathway of conserved residues</article-title><source>Journal of Molecular Biology</source><volume>373</volume><fpage>1361</fpage><lpage>1373</lpage><pub-id pub-id-type="doi">10.1016/j.jmb.2007.08.059</pub-id></element-citation></ref><ref id="bib58"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Taylor</surname><given-names>NR</given-names></name></person-group><year>2013</year><article-title>Small world network strategies for studying protein structures and binding</article-title><source>Computational and Structural Biotechnology Journal</source><volume>5</volume><fpage>e201302006</fpage><pub-id pub-id-type="doi">10.5936/csbj.201302006</pub-id></element-citation></ref><ref id="bib59"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thanassi</surname><given-names>DG</given-names></name><name><surname>Saulino</surname><given-names>ET</given-names></name><name><surname>Hultgren</surname><given-names>SJ</given-names></name></person-group><year>1998</year><article-title>The chaperone/usher pathway: a major terminal branch of the general secretory pathway</article-title><source>Current Opinion in Microbiology</source><volume>1</volume><fpage>223</fpage><lpage>231</lpage><pub-id pub-id-type="doi">10.1016/S1369-5274(98)80015-5</pub-id></element-citation></ref><ref id="bib60"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thanassi</surname><given-names>DG</given-names></name><name><surname>Stathopoulos</surname><given-names>C</given-names></name><name><surname>Dodson</surname><given-names>K</given-names></name><name><surname>Geiger</surname><given-names>D</given-names></name><name><surname>Hultgren</surname><given-names>SJ</given-names></name></person-group><year>2002</year><article-title>Bacterial outer membrane ushers contain distinct targeting and assembly domains for pilus biogenesis</article-title><source>Journal of Bacteriology</source><volume>184</volume><fpage>6260</fpage><lpage>6269</lpage><pub-id pub-id-type="doi">10.1128/JB.184.22.6260-6269.2002</pub-id></element-citation></ref><ref id="bib61"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tsai</surname><given-names>CJ</given-names></name><name><surname>Del Sol</surname><given-names>A</given-names></name><name><surname>Nussinov</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>Protein allostery, signal transmission and dynamics: a classification scheme of allosteric mechanisms</article-title><source>Molecular Biosystems</source><volume>5</volume><fpage>207</fpage><lpage>216</lpage><pub-id pub-id-type="doi">10.1039/b819720b</pub-id></element-citation></ref><ref id="bib62"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Van Der Spoel</surname><given-names>D</given-names></name><name><surname>Lindahl</surname><given-names>E</given-names></name><name><surname>Hess</surname><given-names>B</given-names></name><name><surname>Groenhof</surname><given-names>G</given-names></name><name><surname>Mark</surname><given-names>AE</given-names></name><name><surname>Berendsen</surname><given-names>HJ</given-names></name></person-group><year>2005</year><article-title>GROMACS: fast, flexible, and free</article-title><source>Journal of Computational Chemistry</source><volume>26</volume><fpage>1701</fpage><lpage>1718</lpage><pub-id pub-id-type="doi">10.1002/jcc.20291</pub-id></element-citation></ref><ref id="bib63"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Volkan</surname><given-names>E</given-names></name><name><surname>Kalas</surname><given-names>V</given-names></name><name><surname>Pinkner</surname><given-names>JS</given-names></name><name><surname>Dodson</surname><given-names>KW</given-names></name><name><surname>Henderson</surname><given-names>NS</given-names></name><name><surname>Pham</surname><given-names>T</given-names></name><name><surname>Waksman</surname><given-names>G</given-names></name><name><surname>Delcour</surname><given-names>AH</given-names></name><name><surname>Thanassi</surname><given-names>DG</given-names></name><name><surname>Hultgren</surname><given-names>SJ</given-names></name></person-group><year>2013</year><article-title>Molecular basis of usher pore gating in <italic>Escherichia coli</italic> pilus biogenesis</article-title><source>Proceedings of the National Academy of Sciences of USA</source><volume>110</volume><fpage>20741</fpage><lpage>20746</lpage><pub-id pub-id-type="doi">10.1073/pnas.1320528110</pub-id></element-citation></ref><ref id="bib64"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Waksman</surname><given-names>G</given-names></name><name><surname>Hultgren</surname><given-names>SJ</given-names></name></person-group><year>2009</year><article-title>Structural biology of the chaperone-usher pathway of pilus biogenesis</article-title><source>Nature Reviews Microbiology</source><volume>7</volume><fpage>765</fpage><lpage>774</lpage><pub-id pub-id-type="doi">10.1038/nrmicro2220</pub-id></element-citation></ref></ref-list></back><sub-article article-type="article-commentary" id="SA1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.03532.019</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Dötsch</surname><given-names>Volker</given-names></name><role>Reviewing editor</role><aff><institution>Goethe University</institution>, <country>Germany</country></aff></contrib></contrib-group></front-stub><body><boxed-text><p>eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see <ext-link ext-link-type="uri" xlink:href="http://elifesciences.org/review-process">review process</ext-link>). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.</p></boxed-text><p>Thank you for sending your work entitled "Allosteric Signalling in the Outer Membrane Translocation Domain of PapC Usher" for consideration at <italic>eLife</italic>. Your article has been favorably evaluated by John Kuriyan (Senior editor) and Volker Dötsch (Reviewing editor), and two other reviewers (Nir Ben-Tal and Michael Wiener).</p><p>The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.</p><p>The manuscript describes mechanistic studies of the outer membrane translocation domain of PapC Usher, a protein required for assembly and secretion of P pili in uropathogenic <italic>E. coli</italic>, which is involved in kidney infection. This is an interesting paper that utilizes a hybrid computational approach to identify amino acid residues involved in signaling/allosteric regulation. The use of multiple computational methods to select/filter those residues most likely to be significant adds confidence to their selection, and enabled the verification experiments to be better focused.</p><p>A few points should, however, be addressed:</p><p>1) While the three steps of RIN, evolutionary analysis, and interaction analysis are presented and well-described, the relative contributions of each to the final Hot Spot sub-network are not provided. For example, if evolutionary analysis were not included, how would this have affected the identified Hot Spot sub-network? Please include some information, like a Venn diagram, to illustrate how each step selected/filtered the results. Or, what if order of application of the steps were varied?</p><p>2) For the electrophysiology assay the authors wanted to purify all 10 mutants that had shown an effect in either the hemagglutination or in the permeability assay. Only 7 were stable enough for further analysis. This is to some extent a contradiction to the earlier statement that no difference in stability of the mutants relative to wild type protein has been seen in the OM of the bacteria. Is this only a detergent effect or can it be ruled out that these mutants are also less stable in the OM and therefore the results are influenced by differences in folding/stability?</p><p>3) The goal of this project was to use a combination of in silico tools to identify residues that play key roles in the allosteric regulation of the pore formation. One particularly important aspect that was noted by all reviewers is that this combination of multiple analytical methods probably decreases "false positives", i.e. residues that one or more methods selects as being functionally important, but another method (or two) filters out. Do you have identified residues that one or more methods select and that are removed by the remaining method(s)? A panel of three figures showing the residues that "light up" with each method would be useful.</p><p>In addition, are there residues located close to the allosteric residue communities in the structure that would have been selected by visual inspection of the structure as probably functionally important but that have not scored high in your combinatorial in silico approach? If there would be even experimental data showing that such "negative control" residues exist and can be predicted as functionally not important that would strongly support your method.</p><p>4) Along the same lines: 15 mutations were made but only 7 showed an effect in the hemagglutination assay and 10 if combined with the antibiotic resistance assay. Why do the other mutations not have the expected effect?</p></body></sub-article><sub-article article-type="reply" id="SA2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.03532.020</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><p><italic>1) While the three steps of RIN, evolutionary analysis, and interaction analysis are presented and well-described, the relative contributions of each to the final Hot Spot sub-network are not provided</italic>. <italic>For example, if evolutionary analysis were not included, how would this have affected the identified Hot Spot sub-network? Please include some information, like a Venn diagram, to illustrate how each step selected/filtered the results. Or, what if order of application of the steps were varied?</italic></p><p>We thank the reviewers for raising this issue, which is clearly very important but only partly addressed in the manuscript. We now added some new details about our method to clarify some of these issues.</p><p>First, in the Results we expanded the paragraph entitled “Identifying Allosteric ‘Hot Spots’ From a Hybrid Residue Interaction Network” that refers to the filtering issues (sentences starting “This initial large sub- network is formed by 208 nodes (residues) connected by 456 NCI edges (in the native RIN)”).</p><p>Second, we would like to highlight that the order of the filters does not affect the identified Hot Spot sub-network, which rises from the intersection of the different sets of networks resulting from the application of the above filters. To clarify this point we added a supplementary figure (<xref ref-type="fig" rid="fig4s1">Figure 4–figure supplement 1A</xref>), showing a Venn diagram of each step involved in the method.</p><p>Additionally, we now updated the section entitled “Allosteric ‘hot spots’ sub-network reconstruction” to clarify that our filters are used in combination (sentences starting “Using the hybrid network and starting from the α-helix region (residues 445-468) and β-hairpin residues (residues 230-240) we generated a sub-network of first neighbours residues based only on the NCI score higher than 0.3 in the native RIN”).</p><p>Note that originally we wrongly stated in the Methods that we used a NCI score greater than 0.3 for the selection of the weakened interactions in the section above. Now this point is corrected.</p><p>Finally, we changed <xref ref-type="fig" rid="fig4">Figure 4</xref> to clarify the issue of the order by which the filters are applied (i.e. the two filters are indicated in the same level).</p><p><italic>2) For the electrophysiology assay the authors wanted to purify all 10 mutants that had shown an effect in either the hemagglutination or in the permeability assay. Only 7 were stable enough for further analysis</italic>. <italic>This is to some extent a contradiction to the earlier statement that no difference in stability of the mutants relative to wild type protein has been seen in the OM of the bacteria. Is this only a detergent effect or can it be ruled out that these mutants are also less stable in the OM and therefore the results are influenced by differences in folding/stability?</italic></p><p>We appreciate the concerns of the reviewers. The overall expression levels and folding of the 10 examined PapC mutants in the outer membrane are, indeed, similar to wild-type PapC usher. All PapC variants used in the electrophysiology assay for single channel recording were equally stable during the OM extraction procedure as wild-type. However, some PapC mutants behave well in their physiological environment within the membrane, but are unstable once extracted out of the outer membrane for purification. This is the case in 3 of the 10 purified mutants. Indeed, during the OM extraction procedure, D249A, R442A and Y260A mutants were not stable enough to maintain their native conformation in detergents, likely due to the loss of the membrane bilayer environment.</p><p>In response to the suggestions, we added a sentence to the section entitled “Electrophysiological analysis of selected mutants” in the Results:</p><p>“During the OM extraction procedure, D249A, R442A, and Y260A were not stable enough due to the loss of the membrane bilayer environment and inability to maintain their native conformation in detergents.”</p><p><italic>3) The goal of this project was to use a combination of in silico tools to identify residues that play key roles in the allosteric regulation of the pore formation. One particularly important aspect that was noted by all reviewers is that this combination of multiple analytical methods probably decreases "false positives", i.e. residues that one or more methods selects as being functionally important, but another method (or two) filters out. Do you have identified residues that one or more methods select and that are removed by the remaining method(s)? A panel of three figures showing the residues that "light up" with each method would be useful</italic>.</p><p><italic>In addition, are there residues located close to the allosteric residue communities in the structure that would have been selected by visual inspection of the structure as probably functionally important but that have not scored high in your combinatorial in silico approach? If there would be even experimental data showing that such "negative control" residues exist and can be predicted as functionally not important that would strongly support your method.</italic></p><p>We thank the reviewers for raising this interesting point.</p><p>The residues that “light up” as a result of the evolutionary information are shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>. Additionally, to highlight the residues that light up as a result of the dynamic information in each of the mutants and their combination with the native PapC we have now added new supplemental figures: <xref ref-type="fig" rid="fig2s1">Figure 2–figure supplement 1</xref> and <xref ref-type="fig" rid="fig2s2">Figure 2–figure supplement 2</xref>, showing each of the non-covalent RINs and their combination with the native PapC RIN, respectively.</p><p>As highlighted by the reviewers, the predictive power of our method rises from the combination of multiple analytical methods. These are based both on evolutionary data and on the information regarding the dynamics of the protein (without the need for prior knowledge on large confrontational changes it may undergo). As shown in the new supplementary figure (<xref ref-type="fig" rid="fig4s1">Figure 4–figure supplement 1A</xref>) only 14% of all the residues in the TD were selected as a result of our integrative method, filtering out a large portion of the residues. Thus, the main strength of our method is the ability to explore a large part of the protein. As a result, the method is able to predict not only residues that could in principle be identified by visual inspection (especially if a structure of another conformation is available) but also residues in regions that would not attract our attention at all although they may contain important information in respect to the allosteric network. For example, V327, which was picked by our method and was shown experimentally to have a functional effect on plug displacement, could potentially have been selected for mutagenesis studies by visual inspection. This residue, which is located on P-linker2 (following a proline that is at the border of the mobile part of the linker) as part of the interface with P-linker1 and the barrel wall, has a different interaction pattern in the open state of the usher (based on FimD:FimC:FimH structure (<xref ref-type="bibr" rid="bib49">Phan et al., Nature 2011</xref> Jun 2;474(7349):49-53). On the other hand, the method predicted, for example, two residues to be functionally important, T331 and S444, and it is extremely unlikely that they would have been selected as candidates for mutagenesis studies by visual inspection only (out of the 494 residues comprising the TD). T331 is located on P-linker2 at the interface between P -linker2 and the barrel wall, which stays intact in the open conformation of the usher) S444 β13 is located on base of the α-helix, i.e. on the barrel itself, and it lacks any direct contact with the plug domain. However, both residues were shown to have an effect on plug displacement.</p><p>We have added a dedicated paragraph in a new Discussion section (sentences starting “Interestingly, our community residues in some cases are found to have very different patterns of interactions in an alternative conformation”).</p><p><italic>4) Along the same lines: 15 mutations were made but only 7 showed an effect in the hemagglutination assay and 10 if combined with the antibiotic resistance assay</italic>. <italic>Why do the other mutations not have the expected effect?</italic></p><p>Our method helps to guide the experimental investigation pointing towards regions of interest in the structure (“communities”). Along the line of the answer to question #3, it is clear that the method is useful to pinpointing key candidates while exploring a large part of the investigated protein without prior knowledge of conformational changes. Additionally one should keep in mind that the hemagglutination assay and antibiotic sensitivity assay (used in this study) are screening tools, and as such, lack in sensitivity to pick up smaller changes in the ability to assemble pili or channel activity of the usher. For this reason we further investigated the key candidate mutants using a more sensitive assay such as single channel recording.</p><p>We have now added a sentence in the Results section entitled ‘Experimental Analysis of Residues In The Hot Spot Sub-Network’ to clarify this point:</p><p>“Y329A, R337A, S363A and S420A did not appear to perturb the allosteric signalling within PapC, showing the same antibiotic sensitivity phenotype and ability to assemble pili of the native PapC (<xref ref-type="table" rid="tbl4">Table 4</xref>). However, the hemagglutination assay and antibiotic sensitivity assay are screening tools, and as such, they lack in sensitivity to pick up smaller changes in the ability to assemble pili or channel activity of the usher.”</p></body></sub-article></article> |