Design of symmetric TIM barrel proteins from first principles

Background: Computational protein design is a rapidly maturing field within structural biology, with the goal of designing proteins with custom structures and functions. Such proteins could find widespread medical and industrial applications. Here, we have adapted algorithms from the Rosetta software suite to design much larger proteins, based on ideal geometric and topological criteria. Furthermore, we have developed techniques to incorporate symmetry into designed structures. For our first design attempt, we targeted the (alpha/beta)(8) TIM barrel scaffold. We gained novel insights into TIM barrel folding mechanisms from studying natural TIM barrel structures, and from analyzing previous TIM barrel design attempts. Methods: Computational protein design and analysis was performed using the Rosetta software suite and custom scripts. Genes encoding all designed proteins were synthesized and cloned on the pET20-b vector. Standard circular dichroism and gel chromatographic experiments were performed to determine protein biophysical characteristics. 1D NMR and 2D HSQC experiments were performed to determine protein structural characteristics. Results: Extensive protein design simulations coupled with ab initio modeling yielded several all-atom models of ideal, 4-fold symmetric TIM barrels. Four such models were experimentally characterized. The best designed structure (Symmetrin-1) contained a polar, histidine-rich pore, forming an extensive hydrogen bonding network. Symmetrin-1 was easily expressed and readily soluble. It showed circular dichroism spectra characteristic of well-folded alpha/beta proteins. Temperature melting experiments revealed cooperative and reversible unfolding, with a T-m of 44 degrees C and a Gibbs free energy of unfolding (Delta G degrees) of 8.0 kJ/mol. Urea denaturing experiments confirmed these observations, revealing a C-m of 1.6 M and a Delta G degrees of 8.3 kJ/mol. Symmetrin-1 adopted a monomeric conformation, with an apparent molecular weight of 32.12 kDa, and displayed well resolved 1D-NMR spectra. However, the HSQC spectrum revealed somewhat molten characteristics. Conclusions: Despite the detection of molten characteristics, the creation of a soluble, cooperatively folding protein represents an advancement over previous attempts at TIM barrel design. Strategies to further improve Symmetrin-1 are elaborated. Our techniques may be used to create other large, internally symmetric proteins.

Anisotropic magnetic couplings and structure-driven canted to collinear transitions in Sr2IrO4 by magnetically constrained noncollinear DFT

We study the canted magnetic state in Sr2IrO4 using fully relativistic density functional theory (DFT) including an on-site Hubbard U correction. A complete magnetic phase diagram with respect to the tetragonal distortion and the rotation of IrO6 octahedra is constructed, revealing the presence of two types of canted to collinear magnetic transitions: a spin-flop transition with increasing tetragonal distortion and a complete quenching of the basal weak ferromagnetic moment below a critical octahedral rotation. Moreover, we put forward a scheme to study the anisotropic magnetic couplings by mapping magnetically constrained noncollinear DFT onto a general spin Hamiltonian. This procedure allows for the simultaneous account and direct control of the lattice, spin, and orbital interactions within a fully ab initio scheme. We compute the isotropic, single site anisotropy and Dzyaloshinskii-Moriya (DM) coupling parameters, and clarify that the origin of the canted magnetic state in Sr2IrO4 arises from the structural distortions and the competition between isotropic exchange and DM interactions.

Na2.44Mn1.79(SO4)(3): a new member of the alluaudite family of insertion compounds for sodium ion batteries

Sodium-ion batteries have been extensively pursued as economic alternatives to lithium-ion batteries. Investigating the polyanion chemistry, alluaudite structured Na2Fe2II(SO4)(3) has been recently discovered as a 3.8 V positive electrode material (Barpanda et al., Nature Commun., 5: 4358, 2014). Registering the highest ever Fe-III/Fe-II redox potential (vs. Na/Na+) and formidable energy density, it has opened up a new polyanion family for sodium batteries. Exploring the alluaudite family, here we report isotypical Na2+2xMn2-xII(SO4)(3) (x = 0.22) as a novel high-voltage cathode material for the first time. Following low-temperature (ca. 350 degrees C) solid-state synthesis, the structure of this new alluaudite compound has been solved adopting a monoclinic framework (s.g. C2/c) showing antiferromagnetic ordering at 3.4 K. Synergising experimental and ab initio DFT investigation, Na2+2xMn2-xII(SO4)(3) has been found to be a potential high-voltage (ca. 4.4 V) cathode material for sodium batteries.

X-Ha center dot center dot center dot C hydrogen bonds in n-alkane-HX (X = F, OH) complexes are stronger than C-Ha <-X hydrogen bonds

Computational study of X-Ha <-C and C-Ha <-X hydrogen bonds in n-alkane-HX complexes (X =F,OH, alkane =propane, butane, pentane) has been carried out in this work. Ab initio and density functional theories were used for this study. For n-alkane-H2O complexes both Oa <-H-C and O-Ha <-C hydrogen bonded complex have been found, while for n-alkane-HF complexes, our attempt to optimize Fa <-H-C H-bond was not successful. Like most of the hydrogen bonded systems, strong correlation between binding energy and stretching frequency of H-F and O-H stretching mode was observed. The values of electron density and Laplacian of electron density are within the accepted range for hydrogen bonds. In all these cases, X-Ha <-C hydrogen bonds are found to be stronger than C-Ha <-X hydrogen bonds.

Reconstruction of the Disassembly Pathway of an Icosahedral Viral Capsid and Shape Determination of Two Successive Intermediates

Viral capsids derived from an icosahedral plant virus widely used in physical and nanotechnological investigations were fully dissociated into dimers by a rapid change of pH. The process was probed in vitro at high spatiotemporal resolution by time-resolved small-angle X-ray scattering using a high brilliance synchrotron source. A powerful custom-made global fitting algorithm allowed us to reconstruct the most likely pathway parametrized by a set of stoichiometric coefficients and to determine the shape of two successive intermediates by ab initio calculations. None of these two unexpected intermediates was previously identified in self-assembly experiments, which suggests that the disassembly pathway is not a mirror image of the assembly pathway. These findings shed new light on the mechanisms and the reversibility of the assembly/disassembly of natural and synthetic virus-based systems. They also demonstrate that both the structure and dynamics of an increasing number of intermediate species become accessible to experiments.

Common effect of chemical and external pressures on the magnetic properties of RCoPO (R = La, Pr, Nd, Sm). II.

We investigate the direct correspondence between Co band ferromagnetism and structural parameters in the pnictide oxides RCoPO for different rare-earth ions (R = La, Pr, Nd, Sm) by means of muon-spin spectroscopy and ab initio calculations, complementing our results published previously G. Prando et al., Common effect of chemical and external pressures on the magnetic properties of RCoPO (R = La, Pr), Phys. Rev. B 87, 064401 (2013)]. We find that both the transition temperature to the ferromagnetic phase T-C and the volume of the crystallographic unit cell V are conveniently tuned by the R ionic radius and/or external pressure. We report a linear correlation between T-C and V and our ab initio calculations unambiguously demonstrate a full equivalence of chemical and external pressures. As such, we show that R ions influence the ferromagnetic phase only via the induced structural shrinkage without involving any active role from the electronic f degrees of freedom, which are only giving a sizable magnetic contribution at much lower temperatures.

Rational molecular dynamics scheme for predicting optimum concentration loading of nano-additive in phase change materials

The present study deals with the diffusion and phase transition behaviour of paraffin reinforced with carbon nano-additives namely graphene oxide (GO) and surface functionalized single walled carbon nanotubes (SWCNT). Bulk disordered systems of paraffin hydrocarbons impregnated with carbon nano-additives have been generated in realistic equilibrium conformations for potential application as latent heat storage systems. Ab initio molecular dynamics(MD) in conjugation with COMPASS forcefield has been implemented using periodic boundary conditions. The proposed scheme allows determination of optimum nano-additive loading for improving thermo-physical properties through analysis of mass, thermal and transport properties; and assists in determination of composite behaviour and related performance from microscopic point of view. It was observed that nanocomposites containing 7.8% surface functionalised SWCNT and 55% GO loading corresponds to best latent heat storage system. The propounded methodology could serve as a by-pass route for economically taxing and iterative experimental procedures required to attain the optimum composition for best performance. The results also hint at the large unexplored potential of ab-initio classical MD techniques for predicting performance of new nanocomposites for potential phase change material applications. (C) 2015 Author(s).

Synthesis, and crystal and electronic structure of sodium metal phosphate for use as a hybrid capacitor in non-aqueous electrolyte

Energy storage devices based on sodium have been considered as an alternative to traditional lithium based systems because of the natural abundance, cost effectiveness and low environmental impact of sodium. Their synthesis, and crystal and electronic properties have been discussed, because of the importance of electronic conductivity in supercapacitors for high rate applications. The density of states of a mixed sodium transition metal phosphate (maricite, NaMn1/3Co1/3Ni1/3PO4) has been determined with the ab initio generalized gradient approximation (GGA)+Hubbard term (U) method. The computed results for the mixed maricite are compared with the band gap of the parent NaFePO4 and the electrochemical experimental results are in good agreement. A mixed sodium transition metal phosphate served as an active electrode material for a hybrid supercapacitor. The hybrid device (maricite versus carbon) in a nonaqueous electrolyte shows redox peaks in the cyclic voltammograms and asymmetric profiles in the charge-discharge curves while exhibiting a specific capacitance of 40 F g(-1) and these processes are found to be quasi-reversible. After long term cycling, the device exhibits excellent capacity retention (95%) and coulombic efficiency (92%). The presence of carbon and the nanocomposite morphology, identified through X-ray photoelectron spectroscopy (XPS) and transmission electron microscopy (TEM) studies, ensures the high rate capability while offering possibilities to develop new cathode materials for sodium hybrid devices.

Continuum in the X-Z-Y Weak Bonds: Z= Main Group Elements

The Continuum in the variation of the X-Z bond length change from blue-shifting to red-shifting through zero-shifting in the X-Z---Y complex is inevitable. This has been analyzed by ab-initio molecular orbital calculations using Z= Hydrogen, Halogens, Chalcogens, and Pnicogens as prototypical examples. Our analysis revealed that, the competition between negative hyperconjugation within the donor (X-Z) molecule and Charge Transfer (CT) from the acceptor (Y) molecule is the primary reason for the X-Z bond length change. Here, we report that, the proper tuning of X-and Y-group for a particular Z-can change the blue-shifting nature of X-Z bond to zero-shifting and further to red-shifting. This observation led to the proposal of a continuum in the variation of the X-Z bond length during the formation of X-Z---Y complex. The varying number of orbitals and electrons available around the Z-atom differentiates various classes of weak interactions and leads to interactions dramatically different from the H-Bond. Our explanations based on the model of anti-bonding orbitals can be transferred from one class of weak interactions to another. We further take the idea of continuum to the nature of chemical bonding in general. (C) 2015 Wiley Periodicals, Inc.

Pyrolysis of 3-carene: Experiment, Theory and Modeling

Thermal decomposition studies of 3-carene, a bio-fuel, have been carried out behind the reflected shock wave in a single pulse shock tube for temperature ranging from 920 K to 1220 K. The observed products in thermal decomposition of 3-carene are acetylene, allene, butadiene, isoprene, cyclopentadiene, hexatriene, benzene, toluene and p-xylene. The overall rate constant for 3-carene decomposition was found to be k/s(-1) = 10((9.95 +/- 0.54)) exp(-40.88 +/- 2.71 kcal mol(-1) /RT). Ab-initio theoretical calculations were carried out to find the minimum energy pathway that could explain the formation of the observed products in the thermal decomposition experiments. These calculations were carried out at B3LYP/6-311 + G(d,p) and G3 level of theories. A kinetic mechanism explaining the observed products in the thermal decomposition experiments has been derived. It is concluded that the linear hydrocarbons are the primary products in the pyrolysis of 3-carene.

Weak-coupling superconductivity in a strongly correlated iron pnictide

Iron-based superconductors have been found to exhibit an intimate interplay of orbital, spin, and lattice degrees of freedom, dramatically affecting their low-energy electronic properties, including superconductivity. Albeit the precise pairing mechanism remains unidentified, several candidate interactions have been suggested to mediate the superconducting pairing, both in the orbital and in the spin channel. Here, we employ optical spectroscopy (OS), angle-resolved photoemission spectroscopy (ARPES), ab initio band-structure, and Eliashberg calculations to show that nearly optimally doped NaFe0.978Co0.022As exhibits some of the strongest orbitally selective electronic correlations in the family of iron pnictides. Unexpectedly, we find that the mass enhancement of itinerant charge carriers in the strongly correlated band is dramatically reduced near the Gamma point and attribute this effect to orbital mixing induced by pronounced spin-orbit coupling. Embracing the true band structure allows us to describe all low-energy electronic properties obtained in our experiments with remarkable consistency and demonstrate that superconductivity in this material is rather weak and mediated by spin fluctuations.

Origin of the Spin-Orbital Liquid State in a Nearly J=0 Iridate Ba3ZnIr2O9

We show using detailed magnetic and thermodynamic studies and theoretical calculations that the ground state of Ba3ZnIr2O9 is a realization of a novel spin-orbital liquid state. Our results reveal that Ba3ZnIr2O9 with Ir5+ (5d(4)) ions and strong spin-orbit coupling (SOC) arrives very close to the elusive J = 0 state but each Ir ion still possesses a weak moment. Ab initio density functional calculations indicate that this moment is developed due to superexchange, mediated by a strong intradimer hopping mechanism. While the Ir spins within the structural Ir2O9 dimer are expected to form a spin-orbit singlet state (SOS) with no resultant moment, substantial frustration arising from interdimer exchange interactions induce quantum fluctuations in these possible SOS states favoring a spin-orbital liquid phase down to at least 100 mK.

Molecular dynamics and simulations study on the vibrational and electronic solvatochromism of benzophenone

Solvent plays a key role in diverse physico-chemical and biological processes. Therefore, understanding solute-solvent interactions at the molecular level of detail is of utmost importance. A comprehensive solvatochromic analysis of benzophenone (Bzp) was carried out in various solvents using Raman and electronic spectroscopy, in conjunction with Density Functional Theory (DFT) calculations of supramolecular solute-solvent clusters generated using classical Molecular Dynamics Simulations (c-MDSs). The >C=O stretching frequency undergoes a bathochromic shift with solvent polarity. Interestingly, in protic solvents this peak appears as a doublet: c-MDS and ad hoc explicit solvent ab initio calculations suggest that the lower and higher frequency peaks are associated with the hydrogen bonded and dangling carbonyl group of Bzp, respectively. Additionally, the dangling carbonyl in methanol (MeOH) solvent is 4 cm(-1) blue-shifted relative to acetonitrile solvent, despite their similar dipolarity/polarizability. This suggests that the cybotactic region of the dangling carbonyl group in MeOH is very different from its bulk solvent structure. Therefore, we propose that this blue-shift of the dangling carbonyl originates in the hydrophobic solvation shell around it resulting from extended hydrogen bonding network of the protic solvents. Furthermore, the 1(1)n pi* (band I) and 1(1)pi pi* (band II) electronic transitions show a hypsochromic and bathochromic shift, respectively. In particular, these shifts in protic solvents are due to differences in their excited state-hydrogen bonding mechanisms. Additionally, a linear relationship is obtained for band I and the >C=O stretching frequency (cm(-1)), which suggests that the different excitation wavelengths in band I correspond to different solvation states. Therefore, we hypothesize that the variation in excitation wavelengths in band I could arise from different solvation states leading to varying solvation dynamics. This will have implications for ultrafast processes associated with electron-transfer, charge transfer, and also the photophysical aspects of excited states. (C) 2016 AIP Publishing LLC.

Ionothermal Synthesis of High-Voltage Alluaudite Na2+2xFe2-x(SO4)(3) Sodium Insertion Compound: Structural, Electronic, and Magnetic Insights

Exploring future cathode materials for sodium-ion batteries, alluaudite class of Na2Fe2II(SO4)(3) has been recently unveiled as a 3.8 V positive insertion candidate (Barpanda et al. Nat. Commun. 2014, 5, 4358). It forms an Fe-based polyanionic compound delivering the highest Fe-redox potential along with excellent rate kinetics and reversibility. However, like all known SO4-based insertion materials, its synthesis is cumbersome that warrants careful processing avoiding any aqueous exposure. Here, an alternate low temperature ionothermal synthesis has been described to produce the alluaudite Na2+2xFe2-xII(SO4)(3). It marks the first demonstration of solvothermal synthesis of alluaudite Na2+2xM2-xII(SO4)(3) (M = 3d metals) family of cathodes. Unlike classical solid-state route, this solvothermal route favors sustainable synthesis of homogeneous nanostructured alluaudite products at only 300 degrees C, the lowest temperature value until date. The current work reports the synthetic aspects of pristine and modified ionothermal synthesis of Na2+2xFe2-xII(SO4)(3) having tunable size (300 nm similar to 5 mu m) and morphology. It shows antiferromagnetic ordering below 12 K. A reversible capacity in excess of 80 mAh/g was obtained with good rate kinetics and cycling stability over 50 cycles. Using a synergistic approach combining experimental and ab initio DFT analysis, the structural, magnetic, electronic, and electrochemical properties and the structural limitation to extract full capacity have been described.

Peierls-Nabarro model of interfacial misfit dislocation: An analytic solution

We propose a method to treat the interfacial misfit dislocation array following the original Peierls-Nabarro's ideas. A simple and exact analytic solution is derived in the extended Peierls-Nabarro's model, and this solution reflects the core structure and the energy of misfit dislocation, which depend on misfit and bond strength. We also find that only with beta < 0.2 the structure of interface can be represented by an array of singular Volterra dislocations, which conforms to those of atomic simulation. Interfacial energy and adhesive work can be estimated by inputting ab initio calculation data into the model, and this shows the method can provide a correlation between the ab initio calculations and elastic continuum theory.