Superfast Mass Transport of Na/K Via Mesochannels for Dendrite-Free Metal Batteries

Fast ion diffusion in anode hosts enabling uniform distribution of Li/Na/K is essential for achieving dendrite-free alkali-metal batteries. Common strategies, e.g. expanding the interlayer spacing of anode materials, can enhance bulk diffusion of Li but are less efficient for Na and K due to their larger ionic radius. Herein, a universal strategy to drastically improve the mass-transport efficiency of Na/K by introducing open mesochannels in carbon hosts is proposed. Such pore engineering can increase the accessible surface area by one order of magnitude, thus remarkably accelerating surface diffusion, as visualized by in situ transmission electron microscopy. In particular, once the mesochannels are filled by the Na/K metals, they become the superfast channels for mass transport via the mechanism of interfacial diffusion. Thus-modified carbon hosts enable Na/K filling in their inner cavities and uniform deposition across the whole electrodes with fast kinetics. The resulting Na-metal anodes can exhibit stable dendrite-free cycling with outstanding rate performance at a high current density of up to 30 mA cm^-2 . This work presents an inspiring attempt to address the sluggish transport issue of Na/K, as well as valuable insights into the mass-transport mechanism in porous anodes for high-performance alkali-metal storage.

New Model for Stacking Monomers in Filamentous Actin from Skeletal Muscles of Oryctolagus cuniculus

To date, some scientific evidence (limited proteolysis, mass spectrometry analysis, electron microscopy (EM)) has accumulated, which indicates that the generally accepted model of double-stranded of filamentous actin (F-actin) organization in eukaryotic cells is not the only one. This entails an ambiguous understanding of many of the key cellular processes in which F-actin is involved. For a detailed understanding of the mechanism of F-actin assembly and actin interaction with its partners, it is necessary to take into account the polymorphism of the structural organization of F-actin at the molecular level. Using electron microscopy, limited proteolysis, mass spectrometry, X-ray diffraction, and structural modeling we demonstrated that F-actin presented in the EM images has no double-stranded organization, the regions of protease resistance are accessible for action of proteases in F-actin models. Based on all data, a new spatial model of filamentous actin is proposed, and the F-actin polymorphism is discussed.

Comparative studies on ion-pair energetic, distribution among three domains of life: Archaea, eubacteria, and eukarya

Salt-bridges play a unique role in the structural and functional stability of proteins, especially under harsh environments. How these salt-bridges contribute to the overall thermodynamic stability of protein structure and function across different domains of life is elusive still date. To address the issue, statistical analyses on the energies of salt-bridges, involved in proteins' structure and function, are performed across three domains of life, that is, archaea, eubacteria, and eukarya. Results show that although the majority of salt-bridges are stable and conserved, yet the stability of archaeal proteins (∆∆G_net = -5.06 ± 3.8) is much more than that of eubacteria (∆∆G_net = -3.7 ± 2.9) and eukarya (∆∆G_net = -3.54 ± 3.1). Unlike earlier study with archaea, in eukarya and eubacteria, not all buried salt-bridge in our dataset are stable. Buried salt-bridges play surprising role in protein stability, whose variations are clearly observed among these domains. Greater desolvation penalty of buried salt-bridges is compensated by stable network of salt-bridges apart from equal contribution of bridge and background energy terms. On the basis proteins' secondary structure, topology, and evolution, our observation shows that salt-bridges when present closer to each other in sequence tend to form a greater number. Overall, our comparative study provides insight into the role of specific electrostatic interactions in proteins from different domains of life, which we hope, would be useful for protein engineering and bioinformatics study.

Structural features affecting the enzymatic digestibility of pine wood pretreated with ionic liquids

Pretreating lignocellulosic biomass with certain ionic liquids results in structural and chemical changes that make the biomass more digestible by enzymes. In this study, pine wood was pretreated with 1-ethyl-3-methylimidazolium chloride/acetate ([C2 mim]Cl and [C2 mim][OAc]) at different temperatures to investigate the relative importance of substrate features, such as accessible surface area, cellulose crystallinity, and lignin content, on enzymatic digestibility. The ionic liquid pretreatments resulted in glucan conversions ranging from 23% to 84% on saccharification of the substrates, with [C2 mim][OAc] being more effective than [C2 mim]Cl. The pretreatments resulted in no delignification of the wood, some loss of cellulose crystallinity under certain conditions, and varying levels of increased surface area. Enzymatic digestibility closely correlated with accessible surface area and porosity measurements obtained using Simons' staining and thermoporosimetry techniques. Increased accessible surface area was identified as the principal structural feature responsible for the improved enzymatic digestibility.

Accurate single-sequence prediction of solvent accessible surface area using local and global features

We present a new approach for predicting the Accessible Surface Area (ASA) using a General Neural Network (GENN). The novelty of the new approach lies in not using residue mutation profiles generated by multiple sequence alignments as descriptive inputs. Instead we use solely sequential window information and global features such as single-residue and two-residue compositions of the chain. The resulting predictor is both highly more efficient than sequence alignment-based predictors and of comparable accuracy to them. Introduction of the global inputs significantly helps achieve this comparable accuracy. The predictor, termed ASAquick, is tested on predicting the ASA of globular proteins and found to perform similarly well for so-called easy and hard cases indicating generalizability and possible usability for de-novo protein structure prediction. The source code and a Linux executables for GENN and ASAquick are available from Research and Information Systems at http://mamiris.com, from the SPARKS Lab at http://sparks-lab.org, and from the Battelle Center for Mathematical Medicine at http://mathmed.org.

Sequence-specific transitions of the torsion angle gamma change the polar-hydrophobic profile of the DNA grooves: implication for indirect protein-DNA recognition

Variations of the shape and polarity of the DNA grooves caused by changes of the DNA conformation play an important role in the DNA readout. Despite the fact that non-canonical trans and gauche- conformations of the DNA backbone angle γ (O5'-C5'-C4'-C3') are frequently found in the DNA crystal structures, their possible role in the DNA recognition has not been studied systematically. In order to fill in this gap, we analyze the available high-resolution crystal structures of the naked and complexed DNA. The analysis shows that the non-canonical γ angle conformations are present both in the naked and bound DNA, more often in the bound vs. naked DNA, and in the nucleotides with the A-like vs. the B-like sugar pucker. The alternative angle γ torsions are more frequently observed in the purines with the A-like sugar pucker and in the pyrimidines with the B-like sugar conformation. The minor groove of the nucleotides with non-canonical γ angle conformation is more polar, while the major groove is more hydrophobic than in the nucleotides with the classical γ torsions due to variations in exposure of the polar and hydrophobic groups of the DNA backbone. The propensity of the nucleotides with different γ angle conformations to participate in the protein-nucleic acid contacts in the minor and major grooves is connected with their sugar pucker and sequence-specific. Our findings imply that the angle γ transitions contribute to the process of the protein-DNA recognition due to modification of the polar/hydrophobic profile of the DNA grooves.

The C2B domain is the primary Ca2+ sensor in DOC2B: a structural and functional analysis

DOC2B (double-C2 domain) protein is thought to be a high-affinity Ca(2+) sensor for spontaneous and asynchronous neurotransmitter release. To elucidate the molecular features underlying its physiological role, we determined the crystal structures of its isolated C2A and C2B domains and examined their Ca(2+)-binding properties. We further characterized the solution structure of the tandem domains (C2AB) using small-angle X-ray scattering. In parallel, we tested structure-function correlates with live cell imaging tools. We found that, despite striking structural similarity, C2B binds Ca(2+) with considerably higher affinity than C2A. The C2AB solution structure is best modeled as two domains with a highly flexible orientation and no difference in the presence or absence of Ca(2+). In addition, kinetic studies of C2AB demonstrate that, in the presence of unilamellar vesicles, Ca(2+) binding is stabilized, as reflected by the ~10-fold slower rate of Ca(2+) dissociation than in the absence of vesicles. In cells, isolated C2B translocates to the plasma membrane (PM) with an EC50 of 400 nM while the C2A does not translocate at submicromolar Ca(2+) concentrations, supporting the biochemical observations. Nevertheless, C2AB translocates to the PM with an ~2-fold lower EC50 and to a greater extent than C2B. Our results, together with previous studies, reveal that the C2B is the primary Ca(2+) sensing unit in DOC2B, whereas C2A enhances the interaction of C2AB with the PM.

The Role of Aromatic-Aromatic Interactions in Strand-Strand Stabilization of β-Sheets

Aromatic-aromatic interactions have long been believed to play key roles in protein structure, folding, and binding functions. However, we still lack full understanding of the contributions of aromatic-aromatic interactions to protein stability and the timing of their formation during folding. Here, using an aromatic ladder in the β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), as a case study, we find that aromatic π stacking plays a greater role in the Phe65-Phe71 cross-strand pair, while in another pair, Phe50-Phe65, hydrophobic interactions are dominant. The Phe65-Phe71 pair spans β-strands 4 and 5 in the β-barrel, which lack interstrand hydrogen bonding, and we speculate that it compensates energetically for the absence of strand-strand backbone interactions. Using perturbation analysis, we find that both aromatic-aromatic pairs form after the transition state for folding of CRABP1, thus playing a role in the final stabilization of the β-sheet rather than in its nucleation as had been earlier proposed. The aromatic interaction between strands 4 and 5 in CRABP1 is highly conserved in the intracellular lipid-binding protein (iLBP) family, and several lines of evidence combine to support a model wherein it acts to maintain barrel structure while allowing the dynamic opening that is necessary for ligand entry. Lastly, we carried out a bioinformatics analysis and found 51 examples of aromatic-aromatic interactions across non-hydrogen-bonded β-strands outside the iLBPs, arguing for the generality of the role played by this structural motif.

Insights of interaction between small and large subunits of ADP-glucose pyrophosphorylase from bread wheat (Triticum aestivum L.)

Lack of knowledge of three dimensional structures of small and large subunits of ADP- glucose pyrophosphorylase (AGPase) in wheat has hindered efforts to understand the binding specifities of substrate and catalytic mechanism. Thus, to understand the structure activity relationship, 3D structures were built by homology modelling based on crystal structure of potato tuber ADP-glucose pyrophosphorylase. Selected models were refined by energy minimization and further validated by Procheck and Prosa-web analysis. Ramachandran plot showed that overall main chain and side chain parameters are favourable. Moreover, Z-score of the models from Prosa-web analysis gave the conformation that they are in the range of the template. Interaction analysis depicts the involvement of six amino acids in hydrogen bonding (AGP-SThr422-AGP-LMet138, AGP- SArg420-AGP-LGly47, AGP-SSer259-AGP-LSer306, AGP-SGlu241-AGP-LIle311, AGPSGln113- AGP-LGlu286 and AGP-SGln70-AGP-LLys291). Fifteen amino acids of small subunit were able to make hydrophobic contacts with seventeen amino acids of large subunit. Furthermore, decrease in the solvent accessible surface area in the amino acids involved in interaction were also reported. All the distances were formed in between 2.27 to 3.78Å. The present study focussed on heterodimeric structure of (AGPase). This predicted complex not only enhance our understanding of the interaction mechanism between these subunits (AGP-L and AGP-S) but also enable to further study to obtain better variants of this enzyme for the improvement of the plant yield.

Anisotropic solvent model of the lipid bilayer. 1. Parameterization of long-range electrostatics and first solvation shell effects

A new implicit solvation model was developed for calculating free energies of transfer of molecules from water to any solvent with defined bulk properties. The transfer energy was calculated as a sum of the first solvation shell energy and the long-range electrostatic contribution. The first term was proportional to solvent accessible surface area and solvation parameters (σ(i)) for different atom types. The electrostatic term was computed as a product of group dipole moments and dipolar solvation parameter (η) for neutral molecules or using a modified Born equation for ions. The regression coefficients in linear dependencies of solvation parameters σ(i) and η on dielectric constant, solvatochromic polarizability parameter π*, and hydrogen-bonding donor and acceptor capacities of solvents were optimized using 1269 experimental transfer energies from 19 organic solvents to water. The root-mean-square errors for neutral compounds and ions were 0.82 and 1.61 kcal/mol, respectively. Quantification of energy components demonstrates the dominant roles of hydrophobic effect for nonpolar atoms and of hydrogen-bonding for polar atoms. The estimated first solvation shell energy outweighs the long-range electrostatics for most compounds including ions. The simplicity and computational efficiency of the model allows its application for modeling of macromolecules in anisotropic environments, such as biological membranes.

HIV-1 envelope accessible surface and polarity: clade, blood, and brain

The human immunodeficiency virus type-1 (HIV-1) gp160 (gp120-gp41 complex) trimer envelope (ENV) protein is a potential vaccine candidate for HIV/AIDS. HIV-1 vaccine development has been problematic and charge polarity as well as sequence variation across clades may relate to the difficulties. Further obstacles are caused by sequence variation between blood and brain-derived sequences, since the brain is a separate compartment for HIV-1 infection. We utilize a threedimensional residue measure of solvent exposure, accessible surface area (ASA), which shows that major segments of gp120 and gp41 known structures are solvent exposed across clades. We demonstrate a large percent sequence polarity for solvent exposed residues in gp120 and gp41. The range of sequence polarity varies across clades, blood, and brain from different geographical locations. Regression analysis shows that blood and brain gp120 and gp41 percent sequence polarity range correlate with mean Shannon entropy. These results point to the use of protein modifications to enhance HIV-1 ENV vaccines across multiple clades, blood, and brain. It should be noted that we do not address the issue of protein glycosylation here; however, this is an important issue for vaccine design and development.

ABBREVIATIONS

HIV-1 - human immunodeficiency virus type 1, AIDS - acquired immunodeficiency syndrome, ENV - envelope, gp160 - 160,000d glycoprotein, gp120 - 120,000d glycoprotein, gp41 - 41,000d glycoprotein, LANL - Los Alamos National Laboratories, PDB - Protein Data Bank, HVTN - STEP HIV vaccine trial, AA - amino acids, MSA - multiple sequence alignment, ASA - accessible surface area, SNPs- single nucleotide polymorphisms, HAART - Highly Active Antiretroviral Therapy, CCR5 - C-C chemokine receptor type 5, CNS - central nervous system, HIVE - HIV encephalitis, P - polarity, NP - non-polarity, CTL - cytotoxic T lymphocyte, NIAID - National Institute of Allergy and Infectious Diseases.

Strand 6B deformation and residues exposure towards N-terminal end of helix B during proteinase inhibition by Serpins

Serine Protease inhibitors (Serpins) like antithrombin, antitrypsin, neuroserpin, antichymotrypsin, protein C-inhibitor and plasminogen activator inhibitor is involved in important biological functions like blood coagulation, fibrinolysis, inflammation, cell migration and complement activation. Serpins native state is metastable, which undergoes transformation to a more stable state during the process of protease inhibition. Serpins are prone to conformation defects, however little is known about the factors and mechanisms which promote its conformational change and misfolding. Helix B region in serpins is with several point mutations which result in pathological conditions due to polymerization. Helix B analysis for residue burial and cavity was undertaken to understand its role in serpin structure function. A structural overlap and an accessible surface area analysis showed the deformation of strand 6B and exposure of helix B at N-terminal end in cleaved conformation but not in the native and latent conformation of various inhibitory serpins. A cleaved polymer like conformation of antitrypsin also showed deformation of s6B and helix B exposure. Cavity analysis showed that helix B residues were part of the largest cavity in most of the serpins in the native state which increase in size during the transformation to cleaved and latent states. These data for the first time show the importance of strand 6B deformation and exposure of helix B in smooth insertion of the reactive center loop during serpin inhibition and indicate that helix B exposure due to variants may increase its polymer propensity.

Probing protein conformation with a minimal photochemical reagent

3H-diazirine (3H-DZN), a photoreactive gas similar in size to water, was used to probe the topography of the surface and inner space of proteins. On photolysis 3H-DZN generates 3H-methylene carbene, which reacts unselectively with its molecular cage, inserting even into C-H bonds. Labeling of bovine alpha-lactalbumin (alpha-LA, MW: 14,200) with 1 mM (3)H-DZN yielded 0.0041 mol CH2/mol of protein, in agreement with the expectation for an unspecific surface-labeling phenomenon. The cooperative urea-induced unfolding of alpha-LA, as monitored by the extent of 3H-methylene labeling, agrees with that measured by circular dichroism spectroscopy in the far and near ultraviolet regions. At 8 M urea, the unfolded state U was labeled 25-30% more than the native state N primarily because of the increase in the accessible surface area (ASA) of the protein occurring upon unfolding. However, this result lies below the approximately 100% increment expected from theoretical estimates of ASA of state U. Among other factors, most likely the existence of a residual structure in U, that involves helices H2 and H4 of the alpha subdomain, might account for this fact, as shown by a comparative analysis of peptide labeling patterns of N and U samples. In this paper, we demonstrate the usefulness of the 3H-methylene labeling method to monitor conformational transitions and map solvent accessibility along the polypeptide sequence, thus opening the possibility of outlining structural features of nonnative states (i.e., denatured states, molten globule). We anticipate that this technique also would help to identify ligand binding and oligomerization sites in proteins.

Conformational strictness required for maximum activity and stability of bovine pancreatic ribonuclease A as revealed by crystallographic study of three Phe120 mutants at 1.4 A resolution

The replacement of Phe120 with other hydrophobic residues causes a decrease in the activity and thermal stability in ribonuclease A (RNase A). To explain this, the crystal structures of wild-type RNase A and three mutants--F120A, F120G, and F120W--were analyzed up to a 1.4 A resolution. Although the overall backbone structures of all mutant samples were nearly the same as that of wild-type RNase A, except for the C-terminal region of F120G with a high B-factor, two local conformational changes were observed at His119 in the mutants. First, His119 of the wild-type and F120W RNase A adopted an A position, whereas those of F120A and F120G adopted a B position, but the static crystallographic position did not reflect either the efficiency of transphosphorylation or the hydrolysis reaction. Second, His119 imidazole rings of all mutant enzymes were deviated from that of wild-type RNase A, and those of F120W and F120G appeared to be "inside out" compared with that of wild-type RNase A. Only approximately 1 A change in the distance between N(epsilon2) of His12 and N(delta1) of His119 causes a drastic decrease in k(cat), indicating that the active site requires the strict positioning of the catalytic residues. A good correlation between the change in total accessible surface area of the pockets on the surface of the mutant enzymes and enthalpy change in their thermal denaturation also indicates that the effects caused by the replacements are not localized but extend to remote regions of the protein molecule.