Carmen Domene

Structural basis for binding of hypoxia-inducible factor to the oxygen-sensing prolyl hydroxylases

The oxygen-dependent hydroxylation of proline residues in the alpha subunit of hypoxia-inducible transcription factor (HIFalpha) is central to the hypoxic response in animals. Prolyl hydroxylation of HIFalpha increases its binding to the von Hippel-Lindau protein (pVHL), so signaling for degradation via the ubiquitin-proteasome system. The HIF prolyl hydroxylases (PHDs, prolyl hydroxylase domain enzymes) are related to the collagen prolyl hydroxylases, but form unusually stable complexes with their Fe(II) cofactor and 2-oxoglutarate cosubstrate. We report crystal structures of the catalytic domain of PHD2, the most important of the human PHDs, in complex with the C-terminal oxygen-dependent degradation domain of HIF-1alpha. Together with biochemical analyses, the results reveal that PHD catalysis involves a mobile region that isolates the hydroxylation site and stabilizes the PHD2.Fe(II).2OG complex. The results will be of use in the design of PHD inhibitors aimed at treating anemia and ischemic disease.

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 Å) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (∼1 Å) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge‐bending of the pore‐lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.

Potassium Channel, Ions, and Water: Simulation Studies Based on the High Resolution X-Ray Structure of KcsA

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.

Conformational sampling and dynamics of membrane proteins from 10-nanosecond computer simulations.

In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small β‐barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the Cα‐atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time‐dependence of the average structure was analyzed to obtain single‐domain information. The membrane‐embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extra‐membranous domains lacking secondary structure are often markedly under‐sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B‐factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations. Proteins 2004.

On conduction in a bacterial sodium channel.

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 Ca2+ 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.

Atypical mechanism of conduction in potassium channels

Potassium channels can conduct passively K+ ions with rates of up to ≈108 ions per second at physiological conditions, and they are selective to these species by a factor of 104 over Na+ ions. Ion conduction has been proposed to involve transitions between 2 main states, with 2 or 3 K+ ions occupying the selectivity filter separated by an intervening water molecule. The largest free energy barrier of such a process was reported to be of the order of 2–3 kcal mol−1. Here, we present an alternative mechanism for conduction of K+ in potassium channels where site vacancies are involved, and we propose that coexistence of several ion permeation mechanisms is energetically possible. Conduction can be described as a more anarchic phenomenon than previously characterized by the concerted translocations of K+–water–K+.

Dynamics of a bacterial multidrug ABC transporter in the inward- and outward-facing conformations

The study of membrane proteins remains a challenging task, and approaches to unravel their dynamics are scarce. Here, we applied hydrogen/deuterium exchange (HDX) coupled to mass spectrometry to probe the motions of a bacterial multidrug ATP-binding cassette (ABC) transporter, BmrA, in the inward-facing (resting state) and outward-facing (ATP-bound) conformations. Trypsin digestion and global or local HDX support the transition between inward- and outward-facing conformations during the catalytic cycle of BmrA. However, in the resting state, peptides from the two intracellular domains, especially ICD2, show a much faster HDX than in the closed state. This shows that these two subdomains are very flexible in this conformation. Additionally, molecular dynamics simulations suggest a large fluctuation of the Cα positions from ICD2 residues in the inward-facing conformation of a related transporter, MsbA. These results highlight the unexpected flexibility of ABC exporters in the resting state and underline the power of HDX coupled to mass spectrometry to explore conformational changes and dynamics of large membrane proteins.

Conformational Changes and Gating at the Selectivity Filter of Potassium Channels

The translocation of ions and water across cell membranes is a prerequisite for many of lifes processes. K(+) channels are a diverse family of integral membrane proteins through which K(+) can pass selectively. There is an ongoing debate about the nature of conformational changes associated with the opening and closing and conductive and nonconductive states of potassium (K(+)) channels. These changes depend on the membrane potential, the K(+) concentration gradient, and large scale motions of transmembrane helices and associated residues. Experiments also suggest that local structural changes in the selectivity filter may act as the dominant gate referred to as C-type inactivation. Herein we present an extensive computational study on KirBac, which supports the existence of a physical gate or constriction in the selectivity filter (SF) of K(+) channels. Our computations identify a new selectivity filter structure, which is likely associated with C-type inactivation. Specifically, the four peptide chains that comprise the filter adopt an unusual structure in which their dihedrals alternate between left- and right-handed Ramachandran angles, which also justifies the need for conservation of glycine in the K(+) selectivity filter, since it is the only residue able to play this bifunctional role.

Membrane protein simulations: Ion channels and bacterial outer membrane proteins

Publisher Summary The chapter focuses on two classes that have been investigated by a number of research groups—namely, ion channels and bacterial outer membrane proteins—rather than attempt to cover all simulations of membrane proteins. Ion channels play a key role in the electrical activity of cells of the nervous system. However, channels are also found in the membranes of non-excitable cells and a wide range of organisms from viruses to plants. The chapter focuses on simulation studies of a bacterial K channel, KcsA, and of related K channel species. Gram-negative bacterial cells are enclosed by two membranes. The outer membrane is highly asymmetric, whereas its inner leaflet has a similar phospholipid composition to that of the inner membrane: its outer leaflet contains lipo-polysaccharide molecules. A combination of molecular modeling and simulation helps us to extrapolate from the structure of prokaryotic membrane proteins to the structure and dynamics of their human homologues, which may also aid in experimental structure determination. Ultimately, simulation offers the prospect of understanding the effects of mutations on physiological function in a variety of disease states, as well as aiding the discovery and design of drugs to help combat such diseases.

Many-body effects and simulations of potassium channels

The electronic polarizability of an ion or a molecule is a measure of the relative tendency of its electron cloud to be distorted from its normal shape by an electric field. On the molecular scale, in a condensed phase, any species sits in an electric field due to its neighbours, and the resulting polarization is an important contribution to the total interaction energy. Electrostatic interactions are crucial for determining the majority of chemical–physical properties of the system and electronic polarization is a fundamental component of these interactions. Thus, polarization effects should be taken into account if accurate descriptions are desired. In classical computer simulations, the forces required to drive the system are typically based on interatomic interaction potentials derived in part from electronic structure calculations or from experimental data. Owing to the difficulties in including polarization effects in classical force fields, most of them are based just on pairwise additive interaction potentials. At present, major efforts are underway to develop polarizable interaction potentials for biomolecular simulations. In this review, various ways of introducing explicit polarizability into biomolecular models and force fields are reviewed, and the progress that might be achieved in applying such methods to study potassium channels is described.