Developing and validating Fuzzy-Border continuum solvation model with POlarizable Simulations Second order Interaction Model (POSSIM) force field for proteins

The accurate, fast and low cost computational tools are indispensable for studying the structure and dynamics of biological macromolecules in aqueous solution. The goal of this thesis is development and validation of continuum Fuzzy-Border (FB) solvation model to work with the Polarizable Simulations Second-order Interaction Model (POSSIM) force field for proteins developed by Professor G A Kaminski. The implicit FB model has advantages over the popularly used Poisson Boltzmann (PB) solvation model. The FB continuum model attenuates the noise and convergence issues commonly present in numerical treatments of the PB model by employing fixed position cubic grid to compute interactions. It also uses either second or first-order approximation for the solvent polarization which is similar to the second-order explicit polarization applied in POSSIM force field. The FB model was first developed and parameterized with nonpolarizable OPLS-AA force field for small molecules which are not only important in themselves but also building blocks of proteins and peptide side chains. The hydration parameters are fitted to reproduce the experimental or quantum mechanical hydration energies of the molecules with the overall average unsigned error of ca. 0.076kcal/mol. It was further validated by computing the absolute pKa values of 11 substituted phenols with the average unsigned error of 0.41pH units in comparison with the quantum mechanical error of 0.38pH units for this set of molecules. There was a good transferability of hydration parameters and the results were produced only with fitting of the specific atoms to the hydration energy and pKa targets. This clearly demonstrates the numerical and physical basis of the model is good enough and with proper fitting can reproduce the acidity constants for other systems as well. After the successful development of FB model with the fixed charges OPLS-AA force field, it was expanded to permit simulations with Polarizable Simulations Second-order Interaction Model (POSSIM) force field. The hydration parameters of the small molecules representing analogues of protein side chains were fitted to their solvation energies at 298.15K with an average error of ca.0.136kcal/mol. Second, the resulting parameters were used to reproduce the pKa values of the reference systems and the carboxylic (Asp7, Glu10, Glu19, Asp27 and Glu43) and basic residues (Lys13, Lys29, Lys34, His52 and Lys55) of the turkey ovomucoid third domain (OMTKY3) protein. The overall average unsigned error in the pKa values of the acid residues was found to be 0.37pH units and the basic residues was 0.38 pH units compared to 0.58pH units and 0.72 pH units calculated previously using polarizable force field (PFF) and Poisson Boltzmann formalism (PBF) continuum solvation model. These results are produced with fitting of specific atoms of the reference systems and carboxylic and basic residues of the OMTKY3 protein. Since FB model has produced improved pKa shifts of carboxylic residues and basic protein residues in OMTKY3 protein compared to PBF/PFF, it suggests the methodology of first-order FB continuum solvation model works well in such calculations. In this study the importance of explicit treatment of the electrostatic polarization in calculating pKa of both acid and basic protein residues is also emphasized. Moreover, the presented results demonstrate not only the consistently good degree of accuracy of protein pKa calculations with the second-degree POSSIM approximation of the polarizable calculations and the first-order approximation used in the Fuzzy-Border model for the continuum solvation energy, but also a high degree of transferability of both the POSSIM and continuum solvent Fuzzy Border parameters. Therefore, the FB model of solvation combined with the POSSIM force field can be successfully applied to study the protein and protein-ligand systems in water

Predicting Skin Permeability by means of Computational Approaches : Reliability and Caveats in Pharmaceutical Studies

© 2019 American Chemical Society.The skin is the main barrier between the internal body environment and the external one. The characteristics of this barrier and its properties are able to modify and affect drug delivery and chemical toxicity parameters. Therefore, it is not surprising that permeability of many different compounds has been measured through several in vitro and in vivo techniques. Moreover, many different in silico approaches have been used to identify the correlation between the structure of the permeants and their permeability, to reproduce the skin behavior, and to predict the ability of specific chemicals to permeate this barrier. A significant number of issues, like interlaboratory variability, experimental conditions, data set building rationales, and skin site of origin and hydration, still prevent us from obtaining a definitive predictive skin permeability model. This review wants to show the main advances and the principal approaches in computational methods used to predict this property, to enlighten the main issues that have arisen, and to address the challenges to develop in future research.Peer reviewedFinal Accepted Versio

What Controls the Controller: Structure and Function Characterizations of Transcription Factor PU.1 Uncover Its Regulatory Mechanism

The ETS family transcription factor PU.1/Spi-1 is a master regulator of the self-renewal of hematopoietic stem cells and their differentiation along both major lymphoid and myeloid branches. PU.1 activity is determined in a dosage-dependent manner as a function of both its expression and real-time regulation at the DNA level. While control of PU.1 expression is well established, the molecular mechanisms of its real-time regulation remain elusive. Our work is focused on discovering a complete regulatory mechanism that governs the molecular interactions of PU.1. Structurally, PU.1 exhibits a classic transcription factor architecture in which intrinsically disordered regions (IDR), consisting of 66% of its primary structure, are tethered to a well-structured DNA binding domain. The transcriptionally active form of PU.1 is a monomer that binds target DNA sites as a 1:1 complex. Our investigations show that IDRs of PU.1 reciprocally control two separate inactive dimeric forms, with and without DNA. At high concentrations, PU.1 forms a non-canonical 2:1 complex at a single DNA specific site. In the absence of DNA, PU.1 also forms a dimer, but it is incompatible with DNA binding. The DNA-free PU.1 dimer is further promoted by phosphomimetic mutants of IDR residues that are phosphorylated in B-lymphocytic activation. These results lead us to postulate a model of real-time PU.1 regulation, unknown in the ETS family, where independent dimeric forms antagonize each other to control the dosage of active PU.1 monomer at its target DNA sites. To demonstrate the biological relevance of our model, cellular assays probing PU.1-specific reporters and native target genes show that PU.1 transactivation exhibits a distinct dose response consistent with negative feedback. In summary, we have established the first model for the general real-time regulation of PU.1 at the DNA/protein level, without the need for recruiting specific binding partners. These novel interactions present potential therapeutic targets for correcting de-regulated PU.1 dosage in hematologic disorders, including leukemia, lymphoma, and myeloma

The C Terminus of the Ribosomal-Associated Protein LrtA Is an Intrinsically Disordered Oligomer

The 191-residue-long LrtA protein of Synechocystis sp. PCC 6803 is involved in post-stress survival and in stabilizing 70S ribosomal particles. It belongs to the hibernating promoting factor (HPF) family, intervening in protein synthesis. The protein consists of two domains: The N-terminal region (N-LrtA, residues 1-101), which is common to all the members of the HPF, and seems to be well-folded; and the C-terminal region (C-LrtA, residues 102-191), which is hypothesized to be disordered. In this work, we studied the conformational preferences of isolated C-LrtA in solution. The protein was disordered, as shown by computational modelling, 1D-H-1 NMR, steady-state far-UV circular dichroism (CD) and chemical and thermal denaturations followed by fluorescence and far-UV CD. Moreover, at physiological conditions, as indicated by several biochemical and hydrodynamic techniques, isolated C-LrtA intervened in a self-association equilibrium, involving several oligomerization reactions. Thus, C-LrtA was an oligomeric disordered protein.This research was funded by Spanish Ministry of Economy and Competitiveness [CTQ2015-64445-R (to J.L.N.) and MAT2015-63704-P (to A.A.), with Fondo Social Europeo (ESF)], and by the Basque Government [IT-654-13 (to A.A.)

The C Terminus of the Ribosomal-Associated Protein LrtA Is an Intrinsically Disordered Oligomer

The 191-residue-long LrtA protein of Synechocystis sp. PCC 6803 is involved in post-stress survival and in stabilizing 70S ribosomal particles. It belongs to the hibernating promoting factor (HPF) family, intervening in protein synthesis. The protein consists of two domains: The N-terminal region (N-LrtA, residues 1-101), which is common to all the members of the HPF, and seems to be well-folded; and the C-terminal region (C-LrtA, residues 102-191), which is hypothesized to be disordered. In this work, we studied the conformational preferences of isolated C-LrtA in solution. The protein was disordered, as shown by computational modelling, 1D-H-1 NMR, steady-state far-UV circular dichroism (CD) and chemical and thermal denaturations followed by fluorescence and far-UV CD. Moreover, at physiological conditions, as indicated by several biochemical and hydrodynamic techniques, isolated C-LrtA intervened in a self-association equilibrium, involving several oligomerization reactions. Thus, C-LrtA was an oligomeric disordered protein.This research was funded by Spanish Ministry of Economy and Competitiveness [CTQ2015-64445-R (to J.L.N.) and MAT2015-63704-P (to A.A.), with Fondo Social Europeo (ESF)], and by the Basque Government [IT-654-13 (to A.A.)

The Two-State Prehensile Tail of the Antibacterial Toxin Colicin N

Intrinsically disordered regions within proteins are critical elements in many biomolecular interactions and signaling pathways. Antibacterial toxins of the colicin family, which could provide new antibiotic functions against resistant bacteria, contain disordered N-terminal translocation domains (T-domains) that are essential for receptor binding and the penetration of the Escherichia coli outer membrane. Here we investigate the conformational behavior of the T-domain of colicin N (ColN-T) to understand why such domains are widespread in toxins that target Gram-negative bacteria. Like some other intrinsically disordered proteins in the solution state of the protein, ColN-T shows dual recognition, initially interacting with other domains of the same colicin N molecule and later, during cell killing, binding to two different receptors, OmpF and TolA, in the target bacterium. ColN-T is invisible in the high-resolution x-ray model and yet accounts for 90 of the toxin’s 387 amino acid residues. To reveal its solution structure that underlies such a dynamic and complex system, we carried out mutagenic, biochemical, hydrodynamic and structural studies using analytical ultracentrifugation, NMR, and small-angle x-ray scattering on full-length ColN and its fragments. The structure was accurately modeled from small-angle x-ray scattering data by treating ColN as a flexible system, namely by the ensemble optimization method, which enables a distribution of conformations to be included in the final model. The results reveal, to our knowledge, for the first time the dynamic structure of a colicin T-domain. ColN-T is in dynamic equilibrium between a compact form, showing specific self-recognition and resistance to proteolysis, and an extended form, which most likely allows for effective receptor binding