Electrostatics of nanosystems: application to microtubules and the ribosome

Electrostatic potential of aminoacyl-tRNA synthetase navigates tRNA on its pathway to the binding site

In the first stage of a diffusion-controlled enzymatic reaction, aminoacyl-tRNA synthetases (aaRSs) interact with cognate tRNAs forming non-specific encounters. The aaRSs catalyzing the same overall aminoacylation reaction vary greatly in subunit organization, structural domain composition and amino acid sequence. The diffusional association of aaRS and tRNA was found to be governed by long-range electrostatic interactions when the homogeneous negative potential of tRNA fits to the patches of positive potential produced by aaRS; one patch for each tRNA substrate molecule. Considering aaRS as a molecule with anisotropic reactivity and on the basis of continuum electrostatics and Smoluchowski's theory, the reaction conditions for tRNA-aaRS diffusional encounters were formulated. The domains, categorized as enzymatically relevant, appeared to be non-essential for field sculpturing at long distances. On the other hand, a set of complementary domains exerts primary control on the aaRS isopotential surface formation. Subdividing the aaRS charged residues into native, conservative and non-conservative subsets, we evaluated the contribution of each group to long-range electrostatic potential. Surprisingly, the electrostatic potential landscapes generated by native and non-conservative subsets are fairly similar, thus suggesting the non-conservative subset is developed specifically for efficient tRNA attraction.

Structural basis of ARNT PAS-B dimerization: use of a common beta-sheet interface for hetero- and homodimerization

The aryl hydrocarbon receptor nuclear translocator (ARNT) is a promiscuous bHLH-PAS (Per-ARNT-Sim) protein that forms heterodimeric transcriptional regulator complexes with several other bHLH-PAS subunits to control a variety of biological pathways, some of which are centrally involved in disease initiation and/or progression. One of these is the hypoxia response pathway, which allows eukaryotic cells to respond to low oxygen tension via the formation of a heterodimeric complex between ARNT and another bHLH-PAS protein, the hypoxia-inducible factor alpha (HIF-alpha). We have previously shown that the C-terminal PAS domains of an HIF-alpha isoform (HIF-2alpha) and ARNT interact in vitro, and that mutations in the solvent-exposed beta-sheet surface of the HIF-2alpha domain not only disrupt this interaction, but also greatly attenuate the hypoxia response in living cells. Here, we have solved the solution structure of the corresponding PAS domain of ARNT and show that it utilizes a very similar interface for the interaction with the HIF-2alpha PAS domain. We also show that this domain self-associates in a concentration-dependent manner, and that the interface used in this homodimeric complex is very similar to that used in the formation of heterodimer. In addition, using experimentally derived NMR restraints, we used the program HADDOCK to calculate a low-resolution model of the complex formed in solution by these two PAS domains, and confirm the validity of this model using site-directed spin labeling to obtain long-range distance information in solution. With this information, we propose a model for the mode of multi-PAS domain interaction in bHLH-PAS transcriptional activation complexes.

Molecular analysis of anti-N-propionyl Neisseria meningitidis group B polysaccharide monoclonal antibodies

The capsular polysaccharide of Neisseria meningitidis group B (MBPS) is a polymer of alpha (2-->8) N-acetyl neuraminic acid, which is chemically identical to polysialic acid (PSA) expressed in human tissues. Antibodies from mice immunized with a MBPS-protein conjugate vaccine in which N-acetyl groups have been replaced by propionyl groups (N-Pr MBPS) can be bactericidal and show minimal or no cross-reactivity with human PSA. To investigate the molecular basis for antigen recognition, we cloned and sequenced the variable region (V) genes of five bactericidal anti-N-Pr MBPS murine mAbs and produced computer models of the combining sites. The results were compared to those reported in the literature for two autoreactive anti-MBPS. The V region genes of the anti-N-Pr MBPS mAbs and the anti-MBPS autoreactive mAbs are derived from a limited set of germline V, J, and D genes. However, the anti-N-Pr MBPS mAbs are more mutated than the anti-MBPS mAbs and the former use V-D-J editing that introduces arginine in H-CDR3. Models of the respective combining sites indicate that the anti-MBPS or anti-N-Pr MBPS mAbs that react with host PSA have relatively wide and shallow grooves with a high overall positive charge, consistent with recognition of extended helical polysaccharide structures recognized by the autoreactive mAbs. In contrast, anti-N-Pr MBPS mAbs that do not react with host PSA contain pockets and deep clefts that are consistent with recognition of discrete structural features of individual residues.

Heparin-binding Sites in Granulocyte-Macrophage Colony-stimulating Factor

The biological activity of granulocyte-macrophage colony-stimulating factor (GM-CSF) is modulated by the sulfated glycosaminoglycans (GAGs) heparan sulfate and heparin. However, the molecular mechanisms involved in such interactions are still not completely understood. We have proposed previously that helix C, one of the four alpha-helices of human GM-CSF (hGM-CSF), contains a GAG-binding site in which positively charged residues are spatially positioned for interaction with the sulfate moieties of the GAGs (Wettreich, A., Sebollela, A., Carvalho, M. A., Azevedo, S. P., Borojevic, R., Ferreira, S. T., and Coelho-Sampaio, T. (1999) J. Biol. Chem. 274, 31468-31475). Protonation of two histidine residues (His83 and His87) in helix C of hGM-CSF appears to act as a pH-dependent molecular switch to control the interaction with GAGs. Based on these findings, we have now generated a triple mutant form of murine GM-CSF (mGM-CSF) in which three noncharged residues in helix C of the murine factor (Tyr83, Gln85, and Tyr87) were replaced by the corresponding basic residues present in hGM-CSF (His83, Lys85, and His87). Binding assays on heparin-Sepharose showed that, at acidic pH, the triple mutant mGM-CSF binds to immobilized heparin with significantly higher affinity than wild type (WT) mGM-CSF and that neither protein binds to the column at neutral pH. The fact that even WT mGM-CSF binds to heparin at acidic pH indicates the existence of a distinct, lower affinity heparin-binding site in the protein. Chemical modification of the single histidine residue (His15) located in helix A of WT mGM-CSF with diethyl pyrocarbonate totally abolished binding to immobilized heparin. Moreover, replacement of His15 for an alanine residue significantly reduced the affinity of mGM-CSF for heparin at pH 5.0 and completely blocked heparin binding to a synthetic peptide corresponding to helix A of GM-CSF. These results indicate a major role of histidine residues in the regulation of the binding of GM-CSF to GAGs, supporting the notion that an acidic microenvironment is required for GM-CSF-dependent regulation of target cells. In addition, our results provide insight into the molecular basis of the strict species specificity of the biological activity of GM-CSF.

The structure at 2.5 A resolution of human basophilic leukemia-expressed protein BLES03

The crystal structure of the human basophilic leukemia-expressed protein (BLES03, p5326, Hs.433573) was determined by single-wavelength anomalous diffraction and refined to an R factor of 18.8% (Rfree = 24.5%) at 2.5 A resolution. BLES03 shows no detectable sequence similarity to any functionally characterized proteins using state-of-the-art sequence-comparison tools. The structure of BLES03 adopts a fold similar to that of eukaryotic transcription initiation factor 4E (eIF4E), a protein involved in the recognition of the cap structure of eukaryotic mRNA. In addition to fold similarity, the electrostatic surface potentials of BLES03 and eIF4E show a clear conservation of basic and acidic patches. In the crystal lattice, the acidic amino-terminal helices of BLES03 monomers are bound within the basic cavity of symmetry-related monomers in a manner analogous to the binding of mRNA by eIF4E. Interestingly, the gene locus encoding BLES03 is located between genes encoding the proteins DRAP1 and FOSL1, both of which are involved in transcription initiation. It is hypothesized that BLES03 itself may be involved in a biochemical process that requires recognition of nucleic acids.

Modeling Tubulin at Interfaces. Immobilization of Microtubules on Self-Assembled Monolayers

A theoretical study of protein docking to self-assembled monolayers using a new approach is presented. Docking experiments based on space complementarity implemented in FTDock software were performed for three different proteins: tubulin dimer, cytochrome c, and lysozyme. The proteins were adsorbed on alkanethiol surfaces with different terminating groups and 50,000 best orientations of each protein were analyzed. For all systems three filters based on different chemical and biological approaches were applied. Correctly docked proteins for the cytochrome c and lysozyme systems were found in a list of the first 12 results after applying the geometrical and grouping filter and in a list of the first 3 results after applying the biological filter. We have found that alkanethiol monolayers with odd and even numbers of -CH2- groups have similar properties in terms of interactions with the two proteins. Docking of the tubulin dimer revealed that the orientation favored from the applicational point of view can be found in a list of the first 14 results for monolayers with different terminating groups and that there may be a noticeable difference in tubulin dimer interactions with alkanethiol chains of various length. The results for tubulin dimer docking combined with microtubules ability of reversible assembly suggest that these biological structures may become good candidates to serve as templates for fabrication of nanowires and other nanoscale electronic devices. The new method of theoretical docking presented may be used as a fast and reliable tool complementing other theoretical and experimental techniques of exploring other protein-surface interfaces.

Compressed representations of macromolecular structures and properties

We introduce a new and unified, compressed volumetric representation for macromolecular structures at varying feature resolutions, as well as for many computed associated properties. Important caveats of this compressed representation are fast random data access and decompression operations. Many computational tasks for manipulating large structures, including those requiring interactivity such as real-time visualization, are greatly enhanced by utilizing this compact representation. The compression scheme is obtained by using a custom designed hierarchical wavelet basis construction. Due to the continuity offered by these wavelets, we retain very good accuracy of molecular surfaces, at very high compression ratios, for macromolecular structures at multiple resolutions.

Radial compression of microtubules and the mechanism of action of taxol and associated proteins

Microtubules (MTs) are hollow cylindrical polymers composed of alphabeta-tubulin heterodimers that align head-to-tail in the MT wall, forming linear protofilaments that interact laterally. We introduce a probe of the interprotofilament interactions within MTs and show that this technique gives insight into the mechanisms by which MT-associated proteins (MAPs) and taxol stabilize MTs. In addition, we present further measurements of the mechanical properties of MT walls, MT-MT interactions, and the entry of polymers into the MT lumen. These results are obtained from a synchrotron small angle x-ray diffraction (SAXRD) study of MTs under osmotic stress. Above a critical osmotic pressure, P(cr), we observe rectangular bundles of MTs whose cross sections have buckled to a noncircular shape; further increases in pressure continue to distort MTs elastically. The P(cr) of approximately 600 Pa provides, for the first time, a measure of the bending modulus of the interprotofilament bond within an MT. The presence of neuronal MAPs greatly increases P(cr), whereas surprisingly, the cancer chemotherapeutic drug taxol, which suppresses MT dynamics and inhibits MT depolymerization, does not affect the interprotofilament interactions. This SAXRD-osmotic stress technique, which has enabled measurements of the mechanical properties of MTs, should find broad application for studying interactions between MTs and of MTs with MAPs and MT-associated drugs.

PCE: web tools to compute protein continuum electrostatics

PCE (protein continuum electrostatics) is an online service for protein electrostatic computations presently based on the MEAD (macroscopic electrostatics with atomic detail) package initially developed by D. Bashford [(2004) Front Biosci., 9, 1082–1099]. This computer method uses a macroscopic electrostatic model for the calculation of protein electrostatic properties, such as pK_a values of titratable groups and electrostatic potentials. The MEAD package generates electrostatic energies via finite difference solution to the Poisson–Boltzmann equation. Users submit a PDB file and PCE returns potentials and pK_a values as well as color (static or animated) figures displaying electrostatic potentials mapped on the molecular surface. This service is intended to facilitate electrostatics analyses of proteins and thereby broaden the accessibility to continuum electrostatics to the biological community. PCE can be accessed at .

Structure of the Flavivirus helicase: implications for catalytic activity, protein interactions, and proteolytic processing

Yellow fever virus (YFV), a member of the Flavivirus genus, has a plus-sense RNA genome encoding a single polyprotein. Viral protein NS3 includes a protease and a helicase that are essential to virus replication and to RNA capping. The 1.8-A crystal structure of the helicase region of the YFV NS3 protein includes residues 187 to 623. Two familiar helicase domains bind nucleotide in a triphosphate pocket without base recognition, providing a site for nonspecific hydrolysis of nucleoside triphosphates and RNA triphosphate. The third, C-terminal domain has a unique structure and is proposed to function in RNA and protein recognition. The organization of the three domains indicates that cleavage of the viral polyprotein NS3-NS4A junction occurs in trans.

Structure and function of a potent agonist for the semi-invariant natural killer T cell receptor

Natural killer T cells express a conserved, semi-invariant alphabeta T cell receptor that has specificity for self glycosphingolipids and microbial cell wall alpha-glycuronosylceramide antigens presented by CD1d molecules. Here we report the crystal structure of CD1d in complex with a short-chain synthetic variant of alpha-galactosylceramide at a resolution of 2.2 A. This structure elucidates the basis for the high specificity of these microbial ligands and explains the restriction of the alpha-linkage as a unique pathogen-specific pattern-recognition motif. Comparison of the binding of altered lipid ligands to CD1d and T cell receptors suggested that the differential T helper type 1-like and T helper type 2-like properties of natural killer T cells may originate largely from differences in their 'loading' in different cell types and hence in their tissue distribution in vivo.

Molecular dynamics of a protein surface: ion-residues interactions

Time-resolved measurements indicated that protons could propagate on the surface of a protein or a membrane by a special mechanism that enhanced the shuttle of the proton toward a specific site. It was proposed that a suitable location of residues on the surface contributes to the proton shuttling function. In this study, this notion was further investigated by the use of molecular dynamics simulations, where Na(+) and Cl(-) are the ions under study, thus avoiding the necessity for quantum mechanical calculations. Molecular dynamics simulations were carried out using as a model a few Na(+) and Cl(-) ions enclosed in a fully hydrated simulation box with a small globular protein (the S6 of the bacterial ribosome). Three independent 10-ns-long simulations indicated that the ions and the protein's surface were in equilibrium, with rapid passage of the ions between the protein's surface and the bulk. However, it was noted that close to some domains the ions extended their duration near the surface, thus suggesting that the local electrostatic potential hindered their diffusion to the bulk. During the time frame in which the ions were detained next to the surface, they could rapidly shuttle between various attractor sites located under the electrostatic umbrella. Statistical analysis of the molecular dynamics and electrostatic potential/entropy consideration indicated that the detainment state is an energetic compromise between attractive forces and entropy of dilution. The similarity between the motion of free ions next to a protein and the proton transfer on the protein's surface are discussed.

The Structural Basis for Substrate Specificity in DNA Topoisomerase IV

Most bacteria possess two type IIA topoisomerases, DNA gyrase and topo IV, that together help manage chromosome integrity and topology. Gyrase primarily introduces negative supercoils into DNA, an activity mediated by the C-terminal domain of its DNA binding subunit (GyrA). Although closely related to gyrase, topo IV preferentially decatenates DNA and relaxes positive supercoils. Here we report the structure of the full-length Escherichia coli ParC dimer at 3.0 A resolution. The N-terminal DNA binding region of ParC is highly similar to that of GyrA, but the ParC dimer adopts a markedly different conformation. The C-terminal domain (CTD) of ParC is revealed to be a degenerate form of the homologous GyrA CTD, and is anchored to the top of the N-terminal domains in a configuration different from that thought to occur in gyrase. Biochemical assays show that the ParC CTD controls the substrate specificity of topo IV, likely by capturing DNA segments of certain crossover geometries. This work delineates strong mechanistic parallels between topo IV and gyrase, while explaining how structural differences between the two enzyme families have led to distinct activity profiles. These findings in turn explain how the structures and functions of bacterial type IIA topoisomerases have evolved to meet specific needs of different bacterial families for the control of chromosome superstructure.

Dimerization of doublesex is mediated by a cryptic ubiquitin-associated domain fold: implications for sex-specific gene regulation

Male- and female-specific isoforms of the Doublesex (DSX) transcription factor regulate somatic sexual differentiation in Drosophila. The isoforms (DSX(M) and DSX(F)) share an N-terminal DNA binding domain (the DM motif), broadly conserved among metazoan sex-determining pathways. DM-DNA recognition is enhanced by a C-terminal dimerization domain. The crystal structure of this domain, determined at a resolution of 1.6 A, reveals a novel dimeric arrangement of ubiquitin-associated (UBA) folds. Although this alpha-helical motif is well characterized in pathways of DNA repair and subcellular trafficking, to our knowledge this is its first report in a transcription factor. Dimerization is mediated by a non-canonical hydrophobic interface extrinsic to the putative ubiquitin binding surface. Key side chains at this interface, identified by alanine scanning mutagenesis, are conserved among DSX homologs. The mechanism of dimerization is thus unrelated to the low affinity domain swapping observed among ubiquitin-associated CUE domains. The unexpected observation of a ubiquitin-associated fold in DSX extends the repertoire of alpha-helical dimerization elements in transcription factors. The possibility that the ubiquitination machinery participates in the regulation of sexual dimorphism is discussed.

Autophosphorylation of Archaeoglobus fulgidus Rio2 and crystal structures of its nucleotide-metal ion complexes

The highly conserved, atypical RIO serine protein kinases are found in all organisms, from archaea to man. In yeast, the kinase activity of Rio2 is necessary for the final processing step of maturing the 18S ribosomal rRNA. We have previously shown that the Rio2 protein from Archaeoglobus fulgidus contains both a small kinase domain and an N-terminal winged helix domain. Previously solved structures using crystals soaked in nucleotides and Mg2+ or Mn2+ showed bound nucleotide but no ordered metal ions, leading us to the conclusion that they did not represent an active conformation of the enzyme. To determine the functional form of Rio2, we crystallized it after incubation with ATP or ADP and Mn2+. Co-crystal structures of Rio2-ATP-Mn and Rio2-ADP-Mn were solved at 1.84 and 1.75 angstroms resolution, respectively. The gamma-phosphate of ATP is firmly positioned in a manner clearly distinct from its location in canonical serine kinases. Comparison of the Rio2-ATP-Mn complex with the Rio2 structure with no added nucleotides and with the ADP complex indicates that a flexible portion of the Rio2 molecule becomes ordered through direct interaction between His126 and the gamma-phosphate oxygen of ATP. Phosphopeptide mapping of the autophosphorylation site of Rio2 identified Ser128, within the flexible loop and directly adjacent to the part that becomes ordered in response to ATP, as the target. These results give us further information about the nature of the active site of Rio2 kinase and suggest a mechanism of regulation of its enzymatic activity.

Exploring Assembly Energetics of the 30S Ribosomal Subunit Using an Implicit Solvent Approach

To explore the relationship between the assembly of the 30S ribosomal subunit and interactions among the constituent components, 16S RNA and proteins, relative binding free energies of the T. thermophilus 30S proteins to the 16S RNA were studied based on an implicit solvent model of electrostatic, nonpolar, and entropic contributions. The late binding proteins in our assembly map were found not to bind to the naked 16S RNA. The 5' domain early kinetic class proteins, on average, carry the highest positive charge, get buried the most upon binding to 16S RNA, and show the most favorable binding. Some proteins (S10/S14, S6/S18, S13/S19) have more stabilizing interactions while binding as dimers. Our computed assembly map resembles that of E. coli; however, the central domain path is more similar to that of A. aeolicus, a hyperthermophilic bacteria.

Small-angle X-ray scattering reveals the solution structure of the full-length DNA gyrase a subunit

DNA gyrase is the topoisomerase uniquely able to actively introduce negative supercoils into DNA. Vital in all bacteria, but absent in humans, this enzyme is a successful target for antibacterial drugs. From biophysical experiments in solution, we report the low-resolution structure of the full-length A subunit (GyrA). Analytical ultracentrifugation shows that GyrA is dimeric, but nonglobular. Ab initio modeling from small-angle X-ray scattering allows us to retrieve the molecular envelope of GyrA and thereby the organization of its domains. The available crystallographic structure of the amino-terminal domain (GyrA59) forms a dimeric core, and two additional pear-shaped densities closely flank it in an unexpected position. Each accommodates very well a carboxyl-terminal domain (GyrA-CTD) built from a homologous crystallographic structure. The uniqueness of gyrase is due to the ability of the GyrA-CTDs to wrap DNA. Their position within the GyrA structure strongly suggests a large conformation change of the enzyme upon DNA binding.

Structure of the zinc-binding domain of an essential component of the hepatitis C virus replicase

Hepatitis C virus (HCV) is a human pathogen affecting nearly 3% of the world's population. Chronic infections can lead to cirrhosis and liver cancer. The RNA replication machine of HCV is a multi-subunit membrane-associated complex. The non-structural protein NS5A is an active component of HCV replicase, as well as a pivotal regulator of replication and a modulator of cellular processes ranging from innate immunity to dysregulated cell growth. NS5A is a large phosphoprotein (56-58 kDa) with an amphipathic alpha-helix at its amino terminus that promotes membrane association. After this helix region, NS5A is organized into three domains. The N-terminal domain (domain I) coordinates a single zinc atom per protein molecule. Mutations disrupting either the membrane anchor or zinc binding of NS5A are lethal for RNA replication. However, probing the role of NS5A in replication has been hampered by a lack of structural information about this multifunctional protein. Here we report the structure of NS5A domain I at 2.5-A resolution, which contains a novel fold, a new zinc-coordination motif and a disulphide bond. We use molecular surface analysis to suggest the location of protein-, RNA- and membrane-interaction sites.

Crystal Structure of Dynein Light Chain TcTex-1

TcTex-1, one of three dynein light chains of the dynein motor complex, has been implicated in targeting and binding cargoes to cytoplasmic dynein for retrograde or apical transport. Interactions between TcTex-1 and a diverse set of proteins such as the dynein intermediate chain, Fyn, DOC2, FIP1, the poliovirus receptor, CD155, and the rhodopsin cytoplasmic tail have been reported; yet, despite the broad range of targets, a consensus binding sequence remains uncertain. Consequently, we have solved the crystal structure of the full-length Drosophila homolog of TcTex-1 to 1.7 A resolution using MAD phasing to gain insight into its function and target specificity. The structure is homodimeric with a domain swapping of beta-strand 2 and has a fold similar to the dynein light chain, LC8. Based on structural alignment, the TcTex-1 and LC8 sequences show no identity, although the root mean square deviation between secondary structural elements is less than 1.6 A. Moreover, the N terminus, which is equivalent to beta-strand 1 in LC8, is splayed out and binds to a crystallographic dimer as an anti-parallel beta-strand at the same position as the neuronal nitric-oxide synthase peptide in the LC8 complex. Similarity to LC8 and comparison to the LC8-neuronal nitricoxide synthase complex suggest that TcTex-1 binds its targets in a similar manner as LC8 and provides insight to the lack of strict sequence identity among the targets for TcTex-1.

Universal method for computation of electrostatic potentials

A computational approach to determine electrostatic interaction and gravitational potentials by performing direct numerical integration is presented. The potential is expanded using finite-element functions of arbitrary order. The method does not involve any solutions of systems of linear equations. The potential is instead obtained as a sum of differential contributions. Thus, no boundary conditions for the potential are needed. It is computationally efficient and well suited for parallel computers, since the innermost loops constitute matrix multiplications and the outer ones can be used as parallel indices. Without using prescreening or other computational tricks to speed up the calculation, the algorithm scales as N4/3 where N denotes the grid size.

The Mechanism of Proton Transfer between Adjacent Sites Exposed to Water

The surface of a protein, or a membrane, is spotted with a multitude of proton-binding sites, some of which are only a few angstroms apart. When a proton is released from one site, it propagates through the water by a random walk under the bias of the local electrostatic potential determined by the distribution of the charges on the protein. Some of the released protons disperse into the bulk, but during the first few nanoseconds, the released protons can be trapped by encounter with nearby acceptor sites. This process resembles a scenario which corresponds with the time-dependent Debye-Smoluchowski equation. In the present study, we investigated the mechanism of proton transfer between sites that are only a few angstroms apart, using as a model the proton exchange between sites on a small molecule, fluorescein, having two, spectrally distinguishable, proton-binding sites. The first site is the oxyanion on the chromophore ring structure. The second site is the carboxylate moiety on the benzene ring of the molecule. Through our experiments, we were able to reconstruct the state of protonation at each site and the velocity of proton transfer between them. The fluorescein was protonated by a few nanosecond long proton pulse under specific conditions that ensured that the dye molecules would be protonated only by a single proton. The dynamics of the protonation of the chromophore were measured under varying initial conditions (temperature, ionic strength, and different solvents (H(2)O or D(2)O)), and the velocity of the proton transfer between the two sites was extracted from the overall global analysis of the signals. The dynamics of the proton transfer between the two proton-binding sites of the fluorescein indicated that the efficiency of the site-to-site proton transfer is very sensitive to the presence of the screening electrolyte and has a very high kinetic isotope effect (KIE = 55). These two parameters clearly distinguish the mechanism from proton diffusion in bulk water. The activation energy of the reaction (E(a) = 11 kcal mol(-1)) is also significantly higher than the activation energy for proton dissociation in bulk water (E(a) approximately 2.5 kcal mol(-1)). These observations are discussed with respect to the effect of the solute on the water molecules located within the solvation layer.

Crystal structure of the RNA 2'-phosphotransferase from Aeropyrum pernix K1

In the final step of tRNA splicing, the 2'-phosphotransferase catalyzes the transfer of the extra 2'-phosphate from the precursor-ligated tRNA to NAD. We have determined the crystal structure of the 2'-phosphotransferase protein from Aeropyrum pernix K1 at 2.8 Angstroms resolution. The structure of the 2'-phosphotransferase contains two globular domains (N and C-domains), which form a cleft in the center. The N-domain has the winged helix motif, a subfamily of the helix-turn-helix family, which is shared by many DNA-binding proteins. The C-domain of the 2'-phosphotransferase superimposes well on the NAD-binding fold of bacterial (diphtheria) toxins, which catalyze the transfer of ADP ribose from NAD to target proteins, indicating that the mode of NAD binding by the 2'-phosphotransferase could be similar to that of the bacterial toxins. The conserved basic residues are assembled at the periphery of the cleft and could participate in the enzyme contact with the sugar-phosphate backbones of tRNA. The modes by which the two functional domains recognize the two different substrates are clarified by the present crystal structure of the 2'-phosphotransferase.

Crystal structure of pyruvate dehydrogenase kinase 3 bound to lipoyl domain 2 of human pyruvate dehydrogenase complex

The human pyruvate dehydrogenase complex (PDC) is regulated by reversible phosphorylation by four isoforms of pyruvate dehydrogenase kinase (PDK). PDKs phosphorylate serine residues in the dehydrogenase (E1p) component of PDC, but their amino-acid sequences are unrelated to eukaryotic Ser/Thr/Tyr protein kinases. PDK3 binds to the inner lipoyl domains (L2) from the 60-meric transacetylase (E2p) core of PDC, with concomitant stimulated kinase activity. Here, we present crystal structures of the PDK3-L2 complex with and without bound ADP or ATP. These structures disclose that the C-terminal tail from one subunit of PDK3 dimer constitutes an integral part of the lipoyl-binding pocket in the N-terminal domain of the opposing subunit. The two swapped C-terminal tails promote conformational changes in active-site clefts of both PDK3 subunits, resulting in largely disordered ATP lids in the ADP-bound form. Our structural and biochemical data suggest that L2 binding stimulates PDK3 activity by disrupting the ATP lid, which otherwise traps ADP, to remove product inhibition exerted by this nucleotide. We hypothesize that this allosteric mechanism accounts, in part, for E2p-augmented PDK3 activity.

High resolution crystal structures of human Rab4a in its active and inactive conformations

The Ras-related human GTPase Rab4a is involved in the regulation of endocytosis through the sorting and recycling of early endosomes. Towards further insight, we have determined the three-dimensional crystal structure of human Rab4a in its GppNHp-bound state to 1.6 Angstroms resolution and in its GDP-bound state to 1.8 Angstroms resolution, respectively. Despite the similarity of the overall structure with other Rab proteins, Rab4a displays significant differences. The structures are discussed with respect to the recently determined structure of human Rab5a and its complex with the Rab5-binding domain of the bivalent effector Rabaptin-5. The Rab4 specific residue His39 modulates the nucleotide binding pocket giving rise to a reduced rate for nucleotide hydrolysis and exchange. In comparison to Rab5, Rab4a has a different GDP-bound conformation within switch 1 region and displays shifts in position and orientation of the hydrophobic triad. The observed differences at the S2-L3-S3 region represent a new example of structural plasticity among Rab proteins and may provide a structural basis to understand the differential binding of similar effector proteins.

Structural basis for the interaction of Asf1 with histone H3 and its functional implications

Asf1 is a conserved histone chaperone implicated in nucleosome assembly, transcriptional silencing, and the cellular response to DNA damage. We solved the NMR solution structure of the N-terminal functional domain of the human Asf1a isoform, and we identified by NMR chemical shift mapping a surface of Asf1a that binds the C-terminal helix of histone H3. This binding surface forms a highly conserved hydrophobic groove surrounded by charged residues. Mutations within this binding site decreased the affinity of Asf1a for the histone H3/H4 complex in vitro, and the same mutations in the homologous yeast protein led to transcriptional silencing defects, DNA damage sensitivity, and thermosensitive growth. We have thus obtained direct experimental evidence of the mode of binding between a histone and one of its chaperones and genetic data suggesting that this interaction is important in both the DNA damage response and transcriptional silencing.

Relevant distance between two different instances of the same potential energy in protein folding

In the context of complex systems and, particularly, of protein folding, a physically meaningful distance is defined which allows to make useful statistical statements about the way in which energy differences are modified when two different instances of the same potential-energy function are used. When the two instances arise from the fact that different algorithms or different approximations are used, the distance herein defined may be used to evaluate the relative accuracy of the two methods. When the difference is due to a change in the free parameters of which the potential depends on, the distance can be used to quantify, in each region of parameter space, the robustness of the modeling to such a change and this, in turn, may be used to assess the significance of a parameters' fit. Both cases are illustrated with a practical example: the study of the Poisson-based solvation energy in the Trp-Cage protein (PDB code 1L2Y).

The 1.51-Angstrom structure of the poxvirus L1 protein, a target of potent neutralizing antibodies

Although eradicated from nature more than two decades ago, the threat of smallpox has reemerged because of concerns over its use as a biological weapon. We present the structure of the poxvirus L1 protein, a molecule that is conserved throughout the poxvirus family and is nearly identical in vaccinia virus and in variola virus, which causes smallpox. L1 is a myristoylated envelope protein that is a potent target for neutralizing antibodies and an important component of current experimental vaccines. The L1 structure reveals a hydrophobic cavity located adjacent to its N terminus. The cavity would be capable of shielding the myristate moiety, which is essential for virion assembly. The structure of L1 is a step in the elucidation of molecular mechanisms common to all poxviruses that may stimulate the design of safer vaccines and new antipoxvirus drugs.

Conformational energy landscape of the acyl pocket loop in acetylcholinesterase: a Monte Carlo-generalized Born model study

The X-ray crystal structure of the reaction product of acetylcholinesterase (AChE) with the inhibitor diisopropylphosphorofluoridate (DFP) showed significant structural displacement in a loop segment of residues 287-290. To understand this conformational selection, a Monte Carlo (MC) simulation study was performed of the energy landscape for the loop segment. A computational strategy was applied by using a combined simulated annealing and room temperature Metropolis sampling approach with solvent polarization modeled by a generalized Born (GB) approximation. Results from thermal annealing reveal a landscape topology of broader basin opening and greater distribution of energies for the displaced loop conformation, while the ensemble average of conformations at 298 K favored a shift in populations toward the native by a free-energy difference in good agreement with the estimated experimental value. Residue motions along a reaction profile of loop conformational reorganization are proposed where Arg-289 is critical in determining electrostatic effects of solvent interaction versus Coulombic charging.

Structure of an XPF endonuclease with and without DNA suggests a model for substrate recognition

The XPF/Mus81 structure-specific endonucleases cleave double-stranded DNA (dsDNA) within asymmetric branched DNA substrates and play an essential role in nucleotide excision repair, recombination and genome integrity. We report the structure of an archaeal XPF homodimer alone and bound to dsDNA. Superposition of these structures reveals a large domain movement upon binding DNA, indicating how the (HhH)(2) domain and the nuclease domain are coupled to allow the recognition of double-stranded/single-stranded DNA junctions. We identify two nonequivalent DNA-binding sites and propose a model in which XPF distorts the 3' flap substrate in order to engage both binding sites and promote strand cleavage. The model rationalises published biochemical data and implies a novel role for the ERCC1 subunit of eukaryotic XPF complexes.

Continuum diffusion reaction rate calculations of wild-type and mutant mouse acetylcholinesterase: adaptive finite element analysis

As described previously, continuum models, such as the Smoluchowski equation, offer a scalable framework for studying diffusion in biomolecular systems. This work presents new developments in the efficient solution of the continuum diffusion equation. Specifically, we present methods for adaptively refining finite element solutions of the Smoluchowski equation based on a posteriori error estimates. We also describe new, molecular-surface-based models, for diffusional reaction boundary criteria and compare results obtained from these models with the traditional spherical criteria. The new methods are validated by comparison of the calculated reaction rates with experimental values for wild-type and mutant forms of mouse acetylcholinesterase. The results show good agreement with experiment and help to define optimal reactive boundary conditions.

Tetrameric mouse acetylcholinesterase: continuum diffusion rate calculations by solving the steady-state Smoluchowski equation using finite element methods

The tetramer is the most important form for acetylcholinesterase in physiological conditions, i.e., in the neuromuscular junction and the nervous system. It is important to study the diffusion of acetylcholine to the active sites of the tetrameric enzyme to understand the overall signal transduction process in these cellular components. Crystallographic studies revealed two different forms of tetramers, suggesting a flexible tetramer model for acetylcholinesterase. Using a recently developed finite element solver for the steady-state Smoluchowski equation, we have calculated the reaction rate for three mouse acetylcholinesterase tetramers using these two crystal structures and an intermediate structure as templates. Our results show that the reaction rates differ for different individual active sites in the compact tetramer crystal structure, and the rates are similar for different individual active sites in the other crystal structure and the intermediate structure. In the limit of zero salt, the reaction rates per active site for the tetramers are the same as that for the monomer, whereas at higher ionic strength, the rates per active site for the tetramers are approximately 67%-75% of the rate for the monomer. By analyzing the effect of electrostatic forces on ACh diffusion, we find that electrostatic forces play an even more important role for the tetramers than for the monomer. This study also shows that the finite element solver is well suited for solving the diffusion problem within complicated geometries.

Structure and RNA binding of the third KH domain of poly(C)-binding protein 1

Poly(C)-binding proteins (CPs) are important regulators of mRNA stability and translational regulation. They recognize C-rich RNA through their triple KH (hn RNP K homology) domain structures and are thought to carry out their function though direct protection of mRNA sites as well as through interactions with other RNA-binding proteins. We report the crystallographically derived structure of the third domain of αCP1 to 2.1 Å resolution. αCP1-KH3 assumes a classical type I KH domain fold with a triple-stranded β-sheet held against a three-helix cluster in a βααββα configuration. Its binding affinity to an RNA sequence from the 3′-untranslated region (3′-UTR) of androgen receptor mRNA was determined using surface plasmon resonance, giving a K_d of 4.37 μM, which is indicative of intermediate binding. A model of αCP1-KH3 with poly(C)-RNA was generated by homology to a recently reported RNA-bound KH domain structure and suggests the molecular basis for oligonucleotide binding and poly(C)-RNA specificity.

Recognition of the spliceosomal branch site RNA helix on the basis of surface and electrostatic features

We have investigated electrostatic and surface features of an essential region of the catalytic core of the spliceosome, the eukaryotic precursor messenger (pre-m)RNA splicing apparatus. The nucleophile for the first of two splicing reactions is the 2′-hydroxyl (OH) of the ribose of a specific adenosine within the intron. During assembly of the spliceosome's catalytic core, this adenosine is positioned by pairing with a short region of the U2 small nuclear (sn)RNA to form the pre-mRNA branch site helix. The solution structure of the spliceosomal pre-mRNA branch site [Newby,M.I. and Greenbaum,N.L. (2002) Nature Struct. Biol., 9, 958–965] showed that a phylogenetically conserved pseudouridine (ψ) residue in the segment of U2 snRNA that pairs with the intron induces a markedly different structure compared with that of its unmodified counterpart. In order to achieve a more detailed understanding of the factors that contribute to recognition of the spliceosome's branch site helix and activation of the nucleophile for the first step of pre-mRNA splicing, we have calculated surface areas and electrostatic potentials of ψ-modified and unmodified branch site duplexes. There was no significant difference between the total accessible area or ratio of total polar:nonpolar groups between modified and unmodified duplexes. However, there was substantially greater exposure of nonpolar area of the adenine base, and less exposure of the 2′-OH, in the ψ-modified structure. Electrostatic potentials computed using a hybrid boundary element and finite difference nonlinear Poisson–Boltzmann approach [Boschitsch, A.H. and Fenley, M.O. (2004) J. Comput. Chem., 25, 935–955] revealed a region of exceptionally negative potential in the major groove surrounding the 2′-OH of the branch site adenosine. These surface and electrostatic features may contribute to the overall recognition of the pre-mRNA branch site region by other components of the splicing reaction.

Molecular dynamics simulation of the M2 helices within the nicotinic acetylcholine receptor transmembrane domain: structure and collective motions

Multiple nanosecond duration molecular dynamics simulations were performed on the transmembrane region of the Torpedo nicotinic acetylcholine receptor embedded within a bilayer mimetic octane slab. The M2 helices and M2-M3 loop regions were free to move, whereas the outer (M1, M3, M4) helix bundle was backbone restrained. The M2 helices largely retain their hydrogen-bonding pattern throughout the simulation, with some distortions in the helical end and loop regions. All of the M2 helices exhibit bending motions, with the hinge point in the vicinity of the central hydrophobic gate region (corresponding to residues alphaL251 and alphaV255). The bending motions of the M2 helices lead to a degree of dynamic narrowing of the pore in the region of the proposed hydrophobic gate. Calculations of Born energy profiles for various structures along the simulation trajectory suggest that the conformations of the M2 bundle sampled correspond to a closed conformation of the channel. Principal components analyses of each of the M2 helices, and of the five-helix M2 bundle, reveal concerted motions that may be relevant to channel function. Normal mode analyses using the anisotropic network model reveal collective motions similar to those identified by principal components analyses.

Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets

We predict regulatory targets of vertebrate microRNAs (miRNAs) by identifying mRNAs with conserved complementarity to the seed (nucleotides 2-7) of the miRNA. An overrepresentation of conserved adenosines flanking the seed complementary sites in mRNAs indicates that primary sequence determinants can supplement base pairing to specify miRNA target recognition. In a four-genome analysis of 3' UTRs, approximately 13,000 regulatory relationships were detected above the estimate of false-positive predictions, thereby implicating as miRNA targets more than 5300 human genes, which represented 30% of our gene set. Targeting was also detected in open reading frames. In sum, well over one third of human genes appear to be conserved miRNA targets.

Crystal Structure of AhpE from Mycobacterium tuberculosis, a 1-Cys peroxiredoxin

All living systems require protection against the damaging effects of reactive oxygen species. The genome of Mycobacterium tuberculosis, the cause of TB, encodes a number of peroxidases that are thought to be active against organic and inorganic peroxides, and are likely to play a key role in the ability of this organism to survive within the phagosomes of macrophages. The open reading frame Rv2238c in M.tuberculosis encodes a 153-residue protein AhpE, which is a peroxidase of the 1-Cys peroxiredoxin (Prx) family. The crystal structure of AhpE, determined at 1.87 A resolution (R(cryst)=0.179, R(free)=0.210), reveals a compact single-domain protein with a thioredoxin fold. AhpE forms both dimers and octamers; a tightly-associated dimer and a ring-like octamer, generated by crystallographic 4-fold symmetry. In this native structure, the active site Cys45 is in its oxidized, sulfenic acid (S-O-H) state. A second crystal form of AhpE, obtained after soaking in sodium bromide and refined at 1.90 A resolution (R(cryst)=0.242, R(free)=0.286), reveals the reduced structure. In this structure, a conformational change in an external loop, in two of the four molecules in the asymmetric unit, allows Arg116 to stabilise the Cys45 thiolate ion, and concomitantly closes a surface channel. This channel is identified as the likely binding site for a physiological reductant, and the conformational change is inferred to be important for the reaction cycle of AhpE.

Electrostatic correlations and fluctuations for ion binding to a finite length polyelectrolyte

A statistical mechanical model is presented which explicitly accounts for the fluctuations, the electrostatic, and the excluded volume correlations for ions bound to a polyelectrolyte such as DNA. The method can be employed to treat a wide range of ionic conditions including multivalent ions. The microscopic framework of the theory permits the use of realistic finite length and grooved structural model for the polyelectrolyte and modeling of the finite size of the bound ions. Test against Monte Carlo simulations suggests that the theory can give accurate predictions for the ion distribution and the thermodynamic properties. For multivalent ions, the theory makes improved predictions as compared with the mean-field approach. Moreover, for long polyelectrolyte and dilute salt concentration, the theory predicts ion binding properties that agree with the counterion condensation theory.

Solving the Poisson-Boltzmann equation with the specialized computer chip MD-GRAPE-2

The accurate description of solvation effects is highly desirable in numerous computational chemistry applications. One widely used methodology treats the solvent as a uniform continuum ("implicit solvation"), and describes its net interaction with the solute by solving the Poisson-Boltzmann (PB) equation using the Boundary Element Method (BEM). These calculations are very time consuming using conventional computers. A new, efficient way of solving the PB equation on the special-purpose hardware chip MDGRAPE-2 is presented. The MDGRAPE-2 chip was originally developed to speed up the calculation of pairwise van der Waals and electrostatic interactions. By recasting the BEM equations and applying an iterative solution procedure, a hardware-accelerated PB solver on MDGRAPE-2 could be implemented. The performance and reliability of the method is demonstrated on three examples ranging from small peptides to large proteins, whereby the obtained acceleration factors range from 15-fold to up to 40-fold with no loss of accuracy compared to the conventional approach.

Structural basis for phosphatidylinositol phosphate kinase type Igamma binding to talin at focal adhesions

The cytoskeletal protein talin binds to a short C-terminal sequence in phosphatidylinositol phosphate kinase type Igamma (PIPKIgamma), activating the enzyme and promoting the local production of phosphatidylinositol 4,5 bisphosphate, which regulates focal adhesion dynamics as well as clathrin-mediated endocytosis in neuronal cells. Here we show by crystallographic, NMR, and calorimetric analysis that the phosphotyrosine binding (PTB)-like domain of talin engages the PIPKIgamma C terminus in a mode very similar to that of integrin binding. However, PIPKIgamma binds in the canonical PTB-peptide mode with an SPLH motif replacing the classic NPXY motif. The tighter packing of the SPLH motif against the hydrophobic core of talin may explain the stronger binding of PIPKIgamma. Two tyrosine residues flanking the SPLH motif (Tyr-644 and Tyr-649) have been implicated in the regulation of talin binding. We show that phosphorylation at Tyr-644, a Src phosphorylation site in vivo, has little effect on the binding mode or strength, which is consistent with modeling studies in which the phosphotyrosine makes surface-exposed salt bridges, and we suggest that its strong activating effect arises from the release of autoinhibitory restraints in the full-length PIPKIgamma. Modeling studies suggest that phosphorylation of Tyr-649 will likewise have little effect on talin binding, whereas phosphorylation of the SPLH serine is predicted to be strongly disruptive. Our data are consistent with the proposal that Src activity promotes a switch from integrin binding to PIPKIgamma binding that regulates focal adhesion turnover.

Structure and dynamics of micelle-bound human alpha-synuclein

Misfolding of the protein alpha-synuclein (aS), which associates with presynaptic vesicles, has been implicated in the molecular chain of events leading to Parkinson's disease. Here, the structure and dynamics of micelle-bound aS are reported. Val3-Val37 and Lys45-Thr92 form curved alpha-helices, connected by a well ordered, extended linker in an unexpected anti-parallel arrangement, followed by another short extended region (Gly93-Lys97), overlapping the recently identified chaperone-mediated autophagy recognition motif and a highly mobile tail (Asp98-Ala140). Helix curvature is significantly less than predicted based on the native micelle shape, indicating a deformation of the micelle by aS. Structural and dynamic parameters show a reduced helical content for Ala30-Val37. A dynamic variation in interhelical distance on the microsecond timescale is complemented by enhanced sub-nanosecond timescale dynamics, particularly in the remarkably glycine-rich segments of the helices. These unusually rich dynamics may serve to mitigate the effect of aS binding on membrane fluidity. The well ordered conformation of the helix-helix connector indicates a defined interaction with lipidic surfaces, suggesting that, when bound to larger diameter synaptic vesicles, it can act as a switch between this structure and a previously proposed uninterrupted helix.

A quantitative assessment of models for voltage-dependent gating of ion channels

Voltage-gated ion channels open and close, or "gate," in response to changes in membrane potential. The electric field across the membrane-protein complex exerts forces on charged residues driving the channel into different functional conformations as the membrane potential changes. To act with the greatest sensitivity, charged residues must be positioned at key locations within or near the transmembrane region, which requires desolvating charged groups, a process that can be energetically prohibitive. Although there is good agreement on which residues are involved in this process for voltage-activated potassium channels, several different models of the sensor geometry and gating motions have been proposed. Here we incorporate low-resolution structural information about the channel into a Poisson-Boltzmann calculation to determine solvation barrier energies and gating charge values associated with each model. The principal voltage-sensing helix, S4, is represented explicitly, whereas all other regions are represented as featureless, dielectric media with complex boundaries. From our calculations, we conclude that a pure rotation of the S4 segment within the voltage sensor is incapable of producing the observed gating charge values, although this shortcoming can be partially remedied by first tipping and then minimally translating the S4 helix. Models in which the S4 segment has substantial interaction with the low-dielectric environment of the membrane incur solvation energies of hundreds of k(B)T, and activation times based on these energies are orders of magnitude slower than experimentally observed.

Functional insights from the structure of the multifunctional C345C domain of C5 of complement

The complement protein C5 initiates assembly of the membrane attack complex. This remarkable process results in lysis of target cells and is fundamental to mammalian defense against infection. The 150-amino acid residue domain at the C terminus of C5 (C5-C345C) is pivotal to C5 function. It interacts with enzymes that convert C5 to C5b, the first step in the assembly of the membrane attack complex; it also binds to the membrane attack complex components C6 and C7 with high affinity. Here a recombinant version of this C5-C345C domain is shown to adopt the oligosaccharide/oligonucleotide binding fold, with two helices packed against a five-stranded beta-barrel. The structure is compared with those from the netrin-like module family that have a similar fold. Residues critical to the interaction with C5-convertase cluster on a mobile, hydrophobic inter-strand loop that protrudes from the open face of the beta-barrel. The opposite, helix-dominated face of C5-C345C carries a pair of exposed hydrophobic side chains adjacent to a striking negatively charged patch, consistent with affinity for positively charged factor I modules in C6 and C7. Modeling of homologous domains from complement proteins C3 and C4, which do not participate in membrane attack complex assembly, suggests that this provisionally identified C6/C7-interacting face is indeed specific to C5.

Diffusion of nucleoside triphosphates and role of the entry site to the RNA polymerase II active center

Nucleoside triphosphates (NTPs) diffuse to the active center of RNA polymerase II through a funnel-shaped opening that narrows to a negatively charged pore. Computer simulation shows that the funnel and pore reduce the rate of diffusion by a factor of approximately 2 x 10(-7). The resulting limitation on the rate of RNA synthesis under conditions of low NTP concentration may be overcome by NTP binding to an entry site adjacent to the active center. Binding to the entry site greatly enhances the lifetime of an NTP in the active center region, and it prevents "backtracking" and the consequent occlusion of the active site.

Electrostatic interaction between RNA and protein capsid in cowpea chlorotic mottle virus simulated by a coarse-grain RNA model and a Monte Carlo approach

Although many viruses have been crystallized and the protein capsid structures have been determined by x-ray crystallography, the nucleic acids often cannot be resolved. This is especially true for RNA viruses. The lack of information about the conformation of DNA/RNA greatly hinders our understanding of the assembly mechanism of various viruses. Here we combine a coarse-grain model and a Monte Carlo method to simulate the distribution of viral RNA inside the capsid of cowpea chlorotic mottle virus. Our results show that there is very strong interaction between the N-terminal residues of the capsid proteins, which are highly positive charged, and the viral RNA. Without these residues, the binding energy disfavors the binding of RNA by the capsid. The RNA forms a shell close to the capsid with the highest densities associated with the capsid dimers. These high-density regions are connected to each other in the shape of a continuous net of triangles. The overall icosahedral shape of the net overlaps with the capsid subunit icosahedral organization. Medium density of RNA is found under the pentamers of the capsid. These findings are consistent with experimental observations.

The parameterization and validation of generalized born models using the pairwise descreening approximation

Generalized Born Surface Area (GBSA) models for water using the Pairwise Descreening Approximation (PDA) have been parameterized by two different methods. The first method, similar to that used in previously reported parameterizations, optimizes all parameters against the experimental free energies of hydration of organic molecules. The second method optimizes the PDA parameters to compensate only for systematic errors of the PDA. The best models are compared to Poisson-Boltzmann calculations and applied to the computation of potentials of mean force (PMFs) for the association of various molecules. PMFs present a more rigorous test of the ability of a solvation model to correctly reproduce the screening of intermolecular interactions by the solvent, than its accuracy at predicting free energies of hydration of small molecules. Models derived with the first method are sometimes shown to fail to compute accurate potentials of mean force because of large errors in the computation of Born radii, while no such difficulties are observed with the second method. Furthermore, accurate computation of the Born radii appears to be more important than good agreement with experimental free energies of solvation. We discuss the source of errors in the potentials of mean force and suggest means to reduce them. Our findings suggest that Generalized Born models that use the Pairwise Descreening Approximation and that are derived solely by unconstrained optimization of parameters against free energies of hydration should be applied to the modeling of intermolecular interactions with caution.

Not Ions Alone: Barriers to Ion Permeation in Nanopores and Channels

A hydrophobic pore of subnanometer dimensions can appear impermeable to an ion even though its radius is still much wider than that of the ion. Pores of molecular dimensions can be found, for instance, in carbon nanotubes, zeolites, or ion channel proteins. We quantify this barrier to ion permeation by calculating the potential of mean force from umbrella-sampled molecular dynamics simulations and compare them to continuum-electrostatic Poisson-Boltzmann calculations. The latter fail to describe the ion barrier because they do not account for the properties of water in the pore. The barrier originates from the energetic cost to desolvate the ion. Even in wide pores, which could accommodate an ion and its hydration shell, a barrier of several kT remains because the liquid water phase is not stable in the hydrophobic pore. Thus, the properties of the solvent play a crucial role in determining permeation properties of ions in confinement at the molecular scale.

Structure of a Human A-type Potassium Channel Interacting Protein DPPX, a Member of the Dipeptidyl Aminopeptidase Family

It has recently been reported that dipeptidyl aminopeptidase X (DPPX) interacts with the voltage-gated potassium channel Kv4 and that co-expression of DPPX together with Kv4 pore forming alpha-subunits, and potassium channel interacting proteins (KChIPs), reconstitutes properties of native A-type potassium channels in vitro. Here we report the X-ray crystal structure of the extracellular domain of human DPPX determined at 3.0A resolution. This structure reveals the potential for a surface electrostatic change based on the protonation state of histidine. Subtle changes in extracellular pH might modulate the interaction of DPPX with Kv4.2 and possibly with other proteins. We propose models of DPPX interaction with the voltage-gated potassium channel complex. The dimeric structure of DPPX is highly homologous to the related protein DPP-IV. Comparison of the active sites of DPPX and DPP-IV reveals loss of the catalytic serine residue but the presence of an additional serine near the "active" site. However, the arrangement of residues is inconsistent with that of canonical serine proteases and DPPX is unlikely to function as a protease (dipeptidyl aminopeptidase).

Ribosome motions modulate electrostatic properties

The electrostatic properties of the 70S ribosome of Thermus thermophilus were studied qualitatively by solving the Poisson-Boltzmann (PB) equation in aqueous solution and with physiological ionic strength. The electrostatic potential was calculated for conformations of the ribosome derived by recent normal mode analysis (Tama, F., et al. Proc Natl Acad Sci USA 2003 100, 9319-9323) of the ratchet-like reorganization that occurs during translocation (Frank, J.; Agrawal, R. K. Nature 2000 406, 318-322). To solve the PB equation, effective parameters (charges and radii), applicable to a highly charged backbone model of the ribosome, were developed. Regions of positive potential were found at the binding site of the elongation factors G and Tu, as well as where the release factors bind. Large positive potential areas are especially pronounced around the L11 and L6 proteins. The region around the L1 protein is also positively charged, supporting the idea that L1 may interact with the E-site tRNA during its release from the ribosome after translocation. Functional rearrangement of the ribosome leads to electrostatic changes which may help the translocation of the tRNAs during the elongation stage.