Efficient simulations of the aqueous bio-interface of graphitic nanostructures with a polarisable model.

To fully harness the enormous potential offered by interfaces between graphitic nanostructures and biomolecules, detailed connections between adsorbed conformations and adsorption behaviour are needed. To elucidate these links, a key approach, in partnership with experimental techniques, is molecular simulation. For this, a force-field (FF) that can appropriately capture the relevant physics and chemistry of these complex bio-interfaces, while allowing extensive conformational sampling, and also supporting inter-operability with known biological FFs, is a pivotal requirement. Here, we present and apply such a force-field, GRAPPA, designed to work with the CHARMM FF. GRAPPA is an efficiently implemented polarisable force-field, informed by extensive plane-wave DFT calculations using the revPBE-vdW-DF functional. GRAPPA adequately recovers the spatial and orientational structuring of the aqueous interface of graphene and carbon nanotubes, compared with more sophisticated approaches. We apply GRAPPA to determine the free energy of adsorption for a range of amino acids, identifying Trp, Tyr and Arg to have the strongest binding affinity and Asp to be a weak binder. The GRAPPA FF can be readily incorporated into mainstream simulation packages, and will enable large-scale polarisable biointerfacial simulations at graphitic interfaces, that will aid the development of biomolecule-mediated, solution-based graphene processing and self-assembly strategies.

Structure of the electrical double layer at aqueous gold and silver interfaces for saline solutions

We report the structure of the electrical double layer, determined from molecular dynamics simulations, for a range of saline solutions (NaCl, KCl, MgCl2 and CaCl2) at both 0.16 and 0.60molkg(-1) on different facets of the gold and silver aqueous interfaces. We consider the Au/Ag(111), native Au/Ag(100) and reconstructed Au(100)(5×1) facets. For a given combination of metallic surface and facet, some variations in density profile are apparent across the different cations in solution, with the corresponding chloride counterion profiles remaining broadly invariant. All density profiles at the higher concentration are predicted to be very similar to their low-concentration counterparts. We find that each electrolyte responds differently to the different metallic surface and facets, particularly those of the divalent metal ions. Our findings reveal marked differences in density profiles between facets for a given metallic interface for both Mg(2+) and Ca(2+), with Na(+) and K(+) showing much less distinction. Mg(2+) was the only ion for which we find evidence of materials-dependent differences in interfacial solution structuring between the Ag and Au.

Structure of arginine overlayers at the aqueous gold interface: implications for nanoparticle assembly.

Adsorption of small biomolecules onto the surface of nanoparticles offers a novel route to generation of nanoparticle assemblies with predictable architectures. Previously, ligand-exchange experiments on citrate-capped gold nanoparticles with the amino acid arginine were reported to support linear nanoparticle assemblies. Here, we use a combination of atomistic modeling with experimental characterization to explore aspects of the assembly hypothesis for these systems. Using molecular simulation, we probe the structural and energetic characteristics of arginine overlayers on the Au(111) surface under aqueous conditions at both low- and high-coverage regimes. In the low-density regime, the arginines lie flat on the surface. At constant composition, these overlayers are found to be lower in energy than the densely packed films, although the latter case appears kinetically stable when arginine is adsorbed via the zwitterion group, exposing the charged guanidinium group to the solvent. Our findings suggest that zwitterion-zwitterion hydrogen bonding at the gold surface and minimization of the electrostatic repulsion between adjacent guanidinium groups play key roles in determining arginine overlayer stability at the aqueous gold interface. Ligand-exchange experiments of citrate-capped gold nanoparticles with arginine derivatives agmatine and N-methyl-l-arginine reveal that modification at the guanidinium group significantly diminishes the propensity for linear assembly of the nanoparticles.

What makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces

Investigation of the non-covalent interaction of biomolecules with aqueous graphene interfaces is a rapidly expanding area. However, reliable exploitation of these interfaces in many applications requires that the links between the sequence and binding of the adsorbed peptide structures be clearly established. Molecular dynamics (MD) simulations can play a key role in elucidating the conformational ensemble of peptides adsorbed at graphene interfaces, helping to elucidate these rules in partnership with experimental characterisation. We apply our recently-developed polarisable force-field for biomolecule-graphene interfaces, GRAPPA, in partnership with advanced simulation approaches, to probe the adsorption behaviour of peptides at aqueous graphene. First we determine the free energy of adsorption of all twenty naturally occurring amino acids (AAs) via metadynamics simulations, providing a benchmark for interpreting peptide-graphene adsorption studies. From these free energies, we find that strong-binding amino acids have flat and/or compact side chain groups, and we relate this behaviour to the interfacial solvent structuring. Second, we apply replica exchange with solute tempering simulations to efficiently and widely sample the conformational ensemble of two experimentally-characterised peptide sequences, P1 and its alanine mutant P1A3, in solution and adsorbed on graphene. For P1 we find a significant minority of the conformational ensemble possesses a helical structure, both in solution and when adsorbed, while P1A3 features mostly extended, random-coil conformations. In solution this helical P1 configuration is stabilised through favourable intra-peptide interactions, while the adsorbed structure is stabilised via interaction of four strongly-binding residues, identified from our metadynamics simulations, with the aqueous graphene interface. Our findings rationalise the performance of the P1 sequence as a known graphene binder.

Understanding water and ion transport behaviour and permeability through poly(amide) thin film composite membrane

Molecular dynamics (MD) together with the adaptive biasing force (ABF) and metadynamics free energy calculation methods was used to investigate the permeation properties of salt water through poly(amide) thin film composite reverse osmosis membranes. The thin films were generated by annealing an amorphous cell of poly(amide) chains through an MD method. The MD results showed they have typical structural properties of the active layer of thin film composite membranes and comparable water diffusivity (2.13×10^-5cm²/s for the film with a density of 1.06g/cm³) and permeability (9.27×10^-15cm³cm/cm²sPa) to experimental data. The simulations of water permeation through the films under different transmembrane pressures revealed the behaviours of water molecules in the thin films and the dynamic regimes of water permeation, including Brownian diffusion, flush and jump diffusion regimes. The intermolecular interactions of water and ions with poly(amide) chains showed a strong dependence on the local structure of films. The attraction between water and ploy(amide) molecules can be up to 8.5kcal/mol in dense polymer regions and 5kcal/mol in the pores of about 3nm. The ABF and metadynamics simulations produced the profiles of free energy potential of water and ions along the depth of the thin films, which provided important information for quantitatively determining the barrier energy required for water permeation and rejection of ions. The thin film with a density of 1.06g/cm³ and a thickness of 6nm offers a rejection to Na⁺ but a slight absorption of Cl^- (0.25kcal/mol) at 0.3-0.4nm distance to its surface. Water molecules must overcome 63kcal/mol energy to move to the centre of the film. The dependences of the barrier energy and the water-polymer interaction energy on the local free volume size in the thin film were analysed. The simulations of water permeation under high transmembrane pressures showed a nonlinear response of the concentration and distribution of water molecules in the film to the imposed pressure. Compaction of the film segments close to the porous substrate and water congestion in dense regions significantly influenced the water permeation when the membrane was operated under pressures of more than 3.0MPa.

Tristearin bilayers: structure of the aqueous interface and stability in the presence of surfactants

We report results of atomistic molecular dynamics simulations of an industrially-relevant, exemplar triacylglycerol (TAG), namely tristearin (TS), under aqueous conditions, at different temperatures and in the presence of an anionic surfactant, sodium dodecylbenzene sulphonate (SDBS). We predict the TS bilayers to be stable and in a gel phase at temperatures of 350 K and below. At 370 K the lipid bilayer was able to melt, but does not feature a stable liquid-crystalline phase bilayer at this elevated temperature. We also predict the structural characteristics of TS bilayers in the presence of SDBS molecules under aqueous conditions, where surfactant molecules are found to spontaneously insert into the TS bilayers. We model TS bilayers containing different amounts of SDBS, with the presence of SDBS imparting only a moderate effect on the structure of the system. Our study represents the first step in applying atomistic molecular dynamics simulations to the investigation of TAG-aqueous interfaces. Our results suggest that the CHARMM36 force-field appears suitable for the simulation of such systems, although the phase behaviour of the system may be shifted to lower temperatures than is the case for the actual system. Our findings provide a foundation for further simulation studies of the TS-aqueous interface.

Multipoint haptic guidance for micrograsping systems

The ability to perform accurate micromanipulation offers wide-reaching benefits and is of increasing interest to researchers. Recent research into microgripper, microtweezer, and microforcep systems contributes toward accurate micrograsping and manipulation. Despite these efforts, achieving adequate operator control remains a distinct research challenge. Haptic interfaces interact with the human's haptic modality and offer the ability to enhance the operator's controllability of micromanipulation systems. Our previous work introduced single-point haptic guidance to assist the operator during intracellular microinjection. This paper extends the approach to propose multipoint haptic guidance for micrograsping tasks. Accurate micrograsping is valuable in many applications, including microassembly and biomanipulation. A multipoint haptic gripper facilitates haptic interaction, and haptic guidance assists the operator in controlling systems suitable for micrograsping. Force fields are used to guide the operator to suitable grasp points on micrometer-sized objects and consist of attractive and repulsive forces. The ability of the force field to effectively assist the operator in grasping the cell is evaluated using a virtual environment. Evaluation results demonstrate the ability of the approach to significantly reduce participants' average grasping error.

Improving the description of interactions between Ca2+ and protein carboxylate groups, including γ-carboxyglutamic acid: revised CHARMM22∗ parameters

A reliable description of ion pair interactions for biological systems, particularly those involving polyatomic ions such as carboxylate and divalent ions such as Ca²⁺, using biomolecular force-fields is essential for making useful predictions for a range of protein functions. In particular, the interaction of divalent ions with the double carboxylate group present in γ-carboxyglutamic acid (Gla), relevant to the function of many proteins, is relatively understudied using biomolecular force-fields. Using force-field based metadynamics simulations to predict the free energy of binding between Ca²⁺ and the carboxylate group in liquid water, we show that a widely-used biomolecular force-field, CHARMM22∗, substantially over-estimates the binding strength between Ca²⁺ and the side-chains of both glutamic acid (Glu) and Gla, compared with experimental data obtained for the analogous systems of aqueous calcium-acetate and calcium-malonate. To correct for this, we propose and test a range of modifications to the σ value of the heteroatomic Lennard-Jones interaction between Ca²⁺ and the oxygen of the carboxylate group. Our revised parameter set can recover the same three association modes of this aqueous ion pair as the standard parameter set, and yields free energies of binding for the carboxylate-Ca²⁺ interaction in good agreement with experimental data. The revised parameter set recovers other structural properties of the ion pair in agreement with the standard CHARMM22∗ parameter set.

Testing the transferability of a coarse-grained model to intrinsically disordered proteins

The intermediate-resolution coarse-grained protein model PLUM [T. Bereau and M. Deserno, J. Chem. Phys., 2009, 130, 235106] is used to simulate small systems of intrinsically disordered proteins involved in biomineralisation. With minor adjustments to reduce bias toward stable secondary structure, the model generates conformational ensembles conforming to structural predictions from atomistic simulation. Without additional structural information as input, the model distinguishes regions of the chain by predicted degree of disorder, manifestation of structure, and involvement in chain dimerisation. The model is also able to distinguish dimerisation behaviour between one intrinsically disordered peptide and a closely related mutant. We contrast this against the poor ability of PLUM to model the S1 quartz-binding peptide.

Computations of the flow past a still sphere at moderate reynolds numbers using an immersed boundary method

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Theoretical studies of liquids by computer simulations: the water-acetone mixture.

Monte Carlo simulation results for pure liquid acetone and water-acetone mixtures calculated in the isothermal and isobaric (NPT) ensemble at T=298K and p=1.0atm are presented. The TIP4P model was used for water and optimized potential for liquid simulation (OPLS) force field parameters used for acetone. The results obtained for the average configurational energy as a function of the mole fraction are in good accord with experimental data. Energy partitioning and co-ordination numbers results calculated for equimolar water-acetone solution are compared to similar data obtained for other water-organic liquid mixtures. These results show an increase in water-water interaction energy and co-ordination numbers when the interaction between water and organic liquid molecules decrease. Distribution functions for pure liquid acetone and water-acetone mixtures are presented. Dipole-dipole angular correlation functions obtained for pure liquid acetone show a predominance of dimers with parallel alignment of dipole moments. Radial distribution functions from water-acetone interaction show characteristic features of hydrogen bonded liquids. Radial and angular distribution functions for water-water correlation calculated in pure water and in equimolar water-acetone mixture are compared, showing very similar features in both systems. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Morphology and topography analysis of mesoporous titania templated by micrometric latex sphere arrays

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

On the unzipping mechanisms of carbon nanotubes: Insights from reactive molecular dynamics simulations

Unzipping carbon nanotubes (CNTs) is considered one of the most promising approaches for the controlled and large-scale production of graphene nanoribbons (GNR). These structures are considered of great importance for the development of nanoelectronics because of its dimensions and intrinsic nonzero band gap value. Despite many years of investigations some details on the dynamics of the CNT fracture/unzipping processes remain unclear. In this work we have investigated some of these process through molecular dynamics simulations using reactive force fields (ReaxFF), as implemented in the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. We considered multi-walled CNTs of different dimensions and chiralities and under induced mechanical stretching. Our preliminary results show that the unzipping mechanisms are highly dependent on CNT chirality. Well-defined and distinct fracture patterns were observed for the different chiralities. Armchair CNTs favor the creation of GNRs with well-defined armchair edges, while zigzag and chiral ones produce GNRs with less defined and defective edges. © 2012 Materials Research Society.

Modelagem numérica do escoamento em válvulas automáticas de compressores pelo Método da Fronteira Imersa

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)