Functional evaluation of autism-associated mutations in NHE9.

NHE9 (SLC9A9) is an endosomal cation/proton antiporter with orthologues in yeast and bacteria. Rare, missense substitutions in NHE9 are genetically linked with autism but have not been functionally evaluated. Here we use evolutionary conservation analysis to build a model structure of NHE9 based on the crystal structure of bacterial NhaA and use it to screen autism-associated variants in the human population first by phenotype complementation in yeast, followed by functional analysis in primary cortical astrocytes from mouse. NHE9-GFP localizes to recycling endosomes, where it significantly alkalinizes luminal pH, elevates uptake of transferrin and the neurotransmitter glutamate, and stabilizes surface expression of transferrin receptor and GLAST transporter. In contrast, autism-associated variants L236S, S438P and V176I lack function in astrocytes. Thus, we establish a neurobiological cell model of a candidate gene in autism. Loss-of-function mutations in NHE9 may contribute to autistic phenotype by modulating synaptic membrane protein expression and neurotransmitter clearance.

Structure and flexibility of the C-ring in the electromotor of rotary F(0)F(1)-ATPase of pea chloroplasts.

A ring of 8-15 identical c-subunits is essential for ion-translocation by the rotary electromotor of the ubiquitous F(O)F(1)-ATPase. Here we present the crystal structure at 3.4Šresolution of the c-ring from chloroplasts of a higher plant (Pisum sativum), determined using a native preparation. The crystal structure was found to resemble that of an (ancestral) cyanobacterium. Using elastic network modeling to investigate the ring's eigen-modes, we found five dominant modes of motion that fell into three classes. They revealed the following deformations of the ring: (I) ellipsoidal, (II) opposite twisting of the luminal circular surface of the ring against the stromal surface, and (III) kinking of the hairpin-shaped monomers in the middle, resulting in bending/stretching of the ring. Extension of the elastic network analysis to rings of different c(n)-symmetry revealed the same classes of dominant modes as in P. sativum (c(14)). We suggest the following functional roles for these classes: The first and third classes of modes affect the interaction of the c-ring with its counterparts in F(O), namely subunits a and bb'. These modes are likely to be involved in ion-translocation and torque generation. The second class of deformation, along with deformations of subunits γ and ε might serve to elastically buffer the torque transmission between F(O) and F(1).

New model of cystic fibrosis transmembrane conductance regulator proposes active channel-like conformation.

The cystic fibrosis transmembrane conductance regulator (CFTR) is an unusual ABC transporter, functioning as a chloride channel critical for fluid homeostasis in multiple organs. Disruption of CFTR function is associated with cystic fibrosis making it an attractive therapeutic target. In addition, CFTR blockers are being developed as potential antidiarrheals. CFTR drug discovery is hampered by the lack of high resolution structural data, and considerable efforts have been invested in modeling the channel structure. Although previously published CFTR models that have been made publicly available mostly agree with experimental data relating to the overall structure, they present the channel in an outward-facing conformation that does not agree with expected properties of a "channel-like" structure. Here, we make available a model of CFTR in such a "channel-like" conformation, derived by a unique modeling approach combining restrained homology modeling and ROSETTA refinement. In contrast to others, the present model is in agreement with expected channel properties such as pore shape, dimensions, solvent accessibility, and experimentally derived distances. We have used the model to explore the interaction of open channel blockers within the pore, revealing a common binding mode and ionic interaction with K95, in agreement with experimental data. The binding-site was further validated using a virtual screening enrichment experiment, suggesting the model might be suitable for drug discovery. In addition, we subjected the model to a molecular dynamics simulation, revealing previously unaddressed salt-bridge interactions that may be important for structure stability and pore-lining residues that may take part in Cl(-) conductance.

A structural model of the copper ATPase ATP7B to facilitate analysis of Wilson disease-causing mutations and studies of the transport mechanism.

The copper-transporting ATPase ATP7B has an essential role in human physiology, particularly for the liver and brain function. Inactivation of ATP7B is associated with a severe hepato-neurologic disorder, Wilson disease (WD). Hundreds of WD related mutations have been identified in ATP7B to date. The low frequency and the compound-heterozygous nature of causative mutations complicate the analysis of individual mutants and the establishment of genotype-phenotype correlations. To facilitate studies of disease-causing mutations and mechanistic understanding of WD, we have homology-modelled the ATP7B core (residues 643-1377) using the recent structure of the bacterial copper-ATPase LCopA as a template. The model, supported by evolutionary conservation and hydrophobicity analysis, as well as existing and new mutagenesis data, allows molecular interpretations of experimentally characterized clinical mutations. We also illustrate that structure and conservation can be used to grade potential deleterious effects for many WD mutations, which were clinically detected but have not yet been experimentally characterized. Finally, we compare the structural features of ATP7B and LCopA and discuss specific features of the eukaryotic copper pump.

A model-structure of a periplasm-facing state of the NhaA antiporter suggests the molecular underpinnings of pH-induced conformational changes.

The Escherichia coli NhaA antiporter couples the transport of H(+) and Na(+) (or Li(+)) ions to maintain the proper pH range and Na(+) concentration in cells. A crystal structure of NhaA, solved at pH 4, comprises 12 transmembrane helices (TMs), arranged in two domains, with a large cytoplasm-facing funnel and a smaller periplasm-facing funnel. NhaA undergoes conformational changes, e.g. after pH elevation to alkaline ranges, and we used two computational approaches to explore them. On the basis of pseudo-symmetric features of the crystal structure, we predicted the structural architecture of an alternate, periplasm-facing state. In contrast to the crystal structure, the model presents a closed cytoplasmic funnel, and a periplasmic funnel of greater volume. To examine the transporter functional direction of motion, we conducted elastic network analysis of the crystal structure and detected two main normal modes of motion. Notably, both analyses predicted similar trends of conformational changes, consisting of an overall rotational motion of the two domains around a putative symmetry axis at the funnel centers, perpendicular to the membrane plane. This motion, along with conformational changes within specific helices, resulted in closure at the cytoplasmic end and opening at the periplasmic end. Cross-linking experiments, performed between segments on opposite sides of the cytoplasmic funnel, revealed pH-dependent interactions consistent with the proposed conformational changes. We suggest that the model-structure and predicted motion represent alkaline pH-induced conformational changes, mediated by a cluster of evolutionarily conserved, titratable residues, at the cytoplasmic ends of TMs II, V, and IX.

C(alpha)-trace model of the transmembrane domain of human copper transporter 1, motion and functional implications.

The trimeric human copper transporter 1 (hCTR1) is essential for copper uptake and is implicated in sensitivity to chemotherapy drugs. Using the cryoelectron microscopy (cryoEM) map of hCTR1 and evolutionary data, we constructed a Calpha-trace model of the membrane region. The model structure, supported by mutagenesis data, was used to investigate global dynamics through elastic network models. Identified as dominant hinge regions, hCTR1's MxxxM and GxxxG motifs were shown to have significant roles in functional movements characterized by the two slowest modes of motion. Both modes predicted significant changes at the wide cytoplasmic pore region; the slowest mode introduced a rotational motion around the pore central axis, whereas in the following mode the cytoplasmic parts of the helices approached and moved away from the pore center. In the most cooperative mode, the MxxxM motif in the extracellular narrow region remained static. The second mode of motion, however, predicted a cooperative rotational motion of this copper-binding motif, possibly reflecting activation at the pore's extracellular entrance. We suggest a molecular mechanism of copper transport in which this motif serves both as a gate and as a selectivity filter. We also suggest residues that are responsible for pH activation.

Model-guided mutagenesis drives functional studies of human NHA2, implicated in hypertension.

Human NHA2 is a poorly characterized Na(+)/H(+) antiporter recently implicated in essential hypertension. We used a range of computational tools and evolutionary conservation analysis to build and validate a three-dimensional model of NHA2 based on the crystal structure of a distantly related bacterial transporter, NhaA. The model guided mutagenic evaluation of transport function, ion selectivity, and pH dependence of NHA2 by phenotype screening in yeast. We describe a cluster of essential, highly conserved titratable residues located in an assembly region made of two discontinuous helices of inverted topology, each interrupted by an extended chain. Whereas in NhaA, oppositely charged residues compensate for partial dipoles generated within this assembly, in NHA2, polar but uncharged residues suffice. Our findings led to a model for transport mechanism that was compared to the well-known electroneutral NHE1 and electrogenic NhaA subtypes. This study establishes NHA2 as a prototype for the poorly understood, yet ubiquitous, CPA2 antiporter family recently recognized in plants and metazoans and illustrates a structure-driven approach to derive functional information on a newly discovered transporter.

iDBPs: a web server for the identification of DNA binding proteins.

SUMMARY: The iDBPs server uses the three-dimensional (3D) structure of a query protein to predict whether it binds DNA. First, the algorithm predicts the functional region of the protein based on its evolutionary profile; the assumption is that large clusters of conserved residues are good markers of functional regions. Next, various characteristics of the predicted functional region as well as global features of the protein are calculated, such as the average surface electrostatic potential, the dipole moment and cluster-based amino acid conservation patterns. Finally, a random forests classifier is used to predict whether the query protein is likely to bind DNA and to estimate the prediction confidence. We have trained and tested the classifier on various datasets and shown that it outperformed related methods. On a dataset that reflects the fraction of DNA binding proteins (DBPs) in a proteome, the area under the ROC curve was 0.90. The application of the server to an updated version of the N-Func database, which contains proteins of unknown function with solved 3D-structure, suggested new putative DBPs for experimental studies. AVAILABILITY: http://idbps.tau.ac.il/

Detection of functionally important regions in "hypothetical proteins" of known structure.

Structural genomics initiatives provide ample structures of "hypothetical proteins" (i.e., proteins of unknown function) at an ever increasing rate. However, without function annotation, this structural goldmine is of little use to biologists who are interested in particular molecular systems. To this end, we used (an improved version of) the PatchFinder algorithm for the detection of functional regions on the protein surface, which could mediate its interactions with, e.g., substrates, ligands, and other proteins. Examination, using a data set of annotated proteins, showed that PatchFinder outperforms similar methods. We collected 757 structures of hypothetical proteins and their predicted functional regions in the N-Func database. Inspection of several of these regions demonstrated that they are useful for function prediction. For example, we suggested an interprotein interface and a putative nucleotide-binding site. A web-server implementation of PatchFinder and the N-Func database are available at http://patchfinder.tau.ac.il/.