Dysregulation of transition metal ion homeostasis is the molecular basis for cadmium toxicity in Streptococcus pneumoniae

Cadmium is a transition metal ion that is highly toxic in biological systems. Although relatively rare in the Earth’s crust, anthropogenic release of cadmium since industrialization has increased biogeochemical cycling and the abundance of the ion in the biosphere. Despite this, the molecular basis of its toxicity remains unclear. Here we combine metal-accumulation assays, high-resolution structural data and biochemical analyses to show that cadmium toxicity, in Streptococcus pneumoniae, occurs via perturbation of first row transition metal ion homeostasis. We show that cadmium uptake reduces the millimolar cellular accumulation of manganese and zinc, and thereby increases sensitivity to oxidative stress. Despite this, high cellular concentrations of cadmium (~17 mM) are tolerated, with negligible impact on growth or sensitivity to oxidative stress, when manganese and glutathione are abundant. Collectively, this work provides insight into the molecular basis of cadmium toxicity in prokaryotes, and the connection between cadmium accumulation and oxidative stress.

The molecular basis for the high toxicity of cadmium is unclear. Here, Begg et al. use the bacterium Streptococcus pneumoniae as a model system, and show that cadmium uptake increases sensitivity to oxidative stress by reducing intracellular concentrations of manganese and zinc through different mechanisms.

Mechanism of activation of protein kinase JAK2 by the growth hormone receptor

Signaling from JAK (Janus kinase) protein kinases to STAT (signal transducers and activators of transcription) transcription factors is key to many aspects of biology and medicine, yet the mechanism by which cytokine receptors initiate signaling is enigmatic. We present a complete mechanistic model for activation of receptor-bound JAK2, based on an archetypal cytokine receptor, the growth hormone receptor. For this, we used fluorescence resonance energy transfer to monitor positioning of the JAK2 binding motif in the receptor dimer, substitution of the receptor extracellular domains with Jun zippers to control the position of its transmembrane (TM) helices, atomistic modeling of TM helix movements, and docking of the crystal structures of the JAK2 kinase and its inhibitory pseudokinase domain with an opposing kinase-pseudokinase domain pair. Activation of the receptor dimer induced a separation of its JAK2 binding motifs, driven by a ligand-induced transition from a parallel TM helix pair to a left-handed crossover arrangement. This separation leads to removal of the pseudokinase domain from the kinase domain of the partner JAK2 and pairing of the two kinase domains, facilitating trans-activation. This model may well generalize to other class I cytokine receptors.

The central cavity of ABCB1 undergoes alternating access during ATP hydrolysis

Understanding the process that underlies multidrug recognition and efflux by P-glycoprotein (ABCB1) remains a key biological challenge. Structural data have recently become available for the murine and Caenorhabditis elegans homologues of ABCB1; however all structures were obtained in the absence of nucleotide. A feature of these structures was the presence of a central cavity that is inaccessible from the extracellular face of the protein. To determine the conformational dynamics of this region several residues in transmembrane helices TM6 (331, 343 and 354) and TM12 (980) were mutated to cysteine. Based upon structural predictions, these residues are proposed to line, or reside proximal to, the central cavity. The mutant isoforms were labelled with a paramagnetic probe enabling the application of EPR spectroscopic methods. Power saturation EPR spectra were recorded in the presence of hydrophobic (O2 ) or hydrophilic (NiEDDA) quenching agents to study the local environment of each residue. ABCB1 was trapped in both its nucleotide-bound and post-hydrolytic conformations and EPR spectra were again recorded in the presence and absence of quenching agents. The EPR line shapes provide information on the movements of these residues within TM6 and TM12 during ATP hydrolysis. Rationalization of the data with molecular dynamic simulations indicates that the cavity is converted to a configuration open to the aqueous phase following nucleotide binding, thereby suggesting alternating access to the cavity on opposite sides of the membrane during translocation.

Vancomycin: ligand recognition, dimerization and super-complex formation

The antibiotic vancomycin targets lipid II, blocking cell wall synthesis in Gram-positive bacteria. Despite extensive study, questions remain regarding how it recognizes its primary ligand and what is the most biologically relevant form of vancomycin. In this study, molecular dynamics simulation techniques have been used to examine the process of ligand binding and dimerization of vancomycin. Starting from one or more vancomycin monomers in solution, together with different peptide ligands derived from lipid II, the simulations predict the structures of the ligated monomeric and dimeric complexes to within 0.1 nm rmsd of the structures determined experimentally. The simulations reproduce the conformation transitions observed by NMR and suggest that proposed differences between the crystal structure and the solution structure are an artifact of the way the NMR data has been interpreted in terms of a structural model. The spontaneous formation of both back-to-back and face-to-face dimers was observed in the simulations. This has allowed a detailed analysis of the origin of the cooperatively between ligand binding and dimerization and suggests that the formation of face-to-face dimers could be functionally significant. The work also highlights the possible role of structural water in stabilizing the vancomycin ligand complex and its role in the manifestation of vancomycin resistance.

The Effect of Environment on the Structure of a Membrane Protein: P-Glycoprotein under Physiological Conditions

The stability of the crystal structure of the multidrug transporter P-glycoprotein proposed by Aller et al. (PDBid 3G5U ) has been examined under different environmental conditions using molecular dynamics. We show that in the presence of the detergent cholate, the structure of P-glycoprotein solved at pH 7.5 is stable. However, when incorporated into a cholesterol-enriched POPC membrane in the presence of 150 mM NaCl, the structure rapidly deforms. Only when the simulation conditions closely matched the experimental conditions under which P-glycoprotein is transport active was a stable conformation obtained. Specifically, the presence of Mg(2+), which bound to distinct sites in the nucleotide binding domains (NBDs), and the double protonation of the catalytic histidines (His583 and His1228) and His149 were required. While the structure obtained in a membrane environment under these conditions is very similar to the crystal structure of Aller et al., there are several key differences. The NBDs are in direct contact, reminiscent of the open state of MalK. The angle between the transmembrane domains is also increased, resulting in an outward motion of the intracellular loops. Notably, the structures obtained from the simulations provide a better match to a range of experimental cross-linking data than does the original 3G5U-a crystal structure. This work highlights the effect small changes in environmental conditions can have of the conformation of a membrane protein and the importance of representing the experimental conditions appropriately in modeling studies.

Orientation of μ-conotoxin PIIIA in a sodium channel vestibule, based on voltage dependence of its binding

Mutant cycle analysis has been used in previous studies to constrain possible docking orientations for various toxins. As an independent test of the bound orientation of μ-conotoxin PIIIA, a selectively targeted sodium channel pore blocker, we determined the contributions to binding voltage dependence of specific residues on the surface of the toxin. A change in the "apparent valence" (zδ) of the block, which is associated with a change of a specific toxin charge, reflects a change in the charge movement within the transmembrane electric field as the toxin binds. Toxin derivatives with charge-conserving mutations (R12K, R14K, and K17R) showed zδ values similar to those of wild type (0.61 ± 0.01, mean ± S.E.M.). Charge-changing mutations produced a range of responses. Neutralizing substitutions for Arg14 and Lys17 showed the largest reductions in zδ values, to 0.18 ± 0.06 and 0.20 ± 0.06, respectively, whereas unit charge-changing substitutions for Arg12, Ser13, and Arg20 gave intermediate values (0.24 ± 0.07, 0.33 ± 0.04, and 0.32 ± 0.05), which suggests that each of these residues contributes to the dependence of binding on the transmembrane voltage. Two mutations, R2A and G6K, yielded no significant change in zδ. These observations suggest that the toxin binds with Arg2 and Gly6 facing the extracellular solution, and Arg14 and Lys17 positioned most deeply in the pore. In this study, we used molecular dynamics to simulate toxin docking and performed Poisson-Boltzmann calculations to estimate the changes in local electrostatic potential when individual charges were substituted on the toxin's surface. Consideration of two limiting possibilities suggests that most of the charge movement associated with toxin binding reflects sodium redistribution within the narrow part of the pore.

Effect of poly(ethylene glycol) (PEG) spacers on the conformational properties of small peptides: a molecular dynamics study

Poly(ethylene glycol) (PEG) is used as an inert spacer in a wide range of biotechnological applications such as to display peptides and proteins on surfaces for diagnostic purposes. In such applications it is critical that the peptide is accessible to solvent and that the PEG does not affect the conformational properties of the peptide to which it is attached. Using molecular dynamics (MD) simulation techniques, we have investigated the influence of a commonly used PEG spacer on the conformation properties of a series of five peptides with differing physical-chemical properties (YGSLPQ, VFVVFV, GSGGSG, EEGEEG, and KKGKKG). The conformational properties of the peptides were compared (a) free in solution, (b) attached to a PEG-11 spacer in solution, and (c) constrained to a two-dimensional lattice via a (PEG-11)(3) spacer, mimicking a peptide displayed on a surface as used in microarray techniques. The simulations suggest that the PEG spacer has little effect on the conformational properties of small neutral peptides but has a significant effect on the conformational properties of small highly charged peptides. When constrained to a two-dimensional surface at peptide densities similar to those used experimentally, it was found that the peptides, in particular the polar and nonpolar peptides, aggregated strongly. The peptides also partitioned into the PEG layer. Potentially, this means that at high packing densities only a small fraction of the peptide attached to the surface would in fact be accessible to a potential interaction partner.

Transmembrane helix 12 modulates progression of the ATP catalytic cycle in ABCB1

Multidrug efflux pumps, such as P-glycoprotein (ABCB1), present major barriers to the success of chemotherapy in a number of clinical settings. Molecular details of the multidrug efflux process by ABCB1 remain elusive, in particular, the interdomain communication associated with bioenergetic coupling. The present investigation has focused on the role of transmembrane helix 12 (TM12) in the multidrug efflux process of ABCB1. Cysteine residues were introduced at various positions within TM12, and their effect on ATPase activity, nucleotide binding, and drug interaction were assessed. Mutation of several residues within TM12 perturbed the maximal ATPase activity of ABCB1, and the underlying cause was a reduction in basal (i.e., drug-free) hydrolysis of the nucleotide. Two of the mutations (L976C and F978C) were found to reduce the binding of [gamma-(32)P]-azido-ATP to ABCB1. In contrast, the A980C mutation within TM12 enhanced the rate of ATP hydrolysis; once again, this was due to modified basal activity. Several residues also caused reductions in the potency of stimulation of ATP hydrolysis by nicardipine and vinblastine, although the effects were independent of changes in drug binding per se. Overall, the results indicate that TM12 plays a key role in the progression of the ATP hydrolytic cycle in ABCB1, even in the absence of the transported substrate.

The mechanism of ABC transporters: general lessons from structural and functional studies of an antigenic peptide transporter

The shuttling of substrates across a cellular membrane frequently requires a specialized ATP-binding cassette (ABC) transporter, which couples the energy of ATP binding and hydrolysis to substrate transport. Due to its importance in immunity, the ABC transporter associated with antigen processing (TAP) has been studied extensively and is an excellent model for other ABC transporters. The TAP protein pumps cytosolic peptides into the endoplasmic reticulum for loading onto class I major histocompatibility complex (MHC) for subsequent immune surveillance. Here, we outline a potential mechanism for the TAP protein with supporting evidence from bacterial transporter structures. The similarities and differences between TAP and other transporters support the notion that ABC transporters in general have adapted around a universal transport mechanism.

Non-linear intramolecular interactions and voltage sensitivity of a KV1 family potassium channel from Polyorchis penicillatus (Eschscholtz 1829)

Voltage sensitivity of voltage-gated potassium channels (VKCs) is a primary factor in shaping action potentials in excitable cells. Variation in the amino acid sequence of the channel proteins is responsible for differences in the voltage range over which the channel opens. Thus, understanding how changes in voltage sensitivity are effected by changes in channel protein sequence illuminates the functional evolution of excitability. The K(V)1-family channel jShak1, from the jellyfish Polyorchis penicillatus, differs from most other K(V)1 channels in ways that are useful for studying the problem of how voltage sensitivity is related to channel sequence. We assessed the contributions of changes in sequence of the S4, voltage sensing, helix and changes in one asparagine residue in the S2 helix, to the relative stability of the open and closed states of the channel. Mutation of the neutral S2 residue (Asn227) to glutamate stabilized the open conformation of the channel. Different modifications of charge and length in S4 favoured either the closed conformation or the open conformation. The interactions between pairs of mutations revealed that some of the S4 mutations alter the conformation of the voltage-sensing domain such that the S4 helix is constrained to be closer to the S2 helix than in the wild-type conformation. These results, taken in conjunction with three-dimensional models of the channel, identify intra-molecular interactions that control the balance between open and closed states. These interactions are likely to be relevant to understanding the functional characteristics of members of this channel family from other organisms.

Structure-based interpretation of the mutagenesis database for the nucleotide binding domains of P-glycoprotein

P-glycoprotein (P-gp) is the most intensively studied eukaryotic ATP binding cassette (ABC) transporter, due to its involvement in the multidrug resistance phenotype of a number of cancers. In common with most ABC transporters, P-gp is comprised of two transmembrane domains (TMDs) and two nucleotide binding domains (NBD), the latter coupling ATP hydrolysis with substrate transport (efflux in the case of P-gp). Biochemical investigations over the past twenty years have attempted to unlock mechanistic aspects of P-glycoprotein through scanning and site-directed mutagenesis of both the TMDs and the NBDs. Contemporaneously, crystallographers have elucidated the atomic structure of numerous ABC transporter NBDs, as well as the intact structure (i.e. NBDs and TMDs) of a distantly related ABC-exporter Sav1866. Significantly, the structure of P-gp remains unknown, and only low resolution electron microscopy data exists. Within the current manuscript we employ crystallographic data for homologous proteins, and a molecular model for P-gp, to perform a structural interpretation of the existing "mutagenesis database" for P-gp NBDs. Consequently, this will enable testable predictions to be made that will result in further in-roads into our understanding of this clinically important drug pump.

Residue G346 in transmembrane segment six is involved in inter-domain communication in P-glycoprotein

Multidrug transporters such as P-glycoprotein require considerable inter-domain communication to couple energy utilization with substrate translocation. Elucidation of the regions or residues involved in these communication pathways is a key step in the eventual molecular description of multidrug transport. We used cysteine-scanning mutagenesis to probe the functional involvement of residues along the cytoplasmic half of transmembrane segment 6 (TM6) and its extension toward the nucleotide binding domain. The mutation of one residue (G346C) in this segment adversely affected drug transport in cells. Further investigation using purified protein revealed that the underlying biochemical effect was a reduction in basal ATP hydrolysis. This G346C mutation also affected the stimulation of ATPase activity in a drug dependent manner but had no effect on drug binding, ATP binding, or ADP release. Homology modeling of P-glycoprotein indicated that the G346C mutation caused a steric interaction between TM5 and TM6, thereby precluding a helical movement required to support ATP hydrolysis.

P-glycoprotein models of the apo and ATP-bound states based on homology with Sav1866 and MalK

We exploit the biochemical and sequence similarity between Staphylococcus aureus Sav1866 and P-glycoprotein to develop a homology model of P-glycoprotein representing an ATP-bound state, which captures the major features of the low-resolution EM structure and is consistent with cysteine mutagenesis studies. Using insights from the MalK crystal structures and BtuCD simulations, we model two nucleotide-free conformations. Conformational changes are characterized by pincering rigid-body rotations of the nucleotide-binding domains, inducing transmembrane domain reorganizations which correspond to the two lowest frequency normal modes of the protein. These conformations (see supplementary material) may characterize some of the major steps in the nucleotide catalytic cycle.

Computer simulations of ABC transporter componentsThis paper is one of a selection of papers published in this Special Issue, entitled CSBMCB — Membrane Proteins in Health and Disease

Current computer simulation techniques provide robust tools for studying the detailed structure and functional dynamics of proteins, as well as their interaction with each other and with other biomolecules. In this minireview, we provide an illustration of recent progress and future challenges in computer modeling by discussing computational studies of ATP-binding cassette (ABC) transporters. ABC transporters have multiple components that work in a well coordinated fashion to enable active transport across membranes. The mechanism by which members of this superfamily execute transport remains largely unknown. Molecular dynamics simulations initiated from high-resolution crystal structures of several ABC transporters have proven to be useful in the investigation of the nature of conformational coupling events that may drive transport. In addition, fruitful efforts have been made to predict unknown structures of medically relevant ABC transporters, such as P-glycoprotein, using homology-based computational methods. The various techniques described here are also applicable to gaining an atomically detailed understanding of the functional mechanisms of proteins in general.