MODELLING INTERMOLECULAR FORCES IN BIOMOLECULES: FROM PROTEIN-PROTEIN INTERACTIONS TO HALOGEN BONDS.

The accurate description and evaluation of the intermolecular interactions has a great importance in the molecular modelling of biological systems. Protein-protein interactions, in particular, being involved in virtually every cellular process, are nowadays the object of thorough studies aimed at the understanding and modulation of the underlying mechanism. In the Part I of the present work, it will be shown how, through a computational approach, it is possible to get an in-depth analysis of the network of interactions occurring at the interface between tubulin subunits and how vinblastine, a commonly used anticancer agents, is able to interfere with the correct protein association, so having a therapeutic effect. In the Part II, the self-association of the bacterial protein FtsZ will be studied, showing that the detailed description of the protein-protein interactions can provide key information for the de-novo design of inhibiting molecules. Finally, in the Part II, the attention will be focused on the computational study of halogen bonding, which is found to have a great relevance in the recognition process between biological macromolecules and halogenated agents. In particular, it will be shown that a specific approach is mandatory for its correct description in the framework of the classical force-fields

Molecular dynamics simulations of p97 including covalent, allosteric and ATP-competitive inhibitors

Binary (nucleotide-protein dimer and hexamer complexes) and ternary (nucleotide-protein-inhibitor complexes) p97 complexes were subjected to molecular dynamics simulations in an attempt to further our understanding of the p97 protein oligomer domain stability and, more importantly, of the recently reported diverse molecular mechanisms of inhibition including allosteric, ATP-competitive and covalent inhibitors. Analysis of stable states following equilibration phases indicated a higher intrinsic stability of the homohexamer as opposed to the dimer, and of N-D1 domains as opposed to the D2 domain. The molecular dynamics of the proposed allosteric binding model reproduced important molecular interactions identified experimentally with high frequency throughout the trajectory. Observed conformational changes occurring in the D2 nucleotide binding site provided a novel bind-rearrange-react hypothesis of stepwise molecular events involved in the specific covalent inhibitor mode of actio

Актуальність впровадження систем газового обліку в сучасних умовах

Free energy calculation has long been an important goal for molecular dynamics simulation and force field development, but historically it has been challenged by limited performance, accuracy, and creation of topologies for arbitrary small molecules. This has made it difficult to systematically compare different sets of parameters to improve existing force fields, but in the past few years several authors have developed increasingly automated procedures to generate parameters for force fields such as Amber, CHARMM, and OPLS. Here, we present a new framework that enables fully automated generation of GROMACS topologies for any of these force fields and an automated setup for parallel adaptive optimization of high-throughput free energy calculation by adjusting lambda point placement on the fly. As a small example of this automated pipeline, we have calculated solvation free energies of 50 different small molecules using the GAFF, OPLS-AA, and CGenFF force fields and four different water models, and by including the often neglected polarization costs, we show that the common charge models are somewhat underpolarized.QC 20150505</p

Conjuntos excepcionais e alguns problemas de Mahler

Dissertação (mestrado)—Universidade de Brasília, Instituto de Ciências Exatas, Departamento de Matemática, 2017.Seja f uma função inteira e transcendente. Denotamos por Sf o conjunto de todos os α ∈ ´Q tais que f(α) ∈ ´Q (o conjunto excepcional de f). Nessa dissertação, mostraremos quais subconjuntos de ´Q podem ser o conjunto excepcional de alguma função inteira e transcendente. Além disso, trataremos de dois problemas de Mahler relacionados a propriedades de funções inteiras e transcendentes. Mostraremos que existem funções inteiras e transcendentes que levam um subconjunto dos números de Liouville nele mesmo e daremos uma resposta positiva ao Problema B de Mahler: Problema B: Existe uma função inteira e transcendente f(z) = Σn =0 ∞ a nz n com coeficientes racionais tal que f( ´Q ) ⊆ ´Q e f−1( ´Q ) ⊆ ´Q ? .Let f be an entire transcendental function. We denote by Sf the set of all α ∈ ´Q such that f(α) ∈ ´Q (exceptional set of f). Throughout this dissertation, we will show which subsets of ´Q can be the exceptional set of some entire transcendental function. Moreover, we will deal with two of Mahler’s problems related to properties of entire transcendental functions. We will show that there are entire transcendental functions that map a subset of Liouville numbers in itself and we will give a positive answer for Mahler’s Problem B: Problem B: Is there an entire transcendental function f(z) = Σn =0 ∞ a nz n with rational coefficients such that que f( ´Q ) ⊆ ´Q e f−1( ´Q ) ⊆ ´Q ?

Association of HLA-DRB1 amino acid residues with giant cell arteritis: genetic association study, meta-analysis and geo-epidemiological investigation

Introduction: Giant cell arteritis (GCA) is an autoimmune disease commonest in Northern Europe and Scandinavia. Previous studies report various associations with HLA-DRB1*04 and HLA-DRB1*01; HLA-DRB1 alleles show a gradient in population prevalence within Europe. Our aims were (1) to determine which amino acid residues within HLA-DRB1 best explained HLA-DRB1 allele susceptibility and protective effects in GCA, seen in UK data combined in meta-analysis with previously published data, and (2) to determine whether the incidence of GCA in different countries is associated with the population prevalence of the HLA-DRB1 alleles that we identified in our meta-analysis. Methods: GCA patients from the UK GCA Consortium were genotyped by using single-strand oligonucleotide polymerization, allele-specific polymerase chain reaction, and direct sequencing. Meta-analysis was used to compare and combine our results with published data, and public databases were used to identify amino acid residues that may explain observed susceptibility/protective effects. Finally, we determined the relationship of HLA-DRB1*04 population carrier frequency and latitude to GCA incidence reported in different countries. Results: In our UK data (225 cases and 1378 controls), HLA-DRB1*04 carriage was associated with GCA susceptibility (odds ratio (OR) = 2.69, P = 1.5×10 −11 ), but HLA-DRB1*01 was protective (adjusted OR = 0.55, P = 0.0046). In meta-analysis combined with 14 published studies (an additional 691 cases and 4038 controls), protective effects were seen from HLA-DR2, which comprises HLA-DRB1*15 and HLA-DRB1*16 (OR = 0.65, P = 8.2×10 −6 ) and possibly from HLA-DRB1*01 (OR = 0.73, P = 0.037). GCA incidence (n = 17 countries) was associated with population HLA-DRB1*04 allele frequency (P = 0.008; adjusted R 2 = 0.51 on univariable analysis, adjusted R 2 = 0.62 after also including latitude); latitude also made an independent contribution. Conclusions: We confirm that HLA-DRB1*04 is a GCA susceptibility allele. The susceptibility data are best explained by amino acid risk residues V, H, and H at positions 11, 13, and 33, contrary to previous suggestions of amino acids in the second hypervariable region. Worldwide, GCA incidence was independently associated both with population frequency of HLA-DRB1*04 and with latitude itself. We conclude that variation in population HLA-DRB1*04 frequency may partly explain variations in GCA incidence and that HLA-DRB1*04 may warrant investigation as a potential prognostic or predictive biomarker

Vinblastine perturbation of tubulin protofilament structure: a computational insight

Tubulin is a heterodimeric protein whose self assembly leads to the formation of protofilaments and of more complex structures called microtubules, key components of the cytoskeleton which have a fundamental role in the cell division process. Due to its biological function, tubulin is the target of many antitumoral molecules that exert their action on proliferating tumoral cells. Among these drugs, vinblastine has been widely used in therapy for a long time, albeit its mechanism of interaction with tubulin has remained elusive until recently. Vinblastine acts as a microtubule destabilizing agent and induces the formation of curved or ring-shaped tubulin polymers instead of linear protofilaments in vitro. In this paper we compare, using molecular dynamics simulations and free energy calculations, the network of interactions that allow the assembly of model linear protofilaments with those present in curved tubulin polymers complexed with vinblastine. It is shown that vinblastine, wedging between tubulin heterodimers, actually mediates part of the interactions between them and acts by crosslinking the two proteins, leading to the observed curved polymers rather than to their disassembly

Halogen Bonding in ligand-receptor systems in the framework of classical force fields

Halogen bond is an important non-covalent interaction which is receiving a growing attention in the study of protein-ligand complexes. Many drugs are halogenated molecules and it has been recently shown that many halogenated ligands establish halogen bonds with biomolecules. As the halogen bond nature is due to an anisotropy of the electrostatic potential around halogen atoms, it is not possible to use traditional force fields based on a set of atom-centred charges to study halogen bonds in biomolecules. We show that the introduction of pseudo-atoms on halogens permits us to correctly describe the anisotropy of the electrostatic potential and to perform molecular dynamics simulations on complexes of proteins with halogenated ligands that reproduce experimental values. The results are compared with crystallographic data and with hybrid quantum mechanics/molecular mechanics calculations

Solvent effect on halogen bonding : the case of the I&#183;&#183;&#183;O interaction

The solvent effect on the I efO halogen bonding in complexes of iodobenzene derivatives with formaldehyde has been investigated by systematically varying the substituents on the iodobenzene ring. Calculations have been performed at MP2 and DFT levels of theory, using the aug-cc-pVDZ basis set and the pseudopotential for iodine. Within the DFT approach, a series of the most widely used exchange-correlation functionals have been considered, comprising PBE, PBE0, B3LYP, BH&HLYP, M06-2X and M06-HF. Results obtained in diethylether and in water using the conductor-like polarizable continuum model (CPCM) have been compared with in vacuo results. Though halogen bonding distances were found to systematically shorten when moving from vacuo to diethylether and then to water, the associated interaction energies showed a decrease in absolute value, indicating that solvent has a destabilizing effect on this interaction. By comparison with MP2 results, all the considered functionals, B3LYP excepted, have been found adequate to describe halogen bonding. As far as the interaction energies are concerned, the best performance was obtained with the M06-HF functional in vacuo and the PBE functional in solution. The geometrical parameters characterizing halogen bonds were better reproduced by the M06-2X functional

Halogen Bonds with Benzene : an Assessment of DFT Functionals

The performance of an extensive set of density functional theory functionals has been tested against CCSD(T) and MP2 results, extrapolated to the complete basis set (CBS) limit, for the interaction of either DCl or DBr (D = H, HCC, F, and NC) with the aromatic system of benzene. It was found that double hybrid functionals explicitly including dispersion, that is, B2PLYPD and mPW2PLYPD, provide the better agreement with the CCSD(T)/CBS results on both energies and equilibrium geometry, indicating the importance of dispersive contributions in determining this interaction. Among the less expensive functionals, the better performance is provided by the \u3c9B97X and M062X functionals, while the \u3c9B97XD and B97D functionals are shown to work very well for bromine complexes but not so well for chlorine complexes