Exploring the Effects of Wetting and Free Fatty Acid Deposition on an Atomistic Hair Fiber Surface Model Incorporating Keratin-Associated Protein 5-1.

The complex development of cosmetic and medical formulations relies on an ever-growing accuracy of predictive models of hair surfaces. Hitherto, modeling efforts have focused on the description of 18-methyl eicosanoic acid (18-MEA), the primary fatty acid covalently attached to the hair surface, without explicit modeling of the protein layer. Herein, the molecular details of the outermost surface of the human hair fiber surface, also called the F-layer, were studied using molecular dynamics (MD) simulations. The F-layer is composed primarily of keratin-associated proteins KAP5 and KAP10, which are decorated with 18-MEA on the outer surface of a hair fiber. In our molecular model, we incorporated KAP5-1 and evaluated the surface properties of 18-MEA through MD simulations, resulting in 18-MEA surface density, layer thickness, and tilt angles in agreement with previous experimental and computational studies. Subsequent models with reduced 18-MEA surface density were also generated to mimic damaged hair surfaces. Response to wetting of virgin and damaged hair showed rearrangement of 18-MEA on the surface, allowing for water penetration into the protein layer. To demonstrate a potential use case for these atomistic models, we deposited naturally occurring fatty acids and measured 18-MEA's response in both dry and wet conditions. As fatty acids are often incorporated in shampoo formulations, this work demonstrates the ability to model the adsorption of ingredients on hair surfaces. This study illustrates, for the first time, the complex behavior of a realistic F-layer at the molecular level and opens up the possibility of studying the adsorption behavior of larger, more complex molecules and formulations.

Molecular-Level Examination of Amorphous Solid Dispersion Dissolution.

Amorphous solid dispersions (ASDs) are commonly used to orally deliver small-molecule drugs that are poorly water-soluble. ASDs consist of drug molecules in the amorphous form which are dispersed in a hydrophilic polymer matrix. Producing a high-performance ASD is critical for effective drug delivery and depends on many factors such as solubility of the drug in the matrix and the rate of drug release in aqueous medium (dissolution), which is linked to bioperformance. Often, researchers perform a large number of design iterations to achieve this objective. A detailed molecular-level understanding of the mechanisms behind ASD dissolution behavior would aid in the screening, designing, and optimization of ASD formulations and would minimize the need for testing a wide variety of prototype formulations. Molecular dynamics and related types of simulations, which model the collective behavior of molecules in condensed phase systems, can provide unique insights into these mechanisms. To study the effectiveness of these simulation techniques in ASD formulation dissolution, we carried out dissipative particle dynamics simulations, which are particularly an efficient form of molecular dynamics calculations. We studied two stages of the dissolution process: the early-stage of the dissolution process, which focuses on the dissolution at the ASD/water interface, and the late-stage of the dissolution process, where significant drug release would have occurred and there would be a mixture of drug and polymer molecules in a predominantly aqueous environment. Experimentally, we used Fourier transform infrared spectroscopy to study the interactions between drugs, polymers, and water in the dry and wet states and the chromatographic technique to study the rate of drug and polymer release. Both experiments and simulations provided evidence of polymer microstructures and drug-polymer interactions as important factors for the dissolution behavior of the investigated ASDs, consistent with previous work by Pudlas et al. (Eur. J. Pharm. Sci.2015, 67, 21-31). As experimental and simulation results are consistent and complementary, it is clear that there is significant potential for combined experimental and computational research for a detailed understanding of ASD formulations and, hence, formulation optimization.

Design of Organic Electronic Materials With a Goal-Directed Generative Model Powered by Deep Neural Networks and High-Throughput Molecular Simulations.

In recent years, generative machine learning approaches have attracted significant attention as an enabling approach for designing novel molecular materials with minimal design bias and thereby realizing more directed design for a specific materials property space. Further, data-driven approaches have emerged as a new tool to accelerate the development of novel organic electronic materials for organic light-emitting diode (OLED) applications. We demonstrate and validate a goal-directed generative machine learning framework based on a recurrent neural network (RNN) deep reinforcement learning approach for the design of hole transporting OLED materials. These large-scale molecular simulations also demonstrate a rapid, cost-effective method to identify new materials in OLEDs while also enabling expansion into many other verticals such as catalyst design, aerospace, life science, and petrochemicals.

Characterizing moisture uptake and plasticization effects of water on amorphous amylose starch models using molecular dynamics methods.

Dynamics and thermophysical properties of amorphous starch were explored using molecular dynamics (MD) simulations. Using the OPLS3e force field, simulations of short amylose chains in water were performed to determine force field accuracy. Using well-tempered metadynamics, a free energy map of the two glycosidic angles of an amylose molecule was constructed and compared with other modern force fields. Good agreement of torsional sampling for both solvated and amorphous amylose starch models was observed. Using combined grand canonical Monte Carlo (GCMC)/MD simulations, a moisture sorption isotherm curve is predicted along with temperature dependence. Concentration-dependent activation energies for water transport agree quantitatively with previous experiments. Finally, the plasticization effect of moisture content on amorphous starch was investigated. Predicted glass transition temperature (Tg) depression as a function of moisture content is in line with experimental trends. Further, our calculations provide a value for the dry Tg for amorphous starch, a value which no experimental value is available.

Massive Theoretical Screen of Hole Conducting Organic Materials in the Heteroacene Family by Using a Cloud-Computing Environment.

Materials exhibiting higher mobilities than conventional organic semiconducting materials such as fullerenes and fused thiophenes are in high demand for applications in printed electronics. To discover new molecules in the heteroacene family that might show improved charge mobility, a massive theoretical screen of hole conducting properties of molecules was performed by using a cloud-computing environment. Over 7 000 000 structures of fused furans, thiophenes and selenophenes were generated and 250 000 structures were randomly selected to perform density functional theory (DFT) calculations of hole reorganization energies. The lowest hole reorganization energy calculated was 0.0548 eV for a fused thioacene having 8 aromatics rings. Hole mobilities of compounds with the lowest 130 reorganization energy were further processed by applying combined DFT and molecular dynamics (MD) methods. The highest mobility calculated was 1.02 and 9.65 cm2/(V s) based on percolation and disorder theory, respectively, for compounds containing selenium atoms with 8 aromatic rings. These values are about 20 times higher than those for dinaphthothienothiophene (DNTT).

E-type delayed fluorescence of a phosphine-supported Cu2(mu-NAr2)2 diamond core: harvesting singlet and triplet excitons in OLEDs.

A highly emissive bis(phosphine)diarylamido dinuclear copper(I) complex (quantum yield = 57%) was shown to exhibit E-type delayed fluorescence by variable temperature emission spectroscopy and photoluminescence decay measurement of doped vapor-deposited films. The lowest energy singlet and triplet excited states were assigned as charge transfer states on the basis of theoretical calculations and the small observed S(1)-T(1) energy gap. Vapor-deposited OLEDs doped with the complex in the emissive layer gave a maximum external quantum efficiency of 16.1%, demonstrating that triplet excitons can be harvested very efficiently through the delayed fluorescence channel. The function of the emissive dopant in OLEDs was further probed by several physical methods, including electrically detected EPR, cyclic voltammetry, and photoluminescence in the presence of applied current.

A short yet very weak dative bond: structure, bonding, and energetic properties of N(2)-BH(3).

The structure, bonding, and energetic properties of the N(2)-BH(3) complex are reported as characterized by density functional theory (DFT) and post-Hartree-Fock (HF) calculations. The equilibrium structure of the complex exhibits a short B-N distance near 1.6 A, comparable to that of a strong acid-base complex like H(3)N-BH(3). However, the binding energy is only 5.7 kcal/mol at the CCSD(T)/6-311+G(2df,2dp) level of theory, which is reminiscent of a weak, nonbonded complex. Natural bond orbital (NBO) and atoms in molecules (AIM) analyses of the electron density from both DFT and post-HF calculations do indicate that the extent of charge transfer and covalent character in the B-N dative bond is only somewhat less than in comparable systems with fairly large binding energies (e.g., H(3)N-BH(3) and OC-BH(3)). Energy decomposition analysis indicates key differences between the N(2), CO, and NH(3) complexes, primarily associated with the natures of the lone pairs involved (sp vs sp(3)) and the donor/acceptor characteristics of the relevant occupied and virtual orbitals, both sigma and pi. Also, CCSD/6-311+G(2df,2dp) calculations indicate that the B-N distance potential is rather anharmonic and exhibits a flat, shelf-like region ranging from 2.1 to 2.5 A that lies about 1.5 kcal/mol above the minimum at 1.67 A. However, this region is more sloped and lies about 2.5 kcal/mol above the equilibrium region according to the CCSD(T)/6-311+G(2df,2dp)//CCSD/6-311+G(2df,2dp) potential. A 1D analysis of the vibrational motion along the B-N stretching coordinate in the CCSD/6-311+G(2df,2dp) potential indicates that the average B-N distance in the ground vibrational state is 1.71 A, about 0.04 A longer than the equilibrium distance. Furthermore, the vibrationally averaged distance obtained via an analysis of the CCSD(T)/6-311+G(2df,2dp)//CCSD/6-311+G(2df,2dp) potential was found to be 0.03 A longer than the CCSD(T)/6-311+G(2df,2dp) minimum.

Coumarin-based, electron-trapping iridium complexes as highly efficient and stable phosphorescent emitters for organic light-emitting diodes.

A new class of coumarin-based iridium tris-cyclometalated complexes has been developed. These complexes are highly emissive, with emission colors ranging from green to orange-red. Besides modification of ligand structures, color tuning was realized by incorporation of ligands with different electrochemical properties in a heteroleptic structure. The organic light-emitting diodes (OLEDs) using these compounds as emissive dopants are highly efficient and stable. Unlike other Ir(III) phosphorescent dopants, these coumarin-based Ir(III) dopants can effectively trap and transport electrons in the emissive layer.

Three-dimensional structure determination of N-(p-Tolyl)-dodecylsulfonamide from powder diffraction data and validation of structure using solid-state NMR spectroscopy.

The three-dimensional structure, conformation, and packing of molecules in the solid state are crucial components used in the optimization of many technologically useful materials properties. Single-crystal X-ray diffraction is the traditional and most effective method of determining 3-D structures in the solid state. Obtaining single crystals that are sufficiently large and free of imperfections is often laborious, time-consuming, and, occasionally, impossible. The feasibility of an integrated approach to the determination and verification of a complete three-dimensional structure for a medium-sized organic molecule without using single crystals is demonstrated for the case of an organic stabilizer compound N-(p-tolyl)-dodecylsulfonamide. The approach uses a combination of powder XRD data, several computational packages involving Monte Carlo simulations and ab initio quantum mechanical calculations, and experimental solid-state NMR chemical shifts. Structure elucidation of N-(p-tolyl)-dodecylsulfonamide revealed that the Bravais lattice is monoclinic, with cell dimensions of a = 38.773 A, b = 5.507 A, c = 9.509 A, and beta = 86.35 degrees, and a space group of P21/c.