Molecular Mechanism of Ciprofloxacin Translocation Through the Major Diffusion Channels of the ESKAPE Pathogens Klebsiella pneumoniae and Enterobacter cloacae.

Experimental studies on the translocation and accumulation of antibiotics in Gram-negative bacteria have revealed details of the properties that allow efficient permeation through bacterial outer membrane porins. Among the major outer membrane diffusion channels, OmpF has been extensively studied to understand the antibiotic translocation process. In a few cases, this knowledge has also helped to improve the efficacy of existing antibacterial molecules. However, the extension of these strategies to enhance the efficacy of other existing and novel drugs require comprehensive molecular insight into the permeation process and an understanding of how antibiotic and channel properties influence the effective permeation rates. Previous studies have investigated how differences in antibiotic charge distribution can influence the observed permeation pathways through the OmpF channel, and have shown that the dynamics of the L3 loop can play a dominant role in the permeation process. Here, we perform all-atom simulations of the OmpF orthologs, OmpE35 from Enterobacter cloacae and OmpK35 from Klebsiella pneumoniae. Unbiased simulations of the porins and biased simulations of the ciprofloxacin permeation processes through these channels provide insight into the differences in the permeation pathway and energetics. In addition, we show that similar to the OmpF channel, antibiotic-induced dynamics of the L3 loop are also operative in the orthologs. However, the sequence and structural differences, influence the extent of the L3 loop fluctuations with OmpK35 showing greater stability in unbiased runs and subdued fluctuations in simulations with ciprofloxacin.

Antibiotic Charge Profile Determines the Extent of L3 Dynamics in OmpF: An Expedited Passage for Molecules with a Positive Charge.

Efficient permeation into Gram-negative bacterial cells is a much-desired property in the design of antibacterial agents. The goal is to arrive at one or more chemical modifications of molecules that improve their uptake into the cell while maintaining a good binding affinity to the intracellular target. Previously, we proposed a mechanistic rationale for the fast permeation of bulky antibiotics that involves induced conformational dynamics in the constriction loop L3 of the OmpF channel. This flexibility is caused by the perturbation of a hydrogen bond network stabilizing the L3 loop due to the strong interactions of the positively charged moiety on the antibiotic with the residues of the L3 loop. In the present work, we examine how differences in the charge profile of antibiotic molecules can affect the permeation process, in particular, the L3 dynamics. To this end, we have performed all-atom molecular dynamics simulations to study the permeation process of molecules with differences in the net charge through the Escherichia coli OmpF channel. The results from these simulations suggest that a positively charged moiety on the antibiotic is responsible for strong interactions with the negatively charged residues of the L3 loop, promoting conformational dynamics in the L3 loop. In contrast, antibiotics without a positively charged moiety are unable to initiate such a dynamic response in the L3 loop. This distinct behavior of the L3 loop in the presence of molecules with different charge characteristics provides a plausible mechanism whereby large molecules with an appropriate charge distribution can leverage an L3 dynamic-dependent pathway to permeate efficiently. The results are relevant to the structure-based design of molecules with improved uptake properties achieved through systematic chemical modifications that effectively engage the L3 loop.

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics.

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

Conformational Dynamics of Loop L3 in OmpF: Implications toward Antibiotic Translocation and Voltage Gating.

In the present work, we delineate the molecular mechanism of a bulky antibiotic permeating through a bacterial channel and uncover the role of conformational dynamics of the constriction loop in this process. Using the temperature accelerated sliced sampling approach, we shed light onto the dynamics of the L3 loop, in particular the F118 to S125 segment, at the constriction regions of the OmpF porin. We complement the findings with single channel electrophysiology experiments and applied-field simulations, and we demonstrate the role of hydrogen-bond stabilization in the conformational dynamics of the L3 loop. A molecular mechanism of permeation is put forward wherein charged antibiotics perturb the network of stabilizing hydrogen-bond interactions and induce conformational changes in the L3 segment, thereby aiding the accommodation and permeation of bulky antibiotic molecules across the constriction region. We complement the findings with single channel electrophysiology experiments and demonstrate the importance of the hydrogen-bond stabilization in the conformational dynamics of the L3 loop. The generality of the present observations and experimental results regarding the L3 dynamics enables us to identify this L3 segment as the source of gating. We propose a mechanism of OmpF gating that is in agreement with previous experimental data that showed the noninfluence of cysteine double mutants that tethered the L3 tip to the barrel wall on the OmpF gating behavior. The presence of similar loop stabilization networks in porins of other clinically relevant pathogens suggests that the conformational dynamics of the constriction loop is possibly of general importance in the context of antibiotic permeation through porins.

Atomistic Simulation of Molecules Interacting with Biological Nanopores: From Current Understanding to Future Directions.

Biological nanopores have been at the focus of numerous studies due to their role in many biological processes as well as their (prospective) technological applications. Among many other topics, recent studies on nanopores have addressed two key areas: antibiotic permeation through bacterial channels and sensing of analytes. Although the two areas are quite far apart in terms of their objectives, in both cases atomistic simulations attempt to understand the solute dynamics and the solute-protein interactions within the channel lumen. While decades of studies on various channels have culminated in an improved understanding of the key molecular factors and led to practical applications in some cases, successful utilization is limited. In this Perspective we summarize recent progress in understanding key issues in molecular simulations of antibiotic translocation and in the development of nanopore sensors. Moreover, we comment on possible advancements in computational algorithms that can potentially resolve some of the issues.

Rational design of hyperstable antibacterial peptides for food preservation.

We describe the design of peptides with properties like thermostability, pH stability, and antibacterial activity against a few bacterial food pathogens. Insights obtained from classical structure-function analysis of natural peptides and their mutants through antimicrobial and enzymatic assays are used to rationally develop a set of peptides. pH and thermostability assays were performed to demonstrate robust antimicrobial activity post-treatment with high temperatures and at wide pH ranges. We have also investigated the mode of action of these hyperstable peptides using membrane permeability assays, electron microscopy, and molecular dynamics simulations. Notably, through mutational studies, we show that these peptides elicit their antibacterial action via both membrane destabilization and inhibition of intracellular trypsin-the two functions attributable to separate peptide segments. Finally, toxicity studies and food preservation assays demonstrate the safety and efficacy of the designed peptides for food preservation. Overall, the study provides a general 'blueprint' for the development of stable antimicrobial peptides (AMPs). Insights obtained from this work may also be combined with combinatorial methods in high-throughput studies for future development of antimicrobials for various applications.

Improved Sampling and Free Energy Estimates for Antibiotic Permeation through Bacterial Porins.

Antibiotics enter into bacterial cells via protein channels that serve as low-energy pathways through the outer membrane, which is otherwise impenetrable. Insights into the molecular mechanisms underlying the transport processes are vital for the development of effective antibacterials. A much-desired prerequisite is an accurate and reproducible determination of free energy surfaces for antibiotic translocation, enabling quantitative and meaningful comparisons of permeation mechanisms for different classes of antibiotics. Inefficient sampling along the orthogonal degrees of freedom, for example, in umbrella sampling and metadynamics approaches, is however a key limitation affecting the accuracy and the convergence of free energy estimates. To overcome this limitation, two sampling methods have been employed in the present study that, respectively, combine umbrella sampling and metadynamics-style biasing schemes with temperature acceleration for improved sampling along orthogonal degrees of freedom. As a model for the transport of bulky solutes, the ciprofloxacin-OmpF system has been selected. The well-tempered metadynamics approach with multiple walkers is compared with its "temperature-accelerated" variant in terms of improvements in sampling and convergence of free energy estimates. We find that the inclusion of collective variables governing solute degrees of freedom and solute-water interactions within the sampling scheme largely alleviates sampling issues. Concerning improved sampling and convergence of free energy estimates from independent simulations, the temperature-accelerated sliced sampling approach that combines umbrella sampling with temperature-accelerated molecular dynamics performs even better as shown for the ciprofloxacin-OmpF system.

Molecular docking analysis of Clostridium perfringens beta toxin model with potential inhibitors from the ZINC database.

Beta toxin from Clostridium perfringens after being secreted in gut is capable of causing necrotic enteritis in humans and several other animal species and does not respond to routinely used antibiotics. Therefore, there is a need to design an effective inhibitor for the Clostridium perfringens beta toxin (CPB) using cutting edge drug discovery technologies. Hence, potential CPB inhibitors were identified using computer aided screening of compounds from the ZINC database. Further, we document the molecular docking analysis of Clostridium perfringens beta toxin model (that revealed 4 binding pockets, A-D) with the identified potential inhibitors. We show that ZINC291192 [N-[(1-methylindol-3-yl) methyl eneamino]-7,10-dioxabicyclo[4.4.0]deca-2,4,11-triene-8- carboxamide] has optimal binding features with calculated binding energy of -10.38 kcal/mol and inhibition constant of 24.76 nM for further consideration.

Structural basis of noncanonical polyphenol oxidase activity in DLL-II: A lectin from Dolichos lablab.

Lectins known to possess an additional enzymatic function are called leczymes. Previous studies reported a unique polyphenol oxidase (PPO) activity in DLL-II-a leczyme from Dolichos lablab. DLL-II shares a high sequence and structural homology with DBL-another leczyme from Dolichos biflorus. Incidentally, DBL possesses lipoxygenase activity, but not the PPO activity. Legume lectins usually possess two metal-binding sites A and B. Although these sites are conserved in both DBL and DLL-II, site A in DLL-II is occupied by Mn2+ and site B by Ca2+ . In contrast, DLL-II binds Cu2+ and Ca2+ at sites A and B, respectively. Here, investigating the structural basis of PPO activity in DLL-II, we find that the PPO activity is only dependent on Cu2+ , but not Ca2+ ; and the lectin activity requires only Ca2+ . Further, our analysis suggests that an alternative mechanism of PPO reaction may be operative in DLL-II, which involves a mononuclear Cu2+ metal center; this is in contrast to the bi-nuclear Cu2+ metal center commonly observed in all PPOs. Importantly, structural and computational approaches employed here, we hypothesize possible PPO binding sites and the corresponding migration channels for accessing the active site.

Structural plasticity mediates distinct GAP-dependent GTP hydrolysis mechanisms in Rab33 and Rab5.

The classical GTP hydrolysis mechanism, as seen in Ras, employs a catalytic glutamine provided in cis by the GTPase and an arginine supplied in trans by a GTPase activating protein (GAP). The key idea emergent from a large body of research on small GTPases is that GTPases employ a variety of different hydrolysis mechanisms; evidently, these variations permit diverse rates of GTPase inactivation, crucial for temporal regulation of different biological processes. Recently, we unified these variations and argued that a steric clash between active site residues (corresponding to positions 12 and 61 of Ras) governs whether a GTPase utilizes the cis-Gln or the trans-Gln (from the GAP) for catalysis. As the cis-Gln encounters a steric clash, the Rab GTPases employ the so-called dual finger mechanism where the interacting GAP supplies a trans-Gln for catalysis. Using experimental and computational methods, we demonstrate how the cis-Gln of Rab33 overcomes the steric clash when it is stabilized by a residue in the vicinity. In effect, this demonstrates how both cis-Gln- and trans-Gln-mediated mechanisms could operate in the same GTPase in different contexts, i.e. depending on the GAP that regulates its action. Interestingly, in the case of Rab5, which possesses a higher intrinsic GTP hydrolysis rate, a similar stabilization of the cis-Gln appears to overcome the steric clash. Taken together with the mechanisms seen for Rab1, it is evident that the observed variations in Rab and their GAP partners allow structural plasticity, or in other words, the choice of different catalytic mechanisms.

Disrupting domain-domain interactions is indispensable for EngA-ribosome interactions.

EngA consists of two tandem GTPase-domains-GD1 and GD2-followed by a KH-domain. EngA was considered to be a 50S assembly factor since it was shown to bind 50S and its deletion leads to the accumulation of immature 45S ribosomal subunits. Subsequently, we demonstrated an additional ribosome bound state of EngA bound to 50S, 30S, and 70S. While the former (50S binding) is achieved upon GTP binding at both GD1 and GD2, the latter is formed upon GTP hydrolysis at GD1, which is believed to trigger a large conformational change in the protein. The present study brings out two key aspects of EngA regulation: First, that distinctly stabilized GD1-KH interfaces allows EngA to exist in different ribosome bound states, and second is the importance of these states to ribosome assembly. Our analyses suggest that distinct inter-domain (GD-KH) interfaces are stabilized by interactions arising from unique sets of motifs, conserved across EngA homologues, and seem to be mechanistically linked to GTP/GDP binding. By experimentally measuring binding affinities for several interface mutants, we show that disrupting the interface interactions is necessary to realize EngA-ribosome binding. These findings are also supported by a recent cryo-EM structure of EngA bound to 50S, wherein the GD1-KH interface is completely disrupted leading to an 'extended' or 'open state' of the protein. Overall, it appears that the transition of EngA from a 'closed state' with GD1-KH forming a tight interface, to an 'open state' mediates interaction with ribosomal subunits.

Drug Target Identification and Prioritization for Treatment of Ovine Foot Rot: An In Silico Approach.

Ovine foot rot is an infection of the feet of sheep, mainly caused by Dichelobacter nodosus. In its virulent form, it is highly contagious and debilitating, causing significant losses in the form of decline in wool growth and quality and poor fertility. Current methods of treatment are ineffective in complete eradication. Effective antibiotic treatment of foot rot is hence necessary to ensure better outcomes during control phases by reduction in culling count and the possibility of carriers of the infection. Using computational approaches, we have identified a set of 297 proteins that are essential to the D. nodosus and nonhomologous with sheep proteins. These proteins may be considered as potential vaccine candidates or drug targets for designing antibiotics against the bacterium. This core set of drug targets have been analyzed for pathway annotation to identify 67 proteins involved in unique bacterial pathways. Choke-point analysis on the drug targets identified 138 choke-point proteins, 29 involved in unique bacterial pathways. Subcellular localization was also predicted for each target to identify the ones that are membrane associated or secreted extracellularly. In addition, a total of 13 targets were identified that are common in at least 10 pathogenic bacterial species.