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.

Multiscale molecular dynamics/hydrodynamics implementation of two dimensional “Mercedes Benz” water model

A multiscale Molecular Dynamics/Hydrodynamics implementation of the 2D Mercedes Benz (MB or BN2D) [1] water model is developed and investigated. The concept and the governing equations of multiscale coupling together with the results of the two-way coupling implementation are reported. The sensitivity of the multiscale model for obtaining macroscopic and microscopic parameters of the system, such as macroscopic density and velocity fluctuations, radial distribution and velocity autocorrelation functions of MB particles, is evaluated. Critical issues for extending the current model to large systems are discussed.

Cellular and molecular consequences of S100A4-induced motility in rat breast tumour Rama 37 cells

Since the first discovery of S100 members in 1965, their expressions have been affiliated with numerous biological functions in all cells of the body. However, in the recent years, S100A4, a member of this superfamily has emerged as the central target in generating new avenue for cancer therapy as its overexpression has been correlated with cancer patients’ mortality as well as established roles as motility and metastasis promoter. As it has no catalytic activity, S100A4 has to interact with its target proteins to regulate such effects. Up to date, more than 10 S100A4 target proteins have been identified but the mechanical process regulated by S100A4 to induce motility remains vague. In this work, we demonstrated that S100A4 overexpression resulted in actin filaments disorganisation, reduction in focal adhesions, instability of filopodia as well as exhibiting polarised morphology. However, such effects were not observed in truncated versions of S100A4 possibly highlighting the importance of C terminus of S100A4 target recognition. In order to assess some of the intracellular mechanisms that may be involved in promoting migrations, different strategies were used, including active pharmaceutical agents, inhibitors and knockdown experiments. Treatment of S100A4 overexpressing cells with blebbistatin and Y-27632, non muscle myosin IIA (NMMIIA) inhibitors, as well as knockdown of NMMIIA, resulted in motility enhancement and focal adhesions reduction proposing that NMMIIA assisted S100A4 in regulating cell motility but its presence is not essential. Further work done using Cos 7 cell lines, naturally lacking NMMIIA, further demonstrated that S100A4 is capable of regulating cell motility independent of NMMIIA, possibly through poor maturation of focal adhesion. Given that all these experiments highlighted the independency of NMMIIA towards migration, a protein that has been put at the forefront of S100A4-induced motility, we aimed to gather further understanding regarding the other molecular mechanisms that may be at play for motility. Using high throughput imaging (HCI), 3 compounds were identified to be capable of inhibiting S100A4-mediated migration. Although we have yet to investigate the underlying mechanism for their effects, these compounds have been shown to target membrane proteins and the externalisation of S100 proteins, for at least one of the compounds, leading us to speculate that preventing externalisation of S100A4 could potentially regulate cell motility.

MHC class II binding prediction by molecular docking

Proteins of the Major Histocompatibility Complex (MHC) bind self and nonself peptide antigens or epitopes within the cell and present them at the cell surface for recognition by T cells. All T-cell epitopes are MHC binders but not all MCH binders are T-cell epitopes. The MHC class II proteins are extremely polymorphic. Polymorphic residues cluster in the peptide-binding region and largely determine the MHC's peptide selectivity. The peptide binding site on MHC class II proteins consist of five binding pockets. Using molecular docking, we have modelled the interactions between peptide and MHC class II proteins from locus DRB1. A combinatorial peptide library was generated by mutation of residues at peptide positions which correspond to binding pockets (so called anchor positions). The binding affinities were assessed using different scoring functions. The normalized scoring functions for each amino acid at each anchor position were used to construct quantitative matrices (QM) for MHC class II binding prediction. Models were validated by external test sets comprising 4540 known binders. Eighty percent of the known binders are identified in the best predicted 15% of all overlapping peptides, originating from one protein. © 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.