Charge environments around phosphorylation sites in proteins

Background: Phosphorylation is a central feature in many biological processes. Structural analyses have identified the importance of charge-charge interactions, for example mediating phosphorylation-driven allosteric change and protein binding to phosphopeptides. Here, we examine computationally the prevalence of charge stabilisation around phosphorylated sites in the structural database, through comparison with locations that are not phosphorylated in the same structures. Results: A significant fraction of phosphorylated sites appear to be electrostatically stabilised, largely through interaction with sidechains. Some examples of stabilisation across a subunit interface are evident from calculations with biological units. When considering the immediately surrounding environment, in many cases favourable interactions are only apparent after conformational change that accompanies phosphorylation. A simple calculation of potential interactions at longer-range, applied to non-phosphorylated structures, recovers the separation exhibited by phosphorylated structures. In a study of sites in the Phospho.ELM dataset, for which structural annotation is provided by non-phosphorylated proteins, there is little separation of the known phospho-acceptor sites relative to background, even using the wider interaction radius. However, there are differences in the distributions of patch polarity for acceptor and background sites in the Phospho.ELM dataset. Conclusion: In this study, an easy to implement procedure is developed that could contribute to the identification of phospho-acceptor sites associated with charge-charge interactions and conformational change. Since the method gives information about potential anchoring interactions subsequent to phosphorylation, it could be combined with simulations that probe conformational change. Our analysis of the Phospho.ELM dataset also shows evidence for mediation of phosphorylation effects through (i) conformational change associated with making a solvent inaccessible phospho-acceptor site accessible, and (ii) modulation of protein-protein interactions

Computation of protein geometry and its applications: Packing and function prediction

This chapter discusses geometric models of biomolecules and geometric constructs, including the union of ball model, the weigthed Voronoi diagram, the weighted Delaunay triangulation, and the alpha shapes. These geometric constructs enable fast and analytical computaton of shapes of biomoleculres (including features such as voids and pockets) and metric properties (such as area and volume). The algorithms of Delaunay triangulation, computation of voids and pockets, as well volume/area computation are also described. In addition, applications in packing analysis of protein structures and protein function prediction are also discussed.Comment: 32 pages, 9 figure

Computational structure‐based drug design: Predicting target flexibility

The role of molecular modeling in drug design has experienced a significant revamp in the last decade. The increase in computational resources and molecular models, along with software developments, is finally introducing a competitive advantage in early phases of drug discovery. Medium and small companies with strong focus on computational chemistry are being created, some of them having introduced important leads in drug design pipelines. An important source for this success is the extraordinary development of faster and more efficient techniques for describing flexibility in three‐dimensional structural molecular modeling. At different levels, from docking techniques to atomistic molecular dynamics, conformational sampling between receptor and drug results in improved predictions, such as screening enrichment, discovery of transient cavities, etc. In this review article we perform an extensive analysis of these modeling techniques, dividing them into high and low throughput, and emphasizing in their application to drug design studies. We finalize the review with a section describing our Monte Carlo method, PELE, recently highlighted as an outstanding advance in an international blind competition and industrial benchmarks.We acknowledge the BSC-CRG-IRB Joint Research Program in Computational Biology. This work was supported by a grant from the Spanish Government CTQ2016-79138-R.J.I. acknowledges support from SVP-2014-068797, awarded by the Spanish Government.Peer ReviewedPostprint (author's final draft

Computational approaches to shed light on molecular mechanisms in biological processes

Computational approaches based on Molecular Dynamics simulations, Quantum Mechanical methods and 3D Quantitative Structure-Activity Relationships were employed by computational chemistry groups at the University of Milano-Bicocca to study biological processes at the molecular level. The paper reports the methodologies adopted and the results obtained on Aryl hydrocarbon Receptor and homologous PAS proteins mechanisms, the properties of prion protein peptides, the reaction pathway of hydrogenase and peroxidase enzymes and the defibrillogenic activity of tetracyclines. © Springer-Verlag 2007

Molecular modeling to study dendrimers for biomedical applications

© 2014 by the authors; licensee MDPI; Basel; Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/). Date of Acceptance: 17/11/2014Molecular modeling techniques provide a powerful tool to study the properties of molecules and their interactions at the molecular level. The use of computational techniques to predict interaction patterns and molecular properties can inform the design of drug delivery systems and therapeutic agents. Dendrimers are hyperbranched macromolecular structures that comprise repetitive building blocks and have defined architecture and functionality. Their unique structural features can be exploited to design novel carriers for both therapeutic and diagnostic agents. Many studies have been performed to iteratively optimise the properties of dendrimers in solution as well as their interaction with drugs, nucleic acids, proteins and lipid membranes. Key features including dendrimer size and surface have been revealed that can be modified to increase their performance as drug carriers. Computational studies have supported experimental work by providing valuable insights about dendrimer structure and possible molecular interactions at the molecular level. The progress in computational simulation techniques and models provides a basis to improve our ability to better predict and understand the biological activities and interactions of dendrimers. This review will focus on the use of molecular modeling tools for the study and design of dendrimers, with particular emphasis on the efforts that have been made to improve the efficacy of this class of molecules in biomedical applications.Peer reviewedFinal Published versio

Mechanistic behaviour and molecular interactions of heat shock protein 47 (HSP47)

This project involves the study of heat shock protein 47 (HSP47), which is a molecular chaperone crucial for collagen biosynthesis. It exhibits a high degree of sequence homology with members of the serine protease inhibitor (serpin) superfamily, though HSP47 does not possess the inhibitory activity. It is a single-substrate chaperone, and binds only to collagen. ‘Knock-out’ of the hsp47 gene impairs the secretion of correctly folded collagen triple helix molecules leading to embryonic lethality in mice. Thus the aim of this project was to elucidate the specific mechanism that governs the binding to and release from collagen at the molecular level, known as the ‘pH-switch mechanism’. Emphasis is given on histidine (His) residues as the HSP47-collagen dissociation pH is similar to the pKa of the imidazole side chain of His residues. Site directed mutagenesis was used to mutate surface His residues, based on a mouse HSP47 homology model. The effects of the mutations on the behaviour of HSP47 were then assessed by collagen binding assays and structural analyses with circular dichroism (CD). All mutants were found to have good solubility and retain their binding ability to collagen like wild-type HSP47 in batch assay, but perturbed behaviour was seen in column experiment. Mutation of His residue at position 191 (H191) causes the shift in the collagen dissociation pH, while mutation of H197 and/or 198 disrupt the specific HSP47-collagen interaction. H191, 197 and 198 are predicted to be located in the region near the C-terminus of strand 3 of β-sheet A (s3A) in the homology model, a region specifically known as the ‘breach cluster’ in serpin nomenclature. The extent of conformational rearrangement of this region was further investigated by means of intrinsic tryptophan fluorescence spectroscopy using a series of single tryptophan (Trp) mutants. Results from analyses performed on the mutants did not contradict the observation seen in His mutational work, as Trp residues in the ‘breach’ cluster are likely to be located in the dynamic region of HSP47 pH-triggered conformational change. In conclusion, this study establishes the importance of His residues in the ‘breach cluster’ to HSP47 pH-switch behaviour. Finally, a model for HSP47 pH-switch mechanism was proposed from data obtained via mutagenesis experiments. The model is hoped to assist future research into HSP47 cellular behaviour and will also be of great use in therapeutic applications involving the molecular chaperone

Fluoroketone inhibition of Ca(2+)-independent phospholipase A2 through binding pocket association defined by hydrogen/deuterium exchange and molecular dynamics.

The mechanism of inhibition of group VIA Ca(2+)-independent phospholipase A(2) (iPLA(2)) by fluoroketone (FK) ligands is examined by a combination of deuterium exchange mass spectrometry (DXMS) and molecular dynamics (MD). Models for iPLA(2) were built by homology with the known structure of patatin and equilibrated by extensive MD simulations. Empty pockets were identified during the simulations and studied for their ability to accommodate FK inhibitors. Ligand docking techniques showed that the potent inhibitor 1,1,1,3-tetrafluoro-7-phenylheptan-2-one (PHFK) forms favorable interactions inside an active-site pocket, where it blocks the entrance of phospholipid substrates. The polar fluoroketone headgroup is stabilized by hydrogen bonds with residues Gly486, Gly487, and Ser519. The nonpolar aliphatic chain and aromatic group are stabilized by hydrophobic contacts with Met544, Val548, Phe549, Leu560, and Ala640. The binding mode is supported by DXMS experiments showing an important decrease of deuteration in the contact regions in the presence of the inhibitor. The discovery of the precise binding mode of FK ligands to the iPLA(2) should greatly improve our ability to design new inhibitors with higher potency and selectivity

Identification of the N-terminal Peptide Binding Site of Glucose-regulated Protein 94

Because the stress protein GRP94 can augment presentation of peptides to T cells, it is important to define how it, as well as all other HSP90 family members, binds peptides. Having previously shown that the N-terminal half of GRP94 can account for the peptide binding activity of the full-length protein, we now locate this binding site by testing predictions of a molecular docking model. The best predicted site was on the opposite face of the β sheet from the pan-HSP90 radicicol-binding pocket, in close proximity to a deep hydrophobic pocket. The peptide and radicicol-binding sites are distinct, as shown by the ability of a radicicol-refractive mutant to bind peptide. When the fluorophore acrylodan is attached to Cys(117)within the hydrophobic pocket, its fluorescence is reduced upon peptide binding, consistent with proximity of the two ligands. Substitution of His(125), which contacts the bound peptide, compromises peptide-binding activity. We conclude that peptide binds to the concave face of the β sheet of the N-terminal domain, where binding is regulated during the action cycle of the chaperone