Defluorination Capability of l-2-Haloacid Dehalogenases in the HAD-Like Hydrolase Superfamily Correlates with Active Site Compactness.

l-2-Haloacid dehalogenases, industrially and environmentally important enzymes that catalyse cleavage of the carbon-halogen bond in S-2-halocarboxylic acids, were known to hydrolyse chlorinated, brominated and iodinated substrates but no activity towards fluorinated compounds had been reported. A screen for novel dehalogenase activities revealed four l-2-haloacid dehalogenases capable of defluorination. We now report crystal structures for two of these enzymes, Bpro0530 and Rha0230, as well as for the related proteins PA0810 and RSc1362, which hydrolyse chloroacetate but not fluoroacetate, all at ∼2.2 Šresolution. Overall structure and active sites of these enzymes are highly similar. In molecular dynamics (MD) calculations, only the defluorinating enzymes sample more compact conformations, which in turn allow more effective interactions with the small fluorine atom. Structural constraints, based on X-ray structures and MD calculations, correctly predict the defluorination activity of the homologous enzyme ST2570.

Mechanism of Amyloidogenesis of a Bacterial AAA+ Chaperone.

Amyloids are fibrillar protein superstructures that are commonly associated with diseases in humans and with physiological functions in various organisms. The precise mechanisms of amyloid formation remain to be elucidated. Surprisingly, we discovered that a bacterial Escherichia coli chaperone-like ATPase, regulatory ATPase variant A (RavA), and specifically the LARA domain in RavA, forms amyloids under acidic conditions at elevated temperatures. RavA is involved in modulating the proper assembly of membrane respiratory complexes. LARA contains an N-terminal loop region followed by a ß-sandwich-like folded core. Several approaches, including nuclear magnetic resonance spectroscopy and molecular dynamics simulations, were used to determine the mechanism by which LARA switches to an amyloid state. These studies revealed that the folded core of LARA is amyloidogenic and is protected by its N-terminal loop. At low pH and high temperatures, the interaction of the N-terminal loop with the folded core is disrupted, leading to amyloid formation.

Atomistic picture of conformational exchange in a T4 lysozyme cavity mutant: an experiment-guided molecular dynamics study.

Despite the importance of dynamics to protein function there is little information about the states that are formed as the protein explores its conformational landscape or about the mechanism by which transitions between the different states occur. Here we used a combined NMR spin relaxation and unbiased molecular dynamics (MD) approach to investigate the exchange process by which a cavity in an L99A mutant of T4 lysozyme (T4L 99A) interconverts between an empty and occupied form that involves repositioning of an aromatic residue, Phe114. Although structures of the end-states of the exchange process are available, insight into the mechanism by which the transition takes place cannot be obtained from experiment and the timescales involved are too slow to address using brute force MD. Using spin relaxation NMR methods, we have identified a triple-mutant of T4L that undergoes the same exchange process as T4L L99A but where the minor state lifetime has decreased significantly so that the spontaneous conformational transition to the major state can be studied using all-atom MD simulations. The simulation trajectories were analyzed using Markov state models and the energy landscape so obtained is in good agreement with expectations based on NMR studies. Notably there is no large-scale perturbation of the structure during the transition, multiple intermediates are formed between the two similar exchanging conformations and the free energy barrier between these two well-folded, compact forms is small (6kBT), only slightly larger than for processes considered to be barrierless.

Hydrophobic Gating of Ion Permeation in Magnesium Channel CorA.

Ion channels catalyze ionic permeation across membranes via water-filled pores. To understand how changes in intracellular magnesium concentration regulate the influx of Mg2+ into cells, we examine early events in the relaxation of Mg2+ channel CorA toward its open state using massively-repeated molecular dynamics simulations conducted either with or without regulatory ions. The pore of CorA contains a 2-nm-long hydrophobic bottleneck which remained dehydrated in most simulations. However, rapid hydration or "wetting" events concurrent with small-amplitude fluctuations in pore diameter occurred spontaneously and reversibly. In the absence of regulatory ions, wetting transitions are more likely and include a wet state that is significantly more stable and more hydrated. The free energy profile for Mg2+ permeation presents a barrier whose magnitude is anticorrelated to pore diameter and the extent of hydrophobic hydration. These findings support an allosteric mechanism whereby wetting of a hydrophobic gate couples changes in intracellular magnesium concentration to the onset of ionic conduction.

Catalysis of Na+ permeation in the bacterial sodium channel Na(V)Ab.

Determination of a high-resolution 3D structure of voltage-gated sodium channel Na(V)Ab opens the way to elucidating the mechanism of ion conductance and selectivity. To examine permeation of Na(+) through the selectivity filter of the channel, we performed large-scale molecular dynamics simulations of Na(V)Ab in an explicit, hydrated lipid bilayer at 0 mV in 150 mM NaCl, for a total simulation time of 21.6 µs. Although the cytoplasmic end of the pore is closed, reversible influx and efflux of Na(+) through the selectivity filter occurred spontaneously during simulations, leading to equilibrium movement of Na(+) between the extracellular medium and the central cavity of the channel. Analysis of Na(+) dynamics reveals a knock-on mechanism of ion permeation characterized by alternating occupancy of the channel by 2 and 3 Na(+) ions, with a computed rate of translocation of (6 ± 1) × 10(6) ions⋅s(-1) that is consistent with expectations from electrophysiological studies. The binding of Na(+) is intimately coupled to conformational isomerization of the four E177 side chains lining the extracellular end of the selectivity filter. The reciprocal coordination of variable numbers of Na(+) ions and carboxylate groups leads to their condensation into ionic clusters of variable charge and spatial arrangement. Structural fluctuations of these ionic clusters result in a myriad of ion binding modes and foster a highly degenerate, liquid-like energy landscape propitious to Na(+) diffusion. By stabilizing multiple ionic occupancy states while helping Na(+) ions diffuse within the selectivity filter, the conformational flexibility of E177 side chains underpins the knock-on mechanism of Na(+) permeation.

Construction and validation of a homology model of the human voltage-gated proton channel hHV1.

The topological similarity of voltage-gated proton channels (H(V)1s) to the voltage-sensing domain (VSD) of other voltage-gated ion channels raises the central question of whether H(V)1s have a similar structure. We present the construction and validation of a homology model of the human H(V)1 (hH(V)1). Multiple structural alignment was used to construct structural models of the open (proton-conducting) state of hH(V)1 by exploiting the homology of hH(V)1 with VSDs of K(+) and Na(+) channels of known three-dimensional structure. The comparative assessment of structural stability of the homology models and their VSD templates was performed using massively repeated molecular dynamics simulations in which the proteins were allowed to relax from their initial conformation in an explicit membrane mimetic. The analysis of structural deviations from the initial conformation based on up to 125 repeats of 100-ns simulations for each system reveals structural features consistently retained in the homology models and leads to a consensus structural model for hH(V)1 in which well-defined external and internal salt-bridge networks stabilize the open state. The structural and electrostatic properties of this open-state model are compatible with proton translocation and offer an explanation for the reversal of charge selectivity in neutral mutants of Asp(112). Furthermore, these structural properties are consistent with experimental accessibility data, providing a valuable basis for further structural and functional studies of hH(V)1. Each Arg residue in the S4 helix of hH(V)1 was replaced by His to test accessibility using Zn(2+) as a probe. The two outermost Arg residues in S4 were accessible to external solution, whereas the innermost one was accessible only to the internal solution. Both modeling and experimental data indicate that in the open state, Arg(211), the third Arg residue in the S4 helix in hH(V)1, remains accessible to the internal solution and is located near the charge transfer center, Phe(150).

Structural asymmetry in the magnesium channel CorA points to sequential allosteric regulation.

Magnesium ions (Mg(2+)) are essential for life, but the mechanisms regulating their transport into and out of cells remain poorly understood. The CorA-Mrs2-Alr1 superfamily of Mg(2+) channels represents the most prevalent group of proteins enabling Mg(2+) ions to cross membranes. Thermotoga maritima CorA (TmCorA) is the only member of this protein family whose complete 3D fold is known. Here, we report the crystal structure of a mutant in the presence and absence of divalent ions and compare it with previous divalent ion-bound TmCorA structures. With Mg(2+) present, this structure shows binding of a hydrated Mg(2+) ion to the periplasmic Gly-Met-Asn (GMN) motif, revealing clues of ion selectivity in this unique channel family. In the absence of Mg(2+), TmCorA displays an unexpected asymmetric conformation caused by radial and lateral tilts of protomers that leads to bending of the central, pore-lining helix. Molecular dynamics simulations support these movements, including a bell-like deflection. Mass spectrometric analysis confirms that major proteolytic cleavage occurs within a region that is selectively exposed by such a bell-like bending motion. Our results point to a sequential allosteric model of regulation, where intracellular Mg(2+) binding locks TmCorA in a symmetric, transport-incompetent conformation and loss of intracellular Mg(2+) causes an asymmetric, potentially influx-competent conformation of the channel.

An iris-like mechanism of pore dilation in the CorA magnesium transport system.

Magnesium translocation across cell membranes is essential for numerous physiological processes. Three recently reported crystal structures of the CorA magnesium transport system revealed a surprising architecture, with a bundle of giant alpha-helices forming a 60-A-long pore that extends beyond the membrane before widening into a funnel-shaped cytosolic domain. The presence of divalent cations in putative intracellular regulation sites suggests that these structures correspond to the closed conformation of CorA. To examine the nature of the conduction pathway, we performed 110-ns molecular-dynamics simulations of two of these structures in a lipid bilayer with and without regulatory ions. The results show that a 15-A-long hydrophobic constriction straddling the membrane-cytosol interface constitutes a steric bottleneck whose location coincides with an electrostatic barrier opposing cation translocation. In one of the simulations, structural relaxation after the removal of regulatory ions led to concerted changes in the tilt of the pore helices, resulting in iris-like dilation and spontaneous hydration of the hydrophobic neck. This simple and robust mechanism is consistent with the regulation of pore opening by intracellular magnesium concentration, and explains the unusual architecture of CorA.

Acidity of a Cu-bound histidine in the binuclear center of cytochrome C oxidase.

Cytochrome c oxidase (CcO) is a crucial enzyme in the respiratory chain. Its function is to couple the reduction of molecular oxygen, which takes place in the Fea3-CuB binuclear center, to proton translocation across the mitochondrial membrane. Although several high-resolution structures of the enzyme are known, the molecular basis of proton pumping activation and its mechanism remain to be elucidated. We examine a recently proposed scheme (J. Am. Chem. Soc. 2004, 126, 1858; FEBS Lett. 2004, 566, 126) that involves the deprotonation of the CuB-bound imidazole ring of a histidine (H291 in mammalian CcO) as a key element in the proton pumping mechanism. The central feature of that proposed mechanism is that the pKa values of the imidazole vary significantly depending on the redox state of the metals in the binuclear center. We use density functional theory in combination with continuum electrostatics to calculate the pKa values, successively in bulk water and within the protein, of the Cu-bound imidazole in various Cu- and Cu-Fe complexes. From pKas in bulk water, we derived a value of -266.34 kcal.mol(-1) for the proton solvation free energy (Delta). This estimate is in close agreement with the experimental value of -264.61 kcal.mol(-1) (J. Am. Chem. Soc. 2001, 123, 7314), which reinforces the conclusion that Delta is more negative than previous values used for pKa calculations. Our approach, on the basis of the study of increasingly more detailed models of the CcO binuclear center at different stages of the catalysis, allows us to examine successively the effect of each of the two metals' redox states and of solvation on the acidity of imidazole, whose pKa is approximately 14 in bulk water. This analysis leads to the following conclusions: first, the effect of Cu ligation on the imidazole acidity is negligible regardless of the redox state of the metal. Second, results obtained for Cu-Fe complexes in bulk water indicate that Cu-bound imidazole pKa values lie within the range of 14.8-16.6 throughout binuclear redox states corresponding to the catalytic cycle, demonstrating that the effect of the Fe oxidation states is also negligible. Finally, the low-dielectric CcO proteic environment shifts the acid-base equilibrium toward a neutral imidazole, further increasing the corresponding pKa values. These results are inconsistent with the proposed role of the Cu-bound histidine as a key element in the pumping mechanism. Limitations of continuum solvation models in pKa calculations are discussed.

Structural determinants of proton blockage in aquaporins.

Aquaporins are an important class of membrane channels selective for water and linear polyols but impermeable to ions, including protons. Recent computational studies have revealed that the relay of protons through the water-conduction pathway of aquaporin channels is opposed by a substantial free energy barrier peaking at the signature NPA motifs. Here, free-energy simulations and continuum electrostatic calculations are combined to examine the nature and the magnitude of the contribution of specific structural elements to proton blockage in the bacterial glycerol uptake facilitator, GlpF. Potential of mean-force profiles for both hop and turn steps of structural diffusion in the narrow pore are obtained for artificial variants of the GlpF channel in which coulombic interactions between the pore contents and conserved residues Asn68 and Asn203 at the NPA signature motifs, Arg206 at the selectivity filter, and the peptidic backbone of the two half-helices M3 and M7, which are arranged in head-to-head fashion around the NPA motifs, are turned off selectively. A comparison of these results with electrostatic energy profiles for the translocation of a probe cation throughout the water permeation pathway indicates that the free-energy profile for proton movement inside the narrow pore is dominated by static effects arising from the distribution of charged and polar groups of the channel, whereas dielectric effects contribute primarily to opposing the access of H+ to the pore mouths (desolvation penalty). The single most effective way to abolish the free-energy gradients opposing the movement of H+ around the NPA motif is to turn off the dipole moments of helices M3 and M7. Mutation of either of the two NPA Asn residues to Asp compensates for charge-dipole and dipole-dipole effects opposing the hop and turn steps of structural diffusion, respectively, and dramatically reduces the free energy barrier of proton translocation, suggesting that these single mutants could leak protons.

Molecular basis of proton blockage in aquaporins.

Water transport channels in membrane proteins of the aquaporin superfamily are impermeable to ions, including H+ and OH-. We examine the molecular basis for the blockage of proton translocation through the single-file water chain in the pore of a bacterial aquaporin, GlpF. We compute the reversible thermodynamic work for the two complementary steps of the Grotthuss "hop-and-turn" relay mechanism: consecutive transfers of H+ along the hydrogen-bonded chain (hop) and conformational reorganization of the chain (turn). In the absence of H+, the strong preference for the bipolar orientation of water around the two Asn-Pro-Ala (NPA) motifs lining the pore over both unidirectional polarization states of the chain precludes the reorganization of the hydrogen-bonded network. Inversely, translocation of an excess proton in either direction is opposed by a free-energy barrier centered at the NPA region. Both hop and turn steps of proton translocation are opposed by the electrostatic field of the channel.