Mutation of Tyr137 of the universal Escherichia coli fimbrial adhesin FimH relaxes the tyrosine gate prior to mannose binding

The most prevalent diseases manifested by Escherichia coli are acute and recurrent bladder infections and chronic inflammatory bowel diseases such as Crohn's disease. E. coli clinical isolates express the FimH adhesin, which consists of a mannose-specific lectin domain connected via a pilin domain to the tip of type 1 pili. Although the isolated FimH lectin domain has affinities in the nanomolar range for all high-mannosidic glycans, differentiation between these glycans is based on their capacity to form predominantly hydrophobic interactions within the tyrosine gate at the entrance to the binding pocket. In this study, novel crystal structures of tyrosine-gate mutants of FimH, ligand-free or in complex with heptyl α-d-O-mannopyranoside or 4-biphenyl α-d-O-mannopyranoside, are combined with quantum-mechanical calculations and molecular-dynamics simulations. In the Y48A FimH crystal structure, a large increase in the dynamics of the alkyl chain of heptyl α-d-O-mannopyranoside attempts to compensate for the absence of the aromatic ring; however, the highly energetic and stringent mannose-binding pocket of wild-type FimH is largely maintained. The Y137A mutation, on the other hand, is the most detrimental to FimH affinity and specificity: (i) in the absence of ligand the FimH C-terminal residue Thr158 intrudes into the mannose-binding pocket and (ii) ethylenediaminetetraacetic acid interacts strongly with Glu50, Thr53 and Asn136, in spite of multiple dialysis and purification steps. Upon mutation, pre-ligand-binding relaxation of the backbone dihedral angles at position 137 in the tyrosine gate and their coupling to Tyr48 via the interiorly located Ile52 form the basis of the loss of affinity of the FimH adhesin in the Y137A mutant.

Modeling protein-protein and protein-peptide complexes: CAPRI 6th edition

We present the sixth report evaluating the performance of methods for predicting the atomic resolution structures of protein complexes offered as targets to the community-wide initiative on the Critical Assessment of Predicted Interactions (CAPRI). The evaluation is based on a total of 20,670 predicted models for 8 protein-peptide complexes, a novel category of targets in CAPRI, and 12 protein-protein targets in CAPRI prediction Rounds held during the years 2013-2016. For two of the protein-protein targets, the focus was on the prediction of side-chain conformation and positions of interfacial water molecules. Seven of the protein-protein targets were particularly challenging owing to their multicomponent nature, to conformational changes at the binding site, or to a combination of both. Encouragingly, the very large multiprotein complex with the nucleosome was correctly predicted, and correct models were submitted for the protein-peptide targets, but not for some of the challenging protein-protein targets. Models of acceptable quality or better were obtained for 14 of the 20 targets, including medium quality models for 13 targets and high quality models for 8 targets, indicating tangible progress of present-day computational methods in modeling protein complexes with increased accuracy. Our evaluation suggests that the progress stems from better integration of different modeling tools with docking procedures, as well as the use of more sophisticated evolutionary information to score models. Nonetheless, adequate modeling of conformational flexibility in interacting proteins remains an important area with a crucial need for improvement. Proteins 2017; 85:359-377. © 2016 Wiley Periodicals, Inc.

Introducing a Clustering Step in a Consensus Approach for the Scoring of Protein-Protein Docking Models

Correctly scoring protein-protein docking models to single out native-like ones is an open challenge. It is also an object of assessment in CAPRI (Critical Assessment of PRedicted Interactions), the community-wide blind docking experiment. We introduced in the field the first pure consensus method, CONSRANK, which ranks models based on their ability to match the most conserved contacts in the ensemble they belong to. In CAPRI, scorers are asked to evaluate a set of available models and select the top ten ones, based on their own scoring approach. Scorers’ performance is ranked based on the number of targets/interfaces for which they could provide at least one correct solution. In such terms, blind testing in CAPRI Round 30 (a joint prediction round with CASP11) has shown that critical cases for CONSRANK are represented by targets showing multiple interfaces or for which only a very small number of correct solutions are available. To address these challenging cases, CONSRANK has now been modified to include a contact-based clustering of the models as a preliminary step of the scoring process. We used an agglomerative hierarchical clustering based on the number of common inter-residue contacts within the models. Two criteria, with different thresholds, were explored in the cluster generation, setting either the number of common contacts or of total clusters. For each clustering approach, after selecting the top (most populated) ten clusters, CONSRANK was run on these clusters and the top-ranked model for each cluster was selected, in the limit of 10 models per target. We have applied our modified scoring approach, Clust-CONSRANK, to SCORE_SET, a set of CAPRI scoring models made recently available by CAPRI assessors, and to the subset of homodimeric targets in CAPRI Round 30 for which CONSRANK failed to include a correct solution within the ten selected models. Results show that, for the challenging cases, the clustering step typically enriches the ten top ranked models in native-like solutions. The best performing clustering approaches we tested indeed lead to more than double the number of cases for which at least one correct solution can be included within the top ten ranked models.

Prediction of homoprotein and heteroprotein complexes by protein docking and template-based modeling: A CASP-CAPRI experiment

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein-protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. Proteins 2016; 84(Suppl 1):323-348. © 2016 Wiley Periodicals, Inc.

Molecular docking as a popular tool in drug design, an in silico travel

New molecular modeling approaches, driven by rapidly improving computational platforms, have allowed many success stories for the use of computer-assisted drug design in the discovery of new mechanism-or structure-based drugs. In this overview, we highlight three aspects of the use of molecular docking. First, we discuss the combination of molecular and quantum mechanics to investigate an unusual enzymatic mechanism of a flavoprotein. Second, we present recent advances in anti-infectious agents' synthesis driven by structural insights. At the end, we focus on larger biological complexes made by protein-protein interactions and discuss their relevance in drug design. This review provides information on how these large systems, even in the presence of the solvent, can be investigated with the outlook of drug discovery.

CLUB-MARTINI: Selecting Favourable Interactions amongst Available Candidates, a Coarse-Grained Simulation Approach to Scoring Docking Decoys

Large-scale identification of native binding orientations is crucial for understanding the role of protein-protein interactions in their biological context. Measuring binding free energy is the method of choice to estimate binding strength and reveal the relevance of particular conformations in which proteins interact. In a recent study, we successfully applied coarse-grained molecular dynamics simulations to measure binding free energy for two protein complexes with similar accuracy to full-atomistic simulation, but 500-fold less time consuming. Here, we investigate the efficacy of this approach as a scoring method to identify stable binding conformations from thousands of docking decoys produced by protein docking programs. To test our method, we first applied it to calculate binding free energies of all protein conformations in a CAPRI (Critical Assessment of PRedicted Interactions) benchmark dataset, which included over 19000 protein docking solutions for 15 benchmark targets. Based on the binding free energies, we ranked all docking solutions to select the near-native binding modes under the assumption that the native-solutions have lowest binding free energies. In our top 100 ranked structures, for the ‘easy’ targets that have many near-native conformations, we obtain a strong enrichment of acceptable or better quality structures; for the ‘hard’ targets without near-native decoys, our method is still able to retain structures which have native binding contacts. Moreover, in our top 10 selections, CLUB-MARTINI shows a comparable performance when compared with other state-of-the-art docking scoring functions. As a proof of concept, CLUB-MARTINI performs remarkably well for many targets and is able to pinpoint near-native binding modes in the top selections. To the best of our knowledge, this is the first time interaction free energy calculated from MD simulations have been used to rank docking solutions at a large scale.

Structures of C-mannosylated anti-adhesives bound to the type 1 fimbrial FimH adhesin

Selective inhibitors of the type 1 fimbrial adhesin FimH are recognized as attractive alternatives for antibiotic therapies and prophylaxes against Escherichia coli infections such as urinary-tract infections. To construct these inhibitors, the α-d-mannopyranoside of high-mannose N-glycans, recognized with exclusive specificity on glycoprotein receptors by FimH, forms the basal structure. A hydrophobic aglycon is then linked to the mannose by the O1 oxygen inherently present in the α-anomeric configuration. Substitution of this O atom by a carbon introduces a C-glycosidic bond, which may enhance the therapeutic potential of such compounds owing to the inability of enzymes to degrade C-glycosidic bonds. Here, the first crystal structures of the E. coli FimH adhesin in complex with C-glycosidically linked mannopyranosides are presented. These findings explain the role of the spacer in positioning biphenyl ligands for interactions by means of aromatic stacking in the tyrosine gate of FimH and how the normally hydrated C-glycosidic link is tolerated. As these new compounds can bind FimH, it can be assumed that they have the potential to serve as potent new antagonists of FimH, paving the way for the design of a new family of anti-adhesive compounds against urinary-tract infections.

The Antiadhesive Strategy in Crohn's Disease: Orally Active Mannosides to Decolonize Pathogenic Escherichia coli from the Gut

Blocking the adherence of bacteria to cells is an attractive complementary approach to current antibiotic treatments, which are faced with increasing resistance. This strategy has been particularly studied in the context of urinary tract infections (UTIs), in which the adhesion of pathogenic Escherichia coli strains to uroepithelial cells is prevented by blocking the FimH adhesin expressed at the tips of bacteria organelles called fimbriae. Recently, we extended the antiadhesive concept, showing that potent FimH antagonists can block the attachment of adherent-invasive E. coli (AIEC) colonizing the intestinal mucosa of patients with Crohn's disease (CD). In this work, we designed a small library of analogues of heptyl mannoside (HM), a previously identified nanomolar FimH inhibitor, but one that displays poor antiadhesive effects in vivo. The anomeric oxygen atom was replaced by a sulfur or a methylene group to prevent hydrolysis by intestinal glycosidases, and chemical groups were attached at the end of the alkyl tail. Importantly, a lead compound was shown to reduce AIEC levels in the feces and in the colonic and ileal mucosa after oral administration (10 mg kg(-1) ) in a transgenic mouse model of CD. The compound showed a low bioavailability, preferable in this instance, thus suggesting the possibility of setting up an innovative antiadhesive therapy, based on the water-soluble and non-cytotoxic FimH antagonists developed here, for the CD subpopulation in which AIEC plays a key role.

Balance between Coiled-Coil Stability and Dynamics Regulates Activity of BvgS Sensor Kinase in Bordetella

The two-component system BvgAS controls the expression of the virulence regulon of Bordetella pertussis. BvgS is a prototype of bacterial sensor kinases with extracytoplasmic Venus flytrap perception domains. Following its transmembrane segment, BvgS harbors a cytoplasmic Per-Arnt-Sim (PAS) domain and then a predicted 2-helix coiled coil that precede the dimerization-histidine-phosphotransfer domain of the kinase. BvgS homologs have a similar domain organization, or they harbor only a predicted coiled coil between the transmembrane and the dimerization-histidine-phosphotransfer domains. Here, we show that the 2-helix coiled coil of BvgS regulates the enzymatic activity in a mechanical manner. Its marginally stable hydrophobic interface enables a switch between a state of great rotational dynamics in the kinase mode and a more rigid conformation in the phosphatase mode in response to signal perception by the periplasmic domains. We further show that the activity of BvgS is controlled in the same manner if its PAS domain is replaced with the natural α-helical sequences of PAS-less homologs. Clamshell motions of the Venus flytrap domains trigger the shift of the coiled coil's dynamics. Thus, we have uncovered a general mechanism of regulation for the BvgS family of Venus flytrap-containing two-component sensor kinases.

IMPORTANCE

The two-component system BvgAS of the whooping cough agent Bordetella pertussis regulates the virulence factors necessary for infection in a coordinated manner. BvgS is the prototype of a family of sensor kinase proteins found in major bacterial pathogens. When BvgS functions as a kinase, B. pertussis is virulent, and the bacterium shifts to an avirulent phase after BvgS senses chemicals that make it switch to phosphatase. Our goal is to decipher the signaling mechanisms of BvgS in order to understand virulence regulation in Bordetella, which may lead to new antimicrobial treatments targeting those two-component systems. We discovered that the activity of BvgS is regulated in a mechanical manner. A short region of the protein that precedes the enzymatic domain switches between two states in response to signal perception by other BvgS domains. This switch region is conserved among BvgS homologs, and thus, the regulation uncovered here will likely be relevant for the family.

Signal Transduction by BvgS Sensor Kinase: BINDING OF MODULATOR NICOTINATE AFFECTS THE CONFORMATION AND DYNAMICS OF THE ENTIRE PERIPLASMIC MOIETY

The two-component sensory transduction system BvgAS controls the virulence regulon of the whooping-cough agent Bordetella pertussis. The periplasmic moiety of the homodimeric sensor kinase BvgS is composed of four bilobed Venus flytrap (VFT) perception domains followed by α helices that extend into the cytoplasmic membrane. In the virulent phase, the default state of B. pertussis, the cytoplasmic enzymatic moiety of BvgS acts as kinase by autophosphorylating and transferring the phosphoryl group to the response regulator BvgA. Under laboratory conditions, BvgS shifts to phosphatase activity in response to modulators, notably nicotinate ions. Here we characterized the effects of nicotinate and related modulators on the BvgS periplasmic moiety by using site-directed mutagenesis and in silico and biophysical approaches. Modulators bind with low affinity to BvgS in the VFT2 cavity. Electron paramagnetic resonance shows that their binding globally affects the conformation and dynamics of the periplasmic moiety. Specific amino acid substitutions designed to slacken interactions within and between the VFT lobes prevent BvgS from responding to nicotinate, showing that BvgS shifts from kinase to phosphatase activity in response to this modulator via a tense transition state that involves a large periplasmic structural block. We propose that this transition enables the transmembrane helices to adopt a distinct conformation that sets the cytoplasmic enzymatic moiety in the phosphatase mode. The bona fide, in vivo VFT ligands that remain to be identified are likely to trigger similar effects on the transmembrane and cytoplasmic moieties. This mechanism may be relevant to the other VFT-containing sensor kinases homologous to BvgS.

Virulence regulation with Venus flytrap domains: structure and function of the periplasmic moiety of the sensor-kinase BvgS

Two-component systems (TCS) represent major signal-transduction pathways for adaptation to environmental conditions, and regulate many aspects of bacterial physiology. In the whooping cough agent Bordetella pertussis, the TCS BvgAS controls the virulence regulon, and is therefore critical for pathogenicity. BvgS is a prototypical TCS sensor-kinase with tandem periplasmic Venus flytrap (VFT) domains. VFT are bi-lobed domains that typically close around specific ligands using clamshell motions. We report the X-ray structure of the periplasmic moiety of BvgS, an intricate homodimer with a novel architecture. By combining site-directed mutagenesis, functional analyses and molecular modeling, we show that the conformation of the periplasmic moiety determines the state of BvgS activity. The intertwined structure of the periplasmic portion and the different conformation and dynamics of its mobile, membrane-distal VFT1 domains, and closed, membrane-proximal VFT2 domains, exert a conformational strain onto the transmembrane helices, which sets the cytoplasmic moiety in a kinase-on state by default corresponding to the virulent phase of the bacterium. Signaling the presence of negative signals perceived by the periplasmic domains implies a shift of BvgS to a distinct state of conformation and activity, corresponding to the avirulent phase. The response to negative modulation depends on the integrity of the periplasmic dimer, indicating that the shift to the kinase-off state implies a concerted conformational transition. This work lays the bases to understand virulence regulation in Bordetella. As homologous sensor-kinases control virulence features of diverse bacterial pathogens, the BvgS structure and mechanism may pave the way for new modes of targeted therapeutic interventions.

Regulatory motifs on ISWI chromatin remodelers: molecular mechanisms and kinetic proofreading

Recently, kinetic proofreading scenarios have been proposed for the regulation of chromatin remodeling, first on purely theoretical grounds (Blossey and Schiessel 2008 HFSP J. 2 167-70) and deduced from experiments on the ISWI/ACF system (Narlikar 2010 Curr. Opin. Chem. Biol. 14 660). In the kinetic proofreading scenario of chromatin remodeling, the combination of the recognition of a histone tail state and ATP-hydrolysis in the remodeler motor act together to select (i.e. proofread) a nucleosomal substrate. ISWI remodelers have recently been shown to have an additional level of regulation as they contain auto-inhibitory motifs which need to be inactivated through an interaction with the nucleosome. In this paper we show that the auto-regulatory effect enhances substrate recognition in kinetic proofreading. We further report some suggestive additional insights into the molecular mechanism underlying ISWI-autoregulation.

Blind prediction of interfacial water positions in CAPRI

We report the first assessment of blind predictions of water positions at protein-protein interfaces, performed as part of the critical assessment of predicted interactions (CAPRI) community-wide experiment. Groups submitting docking predictions for the complex of the DNase domain of colicin E2 and Im2 immunity protein (CAPRI Target 47), were invited to predict the positions of interfacial water molecules using the method of their choice. The predictions-20 groups submitted a total of 195 models-were assessed by measuring the recall fraction of water-mediated protein contacts. Of the 176 high- or medium-quality docking models-a very good docking performance per se-only 44% had a recall fraction above 0.3, and a mere 6% above 0.5. The actual water positions were in general predicted to an accuracy level no better than 1.5 Å, and even in good models about half of the contacts represented false positives. This notwithstanding, three hotspot interface water positions were quite well predicted, and so was one of the water positions that is believed to stabilize the loop that confers specificity in these complexes. Overall the best interface water predictions was achieved by groups that also produced high-quality docking models, indicating that accurate modelling of the protein portion is a determinant factor. The use of established molecular mechanics force fields, coupled to sampling and optimization procedures also seemed to confer an advantage. Insights gained from this analysis should help improve the prediction of protein-water interactions and their role in stabilizing protein complexes.

On the molecular basis of D-bifunctional protein deficiency type III

Molecular basis of D-bifunctional protein (D-BP) deficiency was studied with wild type and five disease-causing variants of 3R-hydroxyacyl-CoA dehydrogenase fragment of the human MFE-2 (multifunctional enzyme type 2) protein. Complementation analysis in vivo in yeast and in vitro enzyme kinetic and stability determinants as well as in silico stability and structural fluctuation calculations were correlated with clinical data of known patients. Despite variations not affecting the catalytic residues, enzyme kinetic performance (K(m), V(max) and k(cat)) of the recombinant protein variants were compromised to a varying extent and this can be judged as the direct molecular cause for D-BP deficiency. Protein stability plays an additional role in producing non-functionality of MFE-2 in case structural variations affect cofactor or substrate binding sites. Structure-function considerations of the variant proteins matched well with the available data of the patients.

Oligoarginine vectors for intracellular delivery: role of arginine side-chain orientation in chain length-dependent destabilization of lipid membranes

Arginine-rich peptides receive increased attention due to their capacity to cross different types of membranes and to transport cargo molecules inside cells. Even though peptide-induced destabilization has been investigated extensively, little is known about the peptide side-chain and backbone orientation with respect to the bilayer that may contribute to a molecular understanding of the peptide-induced membrane perturbations. The main objective of this work is to provide a detailed description of the orientation of arginine peptides in the lipid bilayer of PC and negatively charged PG liposomes using ATR-IR spectroscopy and molecular modeling, and to relate these orientational preferences to lipid bilayer destabilization. Molecular modeling showed that above the transition temperature arginine side-chains are preferentially solvent-directed at the PC/water interface whereas several arginine side-chains are pointing towards the PG hydrophobic core. IR dichroic spectra confirmed the orientation of the arginine side chains perpendicular to the lipid-water interface. IR spectra shows an randomly distributed backbone that seems essential to optimize interactions with the lipid membrane. The observed increase of permeation to a fluorescent dye is related to the peptide induced-formation of gauche bonds in the acyl chains. In the absence of hydrophobic residues, insertion of side-chains that favors phosphate/guanidium interaction is another mechanism of membrane permeabilization that has not been further analyzed so far.

Docking and scoring protein interactions: CAPRI 2009

Protein docking algorithms are assessed by evaluating blind predictions performed during 2007-2009 in Rounds 13-19 of the community-wide experiment on critical assessment of predicted interactions (CAPRI). We evaluated the ability of these algorithms to sample docking poses and to single out specific association modes in 14 targets, representing 11 distinct protein complexes. These complexes play important biological roles in RNA maturation, G-protein signal processing, and enzyme inhibition and function. One target involved protein-RNA interactions not previously considered in CAPRI, several others were hetero-oligomers, or featured multiple interfaces between the same protein pair. For most targets, predictions started from the experimentally determined structures of the free (unbound) components, or from models built from known structures of related or similar proteins. To succeed they therefore needed to account for conformational changes and model inaccuracies. In total, 64 groups and 12 web-servers submitted docking predictions of which 4420 were evaluated. Overall our assessment reveals that 67% of the groups, more than ever before, produced acceptable models or better for at least one target, with many groups submitting multiple high- and medium-accuracy models for two to six targets. Forty-one groups including four web-servers participated in the scoring experiment with 1296 evaluated models. Scoring predictions also show signs of progress evidenced from the large proportion of correct models submitted. But singling out the best models remains a challenge, which also adversely affects the ability to correctly rank docking models. With the increased interest in translating abstract protein interaction networks into realistic models of protein assemblies, the growing CAPRI community is actively developing more efficient and reliable docking and scoring methods for everyone to use.

Identification of specific lipid-binding sites in integral membrane proteins

Protein-lipid interactions are increasingly recognized as central to the structure and function of membrane proteins. However, with the exception of simplified models, specific protein-lipid interactions are particularly difficult to highlight experimentally. Here, we used molecular dynamics simulations to identify a specific protein-lipid interaction in lactose permease, a prototypical model for transmembrane proteins. The interactions can be correlated with the functional dependence of the protein to specific lipid species. The technique is simple and widely applicable to other membrane proteins, and a variety of lipid matrices can be used.