Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.     

Muscle electrolyte measurements during and after hypokinesia in determining muscle electrolyte depletion during hypokinesia in the rat

Hypokinesia (diminished movement) induces muscle mineral depletion. However, the mechanism of muscle mineral depletion during hypokinesia (HK) remains unknown. Measuring electrolyte retention and electrolyte values in muscle, plasma, and urine during and after HK, the aim of this study was to discover if HK could depress mineral retention and lead to muscle mineral depletion. Studies were done on 204 13-wk-old male Wistar rats (370–390 g) during 10 d pre-HK period, 98 d HK period, and 15 d post-HK period. Rats were equally divided into two groups: vivarium control rats (VCR) and hypokinetic rats (HKR). All hypokinetic rats were kept for 98 d in small individual cages, which restricted their movements in all directions without hindering food and water intakes. All control rats were housed for 98 d in individual cages under vivarium control conditions. Both groups of rats were pair-fed. During the HK period skeletal muscle sodium (Na), potassium (K), magnesium (Mg), calcium (Ca), and water content and electrolyte retention decreased significantly (p < 0.05), while urinary and plasma electrolyte levels increased significantly (p < 0.05) in HKR compared with their pre-HK values and their respective VCR. During the initial days of the post-HK period, mineral retention increased significantly (p < 0.05), plasma and urinary electrolyte level decreased significantly (p < 0.05), while muscle electrolyte and water content remained significantly (p < 0.05) depressed in HKR compared with VCR. Muscle mineral and water content, electrolyte retention, plasma, and urinary electrolyte values did not change in VCR compared with their pre-HK values. It was concluded that during HK decreased muscle mineral content may suggest muscle mineral depletion, while increased urinary electrolyte loss and muscle mineral depletion may demonstrate reduced mineral retention. Reduced electrolyte excretion and depressed muscle mineral content during post-HK may indicate skeletal muscle mineral depletion during HK. Dissociation between electrolyte retention and muscle mineral depletion may demonstrate the presence of decreased electrolyte retention as the mechanism of muscle electrolyte depletion during prolonged HK.     

Ions and blockers in potassium channels: insights from free energy simulations

Potassium ion channels enable efficient and selective permeation of K+ ions across nonpolar biological membranes. Here we review the results of recent free energy calculations related to the permeation of monovalent cations through K+ channels and to the channel inhibition by blocker compounds. In particular, the progress in computational studies of the bacterial KcsA channel is discussed.     

Computational study of the binding affinity and selectivity of the bacterial ammonium transporter AmtB

We report results from microscopic molecular dynamics and free energy perturbation simulations of substrate binding and selectivity for the Escherichia coli high-affinity ammonium transporter AmtB. The simulation system consists of the protein embedded in a model membrane/water surrounding. The calculated absolute binding free energies for the external NH(4)(+) ions are between -5.8 and -7.3 kcal/mol and are in close agreement with experimental data. The apparent pK(a) of the bound NH(4)(+) increases by more than 4 units, indicating a preference for binding ammonium ion and not neutral ammonia. The external binding site is also selective for NH(4)(+) toward monovalent metal cations by 2.4-4.4 kcal/mol. The externally bound NH(4)(+) shows strong electrostatic interactions with the proximal buried Asp160, stabilized in the anionic form, whereas the interactions with the aromatic rings of Phe107 and Trp148, lining the binding cavity, are less pronounced. Simulated mutation of the highly conserved Asp160 to Asn reduces the pK(a) of the bound ammonium ion by approximately 7 units and causes loss of its binding. The calculations further predict that the substrate affinity of E. coli AmtB depends on the ionization state of external histidines. The computed free energies of hypothetical intermediate states related to transfer of NH(3), NH(4)(+), or H(2)O from the external binding site to the first position inside the internal channel pore favor permeation of the neutral species through the channel interior. However, the predicted change in the apparent pK(a) of NH(4)(+) upon translocation from the external site, Am1, to the first internal site, Am2, indicates that ammonium ion becomes deprotonated only when it enters the channel interior.     

Structure-activity relationship for extracellular block of K+ channels by tetraalkylammonium ions

External tetraalkylammonium ion binding to potassium channels is studied using microscopic molecular modelling methods and the experimental structure of the KcsA channel. Relative binding free energies of the KcsA complexes with Me4N+, Et4N+, and n-Pr4N+ are calculated with the molecular dynamics free energy perturbation approach together with automated ligand docking. The four-fold symmetry of the entrance cavity formed by the Tyr82 residues is found to provide stronger binding for the D2d than for the S4 conformation of the ligands. In agreement with experiment the Et4N+ blocker shows several kcal/mol better binding than the other tetraalkylammonium ions.     

Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping

Background

Natural selection has molded evolution across all taxa. At an arguable date of around 330,000 years ago there were already at least two different types of cattle that became ancestors of nearly all modern cattle, the Bos taurus taurus more adapted to temperate climates and the tropically adapted Bos taurus indicus. After domestication, human selection exponentially intensified these differences. To better understand the genetic differences between these subspecies and detect genomic regions potentially under divergent selection, animals from the International Bovine HapMap Experiment were genotyped for over 770,000 SNP across the genome and compared using smoothed F_ST. The taurine sample was represented by ten breeds and the contrasting zebu cohort by three breeds.

Results

Each cattle group evidenced similar numbers of polymorphic markers well distributed across the genome. Principal components analyses and unsupervised clustering confirmed the well-characterized main division of domestic cattle. The top 1% smoothed F_ST, potentially associated to positive selection, contained 48 genomic regions across 17 chromosomes. Nearly half of the top F_ST signals (n = 22) were previously detected using a lower density SNP assay. Amongst the strongest signals were the BTA7:~50 Mb and BTA14:~25 Mb; both regions harboring candidate genes and different patterns of linkage disequilibrium that potentially represent intrinsic differences between cattle types. The bottom 1% of the smoothed F_ST values, potentially associated to balancing selection, included 24 regions across 13 chromosomes. These regions often overlap with copy number variants, including the highly variable region at BTA23:~24 Mb that harbors a large number of MHC genes. Under these regions, 318 unique Ensembl genes are annotated with a significant overrepresentation of immune related pathways.

Conclusions

Genomic regions that are potentially linked to purifying or balancing selection processes in domestic cattle were identified. These regions are of particular interest to understand the natural and human selective pressures to which these subspecies were exposed to and how the genetic background of these populations evolved in response to environmental challenges and human manipulation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-876) contains supplementary material, which is available to authorized users.     

Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle

Background

Signatures of selection are regions in the genome that have been preferentially increased in frequency and fixed in a population because of their functional importance in specific processes. These regions can be detected because of their lower genetic variability and specific regional linkage disequilibrium (LD) patterns.

Methods

By comparing the differences in regional LD variation between dairy and beef cattle types, and between indicine and taurine subspecies, we aim at finding signatures of selection for production and adaptation in cattle breeds. The VarLD method was applied to compare the LD variation in the autosomal genome between breeds, including Angus and Brown Swiss, representing taurine breeds, and Nelore and Gir, representing indicine breeds. Genomic regions containing the top 0.01 and 0.1 percentile of signals were characterized using the UMD3.1 Bos taurus genome assembly to identify genes in those regions and compared with previously reported selection signatures and regions with copy number variation.

Results

For all comparisons, the top 0.01 and 0.1 percentile included 26 and 165 signals and 17 and 125 genes, respectively, including TECRL, BT.23182 or FPPS, CAST, MYOM1, UVRAG and DNAJA1.

Conclusions

The VarLD method is a powerful tool to identify differences in linkage disequilibrium between cattle populations and putative signatures of selection with potential adaptive and productive importance.     

Structural and Biochemical Analysis of Human Pathogenic Astrovirus Serine Protease at 2.0 Å Resolution

Astroviruses are single-stranded RNA viruses with a replication strategy based on the proteolytic processing of a polyprotein precursor and subsequent release of the viral enzymes of replication. So far, the catalytic properties of the astrovirus protease as well as its structure have remained uncharacterized. In this study, the three-dimensional crystal structure of the predicted protease of human pathogenic astrovirus has been solved to 2.0 Å resolution. The protein displays the typical properties of trypsin-like enzymes but also several characteristic features: (i) a catalytic Asp-His-Ser triad in which the aspartate side chain is oriented away from the histidine, being replaced by a water molecule; (ii) a non-common conformation and composition of the S1 pocket; and (iii) the lack of the typical surface β-ribbons together with a “featureless” shape of the substrate-binding site. Hydrolytic activity assays indicate that the S1 pocket recognises Glu and Asp side chains specifically, which, therefore, are predicted to occupy the P1 position on the substrate cleavage site. The positive electrostatic potential featured by the S1 region underlies this specificity. The comparative structural analysis highlights the peculiarity of the astrovirus protease, and differentiates it from the human and viral serine proteases.     

Impact of temperature on Na2Sr(PO4)F:Eu3+ phosphor

In this article, we report the synthesis of Na₂Sr_1‐x(PO₄)F:Eu_x phosphor via a combustion method. The influence of different annealing temperatures on the photoluminescence properties was investigated. The phosphor was excited at both 254 and 393 nm. Na₂Sr_1‐x(PO₄)F:Eu_x³⁺ phosphors emit strong orange and red color at 593 and 612 nm, respectively, under both excitation wavelengths. Na₂Sr_1‐x(PO₄)F:Eu_x³⁺ phosphors annealed at 1050°C showed stronger emission intensity compared with 600, 900 and 1200°C. Moreover, Na₂Sr_1‐x(PO₄)F:Eu_x³⁺ phosphor was found to be more intense when compared with commercial Y₂O₃:Eu³⁺ phosphor. Copyright © 2013 John Wiley & Sons, Ltd.     

K(+)/Na(+) selectivity of the KcsA potassium channel from microscopic free energy perturbation calculations

Microscopic molecular dynamics free energy perturbation calculations of the K(+)/Na(+) selectivity in the KcsA potassium channel, based on its experimental three-dimensional structure, are reported. The relative binding free energies for K(+) and Na(+) in the most relevant ion occupancy states of the four-site selectivity filter are calculated. The previously proposed mechanism for ion permeation through the KcsA channel is predicted, in agreement with available experimental data, to have a significant selectivity for K(+) over Na(+). The calculations also show that the individual 'binding site' selectivities are generally not additive and the doubly loaded states of the filter thus display cooperative effects. The only site that is not K(+) selective is that which is located at the entrance to the internal water cavity, suggesting the possibility that internal Na(+) could block outward currents.