Folding and self-assembly of the TatA translocation pore based on a charge zipper mechanism.

We propose a concept for the folding and self-assembly of the pore-forming TatA complex from the Twin-arginine translocase and of other membrane proteins based on electrostatic "charge zippers." Each subunit of TatA consists of a transmembrane segment, an amphiphilic helix (APH), and a C-terminal densely charged region (DCR). The sequence of charges in the DCR is complementary to the charge pattern on the APH, suggesting that the protein can be "zipped up" by a ladder of seven salt bridges. The length of the resulting hairpin matches the lipid bilayer thickness, hence a transmembrane pore could self-assemble via intra- and intermolecular salt bridges. The steric feasibility was rationalized by molecular dynamics simulations, and experimental evidence was obtained by monitoring the monomer-oligomer equilibrium of specific charge mutants. Similar "charge zippers" are proposed for other membrane-associated proteins, e.g., the biofilm-inducing peptide TisB, the human antimicrobial peptide dermcidin, and the pestiviral E(RNS) protein.

Tripartite efflux pumps of the RND superfamily: what did we learn from computational studies?

Bacterial resistance to antibiotics has been long recognized as a priority to address for human health. Among all micro-organisms, the so-called multi-drug resistant (MDR) bacteria, which are resistant to most, if not all drugs in our current arsenal, are particularly worrisome. The World Health Organization has prioritized the ESKAPE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa and Enterobacter species) pathogens, which include four Gram-negative bacterial species. In these bacteria, active extrusion of antimicrobial compounds out of the cell by means of 'molecular guns' known as efflux pumps is a main determinant of MDR phenotypes. The resistance-nodulation-cell division (RND) superfamily of efflux pumps connecting the inner and outer membrane in Gram-negative bacteria is crucial to the onset of MDR and virulence, as well as biofilm formation. Thus, understanding the molecular basis of the interaction of antibiotics and inhibitors with these pumps is key to the design of more effective therapeutics. With the aim to contribute to this challenge, and complement and inspire experimental research, in silico studies on RND efflux pumps have flourished in recent decades. Here, we review a selection of such investigations addressing the main determinants behind the polyspecificity of these pumps, the mechanisms of substrate recognition, transport and inhibition, as well as the relevance of their assembly for proper functioning, and the role of protein-lipid interactions. The journey will end with a perspective on the role of computer simulations in addressing the challenges posed by these beautifully complex machineries and in supporting the fight against the spread of MDR bacteria.

Interaction of Radiopharmaceuticals with Somatostatin Receptor 2 Revealed by Molecular Dynamics Simulations.

The development of drugs targeting somatostatin receptor 2 (SSTR2), generally overexpressed in neuroendocrine tumors, is focus of intense research. A few molecules in conjugation with radionuclides are in clinical use for both diagnostic and therapeutic purposes. These radiopharmaceuticals are composed of a somatostatin analogue biovector conjugated to a chelator moiety bearing the radionuclide. To date, despite valuable efforts, a detailed molecular-level description of the interaction of radiopharmaceuticals in complex with SSTR2 has not yet been accomplished. Therefore, in this work, we carefully analyzed the key dynamical features and detailed molecular interactions of SSTR2 in complex with six radiopharmaceutical compounds selected among the few already in use (64Cu/68Ga-DOTATATE, 68Ga-DOTATOC, 64Cu-SARTATE) and some in clinical development (68Ga-DOTANOC, 64Cu-TETATATE). Through molecular dynamics simulations and exploiting recently available structures of SSTR2, we explored the influence of the different portions of the compounds (peptide, radionuclide, and chelator) in the interaction with the receptor. We identified the most stable binding modes and found distinct interaction patterns characterizing the six compounds. We thus unveiled detailed molecular interactions crucial for the recognition of this class of radiopharmaceuticals. The microscopically well-founded analysis presented in this study provides guidelines for the design of new potent ligands targeting SSTR2.

Recognition of quinolone antibiotics by the multidrug efflux transporter MexB of Pseudomonas aeruginosa.

The drug/proton antiporter MexB is the engine of the major efflux pump MexAB-OprM in Pseudomonas aeruginosa. This protein is known to transport a large variety of compounds, including antibiotics, thus conferring a multi-drug resistance phenotype. Due to the difficulty of producing co-crystals, only two X-ray structures of MexB in a complex with ligands are available to date, and mechanistic aspects are largely hypothesized based on the body of data collected for the homologous protein AcrB of Escherichia coli. In particular, a recent study (Ornik-Cha, Wilhelm, Kobylka et al., Nat. Commun., 2021, 12, 6919) reported a co-crystal structure of AcrB in a complex with levofloxacin, an antibiotic belonging to the important class of (fluoro)-quinolones. In this work, we performed a systematic ensemble docking campaign coupled to the cluster analysis and molecular-mechanics optimization of docking poses to study the interaction between 36 quinolone antibiotics and MexB. We additionally investigated surface complementarity between each molecule and the transporter and thoroughly assessed the computational protocol adopted against the known experimental data. Our study reveals different binding preferences of the investigated compounds towards the sub-sites of the large deep binding pocket of MexB, supporting the hypothesis that MexB substrates oscillate between different binding modes with similar affinity. Interestingly, small changes in the molecular structure translate into significant differences in MexB-quinolone interactions. All the predicted binding modes are available for download and visualization at the following link: https://www.dsf.unica.it/dock/mexb/quinolones.

Molecular insights into the Patched1 drug efflux inhibitory activity of panicein A hydroquinone: a computational study.

Human Hedgehog receptor Patched1 (PTCH1) is able to efflux chemotherapeutics of different chemical structure out of cancer cells thus contributing to multidrug resistance phenomena in tumor treatment. A screening of natural compounds purified from marine sponges led to the identification of the first PTCH1 efflux inhibitor, panicein A hydroquinone (PAH), demonstrated to increase doxorubicin toxicity in vitro and vemurafenib toxicity in vitro and in vivo. In this work we combined different computational techniques to gain molecular insights of the inhibitory activity of PAH and some of its active and inactive analogues. We first performed a thorough characterization and druggability analysis of the main putative substrate binding pockets known from available cryo-electron microscopy structures. Further, dynamical descriptors of the active and inactive PAH analogues were extracted from microsecond-long all-atom molecular dynamics simulations in water solution. Finally, a blind ensemble docking methodology coupled with the conformational analysis of compounds enabled rationalization of the interaction between PTCH1 and PAH and derivatives in terms of their intrinsic physico-chemical properties. Our results suggest that the Neck pocket is the preferential binding site for PAH analogues on PTCH1, and that compounds assuming an open cylindric-like shape in solution are most likely to be good binders for PTCH1.

Molecular Interactions of Carbapenem Antibiotics with the Multidrug Efflux Transporter AcrB of Escherichia coli.

The drug/proton antiporter AcrB, engine of the major efflux pump AcrAB(Z)-TolC of Escherichia coli and other bacteria, is characterized by its impressive ability to transport chemically diverse compounds, conferring a multi-drug resistance (MDR) phenotype. Although hundreds of small molecules are known to be AcrB substrates, only a few co-crystal structures are available to date. Computational methods have been therefore intensively employed to provide structural and dynamical fingerprints related to transport and inhibition of AcrB. In this work, we performed a systematic computational investigation to study the interaction between representative carbapenem antibiotics and AcrB. We focused on the interaction of carbapenems with the so-called distal pocket, a region known for its importance in binding inhibitors and substrates of AcrB. Our findings reveal how the different physico-chemical nature of these antibiotics is reflected on their binding preference for AcrB. The molecular-level information provided here could help design new antibiotics less susceptible to the efflux mechanism.

Molecular basis for inhibition of AcrB multidrug efflux pump by novel and powerful pyranopyridine derivatives.

The Escherichia coli AcrAB-TolC efflux pump is the archetype of the resistance nodulation cell division (RND) exporters from Gram-negative bacteria. Overexpression of RND-type efflux pumps is a major factor in multidrug resistance (MDR), which makes these pumps important antibacterial drug discovery targets. We have recently developed novel pyranopyridine-based inhibitors of AcrB, which are orders of magnitude more powerful than the previously known inhibitors. However, further development of such inhibitors has been hindered by the lack of structural information for rational drug design. Although only the soluble, periplasmic part of AcrB binds and exports the ligands, the presence of the membrane-embedded domain in AcrB and its polyspecific binding behavior have made cocrystallization with drugs challenging. To overcome this obstacle, we have engineered and produced a soluble version of AcrB [AcrB periplasmic domain (AcrBper)], which is highly congruent in structure with the periplasmic part of the full-length protein, and is capable of binding substrates and potent inhibitors. Here, we describe the molecular basis for pyranopyridine-based inhibition of AcrB using a combination of cellular, X-ray crystallographic, and molecular dynamics (MD) simulations studies. The pyranopyridines bind within a phenylalanine-rich cage that branches from the deep binding pocket of AcrB, where they form extensive hydrophobic interactions. Moreover, the increasing potency of improved inhibitors correlates with the formation of a delicate protein- and water-mediated hydrogen bond network. These detailed insights provide a molecular platform for the development of novel combinational therapies using efflux pump inhibitors for combating multidrug resistant Gram-negative pathogens.

Water-mediated interactions enable smooth substrate transport in a bacterial efflux pump.

BACKGROUND: Efflux pumps of the Resistance-Nodulation-cell Division superfamily confer multi-drug resistance to Gram-negative bacteria. The most-studied polyspecific transporter belonging to this class is the inner-membrane trimeric antiporter AcrB of Escherichia coli. In previous studies, a functional rotation mechanism was proposed for its functioning, according to which the three monomers undergo concerted conformational changes facilitating the extrusion of substrates. However, the molecular determinants and the energetics of this mechanism still remain unknown, so its feasibility must be proven mechanistically. METHODS: A computational protocol able to mimic the functional rotation mechanism in AcrB was developed. By using multi-bias molecular dynamics simulations we characterized the translocation of the substrate doxorubicin driven by conformational changes of the protein. In addition, we estimated for the first time the free energy profile associated to this process. RESULTS: We provided a molecular view of the process in agreement with experimental data. Moreover, we showed that the conformational changes occurring in AcrB enable the formation of a layer of structured waters on the internal surface of the transport channel. This water layer, in turn, allows for a fairly constant hydration of the substrate, facilitating its diffusion over a smooth free energy profile. CONCLUSIONS: Our findings reveal a new molecular mechanism of polyspecific transport whereby water contributes by screening potentially strong substrate-protein interactions. GENERAL SIGNIFICANCE: We provided a mechanistic understanding of a fundamental process related to multi-drug transport. Our results can help rationalizing the behavior of other polyspecific transporters and designing compounds avoiding extrusion or inhibitors of efflux pumps.

AcrB drug-binding pocket substitution confers clinically relevant resistance and altered substrate specificity.

The incidence of multidrug-resistant bacterial infections is increasing globally and the need to understand the underlying mechanisms is paramount to discover new therapeutics. The efflux pumps of Gram-negative bacteria have a broad substrate range and transport antibiotics out of the bacterium, conferring intrinsic multidrug resistance (MDR). The genomes of pre- and posttherapy MDR clinical isolates of Salmonella Typhimurium from a patient that failed antibacterial therapy and died were sequenced. In the posttherapy isolate we identified a novel G288D substitution in AcrB, the resistance-nodulation division transporter in the AcrAB-TolC tripartite MDR efflux pump system. Computational structural analysis suggested that G288D in AcrB heavily affects the structure, dynamics, and hydration properties of the distal binding pocket altering specificity for antibacterial drugs. Consistent with this hypothesis, recreation of the mutation in standard Escherichia coli and Salmonella strains showed that G288D AcrB altered substrate specificity, conferring decreased susceptibility to the fluoroquinolone antibiotic ciprofloxacin by increased efflux. At the same time, the substitution increased susceptibility to other drugs by decreased efflux. Information about drug transport is vital for the discovery of new antibacterials; the finding that one amino acid change can cause resistance to some drugs, while conferring increased susceptibility to others, could provide a basis for new drug development and treatment strategies.

Analysis of fast channel blockage: revealing substrate binding in the microsecond range.

For an antibiotic to be effective, it needs to cross the outer membrane barrier and reach the target inside the cell. Hydrophilic antibiotics, e.g.ß-lactams, use porin channels to cross the outer membrane and accumulate in the periplasm. Experimental determination of antibiotic interactions with porin is performed by using electrophysiology on a single channel level by noise analysis or single event analysis methods. We report a novel framework for analyzing the ion-current noise, taking into account the corrections due to the analogous filter and the sampling procedure, with the goal of extending the time resolution to a range previously inaccessible by event analysis or by conventional noise analysis. The new method allows one to analyse fast binding events and/or the case when the single channel is not completely blocked by the substrate. We demonstrate the power of this approach by using as an example the interactions of meropenem, an antibiotic of the carbapenem family, with the OmpF porin that is considered to be one of the main pathways for antibiotics to enter Escherichia coli. The presence of meropenem in OmpF is detected by ion current blockages, and the on and off rates are estimated from the concentration dependence of the average ion current and of its power spectral density. The obtained average residence time of the antibiotic inside the channel is in the range of a few microseconds, i.e. more than 50 times smaller than the inverse cut-off frequency of the analogous filter.

A Database of Force-Field Parameters, Dynamics, and Properties of Antimicrobial Compounds.

We present an on-line database of all-atom force-field parameters and molecular properties of compounds with antimicrobial activity (mostly antibiotics and some beta-lactamase inhibitors). For each compound, we provide the General Amber Force Field parameters for the major species at physiological pH, together with an analysis of properties of interest as extracted from µs-long molecular dynamics simulations in explicit water solution. The properties include number and population of structural clusters, molecular flexibility, hydrophobic and hydrophilic molecular surfaces, the statistics of intraand inter-molecular H-bonds, as well as structural and dynamical properties of solvent molecules within first and second solvation shells. In addition, the database contains several key molecular parameters, such as energy of the frontier molecular orbitals, vibrational properties, rotational constants, atomic partial charges and electric dipole moment, computed by Density Functional Theory. The present database (to our knowledge the first extensive one including dynamical properties) is part of a wider project aiming to build-up a database containing structural, physico-chemical and dynamical properties of medicinal compounds using different force-field parameters with increasing level of complexity and reliability. The database is freely accessible at http://www.dsf.unica.it/translocation/db/.

Molecular mechanism of viral resistance to a potent non-nucleoside inhibitor unveiled by molecular simulations.

Recently, we reported on a potent benzimidazole derivative (227G) that inhibits the growth of the bovine viral diarrhea virus (BVDV) in cell-based and enzyme assays at nanomolar concentrations. The target of 227G is the viral RNA-dependent RNA polymerase (RdRp), and the I261M mutation located in motif I of the RdRp finger domain was found to induce drug resistance. Here we propose a molecular mechanism for the retained functionality of the enzyme in the presence of the inhibitor, on the basis of a thorough computational study of the apo and holo forms of the BVDV RdRp either in the wild type (wt) or in the form carrying the I261M mutation. Our study shows that although the mutation affects to some extent the structure of the apoenzyme, the functional dynamics of the protein appear to be largely maintained, which is consistent with the retained functionality of this natural mutant. Despite the binding site of 227G not collapsing or undergoing drastic structural changes upon introduction of the I261M substitution, these alterations reflect crucially on the binding mode of 227G, which is significantly different from that found in wt RdRp. In particular, while in the wt system the four loops lining the template entrance site embrace 227G and close the template passageway, in the I261M variant the template entrance is only marginally occluded, allowing in principle the translocation of the template to the interior of the enzyme. In addition, the mutated enzyme in the presence of 227G retains several characteristics of the wt apoprotein. Our work provides an original molecular picture of a resistance mechanism that is consistent with published experimental data.

Molecular mechanism of MBX2319 inhibition of Escherichia coli AcrB multidrug efflux pump and comparison with other inhibitors.

Efflux pumps of the resistance nodulation division (RND) superfamily, such as AcrB, make a major contribution to multidrug resistance in Gram-negative bacteria. The development of inhibitors of the RND pumps would improve the efficacy of current and next-generation antibiotics. To date, however, only one inhibitor has been cocrystallized with AcrB. Thus, in silico structure-based analysis is essential for elucidating the interaction between other inhibitors and the efflux pumps. In this work, we used computer docking and molecular dynamics simulations to study the interaction between AcrB and the compound MBX2319, a novel pyranopyridine efflux pump inhibitor with potent activity against RND efflux pumps of Enterobacteriaceae species, as well as other known inhibitors (D13-9001, 1-[1-naphthylmethyl]-piperazine, and phenylalanylarginine-ß-naphthylamide) and the binding of doxorubicin to the efflux-defective F610A variant of AcrB. We also analyzed the binding of a substrate, minocycline, for comparison. Our results show that MBX2319 binds very tightly to the lower part of the distal pocket in the B protomer of AcrB, strongly interacting with the phenylalanines lining the hydrophobic trap, where the hydrophobic portion of D13-9001 was found to bind by X-ray crystallography. Additionally, MBX2319 binds to AcrB in a manner that is similar to the way in which doxorubicin binds to the F610A variant of AcrB. In contrast, 1-(1-naphthylmethyl)-piperazine and phenylalanylarginine-ß-naphthylamide appear to bind to somewhat different areas of the distal pocket in the B protomer of AcrB than does MBX2319. However, all inhibitors (except D13-9001) appear to distort the structure of the distal pocket, impairing the proper binding of substrates.

Impact of corpus callosum integrity on functional interhemispheric connectivity and cognition in healthy subjects.

To examine the corpus callosum's (CC) integrity in terms of fractional anisotropy (FA) and how it affects resting-state hemispheric connectivity (rs-IHC) and cognitive function in healthy individuals. Sixty-eight healthy individuals were recruited for the study. The global FA (gFA) and FA values of each CC tract (forceps minor, body, tapetum, and forceps major) were evaluated using diffusion-weighted imaging (DWI) sequences. The homotopic functional connectivity technique was used to quantify the effects of FA in the CC tracts on bilateral functional connectivity, including the confounding effect of gFA. Brain regions with higher or lower rs-IHC were identified using the threshold-free cluster enhancement family-wise error-corrected p-value of 0.05. The null hypothesis was rejected if the p-value was ≤ 0.05 for the nonparametric partial correlation technique. Several clusters of increased rs-IHC were identified in relation to the FA of individual CC tracts, each with a unique topographic distribution and extension. Only forceps minor FA values correlated with cognitive scores. The integrity of CC influences rs-IHC differently in healthy subjects. Specifically, forceps minor anisotropy impacts rs-IHC and cognition more than other CC tracts do.

Predicting binding events in very flexible, allosteric, multi-domain proteins.

Knowledge of the structures formed by proteins and small ligands is of fundamental importance for understanding molecular principles of chemotherapy and for designing new and more effective drugs. Due to the still high costs and to the several limitations of experimental techniques, it is most often desirable to predict these ligand-protein complexes in silico, particularly when screening for new putative drugs from databases of millions of compounds. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology aimed at generating bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites. Validation was performed on the enzyme adenylate kinase (ADK), a paradigmatic example of proteins that undergo very large conformational changes upon ligand binding. By only exploiting the unbound structure and the putative binding sites of the protein, we generated a significant fraction of bound-like structures, which employed in ensemble-docking calculations allowed to find native-like poses of substrates, inhibitors, and catalytically incompetent binders. Our protocol provides a general framework for the generation of bound-like conformations of flexible proteins that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening for difficult targets, prediction of the impact of amino acid mutations on structure and dynamics, and protein engineering.

Predicting permeation of compounds across the outer membrane of P. aeruginosa using molecular descriptors.

The ability Gram-negative pathogens have at adapting and protecting themselves against antibiotics has increasingly become a public health threat. Data-driven models identifying molecular properties that correlate with outer membrane (OM) permeation and growth inhibition while avoiding efflux could guide the discovery of novel classes of antibiotics. Here we evaluate 174 molecular descriptors in 1260 antimicrobial compounds and study their correlations with antibacterial activity in Gram-negative Pseudomonas aeruginosa. The descriptors are derived from traditional approaches quantifying the compounds' intrinsic physicochemical properties, together with, bacterium-specific from ensemble docking of compounds targeting specific MexB binding pockets, and all-atom molecular dynamics simulations in different subregions of the OM model. Using these descriptors and the measured inhibitory concentrations, we design a statistical protocol to identify predictors of OM permeation/inhibition. We find consistent rules across most of our data highlighting the role of the interaction between the compounds and the OM. An implementation of the rules uncovered in our study is shown, and it demonstrates the accuracy of our approach in a set of previously unseen compounds. Our analysis sheds new light on the key properties drug candidates need to effectively permeate/inhibit P. aeruginosa, and opens the gate to similar data-driven studies in other Gram-negative pathogens.

Different molecular mechanisms of inhibition of bovine viral diarrhea virus and hepatitis C virus RNA-dependent RNA polymerases by a novel benzimidazole.

The virus-encoded RNA-dependent RNA polymerase (RdRp) has emerged as a primary target in the search for selective inhibitors of Flaviviridae. Recently, we reported on the selective inhibition, in cell-based assays, of both BVDV (EC50 = 0.80 ± 0.06 µM) and HCV (EC50 = 1.11 ± 0.15 µM) by 2-{1-[2-(2,4-dimethoxyphenyl)-1H-benzimidazol-5-yl]ethylidene}hydrazinecarbothioamide (227G). Here we show that, in enzyme assays with recombinant enzymes, 227G inhibits, in a dose-dependent manner, the RdRp of both BVDV (IC50 = 0.0020 ± 0.0004 µM) and HCV (IC50 = 0.40 ± 0.04 µM). Furthermore, we report on the selection and molecular analysis of a BVDV-resistant mutant, characterized by the presence of the I261M mutation. By applying a multilevel computational approach, we identified different 227G binding sites on the two RdRps. They were further validated by the good agreement between the calculated affinities and those extrapolated from IC50 values. Our findings suggest different molecular mechanisms of inhibition of the HCV and BVDV RdRps by 227G and indicate the importance of understanding ligand-enzyme interactions at the molecular level for the rational design of new and more potent leads.

Molecular simulations of SSTR2 dynamics and interaction with ligands.

The cyclic peptide hormone somatostatin regulates physiological processes involved in growth and metabolism, through its binding to G-protein coupled somatostatin receptors. The isoform 2 (SSTR2) is of particular relevance for the therapy of neuroendocrine tumours for which different analogues to somatostatin are currently in clinical use. We present an extensive and systematic computational study on the dynamics of SSTR2 in three different states: active agonist-bound, inactive antagonist-bound and apo inactive. We exploited the recent burst of SSTR2 experimental structures to perform µs-long multi-copy molecular dynamics simulations to sample conformational changes of the receptor and rationalize its binding to different ligands (the agonists somatostatin and octreotide, and the antagonist CYN154806). Our findings suggest that the apo form is more flexible compared to the holo ones, and confirm that the extracellular loop 2 closes upon the agonist octreotide but not upon the antagonist CYN154806. Based on interaction fingerprint analyses and free energy calculations, we found that all peptides similarly interact with residues buried into the binding pocket. Conversely, specific patterns of interactions are found with residues located in the external portion of the pocket, at the basis of the extracellular loops, particularly distinguishing the agonists from the antagonist. This study will help in the design of new somatostatin-based compounds for theranostics of neuroendocrine tumours.

Molecular determinants of avoidance and inhibition of Pseudomonas aeruginosa MexB efflux pump.

Transporters of the resistance-nodulation-cell division (RND) superfamily of proteins are the dominant multidrug efflux power of Gram-negative bacteria. The major RND efflux pump of Pseudomonas aeruginosa is MexAB-OprM, in which the inner membrane transporter MexB is responsible for the recognition and binding of compounds. The high importance of this pump in clinical antibiotic resistance made it a subject of intense investigations and a promising target for the discovery of efflux pump inhibitors. This study is focused on a series of peptidomimetic compounds developed as effective inhibitors of MexAB-OprM. We performed multi-copy molecular dynamics simulations, machine-learning (ML) analyses, and site-directed mutagenesis of MexB to investigate interactions of MexB with representatives of efflux avoiders, substrates, and inhibitors. The analysis of both direct and water-mediated protein-ligand interactions revealed characteristic patterns for each class, highlighting significant differences between them. We found that efflux avoiders poorly interact with the access binding site of MexB, and inhibition engages amino acid residues that are not directly involved in binding and transport of substrates. In agreement, machine-learning models selected different residues predictive of MexB substrates and inhibitors. The differences in interactions were further validated by site-directed mutagenesis. We conclude that the substrate translocation and inhibition pathways of MexB split at the interface (between the main putative binding sites) and at the deep binding pocket and that interactions outside of the hydrophobic patch contribute to the inhibition of MexB. This molecular-level information could help in the rational design of new inhibitors and antibiotics less susceptible to the efflux mechanism. IMPORTANCE Multidrug transporters recognize and expel from cells a broad range of ligands including their own inhibitors. The difference between the substrate translocation and inhibition routes remains unclear. In this study, machine learning and computational and experimental approaches were used to understand dynamics of MexB interactions with its ligands. Our results show that some ligands engage a certain combination of polar and charged residues in MexB binding sites to be effectively expelled into the exit funnel, whereas others engage aromatic and hydrophobic residues that slow down or hinder the next step in the transporter cycle. These findings suggest that all MexB ligands fit into this substrate-inhibitor spectrum depending on their physico-chemical structures and properties.

A novel dendrimeric peptide with antimicrobial properties: structure-function analysis of SB056.

The novel antimicrobial peptide with a dimeric dendrimer scaffold, SB056, was empirically optimized by high-throughput screening. This procedure produced an intriguing primary sequence whose structure-function analysis is described here. The alternating pattern of hydrophilic and hydrophobic amino acids suggests the possibility that SB056 is a membrane-active peptide that forms amphiphilic ß-strands in a lipid environment. Circular dichroism confirmed that the cationic SB056 folds as ß-sheets in the presence of anionic vesicles. Lipid monolayer surface pressure experiments revealed unusual kinetics of monolayer penetration, which suggest lipid-induced aggregation as a membranolytic mechanism. NMR analyses of the linear monomer and the dendrimeric SB056 in water and in 30% trifluoroethanol, on the other hand, yielded essentially unstructured conformations, supporting the excellent solubility and storage properties of this compound. However, simulated annealing showed that many residues lie in the ß-region of the Ramachandran plot, and molecular-dynamics simulations confirmed the propensity of this peptide to fold as a ß-type conformation. The excellent solubility in water and the lipid-induced oligomerization characteristics of this peptide thus shed light on its mechanism of antimicrobial action, which may also be relevant for systems that can form toxic ß-amyloid fibrils when in contact with cellular membranes. Functionally, SB056 showed high activity against Gram-negative bacteria and some limited activity against Gram-positive bacteria. Its potency against Gram-negative strains was comparable (on a molar basis) to that of colistin and polymyxin B, with an even broader spectrum of activity than numerous other reference compounds.