Evaluating protein cross-linking as a therapeutic strategy to stabilize SOD1 variants in a mouse model of familial ALS.

Mutations in the gene encoding Cu-Zn superoxide dismutase 1 (SOD1) cause a subset of familial amyotrophic lateral sclerosis (fALS) cases. A shared effect of these mutations is that SOD1, which is normally a stable dimer, dissociates into toxic monomers that seed toxic aggregates. Considerable research effort has been devoted to developing compounds that stabilize the dimer of fALS SOD1 variants, but unfortunately, this has not yet resulted in a treatment. We hypothesized that cyclic thiosulfinate cross-linkers, which selectively target a rare, 2 cysteine-containing motif, can stabilize fALS-causing SOD1 variants in vivo. We created a library of chemically diverse cyclic thiosulfinates and determined structure-cross-linking-activity relationships. A pre-lead compound, "S-XL6," was selected based upon its cross-linking rate and drug-like properties. Co-crystallographic structure clearly establishes the binding of S-XL6 at Cys 111 bridging the monomers and stabilizing the SOD1 dimer. Biophysical studies reveal that the degree of stabilization afforded by S-XL6 (up to 24°C) is unprecedented for fALS, and to our knowledge, for any protein target of any kinetic stabilizer. Gene silencing and protein degrading therapeutic approaches require careful dose titration to balance the benefit of diminished fALS SOD1 expression with the toxic loss-of-enzymatic function. We show that S-XL6 does not share this liability because it rescues the activity of fALS SOD1 variants. No pharmacological agent has been proven to bind to SOD1 in vivo. Here, using a fALS mouse model, we demonstrate oral bioavailability; rapid engagement of SOD1G93A by S-XL6 that increases SOD1G93A's in vivo half-life; and that S-XL6 crosses the blood-brain barrier. S-XL6 demonstrated a degree of selectivity by avoiding off-target binding to plasma proteins. Taken together, our results indicate that cyclic thiosulfinate-mediated SOD1 stabilization should receive further attention as a potential therapeutic approach for fALS.

Mechanisms of isoform-specific residue influence on GTP-bound HRas, KRas, and NRas.

HRas, KRas, and NRas are GTPases with a common set of effectors that control many cell-signaling pathways, including proliferation through Raf kinase. Their G-domains are nearly identical in sequence, with a few isoform-specific residues that have an effect on dynamics and biochemical properties. Here, we use accelerated molecular dynamics (aMD) simulations consistent with solution x-ray scattering experiments to elucidate mechanisms through which isoform-specific residues associated with each Ras isoform affects functionally important regions connected to the active site. HRas-specific residues cluster in loop 8 to stabilize the nucleotide-binding pocket, while NRas-specific residues on helix 3 directly affect the conformations of switch I and switch II. KRas, the most globally flexible of the isoforms, shows greatest fluctuations in the switch regions enhanced by a KRas-specific residue in loop 7 and a highly dynamic loop 8 region. The analysis of isoform-specific residue effects on Ras proteins is supported by NMR experiments and is consistent with previously published biochemical data.

A Simplified Adult Expiratory Support Device: First Clinical Cases in a Low-Resource Setting.

PURPOSE: In low and middle-income countries, mechanical ventilators or commercially available devices used to offer continuous positive airway pressure are not readily affordable and available. In Ghana, nearly 10% of critically ill patients presenting to the emergency department require ventilator support. DESCRIPTION: We designed, built, and tested a simple expiratory positive airway pressure (EPAP) device to provide adult respiratory support in low resource environments with or without supplemental oxygen and without the need for electricity. EVALUATION: Laboratory tests demonstrated that the device is capable of delivering EPAP at levels expected to provide significant assistance to some patients. We present the first 2 cases where the use of this simple EPAP device provided critical respiratory support during weaning of patients from mechanical ventilation. CONCLUSIONS: A low-cost 3-dimensional printable adult respiratory support device could provide substantial benefit to patients suffering from respiratory distress through the delivery of appropriate levels of EPAP in a low-resource setting with limited infrastructure. Further clinical validation is needed for broader application in low-resource settings.

Raf promotes dimerization of the Ras G-domain with increased allosteric connections.

Ras dimerization is critical for Raf activation. Here we show that the Ras binding domain of Raf (Raf-RBD) induces robust Ras dimerization at low surface densities on supported lipid bilayers and, to a lesser extent, in solution as observed by size exclusion chromatography and confirmed by SAXS. Community network analysis based on molecular dynamics simulations shows robust allosteric connections linking the two Raf-RBD D113 residues located in the Galectin scaffold protein binding site of each Raf-RBD molecule and 85 Å apart on opposite ends of the dimer complex. Our results suggest that Raf-RBD binding and Ras dimerization are concerted events that lead to a high-affinity signaling complex at the membrane that we propose is an essential unit in the macromolecular assembly of higher order Ras/Raf/Galectin complexes important for signaling through the Ras/Raf/MEK/ERK pathway.

Biological and Clinical Consequences of Integrin Binding via a Rogue RGD Motif in the SARS CoV-2 Spike Protein.

Although ACE2 (angiotensin converting enzyme 2) is considered the primary receptor for CoV-2 cell entry, recent reports suggest that alternative pathways may contribute. This paper considers the hypothesis that viral binding to cell-surface integrins may contribute to the high infectivity and widespread extra-pulmonary impacts of the SARS-CoV-2 virus. This potential is suggested on the basis of the emergence of an RGD (arginine-glycine-aspartate) sequence in the receptor-binding domain of the spike protein. RGD is a motif commonly used by viruses to bind cell-surface integrins. Numerous signaling pathways are mediated by integrins and virion binding could lead to dysregulation of these pathways, with consequent tissue damage. Integrins on the surfaces of pneumocytes, endothelial cells and platelets may be vulnerable to CoV-2 virion binding. For instance, binding of intact virions to integrins on alveolar cells could enhance viral entry. Binding of virions to integrins on endothelial cells could activate angiogenic cell signaling pathways; dysregulate integrin-mediated signaling pathways controlling developmental processes; and precipitate endothelial activation to initiate blood clotting. Such a procoagulant state, perhaps together with enhancement of platelet aggregation through virions binding to integrins on platelets, could amplify the production of microthrombi that pose the threat of pulmonary thrombosis and embolism, strokes and other thrombotic consequences. The susceptibility of different tissues to virion-integrin interactions may be modulated by a host of factors, including the conformation of relevant integrins and the impact of the tissue microenvironment on spike protein conformation. Patient-specific differences in these factors may contribute to the high variability of clinical presentation. There is danger that the emergence of receptor-binding domain mutations that increase infectivity may also enhance access of the RGD motif for integrin binding, resulting in viral strains with ACE2 independent routes of cell entry and novel integrin-mediated biological and clinical impacts. The highly infectious variant, B.1.1.7 (or VUI 202012/01), includes a receptor-binding domain amino acid replacement, N501Y, that could potentially provide the RGD motif with enhanced access to cell-surface integrins, with consequent clinical impacts.

Glycine substitution in SH3-SH2 connector of Hck tyrosine kinase causes population shift from assembled to disassembled state.

A combination of small angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations based on a coarse grained model is used to examine the effect of glycine substitutions in the short connector between the SH3 and SH2 domains of Hck, a member of the Src-family kinases. It has been shown previously that the activity of cSrc kinase is upregulated by substitution of 3 residues by glycine in the SH3-SH2 connector. Here, analysis of SAXS data indicates that the population of Hck in the disassembled state increases from 25% in the wild type kinase to 76% in the glycine mutant. This is consistent with the results of free energy perturbation calculations showing that the mutation in the connector shifts the equilibrium from the assembled to the disassembled state. This study supports the notion that the SH3-SH2 connector helps to regulate the activity of tyrosine kinases by shifting the population of the active state of the multidomain protein independent of C-terminal phosphorylation.

FiXR: a framework to reconstruct fiber cross-sections from X-ray fiber diffraction experiments.

Ab initio reconstruction methods have revolutionized the capabilities of small-angle X-ray scattering (SAXS), allowing the data-driven discovery of previously unknown molecular conformations, exploiting optimization heuristics and assumptions behind the composition of globular molecules. While these methods have been successful for the analysis of small particles, their impact on fibrillar assemblies has been more limited. The micrometre-range size of these assemblies and the complex interaction of their periodicities in their scattering profiles indicate that the discovery of fibril structures from SAXS measurements requires novel approaches beyond extending existing tools for molecular discovery. In this work, it is proposed to use SAXS measurements, together with diffraction theory, to infer the electron distribution of the average cross-section of a fiber. This cross-section is modeled as a discrete electron density with continuous support, allowing representations beyond binary distributions. Additional constraints, such as non-negativity or smoothness/connectedness, can also be added to the framework. The proposed approach is tested using simulated SAXS data from amyloid ß fibril models and using measured data of Tobacco mosaic virus from SAXS experiments, recovering the geometry and density of the cross-sections in all cases. The approach is further tested by analyzing SAXS data from different amyloid ß fibril assemblies, with results that are in agreement with previously proposed models from cryo-EM measurements. The limitations of the proposed method, together with an analysis of the robustness of the method and the combination with different experimental sources, are also discussed.

The Structural Basis of Amyloid Strains in Alzheimer's Disease.

Amyloid fibrils represent one of the defining features of Alzheimer's disease (AD). They are made up of protofilaments composed of amyloid ß (Aß) peptides that are held together with extraordinary stability by a network of tight steric zippers and axial hydrogen bonds. This review explores the hypothesis that the peptide conformation within a protofilament represents the physical embodiment of a "strain" of AD. Evidence suggests that within a single strain the fold of individual peptides is invariant. However, the fibrils are capable of structural polymorphism that includes variation in the arrangement of protofilaments into fibrils, the pitch of the resultant fibrils, and the higher-order organization of the plaques into which they aggregate. These intrastrain polymorphisms are separated by low energy barriers, allowing multiple configurations to coexist within a single preparation or tissue. Clinical presentation of different strains may be determined by variation in the way different protofilament structures generate the relevant toxic species, be they monomers, oligomers, or higher-order structures. Evidence reviewed here is consistent with a model in which disease progression is concomitant with a gradual, progressive annealing of amyloid fibrils from benign, loosely packed structures into dense neurotoxic aggregates. This model challenges the commonly held hypothesis that oligomers of Aß peptides are the only active proximate species in neurodegeneration. However, the data do not implicate fibrils themselves. Rather, they cast suspicion on larger-scale supramolecular aggregates as toxic agents. Electron tomography of amyloid plaques in situ strongly suggests that the formation of amyloid aggregates results in perturbation of the cellular membrane integrity, warranting further investigation of this as a potential mode of neurotoxicity. If dense supramolecular amyloid aggregates prove to be important agents of neurodegeneration in AD, this model may also have relevance to other forms of amyloidoses.

The Structural Basis of the Farnesylated and Methylated KRas4B Interaction with Calmodulin.

Ca2+-calmodulin (CaM) extracts KRas4B from the plasma membrane, suggesting that KRas4B/CaM interaction plays a role in regulating Ras signaling. To gain mechanistic insight, we provide a computational model, supported by experimental structural data, of farnesylated/methylated KRas4B1-185 interacting with CaM in solution and at anionic membranes including signaling lipids. Due to multiple interaction modes, we observe diverse conformational ensembles of the KRas4B-CaM complex. A highly populated conformation reveals the catalytic domain interacting with the N-lobe and the hypervariable region (HVR) wrapping around the linker with the farnesyl docking to the extended CaM's C-lobe pocket. Alternatively, KRas4B can interact with collapsed CaM with the farnesyl penetrating CaM's center. At anionic membranes, CaM interacts with the catalytic domain with large fluctuations, drawing the HVR. Signaling lipids establishing strong salt bridges with CaM prevent membrane departure. Membrane-interacting KRas4B-CaM complex can productively recruit phosphatidylinositol 3-kinase α (PI3Kα) to the plasma membrane, serving as a coagent in activating PI3Kα/Akt signaling.

Effects of Catalytic Action and Ligand Binding on Conformational Ensembles of Adenylate Kinase.

Crystal structures of adenylate kinase (AdK) from Escherichia coli capture two states: an "open" conformation (apo) obtained in the absence of ligands and a "closed" conformation in which ligands are bound. Other AdK crystal structures suggest intermediate conformations that may lie on the transition pathway between these two states. To characterize the transition from open to closed states in solution, X-ray solution scattering data were collected from AdK in the apo form and with progressively increasing concentrations of five different ligands. Scattering data from apo AdK are consistent with scattering predicted from the crystal structure of AdK in the open conformation. In contrast, data from AdK samples saturated with Ap5A do not agree with that calculated from AdK in the closed conformation. Using cluster analysis of available structures, we selected representative structures in five conformational states: open, partially open, intermediate, partially closed, and closed. We used these structures to estimate the relative abundances of these states for each experimental condition. X-ray solution scattering data obtained from AdK with AMP are dominated by scattering from AdK in the open conformation. For AdK in the presence of high concentrations of ATP and ADP, the conformational ensemble shifts to a mixture of partially open and closed states. Even when AdK is saturated with Ap5A, a significant proportion of AdK remains in a partially open conformation. These results are consistent with an induced-fit model in which the transition of AdK from an open state to a closed state is initiated by ATP binding.

Calculation of the cross-sectional shape of a fibril from equatorial scattering.

An alternate formulation of helical diffraction theory is used to generate cross-sectional shapes of fibrous structures from equatorial scattering. We demonstrate this approach with computationally generated scattering intensities and then apply it to scattering data from Tobacco Mosaic Virus (TMV) and in vitro assembled fibrils of Aß40 peptides. Refining the cross-sectional shape of TMV from SAXS data collected on a 26mg/ml solution resulted in a circular shape with outer diameter of ∼180Å and inner diameter of ∼40Å consistent with the known structure of TMV. We also utilized this method to analyze the equatorial scattering from TMV collected by Don Caspar from a concentrated (24% ∼295mg/ml) gel of TMV as reported in his Ph.D. thesis in 1955. This data differs from the SAXS data in having a sharp interference peak at ∼250Å spacing, indicative of strong interparticle interactions in the gel. Analysis of this data required consideration of interatomic vectors as long as 2000Å and resulted in generation of images that were interpreted as representative of local organization of TMV particles in the sample. Peaks in the images were separated, on average by about 250Å with a density consistent with Caspar's original measurements. Analysis of SAXS data from Aß fibrils resulted in a cross-sectional shape that could be interpreted in terms of structural models that have been constructed from ssNMR and cryoEM. These results demonstrate an unexpected use of the small-angle region of fiber diffraction patterns to derive fundamental structural properties of scattering objects.

Amyloid structure exhibits polymorphism on multiple length scales in human brain tissue.

Aggregation of Aß amyloid fibrils into plaques in the brain is a universal hallmark of Alzheimer's Disease (AD), but whether plaques in different individuals are equivalent is unknown. One possibility is that amyloid fibrils exhibit different structures and different structures may contribute differentially to disease, either within an individual brain or between individuals. However, the occurrence and distribution of structural polymorphisms of amyloid in human brain is poorly documented. Here we use X-ray microdiffraction of histological sections of human tissue to map the abundance, orientation and structural heterogeneities of amyloid. Our observations indicate that (i) tissue derived from subjects with different clinical histories may contain different ensembles of fibrillar structures; (ii) plaques harboring distinct amyloid structures can coexist within a single tissue section and (iii) within individual plaques there is a gradient of fibrillar structure from core to margins. These observations have immediate implications for existing theories on the inception and progression of AD.

c-Abl Tyrosine Kinase Adopts Multiple Active Conformational States in Solution.

Protein tyrosine kinases of the Abl family have diverse roles in normal cellular regulation and drive several forms of leukemia as oncogenic fusion proteins. In the crystal structure of the inactive c-Abl kinase core, the SH2 and SH3 domains dock onto the back of the kinase domain, resulting in a compact, assembled state. This inactive conformation is stabilized by the interaction of the myristoylated N-cap with a pocket in the C-lobe of the kinase domain. Mutations that perturb these intramolecular interactions result in kinase activation. Here, we present X-ray scattering solution structures of multidomain c-Abl kinase core proteins modeling diverse active states. Surprisingly, the relative positions of the regulatory N-cap, SH3, and SH2 domains in an active myristic acid binding pocket mutant (A356N) were virtually identical to those of the assembled wild-type kinase core, indicating that Abl kinase activation does not require dramatic reorganization of the downregulated core structure. In contrast, the positions of the SH2 and SH3 domains in a clinically relevant imatinib-resistant gatekeeper mutant (T315I) appear to be reconfigured relative to their positions in the wild-type protein. Our results demonstrate that c-Abl kinase activation can occur either with (T315I) or without (A356N) global allosteric changes in the core, revealing the potential for previously unrecognized signaling diversity.

Cell wall targeted in planta iron accumulation enhances biomass conversion and seed iron concentration in Arabidopsis and rice.

Conversion of nongrain biomass into liquid fuel is a sustainable approach to energy demands as global population increases. Previously, we showed that iron can act as a catalyst to enhance the degradation of lignocellulosic biomass for biofuel production. However, direct addition of iron catalysts to biomass pretreatment is diffusion-limited, would increase the cost and complexity of biorefinery unit operations and may have deleterious environmental impacts. Here, we show a new strategy for in planta accumulation of iron throughout the volume of the cell wall where iron acts as a catalyst in the deconstruction of lignocellulosic biomass. We engineered CBM-IBP fusion polypeptides composed of a carbohydrate-binding module family 11 (CBM11) and an iron-binding peptide (IBP) for secretion into Arabidopsis and rice cell walls. CBM-IBP transformed Arabidopsis and rice plants show significant increases in iron accumulation and biomass conversion compared to respective controls. Further, CBM-IBP rice shows a 35% increase in seed iron concentration and a 40% increase in seed yield in greenhouse experiments. CBM-IBP rice potentially could be used to address iron deficiency, the most common and widespread nutritional disorder according to the World Health Organization.

Large-Scale Interlaboratory Study to Develop, Analytically Validate and Apply Highly Multiplexed, Quantitative Peptide Assays to Measure Cancer-Relevant Proteins in Plasma.

There is an increasing need in biology and clinical medicine to robustly and reliably measure tens to hundreds of peptides and proteins in clinical and biological samples with high sensitivity, specificity, reproducibility, and repeatability. Previously, we demonstrated that LC-MRM-MS with isotope dilution has suitable performance for quantitative measurements of small numbers of relatively abundant proteins in human plasma and that the resulting assays can be transferred across laboratories while maintaining high reproducibility and quantitative precision. Here, we significantly extend that earlier work, demonstrating that 11 laboratories using 14 LC-MS systems can develop, determine analytical figures of merit, and apply highly multiplexed MRM-MS assays targeting 125 peptides derived from 27 cancer-relevant proteins and seven control proteins to precisely and reproducibly measure the analytes in human plasma. To ensure consistent generation of high quality data, we incorporated a system suitability protocol (SSP) into our experimental design. The SSP enabled real-time monitoring of LC-MRM-MS performance during assay development and implementation, facilitating early detection and correction of chromatographic and instrumental problems. Low to subnanogram/ml sensitivity for proteins in plasma was achieved by one-step immunoaffinity depletion of 14 abundant plasma proteins prior to analysis. Median intra- and interlaboratory reproducibility was <20%, sufficient for most biological studies and candidate protein biomarker verification. Digestion recovery of peptides was assessed and quantitative accuracy improved using heavy-isotope-labeled versions of the proteins as internal standards. Using the highly multiplexed assay, participating laboratories were able to precisely and reproducibly determine the levels of a series of analytes in blinded samples used to simulate an interlaboratory clinical study of patient samples. Our study further establishes that LC-MRM-MS using stable isotope dilution, with appropriate attention to analytical validation and appropriate quality control measures, enables sensitive, specific, reproducible, and quantitative measurements of proteins and peptides in complex biological matrices such as plasma.

Multidomain assembled states of Hck tyrosine kinase in solution.

An approach combining small-angle X-ray solution scattering (SAXS) data with coarse-grained (CG) simulations is developed to characterize the assembly states of Hck, a member of the Src-family kinases, under various conditions in solution. First, a basis set comprising a small number of assembly states is generated from extensive CG simulations. Second, a theoretical SAXS profile for each state in the basis set is computed by using the Fast-SAXS method. Finally, the relative population of the different assembly states is determined via a Bayesian-based Monte Carlo procedure seeking to optimize the theoretical scattering profiles against experimental SAXS data. The study establishes the concept of basis-set supported SAXS (BSS-SAXS) reconstruction combining computational and experimental techniques. Here, BSS-SAXS reconstruction is used to reveal the structural organization of Hck in solution and the different shifts in the equilibrium population of assembly states upon the binding of different signaling peptides.

Trinucleotide cassettes increase diversity of T7 phage-displayed peptide library.

BACKGROUND: Amino acid sequence diversity is introduced into a phage-displayed peptide library by randomizing library oligonucleotide DNA. We recently evaluated the diversity of peptide libraries displayed on T7 lytic phage and M13 filamentous phage and showed that T7 phage can display a more diverse amino acid sequence repertoire due to differing processes of viral morphogenesis. METHODS: In this study, we evaluated and compared the diversity of a 12-mer T7 phage-displayed peptide library randomized using codon-corrected trinucleotide cassettes with a T7 and an M13 12-mer phage-displayed peptide library constructed using the degenerate codon randomization method. RESULTS: We herein demonstrate that the combination of trinucleotide cassette amino acid codon randomization and T7 phage display construction methods resulted in a significant enhancement to the functional diversity of a 12-mer peptide library. This novel library exhibited superior amino acid uniformity and order-of-magnitude increases in amino acid sequence diversity as compared to degenerate codon randomized peptide libraries. Comparative analyses of the biophysical characteristics of the 12-mer peptide libraries revealed the trinucleotide cassette-randomized library to be a unique resource. CONCLUSION: The combination of T7 phage display and trinucleotide cassette randomization resulted in a novel resource for the potential isolation of binding peptides for new and previously studied molecular targets.

T7 lytic phage-displayed peptide libraries exhibit less sequence bias than M13 filamentous phage-displayed peptide libraries.

We investigated whether the T7 system of phage display could produce peptide libraries of greater diversity than the M13 system of phage display due to the differing processes of lytic and filamentous phage morphogenesis. Using a bioinformatics-assisted computational approach, collections of random peptide sequences obtained from a T7 12-mer library (X(12)) and a T7 7-mer disulfide-constrained library (CX(7)C) were analyzed and compared with peptide populations obtained from New England BioLabs' M13 Ph.D.-12 and Ph.D.-C7C libraries. Based on this analysis, peptide libraries constructed with the T7 system have fewer amino acid biases, increased peptide diversity, and more normal distributions of peptide net charge and hydropathy than the M13 libraries. The greater diversity of T7-displayed libraries provides a potential resource of novel binding peptides for new as well as previously studied molecular targets. To demonstrate their utility, several of the T7-displayed peptide libraries were screened for streptavidin- and neutravidin-binding phage. Novel binding motifs were identified for each protein.

Subtractive transcriptomics: establishing polarity drives in vitro human endothelial morphogenesis.

Although investigations of mature normal and tumor-derived capillaries have resulted in characterization of these structures at the phenotypic level, less is known regarding the initial molecular cues for cellular assembly of endothelial cells into human capillaries. Here, we employ a novel combination of microenvironmental manipulation and microarray data filtration over narrowly delineated temporal data series to identify the morphogenesis component apart from the proliferation component, as pooled human microvascular-derived endothelial cells are induced to form capillary-like structures in vitro in a murine tumor-derived matrix. The 217 morphogenesis-specific genes identified using this subtractive transcriptomics approach are mostly independent of the angiogenic proteins currently used as therapeutic targets for aberrant angiogenesis. Quantitative real-time PCR was used to validate 20% of these transcripts. Immunofluorescent analysis of proliferating and tube-forming cells validates at the protein level the morphogenesis-specific expression pattern of 16 of the 217 gene products identified. The transcripts that are selectively up-regulated in tube-forming endothelial cells reveal a temporal expression pattern of genes primarily associated with intracellular trafficking, guided migration, cytoskeletal reorganization, cellular adhesion, and proliferation inhibition. These data show that a sequential up-regulation of genes that establish and maintain polarity occurs during migration and morphogenesis of in vitro human endothelial cells undergoing tubulogenesis; some of which may well be effective as novel antiangiogenic drug targets.

Phage display reveals multiple contact sites between FhuA, an outer membrane receptor of Escherichia coli, and TonB.

The ferric hydroxamate uptake receptor FhuA from Escherichia coli transports siderophores across the outer membrane (OM). TonB-ExbB-ExbD transduces energy from the cytoplasmic membrane to the OM by contacts between TonB and OM receptors that contain the Ton box, a consensus sequence near the N terminus. Although the Ton box is a region of known contact between OM receptors and TonB, our biophysical studies established that TonB binds to FhuA through multiple regions of interaction. Panning of phage-displayed random peptide libraries (Ph.D.-12, Ph.D.-C7C) against TonB identified peptide sequences that specifically interact with TonB. Analyses of these sequences using the Receptor Ligand Contacts (RELIC) suite of programs revealed clusters of multiply aligned peptides that mapped to FhuA. These clusters localized to a continuous periplasm-accessible surface: Ton box/switch helix; cork domain/beta1 strand; and periplasmic turn 8. Guided by such matches, synthetic oligonucleotides corresponding to DNA sequences identical to fhuA were fused to malE; peptides corresponding to the above regions were displayed at the N terminus of E.coli maltose-binding protein (MBP). Purified FhuA peptides fused to MBP bound specifically to TonB by ELISA. Furthermore, they competed with ligand-loaded FhuA for binding to TonB. RELIC also identified clusters of multiply aligned peptides corresponding to the Ton box regions in BtuB, FepA, and FecA; to periplasmic turn 8 in BtuB and FecA; and to periplasmic turns 1 and 2 in FepA. These experimental outcomes identify specific molecular contacts made between TonB and OM receptors that extend beyond the well-characterized Ton box.