Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

T-cell epitope prediction and immune complex simulation using molecular dynamics:state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

Molecular dynamics simulations:bring biomolecular structures alive on a computer

The molecular dynamics (MD) simulations play a very important role in science today. They have been used successfully in binding free-energy calculations and rational design of drugs and vaccines. MD simulations can help visualize and understand structures and dynamics at an atomistic level when combined with molecular graphics programs. The molecular and atomistic properties can be displayed on a computer in a time-dependent way, which opens a road toward a better understanding of the relationship of structure, dynamics, and function. In this chapter, the basics of MD are explained, together with a step-by-step description of setup and running an MD simulation.

Molecular dynamics simulations of the adsorption and diffusion behavior of pure and mixed alkanes in silicalite

The adsorption and diffusion of mixed hydrocarbon components in silicalite have been studied using molecular dynamic simulation methods. We have investigated the effect of molecular loadings and temperature on the diffusional behavior of both pure and mixed alkane components. For binary mixtures with components of similar sizes, molecular diffusional behavior in the channels was noticed to be reversed as loading is increased. This behavior was noticeably absent for components of different sizes in the mixture. Methane molecules in the methane/propane mixture have the highest diffusion coefficients across the entire loading range. Binary mixtures containing ethane molecules prove more difficult to separate compared to other binary components. In the ternary mixture, however, ethane molecules diffuse much faster at 400 K in the channel with a tendency to separate out quickly from other components. © 2005 Elsevier Inc. All rights reserved.

Evidence-based nanoscopic and molecular framework for excipient functionality in compressed orally disintegrating tablets

The work investigates the adhesive/cohesive molecular and physical interactions together with nanoscopic features of commonly used orally disintegrating tablet (ODT) excipients microcrystalline cellulose (MCC) and D-mannitol. This helps to elucidate the underlying physico-chemical and mechanical mechanisms responsible for powder densification and optimum product functionality. Atomic force microscopy (AFM) contact mode analysis was performed to measure nano-adhesion forces and surface energies between excipient-drug particles (6-10 different particles per each pair). Moreover, surface topography images (100 nm2-10 μm2) and roughness data were acquired from AFM tapping mode. AFM data were related to ODT macro/microscopic properties obtained from SEM, FTIR, XRD, thermal analysis using DSC and TGA, disintegration testing, Heckel and tabletability profiles. The study results showed a good association between the adhesive molecular and physical forces of paired particles and the resultant densification mechanisms responsible for mechanical strength of tablets. MCC micro roughness was 3 times that of D-mannitol which explains the high hardness of MCC ODTs due to mechanical interlocking. Hydrogen bonding between MCC particles could not be established from both AFM and FTIR solid state investigation. On the contrary, D-mannitol produced fragile ODTs due to fragmentation of surface crystallites during compression attained from its weak crystal structure. Furthermore, AFM analysis has shown the presence of extensive micro fibril structures inhabiting nano pores which further supports the use of MCC as a disintegrant. Overall, excipients (and model drugs) showed mechanistic behaviour on the nano/micro scale that could be related to the functionality of materials on the macro scale. © 2014 Al-khattawi et al.

Bridging large and small scales of water models using hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework

This thesis presents a two-dimensional water model investigation and development of a multiscale method for the modelling of large systems, such as virus in water or peptide immersed in the solvent. We have implemented a two-dimensional ‘Mercedes Benz’ (MB) or BN2D water model using Molecular Dynamics. We have studied its dynamical and structural properties dependence on the model’s parameters. For the first time we derived formulas to calculate thermodynamic properties of the MB model in the microcanonical (NVE) ensemble. We also derived equations of motion in the isothermal–isobaric (NPT) ensemble. We have analysed the rotational degree of freedom of the model in both ensembles. We have developed and implemented a self-consistent multiscale method, which is able to communicate micro- and macro- scales. This multiscale method assumes, that matter consists of the two phases. One phase is related to micro- and the other to macroscale. We simulate the macro scale using Landau Lifshitz-Fluctuating Hydrodynamics, while we describe the microscale using Molecular Dynamics. We have demonstrated that the communication between the disparate scales is possible without introduction of fictitious interface or approximations which reduce the accuracy of the information exchange between the scales. We have investigated control parameters, which were introduced to control the contribution of each phases to the matter behaviour. We have shown, that microscales inherit dynamical properties of the macroscales and vice versa, depending on the concentration of each phase. We have shown, that Radial Distribution Function is not altered and velocity autocorrelation functions are gradually transformed, from Molecular Dynamics to Fluctuating Hydrodynamics description, when phase balance is changed. In this work we test our multiscale method for the liquid argon, BN2D and SPC/E water models. For the SPC/E water model we investigate microscale fluctuations which are computed using advanced mapping technique of the small scales to the large scales, which was developed by Voulgarakisand et. al.

A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single "zoom-in" user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.

The structure and dynamics of saxicolous lichen communities

The saxicolous lichen vegetation on Ordovician slate rock at the mouth of the River Dovey, South Merionethshire, Wales was described in relation to several environmental variables which include aspect, slope angle, light intensity, rock porosity, rock microtopography and rock stability. Each of the measured environmental variables was shown to influence the lichen vegetation. A number of groups of species which were characteristic of certain environments were described. The data from the saxicolous lichen communities were analysed using multivariate analysis. Qualitative and quantitative data were ordinated, the qualitative data being easier to interpret ecologically, and site number (which reflects distance from the sea and altitude), rock porosity and light intensity were shown to be important environmental variables. A classification of the data was also carried out. The results of the ordination and classification were combined together and a model constructed which describes saxicolous lichen vegetation. A method which uses the model as an aid to the design and interpretation of field experiments is described. The model is applied to an experiment which investigates the effect on growth of transplanting four saxicolous lichens to different aspects. Growth was inhibited in Physcia orbicularis and Parmelia conspersa on rock surfaces of northwest aspect compared with growth on rock surfaces of southeast aspect. Growth was inhibited in Parmelia glabratula ssp. fuliginosa on rock surfaces of southeast aspect compared with rock surfaces of northwesr aspect. The growth of Parmelia saxatilis was similar at both southeast and northwesr aspects. Growth inhibition or stimulation in thalli of Physcia orbicularis, Parmelia conspersa and Parmelia glabratula ssp. fuliginosa after transplantation was consistent with the predictions of the model while the results for Parmelia saxatilis were not as expected. There was evidence that the frequency of Parmelia conspersa and Parmelia glabratula at a site is related to an effect of the environment on the growth of the thalli. There was also evidence that the frequency of Physcia orbicularis at a site is related to an effect of the environment on the establishment phase of the thalli and for the competitive exclusion of Parmelia saxatilis thalli from southeast facing rock surfaces. The distribution of lichens in relation to height on nine rock surfaces was investigated. It was suggested that the distribution of the lichens was influenced by microclimatic factors which are related to height on the rock, environmental variables which are associated with the rock substratum (e.g. rock porosity and rock microtopography) and by historical factors. The pattern of one crustose and one foliose lichen on four rock surfaces of different aspect and slope was investigated. On the vertically inclined surface the density of small thalli of Buellia aethalea and Parmelia glabratula ssp fuliginosa was correlated with the microtopography of the surface in transects horizontally across the rock surface but not in transects vertically down the rock surface. there were consitent differences in the scale and intensity of pattern horizontally and vertically and also a decrease in the intensity of pattern vertically as the slope of the rock surface decreased. These results were consistent with the suggestion of a gradient of microclimatic factors up the rock. The differences in the scale and intensity of pattern in different size classes in the population were consistent with the changes in pattern with time which have been shown to occur during succession in sand dune and salt marsh vegetation. The relationship between thallus size and height on a rock surface and between the radial growth rate and location of a thallus on a rock surface were investigated. Thalli of Parmelia glabratula ssp. fuliginosa were larger at the top of the rock surface than at the bottom and the data were consistent with the suggestion that the colonisation of the rock surface began at the top and, in time, spread downwards. The radial growth rate of the thalli could not be related to variation in slope, porosity, microtopography or directly to height on the rock but could be related to the horizontal location of the thalli on the rock. These results were consistent with the suggestion that here is a gradient of microclimatic factors across the rock surface which is also modified by height on the rock surface. The succession of lichen communities was described by relating the vegetation to rock porosity, rock microtopography, species diversity and rock stability. An initial stage dominated by crustose lichens leads to communities dominated by crustose, foliose and fruticose species. In the late stages of the succession on some rock surfaces crustose species again become dominant. The occurrence of the climax state and cyclic vegetation change in lichen communities are discussed. A mthod of estimating the age structure of a lichen population by relating thallus size to growth rate is described. The sources of error in the method are discussed in some detail and several refinements suggested to increase the accuracy of the method. The population dynamics of Parmelia glabratula ssp. fuliginosa was investigated by applying life tables to the age structures of eight different populations. The data were consistent with a period of relatively constant recruitment of thalli into the populations. Mortality in lichen populations was divided into deaths which occur after fragmentation of the thallus and deaths which occur after catastrophic environmental events. THe data suggest that the rate of fragmenting death is dependent on the age of the thallus while the rate of catastrophic death is dependent on the number of thalli established in an age class. A comparison of the numbers of thalli in each age class in the eight populations suggested that population density is controlled firstly, by climate and secondly, by variables related to the local rock surface environment. The rate of fragmenting death is related to the diversity of the community and the influence of diversity together with environmental variables in fluctuating or cyclic changes in population number.