Hydrotalcite intercalated siRNA:computational characterization of the interlayer environment

Using molecular dynamics (MD) simulations, we explore the structural and dynamical properties of siRNA within the intercalated environment of a Mg:Al 2:1 Layered Double Hydroxide (LDH) nanoparticle. An ab initio force field (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies: COMPASS) is used for the MD simulations of the hybrid organic-inorganic systems. The structure, arrangement, mobility, close contacts and hydrogen bonds associated with the intercalated RNA are examined and contrasted with those of the isolated RNA. Computed powder X-ray diffraction patterns are also compared with related LDH-DNA experiments. As a method of probing whether the intercalated environment approximates the crystalline or rather the aqueous state, we explore the stability of the principle parameters (e.g., the major groove width) that differentiate both A- and A'- crystalline forms of siRNA and contrast this with recent findings for the same siRNA simulated in water. We find the crystalline forms remain structurally distinct when intercalated, whereas this is not the case in water. Implications for the stability of hybrid LDH-RNA systems are discussed.

Atomic structure of the two intermediate phase glasses SiSe4 and GeSe4

The microscopic origin of the intermediate phase in two prototypical covalently bonded AxB1-x network glass forming systems, where A=Ge or Si, B=Se, and 0=x=1, was investigated by combining neutron diffraction with first-principles molecular-dynamics methods. Specifically, the structure of glassy GeSe4 and SiSe4 was examined, and the calculated total structure factor and total pair-correlation function for both materials are in good agreement with experiment. The structure of both glasses differs markedly from a simple model comprising undefective AB4 corner-sharing tetrahedra in which all A atoms are linked by B2 dimers. Instead, edge-sharing tetrahedra occur and the twofold coordinated Se atoms form three distinct structural motifs, namely, Se-Se2, Se-SeGe (or Se-SeSi), and Se-Ge2 (or Se-Si2). This identifies several of the conformations that are responsible for the structural variability in GexSe1-x and SixSe1-x glasses, a quantity that is linked to the finite width of the intermediate phase window.

Self-assembly of trehalose molecules on a lysozyme surface: the broken glass hypothesis

To help understand how sugar interactions with proteins stabilise biomolecular structures, we compare the three main hypotheses for the phenomenon with the results of long molecular dynamics simulations on lysozyme in aqueous trehalose solution (0.75 M). We show that the water replacement and water entrapment hypotheses need not be mutually exclusive, because the trehalose molecules assemble in distinctive clusters on the surface of the protein. The flexibility of the protein backbone is reduced under the sugar patches supporting earlier findings that link reduced flexibility of the protein with its higher stability. The results explain the apparent contradiction between different experimental and theoretical results for trehalose effects on proteins.

Toward the atomistic simulation of T cell epitopes automated construction of MHC:peptide structures for free energy calculations

Epitopes mediated by T cells lie at the heart of the adaptive immune response and form the essential nucleus of anti-tumour peptide or epitope-based vaccines. Antigenic T cell epitopes are mediated by major histocompatibility complex (MHC) molecules, which present them to T cell receptors. Calculating the affinity between a given MHC molecule and an antigenic peptide using experimental approaches is both difficult and time consuming, thus various computational methods have been developed for this purpose. A server has been developed to allow a structural approach to the problem by generating specific MHC:peptide complex structures and providing configuration files to run molecular modelling simulations upon them. A system has been produced which allows the automated construction of MHC:peptide structure files and the corresponding configuration files required to execute a molecular dynamics simulation using NAMD. The system has been made available through a web-based front end and stand-alone scripts. Previous attempts at structural prediction of MHC:peptide affinity have been limited due to the paucity of structures and the computational expense in running large scale molecular dynamics simulations. The MHCsim server (http://igrid-ext.cryst.bbk.ac.uk/MHCsim) allows the user to rapidly generate any desired MHC:peptide complex and will facilitate molecular modelling simulation of MHC complexes on an unprecedented scale.

The antimicrobial action of pseudin-2

The amphibian antimicrobial peptide pseudin-2 is a peptide derived from the skin of the South-American frog Pseudis paradoxa (Olson et al., 2001). This peptide possesses tremendous potential as a therapeutic lead since it has been shown to possess both antimicrobial as well insulin-releasing properties (Olson et al., 2001; Abdel-Wahab et al., 2008). This study aimed to develop pseudin-2’s potential by understanding and improving its properties as an antimicrobial agent. The structure-function relationships of pseudin-2 were explored using a combination of in-vitro and in-silico techniques, with an aim to predict how the structure of the peptide may be altered in order to improve its efficacy. A library of pseudin-2 mutants was generated by randomizing codons at positions 10, 14 and 18 of a synthetic gene, using NNK saturation mutagenesis. Analysis of these novel peptides broadly confirmed, in line with literature precedent, that anti-microbial activity increases with increased positive charge. Specifically, 2 positively-charged residues at positions 10 and 14 and a hydrophobic at position 18 are preferred. However, substitution at position 14 with some polar, non-charged residues also created peptides with antimicrobial activity. Interestingly, the pseudin-2 analogue [10-E, 14-Q, 18-L] which is identical to pseudin-2, except that the residues at positions 10 and 14 are switched, showed no anti-microbial activity at all. Molecular dynamics simulations of pseudin-2 showed that the peptide possesses two equilibrium structures in a membrane environment: a linear and a kinked a-helix which both embed into the membrane at an angle. Biophysical characterization using circular dichroism spectroscopy confirmed that the peptide is helical within the membrane environment whilst linear dichroism established that the peptide has no defined orientation within the membrane. Collectively, these data indicate that Pseudin-2 exerts its antimicrobial activity via the carpet model.

Structure of glassy GeO2

The full set of partial structure factors for glassy germania, or GeO2, were accurately measured by using the method of isotopic substitution in neutron diffraction in order to elucidate the nature of the pair correlations for this archetypal strong glass former. The results show that the basic tetrahedral Ge(O-1/2)(4) building blocks share corners with a mean inter-tetrahedral Ge-O-Ge bond angle of 132(2)degrees. The topological and chemical ordering in the resultant network displays two characteristic length scales at distances greater than the nearest neighbour. One of these describes the intermediate range order, and manifests itself by the appearance of a first sharp diffraction peak in the measured diffraction patterns at a scattering vector k(FSDP) approximate to 1.53 angstrom(-1), while the other describes so-called extended range order, and is associated with the principal peak at k(PP) = 2.66( 1) angstrom(-1). We find that there is an interplay between the relative importance of the ordering on these length scales for tetrahedral network forming glasses that is dominated by the extended range ordering with increasing glass fragility. The measured partial structure factors for glassy GeO2 are used to reproduce the total structure factor measured by using high energy x-ray diffraction and the experimental results are also compared to those obtained by using classical and first principles molecular dynamics simulations.

Structure of molten TbCl3 measured by neutron diffraction

The total structure factor of molten TbCl3 at 617ºC was measured by using neutron diffraction. The data are in agreement with results from previous experimental work but the use of a diffractometer having an extended reciprocal-space measurement window leads to improved resolution in real space. Significant discrepancies with the results obtained from recent molecular dynamics simulations carried out using a polarizable ion model, in which the interaction potentials were optimized to enhance agreement with previous diffraction data, are thereby highlighted. It is hence shown that there is considerable scope for the development of this model for TbCl3 and for other trivalent metal halide systems spanning a wide range of ion size ratios.

Visualising and controlling the flow in biomolecular systems at and between multiple scales:from atoms to hydrodynamics at different locations in time and space

A novel framework for modelling biomolecular systems at multiple scales in space and time simultaneously is described. The atomistic molecular dynamics representation is smoothly connected with a statistical continuum hydrodynamics description. The system behaves correctly at the limits of pure molecular dynamics (hydrodynamics) and at the intermediate regimes when the atoms move partly as atomistic particles, and at the same time follow the hydrodynamic flows. The corresponding contributions are controlled by a parameter, which is defined as an arbitrary function of space and time, thus, allowing an effective separation of the atomistic 'core' and continuum 'environment'. To fill the scale gap between the atomistic and the continuum representations our special purpose computer for molecular dynamics, MDGRAPE-4, as well as GPU-based computing were used for developing the framework. These hardware developments also include interactive molecular dynamics simulations that allow intervention of the modelling through force-feedback devices.

Multiscale modelling:approaches and challenges

Multiscale systems that are characterized by a great range of spatial–temporal scales arise widely in many scientific domains. These range from the study of protein conformational dynamics to multiphase processes in, for example, granular media or haemodynamics, and from nuclear reactor physics to astrophysics. Despite the diversity in subject areas and terminology, there are many common challenges in multiscale modelling, including validation and design of tools for programming and executing multiscale simulations. This Theme Issue seeks to establish common frameworks for theoretical modelling, computing and validation, and to help practical applications to benefit from the modelling results. This Theme Issue has been inspired by discussions held during two recent workshops in 2013: ‘Multiscale modelling and simulation’ at the Lorentz Center, Leiden (http://www.lorentzcenter.nl/lc/web/2013/569/info.php3?wsid=569&venue=Snellius), and ‘Multiscale systems: linking quantum chemistry, molecular dynamics and microfluidic hydrodynamics’ at the Royal Society Kavli Centre. The objective of both meetings was to identify common approaches for dealing with multiscale problems across different applications in fluid and soft matter systems. This was achieved by bringing together experts from several diverse communities.

Structural determinants of oligomerization of the aquaporin-4 channel

The aquaporins (AQP) family of integral membrane protein channels mediate cellular water and solute flow. Although qualitative and quantitative differences in channel permeability, selectivity, subcellular localization and trafficking responses have been observed for different members of the AQP family, the signature homotetrameric quaternary structure is conserved. Using a variety of biophysical techniques, we show that mutations to an intracellular loop (loop D) of human AQP4 reduce oligomerization. Non-tetrameric AQP4 mutants are unable to relocalize to the plasma membrane in response to changes in extracellular tonicity, despite equivalent constitutive surface expression levels and water permeability to wild-type AQP4. A network of AQP4 loop D hydrogen bonding interactions, identified using molecular dynamics simulations and based on a comparative mutagenic analysis of AQPs 1, 3 and 4, suggest that loop D interactions may provide a general structural framework for tetrameric assembly within the AQP family.

Two-phase flow analogy as an effective boundary condition for modelling liquids at atomistic resolution

A hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework based on the analogy with two-phase hydrodynamics has been extended to dynamically tracking the feature of interest at all-atom resolution. In the model, the hydrodynamics description is used as an effective boundary condition to close the molecular dynamics solution without resorting to standard periodic boundary conditions. The approach is implemented in a popular Molecular Dynamics package GROMACS and results for two biomolecular systems are reported. A small peptide dialanine and a complete capsid of a virus porcine circovirus 2 in water are considered and shown to reproduce the structural and dynamic properties compared to those obtained in theory, purely atomistic simulations, and experiment.

Bioactive sol-gel glasses at the atomic scale:the complementary use of advanced probe and computer modeling methods

Sol-gel-synthesized bioactive glasses may be formed via a hydrolysis condensation reaction, silica being introduced in the form of tetraethyl orthosilicate (TEOS), and calcium is typically added in the form of calcium nitrate. The synthesis reaction proceeds in an aqueous environment; the resultant gel is dried, before stabilization by heat treatment. These materials, being amorphous, are complex at the level of their atomic-scale structure, but their bulk properties may only be properly understood on the basis of that structural insight. Thus, a full understanding of their structure-property relationship may only be achieved through the application of a coherent suite of leading-edge experimental probes, coupled with the cogent use of advanced computer simulation methods. Using as an exemplar a calcia-silica sol-gel glass of the kind developed by Larry Hench, in the memory of whom this paper is dedicated, we illustrate the successful use of high-energy X-ray and neutron scattering (diffraction) methods, magic-angle spinning solid-state NMR, and molecular dynamics simulation as components to a powerful methodology for the study of amorphous materials.

An allosteric role for receptor activity-modifying proteins in defining GPCR pharmacology

G protein-coupled receptors are allosteric proteins that control transmission of external signals to regulate cellular response. Although agonist binding promotes canonical G protein signalling transmitted through conformational changes, G protein-coupled receptors also interact with other proteins. These include other G protein-coupled receptors, other receptors and channels, regulatory proteins and receptor-modifying proteins, notably receptor activity-modifying proteins (RAMPs). RAMPs have at least 11 G protein-coupled receptor partners, including many class B G protein-coupled receptors. Prototypic is the calcitonin receptor, with altered ligand specificity when co-expressed with RAMPs. To gain molecular insight into the consequences of this protein–protein interaction, we combined molecular modelling with mutagenesis of the calcitonin receptor extracellular domain, assessed in ligand binding and functional assays. Although some calcitonin receptor residues are universally important for peptide interactions (calcitonin, amylin and calcitonin gene-related peptide) in calcitonin receptor alone or with receptor activity-modifying protein, others have RAMP-dependent effects, whereby mutations decreased amylin/calcitonin gene-related peptide potency substantially only when RAMP was present. Remarkably, the key residues were completely conserved between calcitonin receptor and AMY receptors, and between subtypes of AMY receptor that have different ligand preferences. Mutations at the interface between calcitonin receptor and RAMP affected ligand pharmacology in a RAMP-dependent manner, suggesting that RAMP may allosterically influence the calcitonin receptor conformation. Supporting this, molecular dynamics simulations suggested that the calcitonin receptor extracellular N-terminal domain is more flexible in the presence of receptor activity-modifying protein 1. Thus, RAMPs may act in an allosteric manner to generate a spectrum of unique calcitonin receptor conformational states, explaining the pharmacological preferences of calcitonin receptor-RAMP complexes. This provides novel insight into our understanding of G protein-coupled receptor-protein interaction that is likely broadly applicable for this receptor class.