Determining molecular forces that stabilize human aquaporin-1

Atomic force microscopy (AFM) was used to measure the forces stabilizing human aquaporin-1 (hAQP1), a tetrameric transmembrane protein that forms highly specific water channels. To this end, the AFM tip was attached to the C-terminus of hAQP1 and secondary structure elements were extracted from the membrane while the single-molecule force-extension curve was being recorded. Force peaks, reflecting the unfolding of secondary structure elements, could be interpreted in depth using the atomic model of hAQP1. Different classes of force-extension curves indicated the existence of alternative unfolding pathways for individual proteins. In addition, transmembrane helices at the periphery of the hAQP1 tetramer exhibited smaller extraction forces than helices at the interface between hAQP1 monomers. These results represent the first direct assessment of intermolecular forces defining the oligomeric state of a membrane protein.     

Crystal structure of human immunoglobulin fragment Fab New refined at 2.0 A resolution.

The three-dimensional structure of the human immunoglobulin fragment Fab New (IgG1, lambda) has been refined to a crystallographic R-factor of 16.9% to 2 A resolution. Rms deviations of the final model from ideal geometry are 0.014 A for bond distances and 3.03 degrees for bond angles. Refinement was based on a new X-ray data set including 28,301 reflections with F > 2.5 sigma(F) from 6.0 to 2.0 A resolution. The starting model for the refinement procedure reported here is from the Brookhaven Protein Data Bank entry 3FAB (rev. 1981). Differences between the initial and final models include modified polypeptide-chain folding in the third complementarity-determining region (CDR3) and the third framework region (FR3) of VH and in some exposed loops of CL and CH1. Amino acid sequence changes were determined at a number of positions by inspection of difference electron density maps. The incorporation of amino acid sequence changes results in an improved VH framework model for the "humanization" of monoclonal antibodies.     

Calmodulin structure refined at 1.7 A resolution.

We have determined and refined the crystal structure of a recombinant calmodulin at 1.7 A resolution. The structure was determined by molecular replacement, using the 2.2 A published native bovine brain structure as the starting model. The final crystallographic R-factor, using 14,469 reflections in the 10.0 to 1.7 A range with structure factors exceeding 0.5 sigma, is 0.216. Bond lengths and bond angle distances have root-mean-square deviations from ideal values of 0.009 A and 0.032 A, respectively. The final model consists of 1279 non-hydrogen atoms, including four calcium ions, 1130 protein atoms, including three Asp118 side-chain atoms in double conformation, 139 water molecules and one ethanol molecule. The electron densities for residues 1 to 4 and 148 of calmodulin are poorly defined, and not included in our model, except for main-chain atoms of residue 4. The calmodulin structure from our crystals is very similar to the earlier 2.2 A structure described by Babu and coworkers with a root-mean-square deviation of 0.36 A. Calmodulin remains a dumb-bell-shaped molecule, with similar lobes and connected by a central alpha-helix. Each lobe contains three alpha-helices and two Ca2+ binding EF hand loops, with a short antiparallel beta-sheet between adjacent EF hand loops and one non-EF hand loop. There are some differences in the structure of the central helix. The crystal packing is extensively studied, and facile crystal growth along the z-axis of the triclinic crystals is explained. Herein, we describe hydrogen bonding in the various secondary structure elements and hydration of calmodulin.     

The refined crystal structure of ribonuclease A at 2.0 A resolution

This paper describes the structure of bovine pancreatic ribonuclease A, refined by a restrained parameter least squares procedure at 2.0 A resolution, and rebuilt using computer graphics. The final agreement factor (formula see text) is 0.159. The positions of the 951 main chain atoms have been determined with an estimated accuracy of 0.17 A. In addition, the model includes a phosphate group in the active site and 176 waters, many of them with partial occupancy. The bond lengths in the refined structure of RNase A differ from the ideal values by an overall root mean square deviation of 0.022 A; the corresponding value for angle distances is 0.06 A. The root mean square deviation of planar atoms from ideality is 0.017 A, and root mean square deviation of the peptide torsion angles from 180 degrees is 3.4 degrees. The model is in good agreement with the final difference Fourier maps. Two active site histidines, His 12 and His 119, form hydrogen bonds to the phosphate ion. His 119 is also hydrogen bonded to the carboxyl of ASp 121 and His 12 to the carbonyl of Thr 45. The structure of the RNase A is very similar to that of RNase S, particularly in the active site region. The root mean square discrepancy of all atoms from residues 1 to 16 and 24 to 123 is 1.06 A and the root mean square discrepancy for the active site region is 0.6 A.     

Induction of human aquaporin-1 gene by retinoic acid in human erythroleukemia HEL cells

Retinoids have been implicated in the control of cell proliferation, differentiation, and developmental processes. We report here that aquaporin-1 (AQP1) is specifically induced by retinoic acid (RA) in human erythroleukemia HEL cells. Both all-trans-RA (ATRA) and 9-cis-RA (9CRA) strongly induced the AQP1 mRNA and protein in a dose-dependent manner. AQP1 protein was mainly expressed in plasma membrane in cells induced by RAs. To identify the RA response element (RARE) in the human AQP1 promoter, the 5(')-flanking region of AQP1 promoter was isolated and transient transfection experiment in HEL cells was performed. Deletion analysis of the AQP1 promoter revealed that one putative DR5-like RARE with five spaces was located in the region from -2218 to -2202; AGGGCAgggacAGGTGA. Electrophoretic mobility shift assay (EMSA) experiment demonstrated that two slowly migrated complexes (C1 and C2) capable of binding the RARE sequence were present in nuclear extracts prepared from cells and the complex C1 was strongly increased in nuclear extracts by RA stimulation. The complexes C1 and C2 were significantly abolished by an excess unlabeled probe. These results indicate that RAs strongly stimulate the human AQP1 gene expression through the RARE and define a novel role in the regulation of erythropoiesis.     

A refined model for the solution structure of oxidized putidaredoxin

A refined model for the solution structure of oxidized putidaredoxin (Pdxo), a Cys4Fe2S2 ferredoxin, has been determined. A previous structure (Pochapsky et al. (1994) Biochemistry 33, 6424-6432; PDB entry ) was calculated using the results of homonuclear two-dimensional NMR experiments. New data has made it possible to calculate a refinement of the original Pdxo solution structure. First, essentially complete assignments for diamagnetic 15N and 13C resonances of Pdxo have been made using multidimensional NMR methods, and 15N- and 13C-resolved NOESY experiments have permitted the identification of many new NOE restraints for structural calculations. Stereospecific assignments for leucine and valine CH3 resonances were made using biosynthetically directed fractional 13C labeling, improving the precision of NOE restraints involving these residues. Backbone dihedral angle restraints have been obtained using a combination of two-dimensional J-modulated 15N,1H HSQC and 3D (HN)CO(CO)NH experiments. Second, the solution structure of a diamagnetic form of Pdx, that of the C85S variant of gallium putidaredoxin, in which a nonligand Cys is replaced by Ser, has been determined (Pochapsky et al. (1998) J. Biomol. NMR 12, 407-415), providing information concerning structural features not observable in the native ferredoxin due to paramagnetism. Third, a crystal structure of a closely related ferredoxin, bovine adrenodoxin, has been published (Müller et al. (1998) Structure 6, 269-280). This structure has been used to model the metal binding site structure in Pdx. A family of fourteen structures is presented that exhibits an rmsd of 0.51 A for backbone heavy atoms and 0.83 A for all heavy atoms. Exclusion of the modeled metal binding loop region reduces overall the rmsd to 0.30 A for backbone atoms and 0.71 A for all heavy atoms.     

Distribution of sodium transporters and aquaporin-1 in the human choroid plexus

The choroid plexus epithelium secretes electrolytes and fluid in the brain ventricular lumen at high rates. Several channels and ion carriers have been identified as likely mediators of this transport in rodent choroid plexus. This study aimed to map several of these proteins to the human choroid plexus. Immunoperoxidase-histochemistry was employed to determine the cellular and subcellular localization of the proteins. The water channel, aquaporin (AQP) 1, was predominantly situated in the apical plasma membrane domain, although distinct basolateral and endothelial immunoreactivity was also observed. The Na⁺-K⁺-ATPase ₁-subunit was exclusively localized apically in the human choroid plexus epithelial cells. Immunoreactivity for the Na⁺-K⁺-2Cl^– cotransporter, NKCC1, was likewise confined to the apical plasma membrane domain of the epithelium. The Cl^–/HCO₃^– exchanger, AE2, was localized basolaterally, as was the Na⁺-dependent Cl^–/HCO₃^– exchanger, NCBE, and the electroneutral Na⁺-HCO₃^– cotransporter, NBCn1. No immunoreactivity was found toward the Na⁺-dependent acid/base transporters NHE1 or NBCe2. Hence, the human choroid plexus epithelium displays an almost identical distribution pattern of water channels and Na⁺ transporters as the rat and mouse choroid plexus. This general cross species pattern suggests central roles for these transporters in choroid plexus functions such as cerebrospinal fluid production. immunohistochemistry; metabolism; cerebrospinal fluid secretion     

Gene structure of human apolipoprotein A1.

Apolipoprotein A1 is the major polypeptide of the human plasma high density lipoprotein (HDL). The structure and function of the apo A1 gene are of interest because of the inverse correlation shown between HDL levels and coronary heart disease. We have determined the nucleotide sequence of the apo A1 gene previously isolated. Its coding sequence is interrupted by two introns of 185 and 588 base pairs. As there may be one or more introns in the 5' non coding region, the 5' end of the gene could not be precisely located. The human apo A1 has an unusual propeptide segment very similar to its rat counterpart. The data reported here provide an essential basis for future studies of structural and functional alleles of the apo A1 gene.     

The structure of cyclolinopeptide A in chloroform refined by RDC measurements

Three‐dimensional structures of molecules traditionally assigned from nuclear Overhauser effects and vicinal coupling constants are recently complemented by measurements of residual dipolar couplings. Residual dipolar couplings measured in a stretched poly(dimethylsiloxane) gel were used to determine the structure of cyclolinopeptide A in chloroform solution at ?50 °C. After structure refinement, conformational details of main cluster were discussed in relation to crystal and nuclear Overhauser effect derived structures. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

The crystal structure of staphylococcal nuclease refined at 1.7 A resolution

The crystal structure of staphylococcal nuclease has been determined to 1.7 A resolution with a final R-factor of 16.2% using stereochemically restrained Hendrickson-Konnert least-squares refinement. The structure reveals a number of conformational changes relative to the structure of the ternary complex of staphylococcal nuclease 1,2 bound with deoxythymidine-3',5'-diphosphate and Ca2+. Tyr-113 and Tyr-115, which pack against the nucleotide base in the nuclease complex, are rotated outward creating a more open binding pocket in the absence of nucleotide. The side chains of Ca2+ ligands Asp-21 and Asp-40 shift as does Glu-43, the proposed general base in the hydrolysis of the 5'-phosphodiester bond. The significance of some changes in the catalytic site is uncertain due to the intrusion of a symmetry related Lys-70 side chain which hydrogen bonds to both Asp-21 and Glu-43. The position of a flexible loop centered around residue 50 is altered, most likely due to conformational changes propagated from the Ca2+ site. The side chains of Arg-35, Lys-84, Tyr-85, and Arg-87, which hydrogen bond to the 3'- and 5'-phosphates of the nucleotide in the nuclease complex, are unchanged in conformation, with packing interactions with adjacent protein side chains sufficient to fix the geometry in the absence of ligand. The nuclease structure presented here, in combination with the stereochemically restrained refinement of the nuclease complex structure at 1.65 A, provides a wealth of structural information for the increasing number of studies using staphylococcal nuclease as a model system of protein structure and function.     

The refined structure of the epsilon-aminocaproic acid complex of human plasminogen kringle 4.

The crystallographic structure of the plasminogen kringle 4-epsilon-aminocaproic acid (ACA) complex (K4-ACA) has been solved by molecular replacement rotation-translation methods utilizing the refined apo-K4 structure as a search model (Mulichak et al., 1991), and it has been refined to an R value of 0.148 at 2.25-A resolution. The K4-ACA structure consists of two interkringle residues, the kringle along with the ACA ligand, and 106 water molecules. The lysine-binding site has been confirmed to be a relatively open and shallow depression, lined by aromatic rings of Trp62, Phe64, and Trp72, which provide a highly nonpolar environment between doubly charged anionic and cationic centers formed by Asp55/Asp57 and Lys35/Arg71. A zwitterionic ACA ligand molecule is held by hydrogen-bonded ion pair interactions and van der Waals contacts between the charged centers. The lysine-binding site of apo-K4 and K4-ACA have been compared: the rms differences in main-chain and side-chain positions are 0.25 and 0.69 A, respectively, both practically within error of the determinations. The largest deviations in the binding site are due to different crystal packing interactions. Thus, the lysine-binding site appears to be preformed, and lysine binding does not require conformational changes of the host. The results of NMR studies of lysine binding with K4 are correlated with the structure of K4-ACA and agree well.     

The refined crystal structure of subtilisin Carlsberg at 2.5 A resolution

We report here the X-ray crystal structure of native subtilisin Carlsberg, solved at 2.5 A resolution by molecular replacement and refined by restrained least squares to a crystallographic residual (Formula see text): of 0.206. we compare this structure to the crystal structure of subtilisin BPN'. We find that, despite 82 amino acid substitutions and one deletion in subtilisin Carlsberg relative to subtilisin BPN', the structures of these enzymes are remarkably similar. We calculate an r.m.s. difference between equivalent alpha-carbon positions in subtilisin Carlsberg and subtilisin BPN' of only 0.55 A. This confirms previous reports of extensive structural homology between these two subtilisins based on X-ray crystal structures of the complex of eglin-c with subtilisin Carlsberg [McPhalen, C.A., Schnebli, H.P. and James, M.N.G. (1985) FEBS Lett., 188, 55; Bode, W., Papamokos, E. and Musil, D. (1987) Eur. J. Biochem., 166, 673-692]. In addition, we find that the native active sites of subtilisins Carlsberg and BPN' are virtually identical. While conservative substitutions at residues 217 and 156 may have subtle effects on the environments of substrate-binding sites S1' and S1 respectively, we find no obvious structural correlate for reports that subtilisins Carlsberg and BPN' differ in their recognition of model substrates. In particular, we find no evidence that the hydrophobic binding pocket S1 in subtilisin Carlsberg is 'deeper', 'narrower' or 'less polar' than the corresponding binding site in subtilisin BPN'.     

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.     

Crystal structure of human rhinovirus serotype 1A (HRV1A)

The structure of human rhinovirus 1A (HRV1A) has been determined to 3.2 A resolution using phase refinement and extension by symmetry averaging starting with phases at 5 A resolution calculated from the known human rhinovirus 14 (HRV14) structure. The polypeptide backbone structures of HRV1A and HRV14 are similar, but the exposed surfaces are rather different. Differential charge distribution of amino acid residues in the "canyon", the putative receptor binding site, provides a possible explanation for the difference in minor versus major receptor group specificities, represented by HRV1A and HRV14, respectively. The hydrophobic pocket in VP1, into which antiviral compounds bind, is in an "open" conformation similar to that observed in drug-bound HRV14. Drug binding in HRV1A does not induce extensive conformational changes, in contrast to the case of HRV14.     

A refined accuracy index to evaluate algorithms of protein secondary structure prediction

Nowadays even a 1% increase of the accuracy for the secondary structure prediction is considered remarkable progress. In this case, we have to consider the reasonableness of the accuracy index Q3, which is used widely. A refined accuracy index, called Q8, is proposed to evaluate algorithms of secondary structure prediction. It is shown that Q8 is superior to the widely used index Q3 in that the former carries more information of the predictive accuracy matrix than does the latter. Therefore, algorithms are evaluated more objectively by Q8 than Q3. Based on 396 nonhomologous proteins, five currently available algorithms of secondary structure prediction were evaluated and compared using the new index Q8. Of the five algorithms, PHD turned out to be the unique algorithm, with Q8 accuracy better than 70%. It is suggested that Q3 should be replaced by Q8 in evaluating secondary structure prediction in future studies.     

The refined 2.3 A crystal structure of human leukocyte elastase in a complex with a valine chloromethyl ketone inhibitor

The stoichiometric complex formed between human leukocyte elastase and a synthetic MeO-Suc-Ala-Ala-Pro-Val chloromethyl ketone inhibitor was co-crystallized and its X-ray structure determined, using Patterson search methods. Its structure has been crystallographically refined to a final R value of 0.145 (8.0 and 2.3 A). The enzyme structure is very similar to that recently observed in a complex formed with the ovomucoid third domain from turkey [(1986) EMBO J. 5,2453-2458]. The rms deviation of all alpha-carbon atoms is 0.32 A. The peptidic inhibitor is bound in a similar overall conformation as the ovomucoid binding segment. Covalent bonds are formed between Val-P1 of the inhibitor and His-57 NE2 and Ser-195 OG of the enzyme. The carbonyl carbon is tetrahedrally deformed to a hemiketal. The valine side chain is arranged in the S1 pocket in the g-conformation.     

Yeast-expressed human membrane protein aquaporin-1 yields excellent resolution of solid-state MAS NMR spectra

One of the biggest challenges in solid-state NMR studies of membrane proteins is to obtain a homogeneous natively folded sample giving high spectral resolution sufficient for structural studies. Eukaryotic membrane proteins are especially difficult and expensive targets in this respect. Methylotrophic yeast Pichia pastoris is a reliable producer of eukaryotic membrane proteins for crystallography and a promising economical source of isotopically labeled proteins for NMR. We show that eukaryotic membrane protein human aquaporin 1 can be doubly (¹³C/¹⁵N) isotopically labeled in this system and functionally reconstituted into phospholipids, giving excellent resolution of solid-state magic angle spinning NMR spectra.     

Crystal and molecular structure of human plasminogen kringle 4 refined at 1.9-A resolution.

The crystal structure of human plasminogen kringle 4 (PGK4) has been solved by molecular replacement using the bovine prothrombin kringle 1 (PTK1) structure as a model and refined by restrained least-squares methods to an R factor of 14.2% at 1.9-A resolution. The K4 structure is similar to that of PTK1, and an insertion of one residue at position 59 of the latter has minimal effect on the protein folding. The PGK4 structure is highly stabilized by an internal hydrophobic core and an extensive hydrogen-bonding network. Features new to this kringle include a cis peptide bond at Pro30 and the presence of two alternate, perpendicular, and equally occupied orientations for the Cys75 side chain. The K4 lysine-binding site consists of a hydrophobic trough formed by the Trp62 and Trp72 indole rings, with anionic (Asp55/Asp57) and cationic (Lys35/Arg71) charge pairs at either end. With the adjacent Asp5 and Arg32 residues, these result in triply charged anionic and cationic clusters (pH of crystals at 6.0), which, in addition to the unusually high accessibility of the Trp72 side chain, serve as an obvious marker of the binding site on the K4 surface. A complex intermolecular interaction occurs between the binding sites of symmetry-related molecules involving a highly ordered sulfate anion of solvation in which the Arg32 side chain of a neighboring kringle occupies the binding site.     

SEL1L, the human homolog of C. elegans sel-1: refined physical mapping, gene structure and identification of polymorphic markers

We have cloned the human full-length cDNA SEL1L, which is highly similar to the C. elegans sel-1 gene, an important negative regulator of the "notch" pathway which acts as a key regulator of the cellular proliferation and specification processes in both vertebrates and invertebrates. The SEL1L gene maps to 14q24.3-31 and here we report its fine localization by HAPPY mapping, which determines its molecular distance to microsatellite markers isolated in the region. We have found two new polymorphic (CA)n microsatellites located in the gene, and have identified the exon-intron boundaries. The gene is composed of 21 exons spanning 70 kb of genomic DNA. Human SEL1L protein exhibits a high degree of similarity compared to the mouse and nematode homologs.     

Structure of prothrombin fragment 1 refined at 2.8 A resolution

The structure of prothrombin fragment 1, solved at 2.8 A resolution (1 A = 0.1 nm) by a combination of multiple and single isomorphous replacement methods utilizing solvent flattening, has been refined by restrained least-squares methods (R = 0.24), solvent not included, using fairly stringent restraints on the molecular geometry and individual thermal parameters. The inner kringle loop possesses significantly lower B-values than the outer loops even though the former also constitutes a surface of the folded kringle structure. This surface forms the Lys sub-site of the fibrin binding site of other kringles. The hydrogen bonding network and ion pair interactions of fragment 1 appear to maintain a compact folded structure among the various loops of the kringle structure. On the other hand, since there is only one hydrogen bond between the kringle and its preceding 30 residues, considerable flexibility is suggested for the Gla-domain consistent with its disorder in crystals. A chitobiose has been located at the Asn77 glycosylation site, but only a single N-acetyl-glucosamine is ordered at Asn101. The lysine binding site region of other kringles is not properly developed in fragment 1, accounting for its lack of Lys/fibrin affinity. Most of the conserved sequence among 11 different kringles is associated with either: (1) protecting the inner loop disulfides Cys87-127, Cys115-139 upon which the folding is based; or (2) a requirement of the lysine binding site. The remainder of the conservation is generally associated with the ten reverse turns of the folding; of these 40 residues, or about half the sequence, 14 are conserved among eight different turns. The intermolecular packing consists of infinite helical columns of fragment 1 molecules related by a crystallographic 4(3) screw axis, which are held together by van der Waals' interactions of aromatic clusters from different molecules related by a crystallographic 2-fold rotation axis.