Kv7.2 (KCNQ2) is the principal molecular component of the slow voltage gated M-channel, which strongly influences neuronal excitability. Calmodulin (CaM) binds to two intracellular C-terminal segments of Kv7.2 channels, helices A and B, and it is required for exit from the endoplasmic reticulum. However, the molecular mechanisms by which CaM controls channel trafficking are currently unknown. Here we used two complementary approaches to explore the molecular events underlying the association between CaM and Kv7.2 and their regulation by Ca(2+). First, we performed a fluorometric assay using dansylated calmodulin (D-CaM) to characterize the interaction of its individual lobes to the Kv7.2 CaM binding site (Q2AB). Second, we explored the association of Q2AB with CaM by NMR spectroscopy, using (15)N-labeled CaM as a reporter. The combined data highlight the interdependency of the N- and C-lobes of CaM in the interaction with Q2AB, suggesting that when CaM binds Ca(2+) the binding interface pivots between the N-lobe whose interactions are dominated by helix B and the C-lobe where the predominant interaction is with helix A. In addition, Ca(2+) makes CaM binding to Q2AB more difficult and, reciprocally, the channel weakens the association of CaM with Ca(2+).

Bacterial outer membrane lipopolysaccharide (LPS) potently stimulates the mammalian innate immune system, and can lead to sepsis, the primary cause of death from infections. LPS is sensed by Toll-like receptor 4 (TLR4) in complex with its lipid-binding coreceptor MD-2, but subtle structural variations in LPS can profoundly modulate the response. To better understand the mechanism of LPS-induced stimulation and bacterial evasion, we have calculated the binding affinity to MD-2 of agonistic and antagonistic LPS variants including lipid A, lipid IVa, and synthetic antagonist Eritoran, and provide evidence that the coreceptor is a molecular switch that undergoes ligand-induced conformational changes to appropriately activate or inhibit the receptor complex. The plasticity of the coreceptor binding cavity is shown to be essential for distinguishing between ligands, whilst similar calculations for a model bacterial LPS bilayer reveal the "membrane-like" nature of the protein cavity. The ability to predict the activity of LPS variants should facilitate the rational design of TLR4 therapeutics.

The bacterial cell envelope is composed of a mixture of different lipids and proteins, making it an inherently complex organelle. The interactions between integral membrane proteins and lipids are crucial for their respective spatial localization within bacterial cells. We have employed microsecond timescale coarse-grained molecular dynamics simulations of vesicles of varying sizes and with a range of protein and lipid compositions, and used novel approaches to measure both local and global system dynamics, the latter based on spherical harmonics analysis. Our results suggest that both hydrophobic mismatch, enhanced by embedded membrane proteins, and curvature based sorting, due to different modes of undulation, may drive assembly in vesicular systems. Interestingly, the modes of undulation of the vesicles were found to be altered by the specific protein and lipid composition of the vesicle. Strikingly, lipid dynamics were shown to be coupled to proteins up to 6 nm from their surface, a substantially larger distance than has previously been observed, resulting in multi-layered annular rings enriched with particular types of phospholipid. Such large protein-lipid complexes may provide a mechanism for long-range communication. Given the complexity of bacterial membranes, our results suggest that subtle changes in lipid composition may have major implications for lipid and protein sorting under a curvature-based membrane-sorting model.

Metabolism of xenobiotics remains a central challenge for the discovery and development of drugs, cosmetics, nutritional supplements, and agrochemicals. Metabolic transformations are frequently related to the incidence of toxic effects that may result from the emergence of reactive species, the systemic accumulation of metabolites, or by induction of metabolic pathways. Experimental investigation of the metabolism of small organic molecules is particularly resource demanding; hence, computational methods are of considerable interest to complement experimental approaches. This review provides a broad overview of structure- and ligand-based computational methods for the prediction of xenobiotic metabolism. Current computational approaches to address xenobiotic metabolism are discussed from three major perspectives: (i) prediction of sites of metabolism (SOMs), (ii) elucidation of potential metabolites and their chemical structures, and (iii) prediction of direct and indirect effects of xenobiotics on metabolizing enzymes, where the focus is on the cytochrome P450 (CYP) superfamily of enzymes, the cardinal xenobiotics metabolizing enzymes. For each of these domains, a variety of approaches and their applications are systematically reviewed, including expert systems, data mining approaches, quantitative structure-activity relationships (QSARs), and machine learning-based methods, pharmacophore-based algorithms, shape-focused techniques, molecular interaction fields (MIFs), reactivity-focused techniques, protein-ligand docking, molecular dynamics (MD) simulations, and combinations of methods. Predictive metabolism is a developing area, and there is still enormous potential for improvement. However, it is clear that the combination of rapidly increasing amounts of available ligand- and structure-related experimental data (in particular, quantitative data) with novel and diverse simulation and modeling approaches is accelerating the development of effective tools for prediction of in vivo metabolism, which is reflected by the diverse and comprehensive data sources and methods for metabolism prediction reviewed here. This review attempts to survey the range and scope of computational methods applied to metabolism prediction and also to compare and contrast their applicability and performance.

Voltage-gated Na⁺-channels are transmembrane proteins that are responsible for the fast depolarizing phase of the action potential in nerve and muscular cells. Selective permeability of Na⁺ over Ca²⁺ or K⁺ ions is essential for the biological function of Na⁺-channels. After the emergence of the first high-resolution structure of a Na⁺-channel, an anionic coordination site was proposed to confer Na⁺ selectivity through partial dehydration of Na⁺ via its direct interaction with conserved glutamate side chains. By combining molecular dynamics simulations and free-energy calculations, a low-energy permeation pathway for Na⁺ ion translocation through the selectivity filter of the recently determined crystal structure of a prokaryotic sodium channel from Arcobacter butzleri is characterised. The picture that emerges is that of a pore preferentially occupied by two ions, which can switch between different configurations by crossing low free-energy barriers. In contrast to K⁺-channels, the movements of the ions appear to be weakly coupled in Na⁺-channels. When the free-energy maps for Na⁺ and K⁺ ions are compared, a selective site is characterised in the narrowest region of the filter, where a hydrated Na⁺ ion, and not a hydrated K⁺ ion, is energetically stable.

The ability of DENV2 to display different morphologies (hence different antigenic properties) complicates vaccine and therapeutics development. Previous studies showed most strains of laboratory adapted DENV2 particles changed from smooth to "bumpy" surfaced morphology when the temperature is switched from 29°C at 37°C. Here we identified five envelope (E) protein residues different between two alternative passage history DENV2 NGC strains exhibiting smooth or bumpy surface morphologies. Several mutations performed on the smooth DENV2 infectious clone destabilized the surface, as observed by cryoEM. Molecular dynamics simulations demonstrated how chemically subtle substitution at various positions destabilized dimeric interactions between E proteins. In contrast, three out of four DENV2 clinical isolates showed a smooth surface morphology at 37°C, and only at high fever temperature (40°C) did they become "bumpy". These results imply vaccines should contain particles representing both morphologies. For prophylactic and therapeutic treatments, this study also informs on which types of antibodies should be used at different stages of an infection, i.e., those that bind to monomeric E proteins on the bumpy surface or across multiple E proteins on the smooth surfaced virus.

The trimeric periplasmic holdase chaperone Skp binds and stabilizes unfolded outer membrane proteins (OMPs) as part of bacterial OMP biogenesis. Skp binds client proteins in its central cavity, thereby reducing its backbone dynamics, but the molecular mechanisms that govern Skp dynamics and adaptation to differently sized clients remains unknown. Here, we employ a combination of microsecond timescale molecular dynamics simulation, small-angle X-ray scattering, and nuclear magnetic resonance spectroscopy to reveal that Skp is remarkably flexible, and features a molecular spring-loaded mechanism in its "tentacle" arms that enables switching between two distinct conformations on sub-millisecond timescales. The conformational switch is executed around a conserved pivot element within the coiled-coil structures of the tentacles, allowing expansion of the cavity and thus accommodation of differently sized clients. The spring-loaded mechanism shows how a chaperone can efficiently modulate its structure and function in an ATP-independent manner.

The ongoing development of drug resistance in HIV continues to push for the need of alternative drug targets in inhibiting HIV. One such target is the Reverse transcriptase (RT) enzyme which is unique and critical in the viral life cycle-a rational target that is likely to have less off-target effects in humans. Serendipitously, we found two chemical scaffolds from the National Cancer Institute (NCI) Diversity Set V that inhibited HIV-1 RT catalytic activity. Computational structural analyses and subsequent experimental testing demonstrated that one of the two chemical scaffolds binds to a novel location in the HIV-1 RT p51 subunit, interacting with residue Y183, which has no known association with previously reported drug resistance. This finding supports the possibility of a novel druggable site on p51 for a new class of non-nucleoside RT inhibitors that may inhibit HIV-1 RT allosterically. Although inhibitory activity was shown experimentally to only be in the micromolar range, the scaffolds serve as a proof-of-concept of targeting the HIV RT p51 subunit, with the possibility of medical chemistry methods being applied to improve inhibitory activity towards more effective drugs.

Fast C-type inactivation confers distinctive functional properties to the hERG potassium channel, and its association to inherited and acquired cardiac arrythmias makes the study of the inactivation mechanism of hERG at the atomic detail of paramount importance. At present, two models have been proposed to describe C-type inactivation in K+-channels. Experimental data and computational work on the bacterial KcsA channel support the hypothesis that C-type inactivation results from a closure of the selectivity filter that sterically impedes ion conduction. Alternatively, recent experimental structures of a mutated Shaker channel revealed a widening of the extracellular portion of the selectivity filter, which might diminish conductance by interfering with the mechanism of ion permeation. Here, we performed molecular dynamics simulations of the wild-type hERG, a non-inactivating mutant (hERG-N629D), and a mutant that inactivates faster than the wild-type channel (hERG-F627Y) to find out which and if any of the two reported C-type inactivation mechanisms applies to hERG. Closure events of the selectivity filter were not observed in any of the simulated trajectories but instead, the extracellular section of the selectivity filter deviated from the canonical conductive structure of potassium channels. The degree of widening of the potassium binding sites at the extracellular entrance of the channel was directly related to the degree of inactivation with hERG-F627Y > wild-type hERG > hERG-N629D. These findings support the hypothesis that C-type inactivation in hERG entails a widening of the extracellular entrance of the channel rather than a closure of the selectivity filter.

Transient receptor potential (TRP) ion channels are important pharmacological targets because of their role in the perception of pain, and so, understanding their chemical regulation is essential for the development of analgesic drugs. Among the currently known TRP channel chemical agonists, capsaicin, the active compound of chili pepper, is probably the most exhaustively studied. The availability of the three-dimensional structure of the vanilloid receptor 1 (TRPV1) has fueled computational studies revealing the molecular details of capsaicin binding modes. Although this is a significant step, a comprehensible binding mechanism or pathway is invaluable for targeting TRP channels in modern pharmacology. In the present work, free-energy and enhanced sampling techniques have been used to explore a possible membrane-mediated pathway for capsaicin to enter the TRPV1 binding pocket where capsaicin accesses the protein starting at the extracellular milieu through the outer leaflet and into its binding site in the protein. The main states visited along this route have been characterized and include (i) a bound state in agreement with the binding mode "head-down, tail-up" and (ii) an alternative state corresponding to a "head-up, tail-down" binding mode. In agreement with previous reports, binding is mediated by both hydrogen bonds and van der Waals interactions, and residue Y511 is crucial for stabilizing the bound state and during the binding process. Together, these results provide a foundation to further understand TRPV channels, and they could be used to guide therapeutic design of selective inhibitors potentially leading to novel avenues for pharmacological applications targeting the TRPV1 channel.

The bacterial outer membrane protein OmpA is composed of an N-terminal 171-residue beta-barrel domain (OmpA(171)) that spans the bilayer and a periplasmic, C-terminal domain of unknown structure. OmpA has been suggested to primarily serve a structural role, as no continuous pore through the center of the barrel can be discerned in the crystal structure of OmpA(171). However, several groups have recorded ionic conductances for bilayer-reconstituted OmpA(171). To resolve this apparent paradox we have used molecular dynamics (MD) simulations on OmpA(171) to explore the conformational dynamics of the protein, in particular the possibility of transient formation of a central pore. A total of 19 ns of MD simulations of OmpA(171) have been run, and the results were analyzed in terms of 1) comparative behavior of OmpA(171) in different bilayer and bilayer-mimetic environments, 2) solvation states of OmpA(171), and 3) pore characteristics in different MD simulations. Significant mobility was observed for residues and water molecules within the beta-barrel. A simulation in which putative gate region side chains of the barrel interior were held in a non-native conformation led to an open pore, with a predicted conductance similar to experimental measurements. The OmpA(171) pore has been shown to be somewhat more dynamic than suggested by the crystal structure. A gating mechanism is proposed to explain its documented channel properties, involving a flickering isomerization of Arg138, forming alternate salt bridges with Glu52 (closed state) and Glu128 (open state).

Non-coding RNAs are crucial regulators for a vast array of cellular processes and have been implicated in human disease. These biological processes represent a hitherto untapped resource in our fight against disease. In this work we identify small molecule inhibitors of a non-coding RNA uridylylation pathway. The TUTase family of enzymes is important for modulating non-coding RNA pathways in both human cancer and pathogen systems. We demonstrate that this new class of drug target can be accessed with traditional drug discovery techniques. Using the Trypanosoma brucei TUTase, RET1, we identify TUTase inhibitors and lay the groundwork for the use of this new target class as a therapeutic opportunity for the under-served disease area of African Trypanosomiasis. In a broader sense this work demonstrates the therapeutic potential for targeting RNA post-transcriptional modifications with small molecules in human disease.

Antimicrobial peptides are small, cationic proteins that can induce lysis of bacterial cells through interaction with their membranes. Different mechanisms for cell lysis have been proposed, but these models tend to neglect the role of the chemical composition of the membrane, which differs between bacterial species and can be heterogeneous even within a single cell. Moreover, the cell envelope of Gram-negative bacteria such as E. coli contains two membranes with differing compositions. To this end, we report the first molecular dynamics simulation study of the interaction of the antimicrobial peptide, polymyxin B1 with complex models of both the inner and outer membranes of E. coli. The results of >16 microseconds of simulation predict that polymyxin B1 is likely to interact with the membranes via distinct mechanisms. The lipopeptides aggregate in the lipopolysaccharide headgroup region of the outer membrane with limited tendency for insertion within the lipid A tails. In contrast, the lipopeptides readily insert into the inner membrane core, and the concomitant increased hydration may be responsible for bilayer destabilization and antimicrobial function. Given the urgent need to develop novel, potent antibiotics, the results presented here reveal key mechanistic details that may be exploited for future rational drug development.

We have designed, synthesised and characterised the effects of a number of novel anthraquinone derivatives and assessed their effects on large conductance, Ca(2+) activated K(+) (BK) channels recorded from rabbit bladder smooth muscle cells using the excised, inside/out configuration of the patch clamp technique. These compounds are members of the GoSlo-SR family of compounds, which potently open BK channels and shift the voltage required for half maximal activation (V1/2) negatively. The efficacy of the anilinoanthraquinone derivatives was enhanced when the size of ring D was increased, since the cyclopentane and cyclohexane derivatives shifted the V1/2, by -24 ± 6 mV and -54 ± 8 mV, respectively, whereas the cycloheptane and cyclooctane derivatives shifted the V1/2 by -61 ± 6 mV and -106 ± 6 mV. To examine if a combination of hydrophobicity and steric bulking of this region further enhanced their ability to open BK channels, we synthesised a number of naphthalene and tetrahydro-naphthalene derivatives. The tetrahydro-2-naphthalene derivative GoSlo-SR-5-69 was the most potent and efficacious of the series since it was able to shift the activation V1/2 by greater than -100 mV when applied at a concentration of 1 μM and had an EC50 of 251 nM, making it one of the most potent and efficacious BK channel openers synthesised to date.

Tetrameric coiled-coil structures are present in many ion channels, often adjacent to a calmodulin (CaM) binding site, although the relationship between the two is not completely understood. Here we examine the dynamic properties of the ABCD domain located in the intracellular C-terminus of tetrameric, voltage-dependent, potassium selective Kv7.2 channels. This domain encompasses the CaM binding site formed by helices A and B, followed by helix C, which is linked to the helix D coiled-coil. The data reveals that helix D stabilizes CaM binding, promoting trans-binding (CaM embracing neighboring subunits), and they suggest that the ABCD domain can be exchanged between subunits of the tetramer. Exchange is faster when mutations in AB weaken the CaM interaction. The exchange of ABCD domains is slower in the presence of Ca2+, indicating that CaM stabilization of the tetrameric assembly is enhanced when loaded with this cation. Our observations are consistent with a model that involves a dynamic mechanism of helix D assembly, which supports reciprocal allosteric coupling between the A-B module and the coiled-coil formed by the helix D. Thus, formation of the distal helix D tetramer influences CaM binding and CaM-dependent Kv7.2 properties, whereas reciprocally, CaM and Ca2+ influence the dynamic behavior of the helix D coiled-coil.

Thrombin-derived C-terminal peptides (TCPs) of about 2 kDa are present in wounds, where they exert anti-endotoxic functions. Employing a combination of nuclear magnetic resonance spectroscopy (NMR), biophysical, mass spectrometry and cellular studies combined with in silico multiscale modelling, we here determine the bound conformation of HVF18 (HVFRLKKWIQKVIDQFGE), a TCP generated by neutrophil elastase, in complex with bacterial lipopolysaccharide (LPS) and define a previously undisclosed interaction between TCPs and human CD14. Further, we show that TCPs bind to the LPS-binding hydrophobic pocket of CD14 and identify the peptide region crucial for TCP interaction with LPS and CD14. Taken together, our results demonstrate the role of structural transitions in LPS complex formation and CD14 interaction, providing a molecular explanation for the previously observed therapeutic effects of TCPs in experimental models of bacterial sepsis and endotoxin shock.

The lactose repressor protein may bind DNA in two possible configurations: a specific one, if the DNA sequence corresponds to a binding site, and a non-specific one otherwise. To find its target sequences, the lactose repressor first binds non-specifically to DNA, and subsequently, it rapidly searches for a binding site. Atomic structures of non-specific and specific complexes are available from crystallographic and nuclear magnetic resonance experiments. However, what remains unknown is a detailed description of the steps that transform the non-specific complex into the specific one. Here, how the protein first recognizes its binding site has been studied using molecular dynamics simulations. The picture that emerges is that of a protein that is as mobile when interacting with non-specific DNA sequences as when free in solution. This high degree of mobility allows the protein to rapidly sample different DNA sequences. In contrast, when the protein encounters a binding site, the configuration ensemble collapses, and the protein sliding movements along the DNA sequence become scarce. The binding energies in the specific and non-specific complexes were analysed using the Molecular Mechanics Poisson Boltzmann Surface Area approach. These results represent a first step towards a throughout characterization of the DNA-recognition process.

Potassium channels are a diverse family of integral membrane proteins through which K(+) can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K(+) channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K(+) concentration is thought to be representative of a K(+) channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K(+) KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K(+) KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K(+) structure of KcsA and KirBac. The low-K(+) KcsA structure is not permeable by Na(+), K(+), or Rb(+), and the selectivity of its binding sites is different from that of the high-K(+) structure.

Interactions of Na(+), K(+), Rb(+), and Cs(+) ions within the selectivity filter of a potassium channel have been investigated via multiple molecular dynamics simulations (total simulation time, 48 ns) based on the high resolution structure of KcsA, embedded in a phospholipid bilayer. As in simulations based on a lower resolution structure of KcsA, concerted motions of ions and water within the filter are seen. Despite the use of a higher resolution structure and the inclusion of four buried water molecules thought to stabilize the filter, this region exhibits a significant degree of flexibility. In particular, pronounced distortion of filter occurs if no ions are present within it. The two most readily permeant ions, K(+) and Rb(+), are similar in their interactions with the selectivity filter. In contrast, Na(+) ions tend to distort the filter by binding to a ring of four carbonyl oxygens. The larger Cs(+) ions result in a small degree of expansion of the filter relative to the x-ray structure. Cs(+) ions also appear to interact differently with the gate region of the channel, showing some tendency to bind within a predominantly hydrophobic pocket. The four water molecules buried between the back of the selectivity filter and the remainder of the protein show comparable mobility to the surrounding protein and do not exchange with water molecules within the filter or the central cavity. A preliminary comparison of the use of particle mesh Ewald versus cutoff protocols for the treatment of long-range electrostatics suggests some difference in the kinetics of ion translocation within the filter.

Maculatin 1.1 (M1.1) is a membrane-active antimicrobial peptide (AMP) from an Australian tree frog that forms a kinked amphipathic alpha-helix in the presence of a lipid bilayer or bilayer-mimetic environment. To help elucidate its mechanism of membrane-lytic activity, we performed a total of approximately 8 micros of coarse-grained molecular dynamics (CG-MD) simulations of M1.1 in the presence of zwitterionic phospholipid membranes. Several systems were simulated in which the peptide/lipid ratio was varied. At a low peptide/lipid ratio, M1.1 adopted a kinked, membrane-interfacial location, consistent with experiment. At higher peptide/lipid ratios, we observed spontaneous, cooperative membrane insertion of M1.1 peptide aggregates. The minimum size for formation of a transmembrane (TM) aggregate was just four peptides. The absence of a simple and well-defined central channel, along with the exclusion of lipid headgroups from the aggregates, suggests that a pore-like model is an unlikely explanation for the mechanism of membrane lysis by M1.1. We also performed an extended 1.25 micros simulation of the permeabilization of a complete liposome by multiple peptides. Consistent with the simpler bilayer simulations, formation of monomeric interfacial peptides and TM peptide clusters was observed. In contrast, major structural changes were observed in the vesicle membrane, implicating induced membrane curvature in the mechanism of active antimicrobial peptide lysis. This contrasted with the behavior of the nonpore-forming model peptide WALP23, which inserted into the vesicle to form extended clusters of TM alpha-helices with relatively little perturbation of bilayer properties.