PyThinFilm: Automated Molecular Dynamics Simulation Protocols for the Generation of Thin Film Morphologies.

The performance of organic optoelectronic devices, such as organic light-emitting diodes (OLEDs) and organic solar cells (OSCs), is intrinsically related to the molecular-scale morphology of the thin films from which they are composed. However, the experimental characterization of morphology at the molecular level is challenging due to the often amorphous or at best semicrystalline nature of these films. Classical molecular modeling techniques, such as molecular dynamics (MD) simulation, are increasingly used to understand the relationship between morphology and the properties of thin-film devices. PyThinFilm (github.com/ATB-UQ/PyThinFilm) is an open-source Python package which allows fully automated MD simulations of thin film growth to be performed using vacuum and/or solution deposition processes. PyThinFilm utilizes the GROMACS simulation package in combination with interaction parameters from the Automated Topology Builder (atb.uq.edu.au). Here, PyThinFilm is described along with an overview of applications in which PyThinFilm has been used to study the thin films of organic semiconductor materials typically used in OLEDs and OSCs.

OFraMP: a fragment-based tool to facilitate the parametrization of large molecules.

An Online tool for Fragment-based Molecule Parametrization (OFraMP) is described. OFraMP is a web application for assigning atomic interaction parameters to large molecules by matching sub-fragments within the target molecule to equivalent sub-fragments within the Automated Topology Builder (ATB, atb.uq.edu.au) database. OFraMP identifies and compares alternative molecular fragments from the ATB database, which contains over 890,000 pre-parameterized molecules, using a novel hierarchical matching procedure. Atoms are considered within the context of an extended local environment (buffer region) with the degree of similarity between an atom in the target molecule and that in the proposed match controlled by varying the size of the buffer region. Adjacent matching atoms are combined into progressively larger matched sub-structures. The user then selects the most appropriate match. OFraMP also allows users to manually alter interaction parameters and automates the submission of missing substructures to the ATB in order to generate parameters for atoms in environments not represented in the existing database. The utility of OFraMP is illustrated using the anti-cancer agent paclitaxel and a dendrimer used in organic semiconductor devices. OFraMP applied to paclitaxel (ATB ID 35922).

The Molecular Origin of Anisotropic Emission in an Organic Light-Emitting Diode.

Atomistic nonequilibrium molecular dynamics simulations have been used to model the induction of molecular orientation anisotropy within the emission layer of an organic light-emitting diode (OLED) formed by vapor deposition. Two emitter species were compared: racemic fac-tris(2-phenylpyridine)iridium(III) (Ir(ppy)3) and trans-bis(2-phenylpyridine)(acetylacetonate)iridium(III) (Ir(ppy)2(acac)). The simulations show that the molecular symmetry axes of both emitters preferentially align perpendicular to the surface during deposition. The molecular arrangement formed on deposition combined with consideration of the transition dipole moments provides insight into experimental reports that Ir(ppy)3 shows isotropic emission, while Ir(ppy)2(acac) displays improved efficiency due to an apparent preferential alignment of the transition dipole vectors parallel to the substrate. The simulations indicate that this difference is not due to differences in the extent of emitter alignment, but rather differences in the direction of the transition dipoles within the two complexes.

Validating lipid force fields against experimental data: Progress, challenges and perspectives.

Biological membranes display a great diversity in lipid composition and lateral structure that is crucial in a variety of cellular functions. Simulations of membranes have contributed significantly to the understanding of the properties, functions and behaviour of membranes and membrane-protein assemblies. This success relies on the ability of the force field used to describe lipid-lipid and lipid-environment interactions accurately, reproducibly and realistically. In this review, we present some recent progress in lipid force-field development and validation strategies. In particular, we highlight how a range of properties obtained from various experimental techniques on lipid bilayers and membranes, can be used to assess the quality of a force field. We discuss the limitations and assumptions that are inherent to both computational and experimental approaches and how these can influence the comparison between simulations and experimental data. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Elucidating the Spatial Arrangement of Emitter Molecules in Organic Light-Emitting Diode Films.

The effect of varying the emitter concentration on the structural properties of an archetypal phosphorescent blend consisting of 4,4'-bis(N-carbazolyl)biphenyl and tris(2-phenylpyridyl)iridium(III) has been investigated using non-equilibrium molecular dynamics (MD) simulations that mimic the process of vacuum deposition. By comparison with reflectometry measurements, we show that the simulations provide an accurate model of the average density of such films. The emitter molecules were found not to be evenly distributed throughout film, but rather they can form networks that provide charge and/or energy migration pathways, even at emitter concentrations as low as ≈5 weight percent. At slightly higher concentrations, percolated networks form that span the entire system. While such networks would give improved charge transport, they could also lead to more non-radiative pathways for the emissive state and a resultant loss of efficiency.

Curved or linear? Predicting the 3-dimensional structure of α-helical antimicrobial peptides in an amphipathic environment.

α-Helical membrane-active antimicrobial peptides (AMPs) are known to act via a range of mechanisms, including the formation of barrel-stave and toroidal pores and the micellisation of the membrane (carpet mechanism). Different mechanisms imply that the peptides adopt different 3D structures when bound at the water-membrane interface, a highly amphipathic environment. Here, an evolutionary algorithm is used to predict the 3D structure of a range of α-helical membrane-active AMPs at the water-membrane interface by optimising amphipathicity. This amphipathic structure prediction (ASP) is capable of distinguishing between curved and linear peptides solved experimentally, potentially allowing the activity and mechanism of action of different membrane-active AMPs to be predicted. The ASP algorithm is accessible via a web interface at http://atb.uq.edu.au/asp/.

Automated partial atomic charge assignment for drug-like molecules: a fast knapsack approach.

A key factor in computational drug design is the consistency and reliability with which intermolecular interactions between a wide variety of molecules can be described. Here we present a procedure to efficiently, reliably and automatically assign partial atomic charges to atoms based on known distributions. We formally introduce the molecular charge assignment problem, where the task is to select a charge from a set of candidate charges for every atom of a given query molecule. Charges are accompanied by a score that depends on their observed frequency in similar neighbourhoods (chemical environments) in a database of previously parameterised molecules. The aim is to assign the charges such that the total charge equals a known target charge within a margin of error while maximizing the sum of the charge scores. We show that the problem is a variant of the well-studied multiple-choice knapsack problem and thus weakly NP -complete. We propose solutions based on Integer Linear Programming and a pseudo-polynomial time Dynamic Programming algorithm. We demonstrate that the results obtained for novel molecules not included in the database are comparable to the ones obtained performing explicit charge calculations while decreasing the time to determine partial charges for a molecule from hours or even days to below a second. Our software is openly available.

Predicting the Prevalence of Alternative Warfarin Tautomers in Solution.

Warfarin, a widely used oral anticoagulant, is prescribed as a racemic mixture. Each enantiomer of neutral Warfarin can exist in 20 possible tautomeric states leading to complex pharmacokinetics and uncertainty as to the relevant species under different conditions. Here, the ability of alternative computational approaches to predict the preferred tautomeric form(s) of neutral Warfarin in different solvents is examined. It is shown that varying the method used to estimate the heat of formation in vacuum (direct or via homodesmic reactions), whether entropic corrections were included, and the method used to estimate the free enthalpy of solvation (i.e., PCM, COSMO, or SMD implicit models or explicit solvent) lead to large differences in the predicted rank and relative populations of the tautomers. In this case, only a combination of the enthalpy of formation using homodesmic reactions and explicit solvent to estimate the free enthalpy of solvation yielded results compatible with the available experimental data. The work also suggests that a small but significant subset of the possible Warfarin tautomers are likely to be physiologically relevant.

Automated Topology Builder Version 3.0: Prediction of Solvation Free Enthalpies in Water and Hexane.

The ability of atomic interaction parameters generated using the Automated Topology Builder and Repository version 3.0 (ATB3.0) to predict experimental hydration free enthalpies (Δ Gwater) and solvation free enthalpies in the apolar solvent hexane (Δ Ghexane) is presented. For a validation set of 685 molecules the average unsigned error (AUE) between Δ Gwater values calculated using the ATB3.0 and experiment is 3.8 kJ·mol-1. The slope of the line of best fit is 1.00, the intercept -1.0 kJ·mol-1, and the R2 0.90. For the more restricted set of 239 molecules used to validate OPLS3 ( J. Chem. Theory Comput. 2016 , 12 , 281 - 296 , DOI: 10.1021/acs.jctc.5b00864 ) the AUE using the ATB3.0 is just 2.7 kJ·mol-1 and the R2 0.93. A roadmap for further improvement of the ATB parameters is presented together with a discussion of the challenges of validating force fields against the available experimental data.

Effect of methyl-branched fatty acids on the structure of lipid bilayers.

Methyl-branched fatty acids are widespread in prokaryotic membranes. Although anteiso and iso branching (that is on the antepenultimate and penultimate carbons) and the presence of multiple methyl branches in the phytanoyl chain are known to modify the thermotropic behavior and enhance the fluidity of lipid bilayers, little is known about the effect of methyl branching on the structure of lipid bilayers. In this study, molecular dynamics simulations are used to examine systematically the impact of one or more methyl branches at different positions along the sn-1 palmitoyl chain on the structural properties of a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayer. It is found that methyl branching reduces lipid condensation, decreases the bilayer thickness, and lowers chain ordering. Branching also results in the formation of kinks at the branching point, thereby enhancing the fluidity of lipid bilayers. Furthermore, this effect varies in a methyl-position-dependent fashion. In the case of polymethylated chains, the simulations suggest that if the gap between the methyl groups is sufficient (two or three carbons), the effects of the methyl branches are additive and equivalent to the combined effect of the corresponding monomethyl-branched lipids.

Some Like It Hot: The Effect of Sterols and Hopanoids on Lipid Ordering at High Temperature.

Sterols and hopanoids have been suggested to reinforce membranes and protect against unfavorable environmental conditions. In particular, hopanoids are found in high concentrations in membranes of thermotolerant and thermophilic bacteria. However, the mechanism whereby sterols and hopanoids stabilize membranes at elevated temperatures is poorly understood. Here, the effect of temperature on the ordering of lipids in bilayers containing cholesterol or the hopanoids bacteriohopanetetrol and diplopterol was explored using molecular dynamics simulations. It is shown that cholesterol induces a high level of ordering over a wide range of temperatures. Bacteriohopanetetrol promotes order within the lipid tails but enhances fluid-like properties of the head groups at high temperatures. In contrast, diplopterol partitions in the midplane of the bilayer. This suggests that individual hopanoids fulfill distinct functions in membranes, with the ordering properties of bacteriohopanetetrol being particularly well suited to maintain the integrity of membranes at temperatures preferred by thermotolerant and thermophilic bacteria.