The structural pathway from its solvated molecular state to the solution crystallisation of the α- and ß-polymorphic forms of para amino benzoic acid.

Para amino benzoic acid (PABA) has two well-characterised α- and ß-polymorphic forms and, whilst both crystallise in the monoclinic space group P21/n, they have quite different crystal chemistry and crystallisability behaviour. Previous work has shown that the molecular conformation deformation energy in the crystalline state is higher for the ß-form than for the α-form and that the lattice energy for the former converges more slowly than for the latter overall. This suggests that not only is there a higher barrier to crystallisation for the ß-form but also that low solution supersaturations might be needed for it to preferentially nucleate. Additionally, solute cluster propensity and solute solvation energetic analysis highlight the importance of an aqueous solvation environment in inhibiting the α-form's strong OH⋯O carboxylic acid hydrogen bond (H-bond) dimer. Despite this, the detailed molecular-scale pathway from solvated molecules to 3D crystallographic structure still remains unclear, most notably regarding how the nucleation process is activated and how, as a result, this mediates the preferential formation of either of the two polymorphic forms. Molecular dynamics (MD) simulations coupled with FTIR studies and intermolecular synthon analysis address this issue through characterisation of the propensity of the incipient bulk synthons that are important in the crystallisation of the two polymorphic forms within the solution state. MD molecular trajectory analysis within crystallisation solutions reveals a greater propensity for OH⋯O synthons (both single H-bonds and homodimers) typical of the α-form and NH⋯O synthons found in both the α- and ß-forms when compared to aqueous solution but much lower propensities for the ß-form's "fingerprinting" OH⋯N and π-π stacking synthons. In contrast, data from the aqueous solution environment reveals a much greater propensity for the ß-form's π-π interaction synthons. IR dilution studies in acetonitrile in the carbonyl region reveal the presence of two CO vibrational stretching bands, whose relative intensities vary as a function of solution dilution. These were assigned to the solvated PABA monomer and a COOH dimer of PABA. Similar data in ethanol shows a main CO stretching band with a shoulder peak suggesting a similar monomer vs. dimer speciation may exist in this solvent. The IR data is consistent with the organic solvent MD data, albeit the corresponding analysis for the aqueous solution was precluded due to the latter's strong OH vibrational mode which restricted validation in aqueous solutions.

Polyguluronate simulations shed light onto the therapeutic action of OligoG CF-5/20.

Chronic mucoid P. aeruginosa cystic fibrosis (CF) lung infections are associated with the development of a biofilm composed of anionic acetylated exopolysaccharide (EPS) alginate, electrostatically stabilised by extracellular Ca2+ ions. OligoG CF-5/20, a low molecular weight guluronate rich oligomer, is emerging as a novel therapeutic capable of disrupting mature P. aeruginosa biofilms. However, its method of therapeutic action on the mucoid biofilm EPS is not definitively known at a molecular level. This work, utilising molecular dynamics (MD) and Density-Functional Theory (DFT), has revealed that OligoG CF-5/20 interaction with the EPS is facilitated solely through bridging Ca2+ ions, which are not liberated from their native EPS binding sites upon OligoG CF-5/20 dispersal, suggesting that OligoG CF-5/20 does not cause disruptions to mature P. aeruginosa biofilms through breaking EPS-Ca2+-EPS ionic cross-links. Rather it is likely that the therapeutic activity arises from sequestering free Ca2+ ions and preventing further Ca2+ induced EPS aggregation.

A QM/MM Study of Nitrite Binding Modes in a Three-Domain Heme-Cu Nitrite Reductase.

Copper-containing nitrite reductases (CuNiRs) play a key role in the global nitrogen cycle by reducing nitrite (NO2-) to nitric oxide, a reaction that involves one electron and two protons. In typical two-domain CuNiRs, the electron is acquired from an external electron-donating partner. The recently characterised Rastonia picketti (RpNiR) system is a three-domain CuNiR, where the cupredoxin domain is tethered to a heme c domain that can function as the electron donor. The nitrite reduction starts with the binding of NO2- to the T2Cu centre, but very little is known about how NO2- binds to native RpNiR. A recent crystallographic study of an RpNiR mutant suggests that NO2- may bind via nitrogen rather than through the bidentate oxygen mode typically observed in two-domain CuNiRs. In this work we have used combined quantum mechanical/molecular mechanical (QM/MM) methods to model the binding mode of NO2- with native RpNiR in order to determine whether the N-bound or O-bound orientation is preferred. Our results indicate that binding via nitrogen or oxygen is possible for the oxidised Cu(II) state of the T2Cu centre, but in the reduced Cu(I) state the N-binding mode is energetically preferred.

DL_ANALYSER Notation for Atomic Interactions (DANAI): A Natural Annotation System for Molecular Interactions, Using Ethanoic Acid Liquid as a Test Case.

The DL_ANALYSER Notation for Atomic Interactions, DANAI, is the notation syntax to describe interactions between molecules. This notation can annotate precisely the detailed atomistic interactions without having to resolve to diagrammatic illustrations, and yet can be interpreted easily by both human users and computational means. By making use of the DL_F Notation, a universal atom typing scheme for molecular simulations, DANAI contains the expression of atomic species in a natural chemical sense. It is implemented within DL_ANALYSER, a general analysis software program for DL_POLY molecular dynamics simulation software. By making references to the molecular dynamics simulations of pure ethanoic acid liquid, it is shown that DL_ANALYSER can identify and distinguish a variety of hydrogen bond and hydrophobic contact networks through the use of the DANAI expression. It was found that the carboxylic groups preferentially orientated in a "head-to-tail" conformation to form hydrogen bonds between the carbonyl oxygen and hydroxyl hydrogen, resulting in a series of linear structures that intertwined with pockets of methyl clusters.

Descriptions and Implementations of DL_F Notation: A Natural Chemical Expression System of Atom Types for Molecular Simulations.

DL_F Notation is an easy-to-understand, standardized atom typesetting expression for molecular simulations for a range of organic force field (FF) schemes such as OPLSAA, PCFF, and CVFF. It is implemented within DL_FIELD, a software program that facilitates the setting up of molecular FF models for DL_POLY molecular dynamics simulation software. By making use of the Notation, a single core conversion module (the DL_F conversion Engine) implemented within DL_FIELD can be used to analyze a molecular structure and determine the types of atoms for a given FF scheme. Users only need to provide the molecular input structure in a simple xyz format and DL_FIELD can produce the necessary force field file for DL_POLY automatically. In commensurate with the development concept of DL_FIELD, which placed emphasis on robustness and user friendliness, the Engine provides a single-step solution to setup complex FF models. This allows users to switch from one of the above-mentioned FF seamlessly to another while at the same time provides a consistent atom typing that is expressed in a natural chemical sense.

Cellulose Iß microfibril interaction with pristine graphene in water: Effects of amphiphilicity by molecular simulation.

Graphene-cellulose interactions have considerable potential in the development of new materials. In previous computational work (Biomacromolecules2016, 16, 1771), we predicted that the model 100 hydrophobic surface of cellulose interacted favourably with pristine graphene in aqueous solution molecular dynamics simulations; conversely, a model of the hydrophilic 010 surface of cellulose exhibited progressive rearrangement to present a more hydrophobic face with the graphene, with weakened hydrogen bonds between cellulose chains and partial permeation of water. Here, we extend this work by simulating the interaction in aqueous solution of the amphiphilic 110 surface of a cellulose Iß microfibril model, comprising 36 chains of 40 glucosyl residues, with an infinite sheet of pristine graphene. This face of the microfibril is of intermediate hydrophilicity and progressively associates with graphene over replicate simulations. As cellulose chains adhere to the graphene surface, forming interactions via its CH and OH groups, we observe a degree of local and global untwisting of the microfibril. Complementary rippling of the graphene surface is also observed, as it adapts to interaction with the microfibril. This adsorption process is accompanied by increased exclusion of water between cellulose and graphene although some water localises between chains at the immediate interface. The predicted propensity of a cellulose microfibril to adsorb spontaneously on the graphene surface, with mutual structural accommodation, highlights the amphiphilic nature of cellulose and the types of interactions that can be harnessed to design new graphene-carbohydrate biopolymer materials.

Data on the intermolecular interactions of 1,1,1,2-tetrafluoroethane liquids from molecular dynamics simulations.

Detailed atomistic interactions of 1,1,1,2-tetrafluoroethane (HFA-134a) liquid were presented in a data format, namely, DL_ANALYSER Notation for Atomic Interactions (DANAI), that annotates precisely the nature of interactions that is discoverable and searchable without having to resolve to diagrammatic illustrations. The datasets were obtained from raw atomic trajectory files of HFA-134a pure liquid models produced by using DL_POLY molecular dynamics software package. The trajectory datafiles contain expressions of atomic species in a natural chemical sense, and hence, provide localized key interactions, 'at a glance', of the liquid model on otherwise a typically disordered system consists of complex network of intermolecular interactions. The data provide insights to detailed structural behavior of molecules in liquid phase, and can be used as cheminformatics comparative investigations, linking to other molecular system models that contain similar interaction types and chemical species. This can form the foundation of investigations into the role of HFA-134a plays within different applications. For example, it can be used to compare structural and atomic interaction differences with alternative refrigerants, or as liquid propellants in pharmaceutical devices when solvating formulation ingredients.

Interactive molecular dynamics in virtual reality for modelling materials and catalysts.

Interactive molecular dynamics simulation in virtual reality (iMD-VR) is emerging as a promising technique in molecular science. Here, we demonstrate its use in a range of fifteen applications in materials science and heterogeneous catalysis. In this work, the iMD-VR package Narupa is used with the MD package, DL_POLY [1]. We show how iMD-VR can be used to: (i) investigate the mechanism of lithium fast ion conduction by directing the formation of defects showing that vacancy transport is favoured over interstitialcy mechanisms, and (ii) guide a molecule through a zeolite pore to explore diffusion within zeolites, examining in detail the motion of methyl n-hexanoate in H-ZSM-5 zeolite and identifying bottlenecks restricting diffusion. iMD-VR allows users to manipulate these systems intuitively, to drive changes in them and observe the resulting changes in structure and dynamics. We make these simulations available, as a resource for both teaching and research. All simulation files, with videos, can be found online (https://doi.org/10.5281/zenodo.8252314) and are provided as open-source material.

Atomic-scale interactions between quorum sensing autoinducer molecules and the mucoid P. aeruginosa exopolysaccharide matrix.

Mucoid Pseudomonas aeruginosa is a prevalent cystic fibrosis (CF) lung coloniser whose chronicity is associated with the formation of cation cross-linked exopolysaccharide (EPS) matrices, which form a biofilm that acts as a diffusion barrier, sequestering cationic and neutral antimicrobials, and making it extremely resistant to pharmacological challenge. Biofilm chronicity and virulence of the colony is regulated by quorum sensing autoinducers (QSAIs), small signalling metabolites that pass between bacteria, through the biofilm matrix, regulating genetic responses on a population-wide scale. The nature of how these molecules interact with the EPS is poorly understood, despite the fact that they must pass through EPS matrix to reach neighbouring bacteria. Interactions at the atomic-scale between two QSAI molecules, C4-HSL and PQS-both utilised by mucoid P. aeruginosa in the CF lung-and the EPS, have been studied for the first time using a combined molecular dynamics (MD) and density functional theory (DFT) approach. A large-scale, calcium cross-linked, multi-chain EPS molecular model was developed and MD used to sample modes of interaction between QSAI molecules and the EPS that occur at physiological equilibrium. The thermodynamic stability of the QSAI-EPS adducts were calculated using DFT. These simulations provide a thermodynamic rationale for the apparent free movement of C4-HSL, highlight key molecular functionality responsible for EPS binding and, based on its significantly reduced mobility, suggest PQS as a viable target for quorum quenching.

Triglycine (GGG) Adopts a Polyproline II (pPII) Conformation in Its Hydrated Crystal Form: Revealing the Role of Water in Peptide Crystallization.

Polyproline II (pPII) is a left-handed 31-helix conformation, which has been observed to be the most abundant secondary structure in unfolded peptides and proteins compared to α-helix and ß-sheet. Although pPII has been reported as the most stable conformation for several unfolded short chain peptides in aqueous solution, it is rarely observed in their solid state. Here, we show for the first time a glycine homopeptide (gly-gly-gly) adopting the pPII conformation in its crystalline dihydrate structure. The single crystal X-ray structure with molecular dynamic simulation suggests that a network of water and the charged carboxylate group is critical in stabilizing the pPII conformation in solid state, offering an insight into the structures of unfolded regions of proteins and the role of water in peptide crystallization.

Enzyme catalysis captured using multiple structures from one crystal at varying temperatures.

High-resolution crystal structures of enzymes in relevant redox states have transformed our understanding of enzyme catalysis. Recent developments have demonstrated that X-rays can be used, via the generation of solvated electrons, to drive reactions in crystals at cryogenic temperatures (100 K) to generate 'structural movies' of enzyme reactions. However, a serious limitation at these temperatures is that protein conformational motion can be significantly supressed. Here, the recently developed MSOX (multiple serial structures from one crystal) approach has been applied to nitrite-bound copper nitrite reductase at room temperature and at 190 K, close to the glass transition. During both series of multiple structures, nitrite was initially observed in a 'top-hat' geometry, which was rapidly transformed to a 'side-on' configuration before conversion to side-on NO, followed by dissociation of NO and substitution by water to reform the resting state. Density functional theory calculations indicate that the top-hat orientation corresponds to the oxidized type 2 copper site, while the side-on orientation is consistent with the reduced state. It is demonstrated that substrate-to-product conversion within the crystal occurs at a lower radiation dose at 190 K, allowing more of the enzyme catalytic cycle to be captured at high resolution than in the previous 100 K experiment. At room temperature the reaction was very rapid, but it remained possible to generate and characterize several structural states. These experiments open up the possibility of obtaining MSOX structural movies at multiple temperatures (MSOX-VT), providing an unparallelled level of structural information during catalysis for redox enzymes.

Understanding gas capacity, guest selectivity, and diffusion in porous liquids.

Porous liquids are a new class of material that could have applications in areas such as gas separation and homogeneous catalysis. Here we use a combination of measurement techniques, molecular simulations, and control experiments to advance the quantitative understanding of these liquids. In particular, we show that the cage cavities remain unoccupied in the absence of a suitable guest, and that the liquids can adsorb large quantities of gas, with gas occupancy in the cages as high as 72% and 74% for Xe and SF6, respectively. Gases can be reversibly loaded and released by using non-chemical triggers such as sonication, suggesting potential for gas separation schemes. Diffusion NMR experiments show that gases are in dynamic equilibrium between a bound and unbound state in the cage cavities, in agreement with recent simulations for related porous liquids. Comparison with gas adsorption in porous organic cage solids suggests that porous liquids have similar gas binding affinities, and that the physical properties of the cage molecule are translated into the liquid state. By contrast, some physical properties are different: for example, solid homochiral porous cages show enantioselectivity for chiral aromatic alcohols, whereas the equivalent homochiral porous liquids do not. This can be attributed to a loss of supramolecular organisation in the isotropic porous liquid.

Active-site protein dynamics and solvent accessibility in native Achromobacter cycloclastes copper nitrite reductase.

Microbial nitrite reductases are denitrifying enzymes that are a major component of the global nitrogen cycle. Multiple structures measured from one crystal (MSOX data) of copper nitrite reductase at 240 K, together with molecular-dynamics simulations, have revealed protein dynamics at the type 2 copper site that are significant for its catalytic properties and for the entry and exit of solvent or ligands to and from the active site. Molecular-dynamics simulations were performed using different protonation states of the key catalytic residues (AspCAT and HisCAT) involved in the nitrite-reduction mechanism of this enzyme. Taken together, the crystal structures and simulations show that the AspCAT protonation state strongly influences the active-site solvent accessibility, while the dynamics of the active-site 'capping residue' (IleCAT), a determinant of ligand binding, are influenced both by temperature and by the protonation state of AspCAT. A previously unobserved conformation of IleCAT is seen in the elevated temperature series compared with 100 K structures. DFT calculations also show that the loss of a bound water ligand at the active site during the MSOX series is consistent with reduction of the type 2 Cu atom.

Study of interactions between polymer nanoparticles and cell membranes at atomistic levels.

Knowledge of how the structure of nanoparticles and the interactions with biological cell membranes is important not only for understanding nanotoxicological effects on human, animal health and the environment, but also for better understanding of nanoparticle fabrication for biomedical applications. In this work, we use molecular modelling techniques, namely molecular dynamics (MD) simulations, to explore how polymer nanoparticles interact with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid cell membranes. Two different polymers have been considered: 100 monomer units of polyethylene (approx. 2.83 kDa) and polystyrene (approx. 10.4 kDa), both of which have wide industrial applications. We found that, despite the polar lipid head groups acting as an effective barrier to prevent the nanoparticles from interacting with the membrane surface, irreversible adhesion can be initiated by insertion of dangling chain ends from the polymer into the hydrophobic interior of the membrane. In addition, alignment of chain segments from the polymers with that of hydrocarbon chains in the interior of the membrane facilitates the complete immersion of the nanoparticles into the cell membrane. These findings highlight the importance of the surface and the topological structures of the polymer particles that dictate the absorption behaviour into the membrane and, subsequently, induce the possible translocation into the cell.

Molecular dynamics using atomic-resolution structure reveal structural fluctuations that may lead to polymerization of human Cu-Zn superoxide dismutase.

Mutations of the gene encoding Cu-Zn superoxide dismutase (SOD1) cause 20% of the familial cases of the progressive neurodegenerative disease ALS. A growing body of evidence suggests that in familial ALS (FALS) it is the molecular behavior of the metal-depleted SOD1 dimer that leads to a gain of toxic properties by misfolding, unfolding, and aggregation. Structural studies have so far provided static snapshots on the behavior of the wild-type enzyme and some of the FALS mutants. New approaches are required to map out the structural trajectories of the molecule. Here, using our 1.15-A resolution structure of fully metallated human SOD1 and highly parallelized molecular dynamics code on a high-performance capability computer, we have undertaken molecular dynamics calculations to 4,000 ps to reveal the first stages of misfolding caused by metal deletion. Large spatial and temporal fluctuations of the "electrostatic" and "Zn-binding" loops adjacent to the metal-binding sites are observed in the apo-enzyme relative to the fully metallated dimer. These early misfolding events expose the beta-barrels of the dimer to the external environment, allowing close interactions with adjacent molecules. Protection of the beta-edge of the protein can be partially restored by incorporating a single Zn molecule per dimer. These calculations reveal an essential step in the formation of the experimentally observed self-aggregations of metal-depleted FALS mutant SOD1. This result also has implications for the role of demetallated wild-type SOD1 in sporadic cases of ALS, for which the molecular cause still remains undiscovered.

Atomistic studies of surface adhesions using molecular-dynamics simulations.

The existing theoretical descriptions of continuum surface adhesions, such as the JKR (Johnson-Kendall-Roberts) model, have been very useful for the interpretation of particle contacts of sizes down to micrometre length-scales. However, the continuum model is expected to fail at atomic length-scales, where discrete atomistic interactions become significant. The crossover length-scales at which the macroscopic phenomena begin to manifest themselves, and how they occur, are equally baffling and remain poorly understood. This paper explores the issue of atomistic surface contacts between two ionic bodies made of similar materials, such as MgO and NaCl, using molecular dynamics. A range of surface phenomenological behaviour has been identified, from simple fracture to plastic dislocation and neck formation, which are not predicted in continuum models. In addition, the influences of body structure with respect to surface contacts will also be discussed briefly.