Optimum quorum quenching bacteria concentration in the better-quality cell entrapping beads to control biofouling in membrane bioreactor.

Quorum quenching (QQ) by cell entrapping beads (CEBs) is known to inhibit biofouling by its biological and physical cleaning effect. Although there are better QQ media reported, due to the ease of fabrication of QQ-CEBs, this study focused on improving the quality of CEBs by comparing two distinct bead-making methods - polyvinyl alcohol-alginate (PVA-alginate) and phase inversion - and on finding the optimum concentration of QQ bacteria in the CEBs. The evaluation of PVA-alginate bead showed better uniformity, and higher mechanical and chemical strength in comparison with the phase inversion bead. Through the operations of two control membrane bioreactors (MBRs) (no bead, vacant bead) and four QQ-MBRs with different Rhodococcus sp. BH4 concentrations (2.5-15 mg cell ml-1) in PVA-alginate CEBs, the maximum QQ effect was observed by 5 mg ml-1 BH4 concentration beads. This implies that an optimum cell concentration of QQ-CEBs is crucial to economically improve MBR performance using QQ.

Effect of Aeration Intensity on Performance of Lab-Scale Quorum-Quenching Membrane Bioreactor.

Biofouling is one of the main drawbacks of membrane bioreactors (MBRs). Among the different methods, the quorum-quenching (QQ) technique is a novel method as it delays biofilm formation on the membrane surface through disruption of bacterial cell-to-cell communication and thus effectively mitigates membrane biofouling. QQ bacteria require a certain concentration of dissolved oxygen to show their best activities. Despite the importance of the amount of aeration, there have not been enough studies on aeration condition utilizing the separate determination of pure QQ effect and physical cleaning effect. This research aimed to find the optimum aeration intensity by separation of the two effects from QQ and physical cleaning. Three bead type conditions (no bead, vacant bead, and QQ beads) at three aeration intensities (1.5, 2.5, and 3.5 L/min representing low, medium, and high aeration intensity) were applied. From the results, no QQ effect and small QQ effect were observed at low and high aeration, while the greatest QQ effect (48.2% of 737 h improvement) was observed at medium aeration. The best performance was observed at high aeration with QQ beads having a 1536 h operational duration (303% improvement compared to the no bead condition); however, this excellent performance was attributed more to the physical cleaning effect than to the QQ effect.

Inhibition of Voltage-Gated K+ Channel Kv1.5 by Antiarrhythmic Drugs.

Molecular dynamics simulations are employed to determine the inhibitory mechanisms of three drugs, 5-(4-phenoxybutoxy)psoralen (PAP-1), vernakalant, and flecainide, on the voltage-gated K+ channel Kv1.5, a target for the treatment of cardiac arrhythmia. At neutral pH, PAP-1 is neutral, whereas the other two molecules carry one positive charge. We show that PAP-1 forms stable dimers in water, primarily through hydrophobic interactions between aromatic rings. All three molecules bind to the cavity between the Ile508 and Val512 residues from the four subunits of the channel. Once bound, the drug molecules are flexible, with the average root-mean-square fluctuation being between 2 and 3 Å, which is larger than the radius of gyration of a bulky amino acid. The presence of a monomeric PAP-1 causes the permeating K+ ion to dehydrate, thereby creating a significant energy barrier. In contrast, vernakalant blocks the ion permeation primarily via an electrostatic mechanism and, therefore, must be in the protonated and charged form to be effective.

Functionalized Fullerene Targeting Human Voltage-Gated Sodium Channel, hNav1.7.

Mutations of hNav1.7 that cause its activities to be enhanced contribute to severe neuropathic pain. Only a small number of hNav1.7 specific inhibitors have been identified, most of which interact with the voltage-sensing domain of the voltage-activated sodium ion channel. In our previous computational study, we demonstrated that a [Lys6]-C84 fullerene binds tightly (affinity of 46 nM) to NavAb, the voltage-gated sodium channel from the bacterium Arcobacter butzleri. Here, we extend this work and, using molecular dynamics simulations, demonstrate that the same [Lys6]-C84 fullerene binds strongly (2.7 nM) to the pore of a modeled human sodium ion channel hNav1.7. In contrast, the fullerene binds only weakly to a mutated model of hNav1.7 (I1399D) (14.5 mM) and a model of the skeletal muscle hNav1.4 (3.7 mM). Comparison of one representative sequence from each of the nine human sodium channel isoforms shows that only hNav1.7 possesses residues that are critical for binding the fullerene derivative and blocking the channel pore.

Molecular dynamics of the honey bee toxin tertiapin binding to Kir3.2.

Tertiapin (TPN), a short peptide isolated from the venom of the honey bee, is a potent and selective blocker of the inward rectifier K+ (Kir) channel Kir3.2. Here we examine in atomic detail the binding mode of TPN to Kir3.2 using molecular dynamics, and deduce the key residue in Kir3.2 responsible for TPN selectivity. The binding of TPN to Kir3.2 is stable when the side chain of either Lys16 (TPNK16-Kir3.2) or Lys17 (TPNK17-Kir3.2) of the toxin protrudes into the channel pore. However, the binding affinity calculated from only TPNK17-Kir3.2 and not TPNK16-Kir3.2 is consistent with experiment, suggesting that Lys17 is the most plausible pore-blocking residue. The alanine mutation of Kir3.2-Glu127, which is not present in TPN-resistant channels, reduces the inhibitory ability of TPN by over 50 fold in TPNK17-Kir3.2, indicating that Kir3.2-Glu127 is important for the selectivity of TPN.

Peptidomimetic Star Polymers for Targeting Biological Ion Channels.

Four end-functionalized star polymers that could attenuate the flow of ionic currents across biological ion channels were first de novo designed computationally, then synthesized and tested experimentally on mammalian K+ channels. The 4-arm ethylene glycol conjugate star polymers with lysine or a tripeptide attached to the end of each arm were specifically designed to mimic the action of scorpion toxins on K+ channels. Molecular dynamics simulations showed that the lysine side chain of the polymers physically occludes the pore of Kv1.3, a target for immuno-suppression therapy. Two of the compounds tested were potent inhibitors of Kv1.3. The dissociation constants of these two compounds were computed to be 0.1 µM and 0.7 µM, respectively, within 3-fold to the values derived from subsequent experiments. These results demonstrate the power of computational methods in molecular design and the potential of star polymers as a new infinitely modifiable platform for ion channel drug discovery.

Computational Studies of Venom Peptides Targeting Potassium Channels.

Small peptides isolated from the venom of animals are potential scaffolds for ion channel drug discovery. This review article mainly focuses on the computational studies that have advanced our understanding of how various toxins interfere with the function of K⁺ channels. We introduce the computational tools available for the study of toxin-channel interactions. We then discuss how these computational tools have been fruitfully applied to elucidate the mechanisms of action of a wide range of venom peptides from scorpions, spiders, and sea anemone.

Molecular dynamics simulations of Na(+) and leucine transport by LeuT.

Molecular dynamics simulations are used to gain insight into the binding of Na(+) and leucine substrate to the bacterial amino acid transporter LeuT, focusing on the crystal structures of LeuT in the outward-open and inward-open states. For both conformations of LeuT, a third Na(+) binding site involving Glu290 in addition to the two sites identified from the crystal structures is observed. Once the negative charge from Glu290 in the inward-open LeuT is removed, the ion bound to the third site is ejected from LeuT rapidly, suggesting that the protonation state of Glu290 regulates Na(+) binding and release. In Cl(-)-dependent transporters where Glu290 is replaced by a neutral serine, a Cl(-) ion would be required to replace the role of Glu290. Thus, the simulations provide insights into understanding Na(+) and substrate transport as well as Cl(-)-independence of LeuT.

Mechanisms and energetics of potassium channel block by local anesthetics and antifungal agents.

Many drug molecules inhibit the conduction of several families of cation channels by binding to a small cavity just below the selectivity filter of the channel protein. The exact mechanisms governing drug-channel binding and the subsequent inhibition of conduction are not well understood. Here the inhibition of two K(+) channel isoforms, Kv1.2 and KCa3.1, by two drug molecules, lidocaine and TRAM-34, is examined in atomic detail using molecular dynamics simulations. A conserved valine-alanine-valine motif in the inner cavity is found to be crucial for drug binding in both channels, consistent with previous studies of similar systems. Potential of mean force calculations show that lidocaine in its charged form creates an energy barrier of ∼6 kT for a permeating K(+) ion when the ion is crossing over the drug, while the neutral form of lidocaine has no significant effect on the energetics of ion permeation. On the other hand, TRAM-34 in the neutral form is able to create a large energy barrier of ∼10 kT by causing the permeating ion to dehydrate. Our results suggest that TRAM-34 analogues that remain neutral and permeable to membranes under acidic conditions common to inflammation may act as possible drug scaffolds for combating local anesthetic failure in inflammation.

Binding modes of two scorpion toxins to the voltage-gated potassium channel kv1.3 revealed from molecular dynamics.

Molecular dynamics (MD) simulations are used to examine the binding modes of two scorpion toxins, margatoxin (MgTx) and hongotoxin (HgTx), to the voltage gated K+ channel, Kv1.3. Using steered MD simulations, we insert either Lys28 or Lys35 of the toxins into the selectivity filter of the channel. The MgTx-Kv1.3 complex is stable when the side chain of Lys35 from the toxin occludes the channel filter, suggesting that Lys35 is the pore-blocking residue for Kv1.3. In this complex, Lys28 of the toxin forms one additional salt bridge with Asp449 just outside the filter of the channel. On the other hand, HgTx forms a stable complex with Kv1.3 when the side chain of Lys28 but not Lys35 protrudes into the filter of the channel. A survey of all the possible favorable binding modes of HgTx-Kv1.3 is carried out by rotating the toxin at 3° intervals around the channel axis while the position of HgTx-Lys28 relative to the filter is maintained. We identify two possible favorable binding modes: HgTx-Arg24 can interact with either Asp433 or Glu420 on the vestibular wall of the channel. The dissociation constants calculated from the two binding modes of HgTx-Kv1.3 differ by approximately 20 fold, suggesting that the two modes are of similar energetics.

Binding of fullerenes and nanotubes to MscL.

Multi-drug resistance is becoming an increasing problem in the treatment of bacterial infections and diseases. The mechanosensitive channel of large conductance (MscL) is highly conserved among prokaryotes. Evidence suggests that a pharmacological agent that can affect the gating of, or block the current through, MscL has significant potential as a new class of antimicrobial compound capable of targeting a range of pathogenic bacteria with minimal side-effects to infected patients. Using molecular dynamics we examine the binding of fullerenes and nanotubes to MscL and demonstrate that both are stable within the MscL pore. We predict that fullerenes will attenuate the flow of ions through MscL by reducing the pore volume available to water and ions, but nanotubes will prevent pore closure resulting in a permanently open pore. Moreover, we confirm experimentally that it is possible to attenuate the flow of ions through MscL using a C60-γ cyclodextrin complex.

Mechanism of µ-conotoxin PIIIA binding to the voltage-gated Na+ channel NaV1.4.

Several subtypes of voltage-gated Na+ (NaV) channels are important targets for pain management. µ-Conotoxins isolated from venoms of cone snails are potent and specific blockers of different NaV channel isoforms. The inhibitory effect of µ-conotoxins on NaV channels has been examined extensively, but the mechanism of toxin specificity has not been understood in detail. Here the known structure of µ-conotoxin PIIIA and a model of the skeletal muscle channel NaV1.4 are used to elucidate elements that contribute to the structural basis of µ-conotoxin binding and specificity. The model of NaV1.4 is constructed based on the crystal structure of the bacterial NaV channel, NaVAb. Six different binding modes, in which the side chain of each of the basic residues carried by the toxin protrudes into the selectivity filter of NaV1.4, are examined in atomic detail using molecular dynamics simulations with explicit solvent. The dissociation constants (Kd) computed for two selected binding modes in which Lys9 or Arg14 from the toxin protrudes into the filter of the channel are within 2 fold; both values in close proximity to those determined from dose response data for the block of NaV currents. To explore the mechanism of PIIIA specificity, a double mutant of NaV1.4 mimicking NaV channels resistant to µ-conotoxins and tetrodotoxin is constructed and the binding of PIIIA to this mutant channel examined. The double mutation causes the affinity of PIIIA to reduce by two orders of magnitude.

Mechanism of tetrodotoxin block and resistance in sodium channels.

Tetrodotoxin (TTX) has been used for many decades to characterize the structure and function of biological ion channels. Yet, the precise mechanism by which TTX blocks voltage-gated sodium (NaV) channels is not fully understood. Here molecular dynamics simulations are used to elucidate how TTX blocks mammalian voltage-gated sodium (Nav) channels and why it fails to be effective for the bacterial sodium channel, NaVAb. We find that, in NaVAb, a sodium ion competes with TTX for the binding site at the extracellular end of the filter, thus reducing the blocking efficacy of TTX. Using a model of the skeletal muscle channel, NaV1.4, we show that the conduction properties of the channel observed experimentally are faithfully reproduced. We find that TTX occludes the entrance of NaV1.4 by forming a network of hydrogen-bonds at the outer lumen of the selectivity filter. The guanidine group of TTX adopts a lateral orientation, rather than pointing into the filter as proposed previously. The acidic residues just above the selectivity filter are important in stabilizing the hydrogen-bond network between TTX and NaV1.4. The effect of two single mutations of a critical tyrosine residue in the filter of NaV1.4 on TTX binding observed experimentally is reproduced using computational mutagenesis.

Multi-ion versus single-ion conduction mechanisms can yield current rectification in biological ion channels.

There is clear evidence that the net magnitude of negative charge at the intracellular end of inwardly rectifying potassium channels helps to generate an asymmetry in the magnitude of the current that will pass in each direction. However, a complete understanding of the physical mechanism that links these charges to current rectification has yet to be obtained. Using Brownian dynamics, we compare the conduction mechanism and binding sites in rectifying and non-rectifying channel models. We find that in our models, rectification is a consequence of asymmetry in the hydrophobicity and charge of the pore lining. As a consequence, inward conduction can occur by a multi-ion conduction mechanism. However, outward conduction is restricted, since there are fewer ions at the intracellular entrance and outwardly moving ions must cross the pore on their own. We pose the question as to whether the same mechanism could be at play in inwardly rectifying potassium channels.

Molecular dynamics simulations of scorpion toxin recognition by the Ca(2+)-activated potassium channel KCa3.1.

The Ca(2+)-activated channel of intermediate-conductance (KCa3.1) is a target for antisickling and immunosuppressant agents. Many small peptides isolated from animal venoms inhibit KCa3.1 with nanomolar affinities and are promising drug scaffolds. Although the inhibitory effect of peptide toxins on KCa3.1 has been examined extensively, the structural basis of toxin-channel recognition has not been understood in detail. Here, the binding modes of two selected scorpion toxins, charybdotoxin (ChTx) and OSK1, to human KCa3.1 are examined in atomic detail using molecular dynamics (MD) simulations. Employing a homology model of KCa3.1, we first determine conduction properties of the channel using Brownian dynamics and ascertain that the simulated results are in accord with experiment. The model structures of ChTx-KCa3.1 and OSK1-KCa3.1 complexes are then constructed using MD simulations biased with distance restraints. The ChTx-KCa3.1 complex predicted from biased MD is consistent with the crystal structure of ChTx bound to a voltage-gated K(+) channel. The dissociation constants (Kd) for the binding of both ChTx and OSK1 to KCa3.1 determined experimentally are reproduced within fivefold using potential of mean force calculations. Making use of the knowledge we gained by studying the ChTx-KCa3.1 complex, we attempt to enhance the binding affinity of the toxin by carrying out a theoretical mutagenesis. A mutant toxin, in which the positions of two amino acid residues are interchanged, exhibits a 35-fold lower Kd value for KCa3.1 than that of the wild-type. This study provides insight into the key molecular determinants for the high-affinity binding of peptide toxins to KCa3.1, and demonstrates the power of computational methods in the design of novel toxins.

Designing a C84 fullerene as a specific voltage-gated sodium channel blocker.

Fullerene derivatives demonstrate considerable potential for numerous biological applications, such as the effective inhibition of HIV protease. Recently, they were identified for their ability to indiscriminately block biological ion channels. A fullerene derivative which specifically blocks a particular ion channel could lead to a new set of drug leads for the treatment of various ion channel-related diseases. Here, we demonstrate their extraordinary potential by designing a fullerene which mimics some of the functions of µ-conotoxin, a peptide derived from cone snail venom which potently binds to the bacterial voltage-gated sodium channel (NavAb). We show, using molecular dynamics simulations, that the C84 fullerene with six lysine derivatives uniformly attached to its surface is selective to NavAb over a voltage-gated potassium channel (Kv1.3). The side chain of one of the lysine residues protrudes into the selectivity filter of the channel, while the methionine residues located just outside of the channel form hydrophobic contacts with the carbon atoms of the fullerene. The modified C84 fullerene strongly binds to the NavAb channel with an affinity of 46 nM but binds weakly to Kv1.3 with an affinity of 3 mM. This potent blocker of NavAb may serve as a structural template from which potent compounds can be designed for the targeting of mammalian Nav channels. There is a genuine need to target mammalian Nav channels as a form of treatment of various diseases which have been linked to their malfunction, such as epilepsy and chronic pain.

Complex structures between the N-type calcium channel (CaV2.2) and ω-conotoxin GVIA predicted via molecular dynamics.

The N-type voltage-gated Ca(2+) channel CaV2.2 is one of the important targets for pain management. ω-Conotoxins isolated from venoms of cone snails, which specifically inhibit CaV2.2, are promising scaffolds for novel analgesics. The inhibitory action of ω-conotoxins on CaV2.2 has been examined experimentally, but the modes of binding of the toxins to this and other related subfamilies of Ca(2+) channels are not understood in detail. Here molecular dynamics simulations are used to construct models of ω-conotoxin GVIA in complex with a homology model of the pore domain of CaV2.2. Three different binding modes in which the side chain of Lys2, Arg17, or Lys24 from the toxin protrudes into the selectivity filter of CaV2.2 are considered. In all the modes, the toxin forms a salt bridge with an aspartate residue of subunit II just above the EEEE ring of the selectivity filter. Using the umbrella sampling technique and potential of mean force calculations, the half-maximal inhibitory concentration (IC50) values are calculated to be 1.5 and 0.7 nM for the modes in which Lys2 and Arg17 occlude the ion conduction pathway, respectively. Both IC50 values compare favorably with the values of 0.04-1.0 nM determined experimentally. The similar IC50 values calculated for the different binding modes demonstrate that GVIA can inhibit CaV2.2 with alternative binding modes. Such a multiple-binding mode mechanism may be common for ω-conotoxins.

Computational methods of studying the binding of toxins from venomous animals to biological ion channels: theory and applications.

The discovery of new drugs that selectively block or modulate ion channels has great potential to provide new treatments for a host of conditions. One promising avenue revolves around modifying or mimicking certain naturally occurring ion channel modulator toxins. This strategy appears to offer the prospect of designing drugs that are both potent and specific. The use of computational modeling is crucial to this endeavor, as it has the potential to provide lower cost alternatives for exploring the effects of new compounds on ion channels. In addition, computational modeling can provide structural information and theoretical understanding that is not easily derivable from experimental results. In this review, we look at the theory and computational methods that are applicable to the study of ion channel modulators. The first section provides an introduction to various theoretical concepts, including force-fields and the statistical mechanics of binding. We then look at various computational techniques available to the researcher, including molecular dynamics, brownian dynamics, and molecular docking systems. The latter section of the review explores applications of these techniques, concentrating on pore blocker and gating modifier toxins of potassium and sodium channels. After first discussing the structural features of these channels, and their modes of block, we provide an in-depth review of past computational work that has been carried out. Finally, we discuss prospects for future developments in the field.

Effect of gating modifier toxins on membrane thickness: implications for toxin effect on gramicidin and mechanosensitive channels.

Various gating modifier toxins partition into membranes and interfere with the gating mechanisms of biological ion channels. For example, GsMTx4 potentiates gramicidin and several bacterial mechanosensitive channels whose gating kinetics are sensitive to mechanical properties of the membrane, whereas binding of HpTx2 shifts the voltage-activity curve of the voltage-gated potassium channel Kv4.2 to the right. The detailed process by which the toxin partitions into membranes has been difficult to probe using molecular dynamics due to the limited time scale accessible. Here we develop a protocol that allows the spontaneous assembly of a polypeptide toxin into membranes in atomistic molecular dynamics simulations of tens of nanoseconds. The protocol is applied to GsMTx4 and HpTx2. Both toxins, released in water at the start of the simulation, spontaneously bind into the lipid bilayer within 50 ns, with their hydrophobic patch penetrated into the bilayer beyond the phosphate groups of the lipids. It is found that the bilayer is about 2 Å thinner upon the binding of a GsMTx4 monomer. Such a thinning effect of GsMTx4 on membranes may explain its potentiation effect on gramicidin and mechanosensitive channels.

Conduction and block of inward rectifier K+ channels: predicted structure of a potent blocker of Kir2.1.

Dysfunction of Kir2.1, thought to be the major component of inward currents, I(K1), in the heart, has been linked to various channelopathies, such as short Q-T syndrome. Unfortunately, currently no known blockers of Kir2.x channels exist. In contrast, Kir1.1b, predominantly expressed in the kidney, is potently blocked by an oxidation-resistant mutant of the honey bee toxin tertiapin (tertiapin-Q). Using various computational tools, we show that both channels are closed by a hydrophobic gating mechanism and inward rectification occurs in the absence of divalent cations and polyamines. We then demonstrate that tertiapin-Q binds to the external vestibule of Kir1.1b and Kir2.1 with K(d) values of 11.6 nM and 131 µM, respectively. We find that a single mutation of tertiapin-Q increases the binding affinity for Kir2.1 by 5 orders of magnitude (K(d) = 0.7 nM). This potent blocker of Kir2.1 may serve as a structural template from which potent compounds for the treatment of various diseases mediated by this channel subfamily, such as cardiac arrhythmia, can be developed.