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.

Supramolecular Hydrogels and Water Channels of Differing Diameters from Dipeptide Isomers

Dipeptides stereoisomers and regioisomers composed of norleucine (Nle) and phenylalanine (Phe) self-assemble into hydrogels under physiological conditions that are suitable for cell culture. The supramolecular behavior, however, differs as the packing modes comprise amphipathic layers or water channels, whose diameter is defined by either four or six dipeptide molecules. A variety of spectroscopy, microscopy, and synchrotron-radiation-based techniques unveil fine details of intermolecular interactions that pinpoint the relationship between the chemical structure and ability to form supramolecular architectures that define soft biomaterials.

Peptide Stereochemistry Effects from pKa-Shift to Gold Nanoparticle Templating in a Supramolecular Hydrogel

The divergent supramolecular behavior of a series of tripeptide stereoisomers was elucidated through spectroscopic, microscopic, crystallographic, and computational techniques. Only two epimers were able to effectively self-organize into amphipathic structures, leading to supramolecular hydrogels or crystals, respectively. Despite the similarity between the two peptides' turn conformations, stereoconfiguration led to different abilities to engage in intramolecular hydrogen bonding. Self-assembly further shifted the pKa value of the C-terminal side chain. As a result, across the pH range 4-6, only one epimer predominated sufficiently as a zwitterion to reach the critical molar fraction, allowing gelation. By contrast, the differing pKa values and higher dipole moment of the other epimer favored crystallization. The four stereoisomers were further tested for gold nanoparticle (AuNP) formation, with the supramolecular hydrogel being the key to control and stabilize AuNPs, yielding a nanocomposite that catalyzed the photodegradation of a dye. Importantly, the AuNP formation occurred without the use of reductants other than the peptide, and the redox chemistry was investigated by LC-MS, NMR, and infrared scattering-type near field optical microscopy (IR s-SNOM). This study provides important insights for the rational design of simple peptides as minimalistic and green building blocks for functional nanocomposites.

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.

Functionally distinct mutations within AcrB underpin antibiotic resistance in different lifestyles

Antibiotic resistance is a pressing healthcare challenge and is mediated by various mechanisms, including the active export of drugs via multidrug efflux systems, which prevent drug accumulation within the cell. Here, we studied how Salmonella evolved resistance to two key antibiotics, cefotaxime and azithromycin, when grown planktonically or as a biofilm. Resistance to both drugs emerged in both conditions and was associated with different substitutions within the efflux-associated transporter, AcrB. Azithromycin exposure selected for an R717L substitution, while cefotaxime for Q176K. Additional mutations in ramR or envZ accumulated concurrently with the R717L or Q176K substitutions respectively, resulting in clinical resistance to the selective antibiotics and cross-resistance to other drugs. Structural, genetic, and phenotypic analysis showed the two AcrB substitutions confer their benefits in profoundly different ways. R717L reduces steric barriers associated with transit through the substrate channel 2 of AcrB. Q176K increases binding energy for cefotaxime, improving recognition in the distal binding pocket, resulting in increased efflux efficiency. Finally, we show the R717 substitution is present in isolates recovered around the world.

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.

Putative binding modes (BMs) of quinolones to the bacterial efflux transporter MexB were identified. Multiple interaction patterns are possible, supporting the hypothesis that substrates oscillate between different BMs with similar affinity.

AB-DB: Force-Field parameters, MD trajectories, QM-based data, and Descriptors of Antimicrobials

Antibiotic resistance is a major threat to public health. The development of chemo-informatic tools to guide medicinal chemistry campaigns in the efficint design of antibacterial libraries is urgently needed. We present AB-DB, an open database of all-atom force-field parameters, molecular dynamics trajectories, quantum-mechanical properties, and curated physico-chemical descriptors of antimicrobial compounds. We considered more than 300 molecules belonging to 25 families that include the most relevant antibiotic classes in clinical use, such as β-lactams and (fluoro)quinolones, as well as inhibitors of key bacterial proteins. We provide traditional descriptors together with properties obtained with Density Functional Theory calculations. Noteworthy, AB-DB contains less conventional descriptors extracted from μs-long molecular dynamics simulations in explicit solvent. In addition, for each compound we make available force-field parameters for the major micro-species at physiological pH. With the rise of multi-drug-resistant pathogens and the consequent need for novel antibiotics, inhibitors, and drug re-purposing strategies, curated databases containing reliable and not straightforward properties facilitate the integration of data mining and statistics into the discovery of new antimicrobials.

Measurement(s)	molecular physical property analysis objective
Technology Type(s)	Computer Modeling

Single-atom substitution enables supramolecular diversity from dipeptide building blocks

Dipeptides are popular building blocks for supramolecular gels that do not persist in the environment and may find various applications. In this work, we show that a simple substitution on the aromatic side-chain of phenylalanine with either fluorine or iodine enables supramolecular diversity upon self-assembly at neutral pH, leading to hydrogels or crystals. Each building block is characterized by ¹H- and ¹³C-NMR spectroscopy, LC-MS, circular dichroism, and molecular models. The supramolecular behaviour is monitored with a variety of techniques, including circular dichroism, oscillatory rheology, transmission electron microscopy, attenuated total reflectance Fourier-transformed infrared spectroscopy, visible Raman spectroscopy, synchrotron-radiation single-crystal X-ray diffraction and UV Resonance Raman spectroscopy, allowing key differences to be pinpointed amongst the halogenated analogues.

Molecular Insights Into Binding and Activation of the Human KCNQ2 Channel by Retigabine

Voltage-gated potassium channels of the Kv7.x family are involved in a plethora of biological processes across many tissues in animals, and their misfunctioning could lead to several pathologies ranging from diseases caused by neuronal hyperexcitability, such as epilepsy, or traumatic injuries and painful diabetic neuropathy to autoimmune disorders. Among the members of this family, the Kv7.2 channel can form hetero-tetramers together with Kv7.3, forming the so-called M-channels, which are primary regulators of intrinsic electrical properties of neurons and of their responsiveness to synaptic inputs. Here, prompted by the similarity between the M-current and that in Kv7.2 alone, we perform a computational-based characterization of this channel in its different conformational states and in complex with the modulator retigabine. After validation of the structural models of the channel by comparison with experimental data, we investigate the effect of retigabine binding on the two extreme states of Kv7.2 (resting-closed and activated-open). Our results suggest that binding, so far structurally characterized only in the intermediate activated-closed state, is possible also in the other two functional states. Moreover, we show that some effects of this binding, such as increased flexibility of voltage sensing domains and propensity of the pore for open conformations, are virtually independent on the conformational state of the protein. Overall, our results provide new structural and dynamic insights into the functioning and the modulation of Kv7.2 and related channels.

Molecular rationale for the impairment of the MexAB-OprM efflux pump by a single mutation in MexA

Efflux pumps of the Resistance-Nodulation-cell Division (RND) superfamily contribute to intrinsic and acquired resistance in Gram-negative pathogens by expelling chemically unrelated antibiotics with high efficiency. They are tripartite systems constituted by an inner-membrane-anchored transporter, an outer membrane factor protein, and a membrane fusion protein. Multimerization of the membrane fusion protein is an essential prerequisite for full functionality of these efflux pumps. In this work, we employed complementary computational techniques to investigate the stability of a dimeric unit of MexA (the membrane fusion protein of the MexAB-OprM RND efflux pump of Pseudomonas aeruginosa), and to provide a molecular rationale for the effect of the G72S substitution, which affects MexAB-OprM functionality by impairing the assembly of MexA. Our findings indicate that: i) dimers of this protein are stable in multiple µs-long molecular dynamics simulations; ii) the mutation drastically alters the conformational equilibrium of MexA, favouring a collapsed conformation that is unlikely to form dimers or higher order assemblies. Unveiling the mechanistic aspects underlying large conformational distortions induced by minor sequence changes is informative to efforts at interfering with the activity of this elusive bacterial weapon. In this respect, our work further confirms how molecular simulations can give important contribution and useful insights to characterize the mechanism of highly complex biological systems.

Mechanistic Duality of Bacterial Efflux Substrates and Inhibitors: Example of Simple Substituted Cinnamoyl and Naphthyl Amides

Antibiotic resistance poses an immediate and growing threat to human health. Multidrug efflux pumps are promising targets for overcoming antibiotic resistance with small-molecule therapeutics. Previously, we identified a diaminoquinoline acrylamide, NSC-33353, as a potent inhibitor of the AcrAB-TolC efflux pump in Escherichia coli. This inhibitor potentiates the antibacterial activities of novobiocin and erythromycin upon binding to the membrane fusion protein AcrA. It is also a substrate for efflux and lacks appreciable intrinsic antibacterial activity of its own in wild-type cells. Here, we have modified the substituents of the cinnamoyl group of NSC-33353, giving rise to analogs that retain the ability to inhibit efflux, lost the features of the efflux substrates, and gained antibacterial activity in wild-type cells. The replacement of the cinnamoyl group with naphthyl isosteres generated compounds that lack antibacterial activity but are both excellent efflux pump inhibitors and substrates. Surprisingly, these inhibitors potentiate the antibacterial activity of novobiocin but not erythromycin. Surface plasmon resonance experiments and molecular docking suggest that the replacement of the cinnamoyl group with naphthyl shifts the affinity of the compounds away from AcrA to the AcrB transporter, making them better efflux substrates and changing their mechanism of inhibition. These results provide new insights into the duality of efflux substrate/inhibitor features in chemical scaffolds that will facilitate the development of new efflux pump inhibitors.

Bacterial efflux transporters' polyspecificity - a gift and a curse?

All mechanisms of clinical antibiotic resistance benefit from activities of polyspecific efflux pumps acting to reduce intracellular accumulation of toxins and antibiotics. In Gram-negative bacteria, the major polyspecific efflux transporters belong to the Resistance-Nodulation-cell Division (RND) superfamily of proteins, which are capable of expelling thousands of structurally diverse compounds. Recent structural and functional advances generated novel insights into mechanisms underlying the biochemical versatility of RND transporters. This opinion article reviews these mechanisms and discusses implications of the polyspecificity of RND transporters for bacterial survival and for the development of efflux pump inhibitors effective in clinics.

Nanoscale Assembly of Functional Peptides with Divergent Programming Elements

Self-assembling peptides are being applied both in the biomedical area and as building blocks in nanotechnology. Their applications are closely linked to their modes of self-assembly, which determine the functional nanostructures that they form. This work brings together two structural elements that direct nanoscale self-association in divergent directions: proline as a β-breaker and the β-structure-associated diphenylalanine motif, into a single tripeptide sequence. Amino acid chirality was found to resolve the tension inherent to these conflicting self-assembly instructions. Stereoconfiguration determined the ability of each of the eight possible Pro-Phe-Phe stereoisomers to self-associate into diverse nanostructures, including nanoparticles, nanotapes, or fibrils, which yielded hydrogels with gel-to-sol transition at a physiologically relevant temperature. Three single-crystal structures and all-atom molecular dynamics simulations elucidated the ability of each peptide to establish key interactions to form long-range assemblies (i,e., stacks leading to gelling fibrils), medium-range assemblies (i.e., stacks yielding nanotapes), or short-range assemblies (i.e., dimers or trimers that further associated into nanoparticles). Importantly, diphenylalanine is known to serve as a binding site for pathological amyloids, potentially allowing these heterochiral systems to influence the fibrillization of other biologically relevant peptides. To probe this hypothesis, all eight Pro-Phe-Phe stereoisomers were tested in vitro on the Alzheimer's disease-associated Aβ(1-42) peptide. Indeed, one nonfibril-forming stereoisomer effectively inhibited Aβ fibrillization through multivalent binding between diphenylalanine motifs. This work thus defined heterochirality as a useful feature to strategically develop future therapeutics to interfere with pathological processes, with the additional value of resistance to protease-mediated degradation and biocompatibility.

Heterochirality and Halogenation Control Phe-Phe Hierarchical Assembly

Diphenylalanine is an amyloidogenic building block that can form a versatile array of supramolecular materials. Its shortcomings, however, include the uncontrolled hierarchical assembly into microtubes of heterogeneous size distribution and well-known cytotoxicity. This study rationalized heterochirality as a successful strategy to address both of these pitfalls and it provided an unprotected heterochiral dipeptide that self-organized into a homogeneous and optically clear hydrogel with excellent ability to sustain fibroblast cell proliferation and viability. Substitution of one l-amino acid with its d-enantiomer preserved the ability of the dipeptide to self-organize into nanotubes, as shown by single-crystal XRD analysis, whereby the pattern of electrostatic and hydrogen bonding interactions of the backbone was unaltered. The effect of heterochirality was manifested in subtle changes in the positioning of the aromatic side chains, which resulted in weaker intermolecular interactions between nanotubes. As a result, d-Phe-l-Phe self-organized into homogeneous nanofibrils with a diameter of 4 nm, corresponding to two layers of peptides around a water channel, and yielded a transparent hydrogel. In contrast with homochiral Phe-Phe stereoisomer, it formed stable hydrogels thermoreversibly. d-Phe-l-Phe displayed no amyloid toxicity in cell cultures with fibroblast cells proliferating in high numbers and viability on this biomaterial, marking it as a preferred substrate over tissue-culture plastic. Halogenation also enabled the tailoring of d-Phe-l-Phe self-organization. Fluorination allowed analogous supramolecular packing as confirmed by XRD, thus nanotube formation, and gave intermediate levels of bundling. In contrast, iodination was the most effective strategy to augment the stability of the resulting hydrogel, although at the expense of optical transparency and biocompatibility. Interestingly, iodine presence hindered the supramolecular packing into nanotubes, resulting instead into amphipathic layers of stacked peptides without the occurrence of halogen bonding. By unravelling fine details to control these materials at the meso- and macro-scale, this study significantly advanced our understanding of these systems.

Cryo-EM Structure and Molecular Dynamics Analysis of the Fluoroquinolone Resistant Mutant of the AcrB Transporter from Salmonella

Salmonella is an important genus of Gram-negative pathogens, treatment of which has become problematic due to increases in antimicrobial resistance. This is partly attributable to the overexpression of tripartite efflux pumps, particularly the constitutively expressed AcrAB-TolC. Despite its clinical importance, the structure of the Salmonella AcrB transporter remained unknown to-date, with much of our structural understanding coming from the Escherichia coli orthologue. Here, by taking advantage of the styrene maleic acid (SMA) technology to isolate membrane proteins with closely associated lipids, we report the very first experimental structure of Salmonella AcrB transporter. Furthermore, this novel structure provides additional insight into mechanisms of drug efflux as it bears the mutation (G288D), originating from a clinical isolate of Salmonella Typhimurium presenting an increased resistance to fluoroquinolones. Experimental data are complemented by state-of-the-art molecular dynamics (MD) simulations on both the wild type and G288D variant of Salmonella AcrB. Together, these reveal several important differences with respect to the E. coli protein, providing insights into the role of the G288D mutation in increasing drug efflux and extending our understanding of the mechanisms underlying antibiotic resistance.

Self-Assembling l-d-l-Tripeptides Dance the Twist

Minimalistic peptides composed of d- and l-amino acids are attractive building blocks for functional supramolecular materials, including catalysts. d-Amino acids have long been known to promote turn conformations in peptides, yet unexpected twists continue to emerge on their effects on self-assembly. The combination of single-crystal X-ray diffraction and full-atom molecular dynamics have finally unraveled fine details of how l-d-l-tripeptides visit different conformations in solution and establish key interactions in supramolecular structures.

Binding and Transport of Carboxylated Drugs by the Multidrug Transporter AcrB

AcrAB(Z)-TolC is the main drug efflux transporter complex in Escherichia coli. The extrusion of various toxic compounds depends on several drug binding sites within the trimeric AcrB transporter. Membrane-localized carboxylated substrates, such as fusidic acid and hydrophobic β-lactams, access the pump via a groove between the transmembrane helices TM1 and TM2. In this article, the transport route from the initial TM1/TM2 groove binding site toward the deep binding pocket located in the periplasmic part has been addressed via molecular modeling studies followed by functional and structural characterization of several AcrB variants. We propose that membrane-embedded drugs bind initially to the TM1/TM2 groove, are oriented by the AcrB PN2 subdomain, and are subsequently transported via a PN2/PC1 interface pathway directly toward the deep binding pocket. Our work emphasizes the exploitation of multiple transport pathways by AcrB tuned to substrate physicochemical properties related to the polyspecificity of the pump.

Molecular Interactions of Cephalosporins with the Deep Binding Pocket of the RND Transporter AcrB

The drug/proton antiporter AcrB, part of the major efflux pump AcrABZ-TolC in Escherichia coli, is characterized by its impressive ability to transport chemically diverse compounds, conferring a multidrug resistance phenotype. However, the molecular features differentiating between good and poor substrates of the pump have yet to be identified. In this work, we combined molecular docking with molecular dynamics simulations to study the interactions between AcrB and two representative cephalosporins, cefepime and ceftazidime (a good and poor substrate of AcrB, respectively). Our analysis revealed different binding preferences of the two compounds toward the subsites of the large deep binding pocket of AcrB. Cefepime, although less hydrophobic than ceftazidime, showed a higher affinity than ceftazidime for the so-called hydrophobic trap, a region known for binding inhibitors and substrates. This supports the hypothesis that surface complementarity between the molecule and AcrB, more than the intrinsic hydrophobicity of the antibiotic, is a feature required for the interaction within this region. Oppositely, the preference of ceftazidime for binding outside the hydrophobic trap might not be optimal for triggering allosteric conformational changes needed to the transporter to accomplish its function. Altogether, our findings could provide valuable information for the design of new antibiotics less susceptible to the efflux mechanism.

Holo-like and Druggable Protein Conformations from Enhanced Sampling of Binding Pocket Volume and Shape

Understanding molecular recognition of small molecules by proteins in atomistic detail is key for drug design. Molecular docking is a widely used computational method to mimic ligand-protein association in silico. However, predicting conformational changes occurring in proteins upon ligand binding is still a major challenge. Ensemble docking approaches address this issue by considering a set of different conformations of the protein obtained either experimentally or from computer simulations, e.g., molecular dynamics. However, holo structures prone to host (the correct) ligands are generally poorly sampled by standard molecular dynamics simulations of the apo protein. In order to address this limitation, we introduce a computational approach based on metadynamics simulations called ensemble docking with enhanced sampling of pocket shape (EDES) that allows holo-like conformations of proteins to be generated by exploiting only their apo structures. This is achieved by defining a set of collective variables that effectively sample different shapes of the binding site, ultimately mimicking the steric effect due to the ligand. We assessed the method on three challenging proteins undergoing different extents of conformational changes upon ligand binding. In all cases our protocol generates a significant fraction of structures featuring a low RMSD from the experimental holo geometry. Moreover, ensemble docking calculations using those conformations yielded in all cases native-like poses among the top-ranked ones.

Relevance of Ebola virus VP35 homo-dimerization on the type I interferon cascade inhibition

Ebola virus high lethality relies on its ability to efficiently bypass the host innate antiviral response, which senses the viral dsRNA through the RIG-I receptor and induces type I interferon α/β production. In the bypassing action, the Ebola virus protein VP35 plays a pivotal role at multiple levels of the RIG-I cascade, masking the viral 5′-triphosphorylated dsRNA from RIG-I, and interacting with other cascade components. The VP35 type I interferon inhibition is exerted by the C-terminal domain, while the N-terminal domain, containing a coiled-coil region, is primarily required for oligomerization. However, mutations at key VP35 residues L90/93/107A (VP35-3m) in the coiled-coil region were reported to affect oligomerization and reduce type I interferon antagonism, indicating a possible but unclear role of homo-oligomerization on VP35 interaction with the RIG-I pathway components. In this work, we investigated the VP35 dimerization thermodynamics and its contribution to type I interferon antagonism by computational and biological methods. Focusing on the coiled-coil region, we combined coarse-grained and all-atom simulations on wild type VP35 and VP35-3m homo-dimerization. According to our results, wild type VP35 coiled-coil is able to self-assemble into dimers, while VP35-3m coiled-coil shows poor propensity to even dimerize. Free-energy calculations confirmed the key role of L90, L93 and L107 in stabilizing the coiled-coil homo-dimeric structure. In vitro type I interferon antagonism studies, using full-length wild type VP35 and VP35-3m, revealed that VP35 homo-dimerization is an essential preliminary step for dsRNA binding, which appears to be the main factor of the VP35 RIG-I cascade inhibition, while it is not essential to block the other steps.

Identification and characterization of carbapenem binding sites within the RND-transporter AcrB

Understanding the molecular determinants for recognition, binding and transport of antibiotics by multidrug efflux systems is important for basic research and useful for the design of more effective antimicrobial compounds. Imipenem and meropenem are two carbapenems whose antibacterial activity is known to be poorly and strongly affected by MexAB-OprM, the major efflux pump transporter in Pseudomonas aeruginosa. However, not much is known regarding recognition and transport of these compounds by AcrAB-TolC, which is the MexAB-OprM homologue in Escherichia coli and by definition the paradigm model for structural studies on efflux pumps. Prompted by this motivation, we unveiled the molecular details of the interaction of imipenem and meropenem with the transporter AcrB by combining computer simulations with biophysical experiments. Regarding the interaction with the two main substrate binding regions of AcrB, the so-called access and deep binding pockets, molecular dynamics simulations revealed imipenem to be more mobile than meropenem in the former, while comparable mobilities were observed in the latter. This result is in line with isothermal titration calorimetry, differential scanning experiments, and binding free energy calculations, indicating a higher affinity for meropenem than imipenem at the deep binding pocket, while both sharing similar affinities at the access pocket. Our findings rationalize how different physico-chemical properties of compounds reflect on their interactions with AcrB. As such, they constitute precious information to be exploited for the rational design of antibiotics able to evade efflux pumps.

Molecular Determinants of the Promiscuity of MexB and MexY Multidrug Transporters of Pseudomonas aeruginosa

Secondary multidrug transporters of the resistance-nodulation-cell division (RND) superfamily contribute crucially to antibiotic resistance in Gram-negative bacteria. Compared to the most studied transporter AcrB of Escherichia coli, little is known about the molecular determinants of distinct polyspecificities of the most important RND transporters MexB and MexY of Pseudomonas aeruginosa. In an effort to add knowledge on this topic, we performed an exhaustive atomic-level comparison of the main putative recognition sites (access and deep binding pockets) in these two Mex transporters. We identified an underlying link between some structural, chemical and dynamical features of the binding pockets and the physicochemical nature of the corresponding substrates recognized by either one or both pumps. In particular, mosaic-like lipophilic and electrostatic surfaces of the binding pockets provide for both proteins several multifunctional sites for diffuse binding of diverse substrates. Specific lipophilicity signatures of the weakly conserved deep pocket suggest a key role of this site as a selectivity filter as in Acr transporters. Finally, the different dynamics of the bottom-loop in MexB and MexY support its possible role in binding of large substrates. Our work represents the first comparative study of the major RND transporters in P. aeruginosa and also the first structure-based study of MexY, for which no experimental structure is available yet.

Molecular Modeling of Multidrug Properties of Resistance Nodulation Division (RND) Transporters

Efflux pumps of the resistance nodulation division (RND) superfamily are among the major contributors to intrinsic and acquired multidrug resistance in Gram-negative bacteria. Structural information on AcrAB-TolC and MexAB-OprM, major efflux pumps of Escherichia coli and Pseudomonas aeruginosa respectively, boosted intensive research aimed at understanding the molecular mechanisms ruling the active extrusion processes. In particular, several studies were devoted to the understanding of the determinants behind the extraordinary broad specificity of the RND transporters AcrB and MexB. In this chapter, we discuss the ever-growing role computational methods have been playing in deciphering key structural and dynamical features of these transporters and of their interaction with substrates and inhibitors. We further discuss and illustrate examples from our lab of how molecular docking, homology modeling, all-atom molecular dynamics simulations and in silico free energy estimations can all together give precious insights into the processes of recognition and extrusion of substrates, as well as on the possible inhibition strategies.

Molecular Rationale behind the Differential Substrate Specificity of Bacterial RND Multi-Drug Transporters

Resistance-Nodulation-cell Division (RND) transporters AcrB and AcrD of Escherichia coli expel a wide range of substrates out of the cell in conjunction with AcrA and TolC, contributing to the onset of bacterial multidrug resistance. Despite sharing an overall sequence identity of ~66% (similarity ~80%), these RND transporters feature distinct substrate specificity patterns whose underlying basis remains elusive. We performed exhaustive comparative analyses of the putative substrate binding pockets considering crystal structures, homology models and conformations extracted from multi-copy μs-long molecular dynamics simulations of both AcrB and AcrD. The impact of physicochemical and topographical properties (volume, shape, lipophilicity, electrostatic potential, hydration and distribution of multi-functional sites) within the pockets on their substrate specificities was quantitatively assessed. Differences in the lipophilic and electrostatic potentials among the pockets were identified. In particular, the deep pocket of AcrB showed the largest lipophilicity convincingly pointing out its possible role as a lipophilicity-based selectivity filter. Furthermore, we identified dynamic features (not inferable from sequence analysis or static structures) such as different flexibilities of specific protein loops that could potentially influence the substrate recognition and transport profile. Our findings can be valuable for drawing structure (dynamics)-activity relationship to be employed in drug design.

Effect of site-directed mutations in multidrug efflux pump AcrB examined by quantitative efflux assays

BACKGROUND

The Resistance-Nodulation-Division (RND) family transporter AcrB plays a major role in the intrinsic and increased resistance of Escherichia coli to a large number of antibiotics. The distal binding pocket within this multidrug efflux transporter is very large, but the effort to define the roles of various residues facing this pocket through site-directed mutagenesis so far involved only the determination of minimal inhibitory concentrations of drugs in mutants.

METHODS

We measured in intact E. coli cells the kinetics of efflux of two substrates, nitrocefin (a cephalosporin) that is predicted mainly to bind to the upper, "groove" domain of the pocket, and L-alanyl-β-naphthylamide (Ala-Naph) that is likely to bind to the lower, "cave" domain, in a number of site-directed mutants of AcrB, where a hydrophobic or aromatic residue was changed into alanine.

RESULTS

The efflux of nitrocefin became attenuated by some mutations in the groove domain, such as I278A and F178A, but in some experiments a mutation in the cave domain, F628A produced a similar result. In some cases an increased value of K_M was detected. The efflux of Ala-Naph was increased by mutations in the cave domain, such as F136A and I626A, but also by those in the groove domain (I277A, I278A, F178A). In most cases the increased V_max values appeared to be responsible. F610A mutation had a profound effect on the efflux of both substrates, as reported earlier.

CONCLUSIONS

Our data show for the first time effects of various substrate-binding pocket mutations on the kinetics of efflux of two substrates by the AcrB pump. They also confirm interactions between substrates and drugs predicted by MD simulation studies, and also reveal areas that need future research.

The Phe-Phe Motif for Peptide Self-Assembly in Nanomedicine

Since its discovery, the Phe-Phe motif has gained in popularity as a minimalist building block to drive the self-assembly of short peptides and their analogues into nanostructures and hydrogels. Molecules based on the Phe-Phe motif have found a range of applications in nanomedicine, from drug delivery and biomaterials to new therapeutic paradigms. Here we discuss the various production methods for this class of compounds, and the characterization, nanomorphologies, and application of their self-assembled nanostructures. We include the most recent findings on their remarkable properties, which hold substantial promise for the creation of the next generation nanomedicines.

Structure-Activity Relationship Study of Hydroxycoumarins and Mushroom Tyrosinase

The structure-activity relationships of four hydroxycoumarins, two with the hydroxyl group on the aromatic ring of the molecule and two with the hydroxyl group replacing hydrogen of the pyrone ring, and their interactions with mushroom tyrosinase were studied. These compounds displayed different behaviors upon action of the enzyme. The two compounds, ar-hydroxylated 6-hydroxycoumarin and 7-hydroxycoumarin, were both weak substrates of the enzyme. Interestingly, in both cases, the product of the catalysis was the 6,7-hydroxycoumarin, although 5,6- and 7,8-isomers could also theoretically be formed. Additionally, both were able to reduce the formation of dopachrome when tyrosinase acted on its typical substrate, L-tyrosine. Although none of the compounds that contained a hydroxyl group on the pyrone ring were substrates of tyrosinase, the 3-hydroxycoumarin was a potent inhibitor of the enzyme, and the 4-hydroxycoumarin was not an inhibitor. These results were compared with those obtained by in silico molecular docking predictions to obtain potentially useful information for the synthesis of new coumarin-based inhibitors that resemble the structure of the 3-hydroxycoumarin.

RND efflux pumps: structural information translated into function and inhibition mechanisms

Efflux pumps of the Resistance Nodulation Division (RND) superfamily play a major role in the intrinsic and acquired resistance of Gram-negative pathogens to antibiotics. Moreover, they are largely responsible for multi-drug resistance (MDR) phenomena in these bacteria. The last decade has seen a sharp increase in the number of experimental and computational studies aimed at understanding their functional mechanisms. Most of these studies focused on the RND drug/proton antiporter AcrB, part of the AcrAB-TolC efflux pump actively recognizing and expelling noxious agents from the interior of bacteria. These studies have been focused on the dynamical interactions between AcrB and its substrates and inhibitors, on the details of the proton translocation mechanisms, and on the way AcrB assembles with protein partners to build up a functional pump. In this review we summarize these advances focusing on the role of AcrB.

Some ligands enhance the efflux of other ligands by the Escherichia coli multidrug pump AcrB

By measuring quantitatively the active efflux of cephalosporins by the RND (resistance-nodulation-division) family efflux pump AcrB in intact cells of Escherichia coli, we found that the simultaneous presence of another substrate, such as chloramphenicol, benzene, cyclohexane, or Arg β-naphthilamide, significantly enhanced the extrusion of cephalosporins. The stimulation occurred also in a strain expressing the covalently linked trimer of AcrB, and thus cannot be ascribed to the enhanced assembly of the trimer from AcrB monomers. When Val139 of AcrB was changed into Phe, the stimulation by benzene was found to occur at much lower concentration of the solvent. A plausible explanation of these observations is that the AcrB pump is constructed to pump out very rapidly the solvent or chloramphenicol molecules, and thus the efflux of cephalosporins, which presumably bind to a different subsite within the large binding pocket of AcrB, can become facilitated. Computer simulations of ligand binding to AcrB, both by docking and by molecular dynamics simulations, produced results supporting and extending this hypothesis. Benzene and the cephalosporin nitrocefin can bind simultaneously to the distal binding pocket of AcrB, both in the wild type and in the V139F variant. Interestingly, while the binding position and strength of benzene are almost unaffected by the presence of nitrocefin, this latter substrate is significantly displaced toward the exit gate in both wild type and mutant transporter in the presence of benzene. Additionally, the cephalosporin efflux may be enhanced by the binding of solvents (sometimes to the cephalosporin-free protomer), which could accelerate AcrB conformational changes necessary for substrate extrusion.

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.

Molecular Dynamics Computer Simulations of Multidrug RND Efflux Pumps

Over-expression of multidrug efflux pumps of the Resistance Nodulation Division (RND) protein super family counts among the main causes for microbial resistance against pharmaceuticals. Understanding the molecular basis of this process is one of the major challenges of modern biomedical research, involving a broad range of experimental and computational techniques. Here we review the current state of RND transporter investigation employing molecular dynamics simulations providing conformational samples of transporter components to obtain insights into the functional mechanism underlying efflux pump-mediated antibiotics resistance in Escherichia coli and Pseudomonas aeruginosa.

Recognition of imipenem and meropenem by the RND-transporter MexB studied by computer simulations

Basic understanding of the means by which multidrug efflux systems can efficiently recognize and transport drugs constitutes a fundamental step toward development of compounds able to tackle the continuous outbreak of new bacterial strains resistant to traditional antibiotics. We applied a series of computational techniques, from molecular docking to molecular dynamics simulations and free energy estimate methods, to determine the differences in the binding properties of imipenem and meropenem, two potent antibiotics of the carbapenem family, to MexB, the RND transporter of the major efflux system of Pseudomonas aeruginosa. We identified and characterized two affinity sites in the periplasmic domain of the transporter, sharing strong similarities with the distal and proximal binding pockets identified in AcrB, the homologue of MexB in Escherichia coli. According to our results, meropenem has a higher affinity to the distal binding pocket than imipenem while both compounds are weakly bound to the proximal pocket. This different behavior is mainly due to the hydration properties of the nonpharmacophore part of the two compounds, being that of imipenem less bulky and hydrophobic. Our data provide for the first time a rationale at molecular level for the experimental evidence indicating meropenem as a compound strongly affected by MexB contrary to imipenem, which is apparently poorly transported by the same pump.

A kinetic Monte Carlo approach to investigate antibiotic translocation through bacterial porins

Many relevant biological processes take place on time scales not reachable by standard all-atom computer simulations. The translocation of antibiotics through non-specific bacterial porins is an example. Microscopic effects compete to determine penetration routes and, consequently, free energy barriers to be overcome. Since bacteria can develop resistance to treatment also by reducing their antibiotic permeability, to understand the microscopic aspects of antibiotic translocation is an important step to rationalize drug design. Here, to investigate the translocation we propose a complete numerical model that combines the diffusion-controlled rate theory and a kinetic Monte Carlo scheme based on both experimental data and microscopically well-founded all-atom simulations. Within our model, an antibiotic translocating through an hour-glass-shaped channel can be described as a molecule moving on a potential of mean force featuring several affinity sites and a high central barrier. The implications of our results for the characterization of antibiotic translocation at in vivo concentrations are discussed. The presence of an affinity site close to the mouth of the channel seems to favor the translocation of antibiotics, the affinity site acting as a particle reservoir. Possible connections between results and the appearance of mutations in clinical strains are also outlined.

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.

Effect of the F610A mutation on substrate extrusion in the AcrB transporter: explanation and rationale by molecular dynamics simulations

The tripartite efflux pump AcrAB-TolC is responsible for the intrinsic and acquired multidrug resistance in Escherichia coli. Its active part, the homotrimeric transporter AcrB, is in charge of the selective binding of substrates and energy transduction. The mutation F610A has been shown to significantly reduce the minimum inhibitory concentration of doxorubicin and many other substrates, although F610 does not appear to interact strongly with them. Biochemical study of transport kinetics in AcrB is not yet possible, except for some β-lactams, and other techniques should supply this important information. Therefore, in this work, we assess the impact of the F610A mutation on the functionality of AcrB by means of computational techniques, using doxorubicin as substrate. We found that the compound slides deeply inside the binding pocket after mutation, increasing the strength of the interaction. During subsequent conformational alterations of the transporter, doxorubicin was either not extruded from the binding site or displaced along a direction other than the one associated with extrusion. Our study indicates how subtle interactions determine the functionality of multidrug transporters, since decreased transport might not be simplistically correlated to decreased substrate binding affinity.

Functional rotation of the transporter AcrB: insights into drug extrusion from simulations

The tripartite complex AcrAB-TolC is the major efflux system in Escherichia coli. It extrudes a wide spectrum of noxious compounds out of the bacterium, including many antibiotics. Its active part, the homotrimeric transporter AcrB, is responsible for the selective binding of substrates and energy transduction. Based on available crystal structures and biochemical data, the transport of substrates by AcrB has been proposed to take place via a functional rotation, in which each monomer assumes a particular conformation. However, there is no molecular-level description of the conformational changes associated with the rotation and their connection to drug extrusion. To obtain insights thereon, we have performed extensive targeted molecular dynamics simulations mimicking the functional rotation of AcrB containing doxorubicin, one of the two substrates that were co-crystallized so far. The simulations, including almost half a million atoms, have been used to test several hypotheses concerning the structure-dynamics-function relationship of this transporter. Our results indicate that, upon induction of conformational changes, the substrate detaches from the binding pocket and approaches the gate to the central funnel. Furthermore, we provide strong evidence for the proposed peristaltic transport involving a zipper-like closure of the binding pocket, responsible for the displacement of the drug. A concerted opening of the channel between the binding pocket and the gate further favors the displacement of the drug. This microscopically well-funded information allows one to identify the role of specific amino acids during the transitions and to shed light on the functioning of AcrB.

In nature, bacteria have to resist several toxic threats to be able to survive, from bile acids in intestines up to antibiotics. The Escherichia coli bacterium, which usually is a commensal inhabitant of human intestines, can also acquire pathogenic properties which would harm the human body. To dispose of toxic compounds, E. coli has developed a protein machinery which is called “efflux pump”. Here, we studied the dynamics of the transporter protein AcrB, a component of the E. coli major efflux system, in complex with an antibiotic (doxorubicin). We used computer simulations to complement the existing experimental data. Our purpose was to gain more detailed insights into the pumping mechanism at the molecular level. In our simulations the drug leaves the binding pocket upon induction of functional rotation in the protein, although a complete extrusion was never observed. A peristaltic motion, which starts with a zipper-like closure of the interior of the protein, is an important step for the extrusion of the drug. Interestingly, such a peristaltic mechanism of pumping has been suggested before on the basis of structural data. The molecular details obtained in this study shall deepen the understanding of the functioning of the efflux pump.