Characterization of a human thyroid microtissue model for testing thyroid disrupting chemicals.

Perturbation of thyroid hormone (T4) synthesis is known to cause numerous developmental, metabolic, and cognitive disorders in humans. Due to species differences in sensitivity to chemical exposures, there is a need for human-based in vitro approaches that recapitulate thyroid cellular architecture and T4 production when screening. To address these limitations, primary human thyrocytes, isolated from healthy adult donor tissues and cryopreserved at passage one (p'1) were characterized for cellular composition, 3D follicular architecture, and thyroglobulin (TG)/T4 expression and inhibition by prototype thyroid disrupting chemicals (TDC). Flow analysis of the post-thaw cell suspension showed >80% EpCAM-positive cells with 10%-50% CD90-positive cells. When seeded onto 96-well Matrigel®-coated plates and treated with bovine thyroid stimulating hormone (TSH), thyrocytes formed 3D microtissues during the initial 4-5 days of culture. The microtissues exhibited a stable morphology and size over a 14-day culture period. TG and T4 production were highest in microtissues when the proportion of CD90-positive cells, seeding density and thyroid stimulating hormone concentrations were between 10%-30%, 6K-12K cells per well, and 0.03-1 mIU/mL, respectively. At maximal TG and T4 production levels, average microtissue diameters ranged between 50 and 200 µm. The T4 IC50 values for two prototype TPO inhibitors, 6-propyl-2-thiouracil and methimazole, were ∼0.7 µM and ∼0.5 µM, respectively, in microtissue cultures treated between days 9 and 14. Overall, p'1 cryopreserved primary human thyrocytes in 3D microtissue culture represent a promising new model system to prioritize potential TDC acting directly on the thyroid as part of a weight-of-evidence hazard characterization.

Challenges of modelling real nanoparticles: Ni@Pt electrocatalysts for the oxygen reduction reaction.

Theoretical/computational methods have been extensively applied to screen possible nano-structures attempting to maximize catalytic and stability properties for applications in electrochemical devices. This work shows that the method used to model core@shell structures is of fundamental importance in order to truly represent the physicochemical changes arising from the formation of a core-shell structure. We demonstrate that using a slab approach for modelling nanoparticles the oxygen adsorption energies are qualitatively well represented. Although this is a good descriptor for the catalytic activity, huge differences are found for the calculated surface stability between the results of a nano-cluster and those of a slab approach. Moreover, for the slab method depending on the geometric properties of the core and their similarity to the elements of the core or shell, contradictory effects are obtained. In order to determine the changes occurring as the number of layers and nano particle size are increased, clusters of Ni@Pt from 13 to 260 atoms were constructed and analyzed in terms of geometric parameters, oxygen adsorption, and dissolution potential shift. It is shown that the results of modelling the Ni@Pt nanoparticles with a cluster approach are in good agreement with experimental geometrical parameters, catalytic activity, and stability of a carefully prepared series of Ni@Pt nanostructures where the shell thickness is systematically changed. The maximum catalytic activity and stability are found for a monolayer of Pt whereas adding a second and third layer the behavior is almost the same than that in pure Pt nanoparticles.

DFT analysis of Li intercalation mechanisms in the Fe-phthalocyanine cathode of Li-ion batteries.

New materials with high intercalation capacity are needed for cathodic materials in order to overcome small capacities at high discharge rates in Li-ion batteries. High intercalation capacities have been reported in the experimental setup using iron phthalocyanine (FePc) as cathodic material; however the real intercalation capacity and the chemistry occurring during the intercalation process are still being debated. In this work we analyze the intercalation of Li atoms in FePc periodic structures using density functional theory including a semi-empirical approach to represent van der Waals (vdW) forces. Within this approach we find intercalation capacities of about 20 Li atoms per FePc molecule at a discharge voltage of ~0.5 V (with respect to Li/Li(+)), and up to 37 Li atoms at lower voltages. The intercalation process is driven mainly by electrostatic interactions between positively charged Li ions and negatively charged FePc molecules, with vdW interactions playing an essential role in reaching the high number of intercalated Li atoms. The reduction of the central Fe atom leading to charges evolving from +1.2 to -0.2 is responsible for the high intercalation voltage; however the further reduction contributions of N, C, and even H atoms make FePc a suitable cathode for Li-ion battery applications.

Interactions of platinum clusters with a graphite substrate.

Density functional theory calculations are used to elucidate the interactions of platinum clusters with graphite; the results are analyzed in terms of geometry, and energetic and electronic properties. Adsorption of platinum clusters from 1 to 38 atoms is evaluated on a 3-layer graphite model structure. The approach includes van der Waals interactions, which have proved to be essential to describe relatively weak interactions. The results show that when interacting with graphite, the clusters tend to slightly wetting the surface. Although the effect is more pronounced in the larger clusters investigated, the energy difference among total, partial, and non-wetting structures in small clusters is very low and may be easily overcome by thermal effects. The van der Waals energy contributes to enhance the graphite-cluster strength and is proportional to the number of interacting atoms at the interface. Small charge transfer takes place from the metallic cluster to the graphite surface and the cluster becomes polarized, with positive values at the interface, and negative values in the top. The interaction with graphite enhances the metallic character of the cluster as shown by density of states analyses. New states resulting from the interaction between graphite and the metal cluster may modify its catalytic behavior.

Hydrogen storage based on physisorption.

Physisorption of molecular hydrogen based on neutral and negatively charged aromatic molecular systems has been evaluated using ab initio calculations to estimate the binding energy, DeltaH, and DeltaG at 298 ( approximately 77 bar) and 77 K (45 bar) in order to compare calculated results with experimental measurements of hydrogen adsorption. The molecular systems used in this study were corannulene (C(20)H(10)), dicyclopenta[def,jkl]triphenylene (C(20)H(10)), 5,8-dioxo-5,8-dihydroindeno[2,1-c]fluorene (C(20)H(10)O(2)), 6-hexyl-5,8-dioxo-5,8-dihydroindeno[2,1-c]fluorene (C(26)H(22)O(2)), coronene (C(24)H(12)), dilithium phthalocyanine (Li(2)Pc, C(32)H(16)Li(2)N(8)), tetrabutylammonium lithium phthalocyanine (TBA-LiPc, C(48)H(52)LiN(9)), and tetramethylammonium lithium phthalocyanine (TMA-LiPc, C(36)H(28)LiN(9)). It was found (a) that the calculated term that corrects 0 K electronic energies to give Gibbs energies (thermal correction to Gibbs energy, TCGE) serves as a good approximation of the adsorbent binding energy required in order for a physisorption process to be thermodynamically allowed and (b) that the binding energy for neutral aromatic molecules varies as a function of curvature (e.g., corannulene versus coronene) or if electron-withdrawing or -donating groups are part of the adsorbent. A negatively charged aromatic ring, the lithium phthalocyanine complex anion, [LiPc](-), introduces charge-induced dipole interactions into the adsorption process, resulting in a doubling of the binding energy of Li(2)Pc relative to corannulene. Experimental hydrogen adsorption results for Li(2)Pc, which are consistent with MD simulation results using chi-Li(2)Pc to simulate the adsorbent, suggest that only one side of the phthalocyanine ring is used in the adsorption process. The introduction of a tetrabutylammonium cation as a replacement for one lithium ion in Li(2)Pc has the effect of increasing the number of hydrogen molecules adsorbed from 10 (3.80 wt %) for Li(2)Pc to 24 (5.93 wt %) at 77 K and 45 bar, suggesting that both sides of the phthalocyanine ring are available for hydrogen adsorption. MD simulations of layered tetramethylammonium lithium phthalocyanine molecular systems illustrate that doubling the wt % H(2) adsorbed is possible via such a system. Ab initio calculations also suggest that layered or sandwich structures can result in significant reductions in the pressure required for hydrogen adsorption.

Addition reactions of alkyl and carboxyl radicals to vinylidene fluoride.

The addition reactions of alkyl radicals CF3* and CH3* and carboxyl radicals C2H5O*, C2H5OCOO*, CF3COO*, and CH3COO* to a vinylidene fluoride (VDF) molecule are studied using ab initio calculations. These radicals were selected because they are intermediate or final products of diacyl peroxides decomposition in the initiation reactions of VDF polymerization. Two combinations of methods for energetics and structure optimization are applied: QCISD/6-311G(d,p)//HF/6-31G(d) and B3LYP/6-311G+(3df, 2p)//B3LYP/6-31G(d). It is found that the formed bond length of the product, the forming bond length of the transition state, and the attack angle of the product structures are not sensitive to the level of theory even though the attack angle of the transition state structures is. Early transition states are obtained upon attack at both high-substituted and nonsubstituted carbon atom VDF ends. Kinetic and thermodynamic control rules play different roles on governing the reactivity of the addition with the studied radicals. Both theoretical methods yield the same trends for the preferential attack site in terms of regioselectivity, barrier energies, and reaction enthalpies. It is shown that the addition reactions of the intermediate radicals C2H5OCOO*, CF3COO*, and CH3COO* of the decomposition of diethyl peroxydicarbonate, trifluoroacetyl peroxide, and diacetyl peroxide initiators yield smaller energy barriers than the additions of the corresponding final radicals, C2H5O*, CF3*, and CH3*; therefore, the reactions of the intermediate radicals should not be ignored when analyzing the initiation process of the VDF polymerization using those initiators.

Effects of confinement on small water clusters structure and proton transport.

Analyses of the structure of two to four water molecule clusters confined between two benzene and between two naphthalene molecules have been performed using ab initio methods. The water clusters tend to maximize the number of hydrogen bonds via formation of a cyclic network. The oxygen atoms locate approximately in the middle of the confined geometry, and the dipole vectors arrange either parallel or pointing to the surfaces. Energy barriers for proton transfer calculated for H3O+-(H2O) complexes in the same confined geometries suggest that there is a specific range of confinement that helps to lower the energy barriers of the proton transfer. When the walls are too close to each other, at a separation of 4 A, the energy barriers are extremely high. Confinement does not lower the barrier energies of proton transfer when the H3O+-(H2O) complexes are located further from each of the surfaces by more than approximately 8 A.

Computational investigation of adsorption of molecular hydrogen on lithium-doped corannulene.

Density functional theory and classical molecular dynamics simulations are used to investigate the prospect of lithium-doped corannulene as adsorbent material for H(2) gas. Potential energy surface scans at the level of B3LYP/6-311G(d,p) show an enhanced interaction of molecular hydrogen with lithium-atom-doped corannulene complexes with respect to that found in undoped corannulene. MP2(FC)/6-31G(d,p) optimizations of 4H(2)-(Li(2)-C(20)H(10)) yield H(2) binding energies of -1.48 kcal/mol for the H(2)-Li interaction and -0.92 kcal/mol for the H(2)-C interaction, whereas values of -0.94 and -0.83 kcal/mol were reported (J. Phys. Chem. B 2006, 110, 7688-7694) for physisorption of H(2) on the concave and the convex side of corannulene using MP2(full)/6-31G(d), respectively. Classical molecular dynamics simulations predict hydrogen uptakes in Li-doped corannulene assemblies that are significantly enhanced with respect to that found in undoped molecules, and the hydrogen uptake ability is dependent on the concentration of lithium dopant. For the Li(6)-C(20)H(10) complex, a hydrogen uptake of 4.58 wt % at 300 K and 230 bar is obtained when the adsorbent molecules are arranged in stack configurations separated by 6.5 A, and with interlayer distances of 10 A, hydrogen uptake reaches 6.5 wt % at 300 K and 215 bar.

Theoretical study on the properties of linear and cyclic amides in gas phase and water solution.

The structural and energetic properties of a group of selected amides, of well-known importance for the design of efficient clathrate inhibitors, are calculated with Hartree-Fock and density functional theory, B3LYP, theoretical levels, and a 6-311++g** basis set in the gas phase and a water solution. The conformational behavior of the molecules is studied through the scanning of the torsional potential energy surfaces and by the analysis of the differences in the energetic and structural properties between the isomers. The properties of the amides in water solution are determined within a self-consistent reaction field approach with a polarizable continuum model that allows the calculation of the different contributions to the free energy of solvation. The calculated barriers to rotation are in good agreement with the available experimental data and the comparison of the gas and water results shows the strong effect of the solute polarization. The properties of different amide-water complexes are calculated and compared with available experimental information.

Investigation of corannulene for molecular hydrogen storage via computational chemistry and experimentation.

Molecular simulations for hydrogen physisorption with corannulene molecules arranged according to their crystal structure result in good agreement with the weight-percent hydrogen stored as determined experimentally employing a 3-g sample of highly crystalline corannulene at ambient temperatures and 72 bar of pressure. Calculated enthalpies of adsorption for corannulene/hydrogen molecular systems obtained from ab initio calculations which take into account electron correlation via second-order Möller-Plesset perturbation theory are in good agreement with literature experimental enthalpies of adsorption for activated carbons interacting with molecular hydrogen. Ab initio results also show that corannulene molecules arranged in a sandwich structure are important for approximately doubling the binding energy of corannulene interacting with molecular hydrogen through a cooperative interaction. To test the effects of finite temperatures and pressures, stack arrays were used as input for molecular dynamics simulations and indicate that physisorption mechanisms including van der Waals forces and dipole-induced dipole interactions may yield enhanced adsorption capacity in relation to other carbon-based materials. These results will be instrumental in identifying interlayer separations of an array of corannulene or related molecules that may provide a high weight percent of physisorbed hydrogen.

Understanding catalysed growth of single-wall carbon nanotubes.

Classical molecular dynamics simulations using a reactive force field, which allows simulation of bond-breaking and bond-forming, are carried out to investigate the several stages of a catalysed synthesis process of single-wall carbon nanotubes. The simulations assume instantaneous catalysis of a precursor gas on the surface of metallic nanoclusters, illustrating how carbon atoms dissolve in the metal cluster and then precipitate on its surface, evolving into various carbon structures, finally forming a cap which eventually grows to a single-wall nanotube. The results are discussed in the context of experimental synthesis results.

Theoretical studies to understand surface chemistry on carbon anodes for lithium-ion batteries: reduction mechanisms of ethylene carbonate.

Reductive decomposition mechanisms for ethylene carbonate (EC) molecule in electrolyte solutions for lithium-ion batteries are comprehensively investigated using density functional theory. In gas phase the reduction of EC is thermodynamically forbidden, whereas in bulk solvent it is likely to undergo one- as well as two-electron reduction processes. The presence of Li cation considerably stabilizes the EC reduction intermediates. The adiabatic electron affinities of the supermolecule Li(+)(EC)n (n = 1-4) successively decrease with the number of EC molecules, independently of EC or Li(+) being reduced. Regarding the reductive decomposition mechanism, Li(+)(EC)n is initially reduced to an ion-pair intermediate that will undergo homolytic C-O bond cleavage via an approximately 11.0 kcal/mol barrier, bringing up a radical anion coordinated with Li(+). Among the possible termination pathways of the radical anion, thermodynamically the most favorable is the formation of lithium butylene dicarbonate, (CH2CH2OCO2Li)2, followed by the formation of one O-Li bond compound containing an ester group, LiO(CH2)2CO2(CH2)2OCO2Li, then two very competitive reactions of the further reduction of the radical anion and the formation of lithium ethylene dicarbonate, (CH2OCO2Li)2, and the least favorable is the formation of a C-Li bond compound (Li carbides), Li(CH2)2OCO2Li. The products show a weak EC concentration dependence as has also been revealed for the reactions of LiCO3(-) with Li(+)(EC)n; that is, the formation of Li2CO3 is slightly more favorable at low EC concentrations, whereas (CH2OCO2Li)2 is favored at high EC concentrations. On the basis of the results presented here, in line with some experimental findings, we find that a two-electron reduction process indeed takes place by a stepwise path. Regarding the composition of the surface films resulting from solvent reduction, for which experiments usually indicate that (CH2OCO2Li)2 is a dominant component, we conclude that they comprise two leading lithium alkyl bicarbonates, (CH2CH2OCO2Li)2 and (CH2OCO2Li)2, together with LiO(CH2)2CO2(CH2)2OCO2Li, Li(CH2)2OCO2Li and Li2CO3.