The structure-dynamics-function relation in proteins: bridging all-atom molecular dynamics, experiments, and simplified models.
Creators: name: Rigoli, Marta
Publisher: name: Universit� degli studi di Trento
Year: 2022

description: Proteins are one of the most studied biological molecules of the last decades. A great amount of experimental techniques provide to researchers direct or indirect informations on proteins structure and function. In silico simulations can be used as a ?computational microscope? giving the possibility to observe protein dynamic properties at atomistic resolution. In this work, various applications of computational methods to biological systems are presented. In particular, all-atom Molecular Dynamics (MD) simulations were employed to investigate the behaviour of proteins at atomstic resolution. The term ?Molecular Dynamics? is usually referred to computational methods used for the simulation of classical many-body systems. These techniques are applied to microscopic systems and they represent a powerful approach for the study of physical processes, providing a tool for their interpretation. They have been widely used in the past decades to elucidate a large variety of molecular processes in different fields such as solid state physics, material science, chemistry, biochemistry and biophysics. Here, all-atom MD simulations were employed to observe equilibrium properties of several biologically relevant proteins. This allowed us to direct perform a comparison of molecular mechanisms occurring at the atomistic level as obtained from in silico studies with experimental data, which usually describe processes at larger length and time scales. These MD simulations were also meant as a starting point for the construction of simplified models, as they were processed through coarse-graining procedures to extrapolate crucial systems features, such as informative protein sites, on the basis of information theory approaches. Specifically we studied the dynamics of pembrolizumab, a humanized immunoglobulin of type G4 (IgG4) used as a therapeutic antibody. It is employed for the treatment of lung cancer, melanoma, stomach and head cancer and Hodgkin?s lymphoma. This antibody interacts with the programmed cell death protein 1 (PD-1) receptor, blocking the suppression of the immune response during cancer development. The studied systems are three: the apo state of pembrolizumab, the holo state (i.e. pembrolizumab bound to PD-1) and the glycosylated apo configuration. Each configuration was simulated for 2?s, for a total of 6?s. The analysis of the trajectories was carried out by combining standard structural analysis techniques and information theory-based measures of correlation. From MD trajectories we could extract valuable informations on the connectivity that exists among the structural domains that compose the antibody structure. Moreover, it was possible to infer which regions are involved in the structural rearrangement in the case of the antigen binding. We could observe that the presence of the antigen reduces the conformational variability of the molecule giving a greater stability to it. The second studied system is the P53 protein complex. In this case we focused on the tetramerization domain (TD) region that is composed by 2 identical dimers and has the function of bringing together the four monomers of the p53 complex. Starting from the observation that in case of the mutation of residue R337 several pathologies are developed in humans, we constructed computational models to reproduce the dynamics of the mutants and investigate their behaviour in silico. We performed simulations for a total of 16 ?s divided in 8 different cases. In the first part of the study the wild type (WT) protein was compared to the R337C and the R337H mutant in three different protonation states: delta protonated Histidine, epsilon protonated Histidine ad double protonated Histidine. In the second part of the study we highlighted the differences between the WT configuration and three rationally designed mutants: R337D-352D, 337R-D352R, R337D-D352R. In this part of the investigation, the importance of the electrostatic interaction between residues R337 and D352 in the stability of the tetramerization do- main was discussed. Furthermore, we matched the obtained computational results of p53 tetramerization domain with functional experiments in yeasts (performed in collaboration with the CIBIO department) of all the simulated forms. The third simulated protein is the zinc sensing transcriptional repressor (CzrA), an homodimeric protein that binds DNA in Staphylococcus aureus. All-atom MD simulations of two different configurations were performed for a total of 4?s, the first one is the WT apo protein while the second is the WT holo system, where the protein is complexed with two Zn ions. In this case, in addition to standard analysis techniques, we applied the mapping entropy minimization protocol to highlight the most informative protein regions, from the perspective of information theory. Finally, our in silico results were compared to available NMR data of the protein itself.