Thermodynamics and stability of the PAAD/DAPIN/PYRIN domain of IFI-16

The PAAD domain is a conserved domain recently identified in more than 35 human proteins that are involved in apoptosis and inflammatory signaling pathways. Structural studies have confirmed that this domain belongs to the death domain superfamily which includes PAAD/CARD/DED/DD families. Recently, the 3D structures determined by NMR of NALP1 and ASC PAAD domain, members of the PAAD family, have shown that it is composed of a 6 helix bundle as with other death domain family members. However, helix-3 in the solved structures is unordered in solution. In this study we compare the thermodynamic, folding and stability properties of different members of the PAAD and CARD families and investigate structural conformational changes induced by the helix inducers trifluoroethanol and SDS on the PAAD domain of IFI16 and on the CARD domain of RAIDD. We show that inside the PAAD and CARD families, members have similar thermodynamic properties, however, the DeltaG of folding for PAAD and CARD members are, respectively, -1.4 and -5.5 kcal mol(-1). This difference is attributed to less alpha helical content for PAAD due to the unfolding of helix-3 that lowers bonded energy and increases disorder when compared to CARD members. Despite identical fold between PAAD and CARD families but limited sequence identity, there are striking differences in the thermodynamics of both families.     

Homology modeling of the Abl-SH3 domain

A tertiary structure model of the Abl-SH3 domain is predicted by using homology modeling techniques coupled to molecular dynamics simulations. Two template proteins were used, Fyn-SH3 and Spc-SH3. The refined model was extensively checked for errors using criteria based on stereochemistry, packing, solvation free-energy, accessible surface areas, and contact analyses. The different checking methods do not totally agree, as each one evaluates a different characteristic of protein structures. Several zones of the protein are more susceptible to incorporating errors. These include residues 13, 15, 35, 39, 45, 46, 50, and 60. An interesting finding is that the measurement of the Cα chirality correlated well with the rest of the criteria, suggesting that this parameter might be a good indicator of correct local conformation. Deviations of more than 4 degrees may be indicative of poor local structure. © 1994 Wiley-Liss, Inc.     

Homology modeling of an immunoglobulin-like domain in the Saccharomyces cerevisiae adhesion protein alpha-agglutinin.

The Saccharomyces cerevisiae adhesion protein alpha-agglutinin is expressed by cells of alpha mating type. On the basis of sequence similarities, alpha-agglutinin has been proposed to contain variable-type immunoglobulin-like (IgV) domains. The low level of sequence similarity to IgV domains of known structure made homology modeling using standard sequence-based alignment algorithms impossible. We have therefore developed a secondary structure-based method that allowed homology modeling of alpha-aggulutinin domain III, the domain most similar to IgV domains. The model was assessed and where necessary refined to accommodate information obtained by biochemical and molecular genetic approaches, including the positions of a disulfide bond, glycosylation sites, and proteolytic sites. The model successfully predicted surface exposure of glycosylation and proteolytic sites, as well as identifying residues essential for binding activity. One side of the domain was predicted to be covered by carbohydrate residues. Surface accessibility and volume packing analyses showed that the regions of the model that have greatest sequence dissimilarity from the IgV consensus sequence are poorly structured in the biophysical sense. Nonetheless, the utility of the model suggests that these alignment and testing techniques should be of general use for building and testing of models of proteins that share limited sequence similarity with known structures.     

Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.     

The PYRIN domain: a member of the death domain-fold superfamily

PYRIN domains were identified recently as putative protein-protein interaction domains at the N-termini of several proteins thought to function in apoptotic and inflammatory signaling pathways. The approximately 95 residue PYRIN domains have no statistically significant sequence homology to proteins with known three-dimensional structure. Using secondary structure prediction and potential-based fold recognition methods, however, the PYRIN domain is predicted to be a member of the six-helix bundle death domain-fold superfamily that includes death domains (DDs), death effector domains (DEDs), and caspase recruitment domains (CARDs). Members of the death domain-fold superfamily are well established mediators of protein-protein interactions found in many proteins involved in apoptosis and inflammation, indicating further that the PYRIN domains serve a similar function. An homology model of the PYRIN domain of CARD7/DEFCAP/NAC/NALP1, a member of the Apaf-1/Ced-4 family of proteins, was constructed using the three-dimensional structures of the FADD and p75 neurotrophin receptor DDs, and of the Apaf-1 and caspase-9 CARDs, as templates. Validation of the model using a variety of computational techniques indicates that the fold prediction is consistent with the sequence. Comparison of a circular dichroism spectrum of the PYRIN domain of CARD7/DEFCAP/NAC/NALP1 with spectra of several proteins known to adopt the death domain-fold provides experimental support for the structure prediction.     

Computational exploration of the binding mode of heme‐dependent stimulators into the active catalytic domain of soluble guanylate cyclase

Soluble guanylate cyclase (sGC), the main target of nitric oxide (NO), has been proven to have a significant role in coronary artery disease, pulmonary hypertension, erectile dysfunction, and myocardial infarction. One of its agonists, BAY 41‐2272 (Riociguat), has been recently approved for treatment of pulmonary arterial hypertension (PHA), while some others are in clinical phases of development. However, the location of the binding sites for the two known types of agonists, heme‐dependent stimulators and heme‐independent activators, is a matter of debate, particularly for the first group where both a location on the regulatory (H‐NOX) and on the catalytic domain have been suggested by different authors. Here, we address its potential location on the catalytic domain, the unique well characterized at the structural level, by an “in silico” approach. Homology models of the catalytic domain of sGC in “inactive” or “active” conformations were constructed using the structure of previously described crystals of the catalytic domains of “inactive” sGCs (2WZ1, 3ET6) and of “active” adenylate cyclase (1CJU). Each model was submitted to six independent molecular dynamics simulations of about 1 μs. Docking of YC‐1, a classic heme‐dependent stimulator, to all frames of representative trajectories of “inactive” and “active” conformations, followed by calculation of absolute binding free energies with the linear interaction energy (LIE) method, revealed a potential high‐affinity binding site on the “active” structure. The site, located between the pseudo‐symmetric and the catalytic site just over the loop β₂–β₃, does not overlap with the forskolin binding site on adenylate cyclases. Proteins 2016; 84:1534–1548. © 2016 Wiley Periodicals, Inc.     

Molecular modeling study of the editing active site of Escherichia coli leucyl-tRNA synthetase: two amino acid binding sites in the editing domain

Aminoacyl-tRNA synthetases (aaRSs) strictly discriminate their cognate amino acids. Some aaRSs accomplish this via proofreading and editing mechanisms. Mursinna and coworkers recently reported that substituting a highly conserved threonine (T252) with an alanine within the editing domain of Escherichia coli leucyl-tRNA synthetase (LeuRS) caused LeuRS to cleave its cognate aminoacylated leucine from tRNA(Leu) (Mursinna et al., Biochemistry 2001;40:5376-5381). To achieve atomic level insight into the role of T252 in LeuRS and the editing reaction of aaRSs, a series of molecular modeling studies including homology modeling and automated docking simulations were carried out. A 3D structure of E. coli LeuRS was constructed via homology modeling using the X-ray structure of Thermus thermophilus LeuRS as a template because the E. coli LeuRS structure is not available from X-ray or NMR studies. However, both the X-ray T. thermophilus and homology-modeled E. coli structures were used in our studies. Amino acid binding sites in the proposed editing domain, which is also called the connective polypeptide 1 (CP1) domain, were investigated by automated docking studies. The root mean square deviation (RMSD) for backbone atoms between the X-ray and homology-modeled structures was 1.18 A overall and 0.60 A for the editing (CP1) domain. Automated docking studies of a leucine ligand into the editing domain were performed for both structures: homology structure of E. coli LeuRS and X-ray structure of T. thermophilus LeuRS for comparison. The results of the docking studies suggested that there are two possible amino acid binding sites in the CP1 domain for both proteins. The first site lies near a threonine-rich region that includes the highly conserved T252 residue, which is important for amino acid discrimination. The second site is located in a flexible loop region surrounded by residues E292, A293, M295, A296, and M298. The important T252 residue is at the bottom of the first binding pocket.     

Computational insights into the binding pattern of mitochondrial calcium uniporter inhibitor through homology modeling,molecular dynamics simulation,binding free energy prediction and density functional theory calculation

Abstract

The mitochondrial calcium uniporter (MCU) is the critical protein of the inner mitochondrial membrane that is the primary mediator for calcium uptake into the mitochondrial matrix. Herein we built the optimal homology model of human MCU which was refined through all-atom molecular dynamics simulation. Then, the binding mode of known inhibitor was predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the docking result and stability of the protein-inhibitor complex. Finally, density functional theory (DFT) calculation enhanced our understanding of the molecular interaction of MCU inhibitor. Our research would provide a deeper insight into the interactions between human MCU and its inhibitor, which boosts to develop novel therapy against MCU related disease.

Communicated by Ramaswamy H. Sarma     

Hsp70 family member, mot-2/mthsp70/GRP75, binds to the cytoplasmic sequestration domain of the p53 protein

Hsp70 family member mot-2/mthsp70/GRP75/PBP74 was shown to bind to the tumor suppressor protein p53. In this study, by in vivo coimmunoprecipitation of mot-2 with p53 and its deletion mutants, the mot-2 binding site of p53 was mapped to its C-terminal amino acid residues 312-352, a region of p53 that includes its cytoplasmic sequestration domain. These data demonstrate that cytoplasmic sequestration and inactivation of p53 by mot-2 occurs by its binding to the cytoplasmic sequestration domain. Therefore, perturbation of mot-p53 interactions can be employed to abrogate cytoplasmic retention of wild-type p53 in tumors.     

Molecular modeling and dynamic simulations of agglutinin-like family members from Candida albicans: New insights into potential targets for the treatment of candidiasis

Infections by Candida albicans in immune compromised patients cause significant morbidity and mortality. In the search for potential molecular targets for drug development, the family of agglutinin-like proteins (Als) in C. albicans have been identified due to numerous attributes associated with high virulence, most prominently due to their role in adherence. Here, molecular models of individual members of the Als family illustrated common and unique structure features. Additionally, dynamic simulations were performed to display regions of high mobility. The results showed variations between Als members in the fluctuation of the A1B1 protein loop, which is located at the entrance to the peptide binding cavity, suggesting that this feature may be a factor contributing to observed differences in affinities to ligands and adhesion properties. Molecular docking results further suggested that ligand affinity could be influenced by movements in the A1B1 loop. In addition, a new site was identified in Als in an area adjacent to the peptide binding cavity that could serve as a new binding site for the design of future anti-adhesion ligands that provide increased specificity inhibiting Als proteins from C. albicans.     

The structure and ligand binding properties of the B. subtilis YkoF gene product, a member of a novel family of thiamin/HMP-binding proteins

The crystal structure of the Bacillus subtilis YkoF gene product, a protein involved in the hydroxymethyl pyrimidine (HMP) salvage pathway, was solved by the multiwavelength anomalous dispersion (MAD) method and refined with data extending to 1.65 A resolution. The atomic model of the protein shows a homodimeric association of two polypeptide chains, each containing an internal repeat of a ferredoxin-like betaalphabetabetaalphabeta fold, as seen in the ACT and RAM-domains. Each repeat shows a remarkable similarity to two members of the COG0011 domain family, the MTH1187 and YBL001c proteins, the crystal structures of which were recently solved by the Northeast Structural Genomics Consortium. Two YkoF monomers form a tightly associated dimer, in which the amino acid residues forming the interface are conserved among family members. A putative small-ligand binding site was located within each repeat in a position analogous to the serine-binding site of the ACT-domain of the Escherichia coli phosphoglycerate dehydrogenase. Genetic data suggested that this could be a thiamin or HMP-binding site. Calorimetric data confirmed that YkoF binds two thiamin molecules with varying affinities and a thiamine-YkoF complex was obtained by co-crystallization. The atomic model of the complex was refined using data to 2.3 A resolution and revealed a unique H-bonding pattern that constitutes the molecular basis of specificity for the HMP moiety of thiamin.     

Structure and dynamics analysis of a new member heparinase II/III of family 12 polysaccharide lyase from Pseudopedobacter saltans by computational modeling and small-angle X-ray scattering

Abstract

The structure of heparinase II/III belonging to family 12 polysaccharide lyase (PsPL12a) from Pseudopedobacter saltans was generated by homology modeling. Multiple sequence alignment showed conserved (Asn216, Tyr270 and His400) and semi-conserved active site amino acid residues. The modeled structure of PsPL12a displayed α/α toroid domain at N-terminal and antiparallel β sheets at C-terminal domain. The modeled structure was similar to those of heparinases from polysaccharide lyase 12 and 21 families. Validation of PsPL12a model by Ramachandran plot showed 94.6% of residues in the favored region, 5.2% of residues in the allowed region and only 0.2% of residues in the outlier region. The area and volume computed for PsPL12a displayed nearly a closed conformation of the active site, similar to HepIII from Bacteroides thetaiotaomicron. The charge calculation on the surface of the PsPL12a structure showed the higher distribution of positive charge in the active site cleft as compared with other homologous structures. Molecular docking study of MD-simulated PsPL12a structure with heparin oligosaccharide showed high binding affinity as compared with heparan sulfate oligosaccharides. Comparison of the active site of modeled PsPL12a with other homologous heparinases revealed putative catalytic triad involving the residues Asn216, His400 and Tyr270. Small-angle X-ray scattering analysis of PsPL12a displayed a fully folded and boxing glove-like envelop.

Communicated by Ramaswamy H. Sarma     

Ligand-supported homology modeling of the human angiotensin II type 1 (AT(1)) receptor: insights into the molecular determinants of telmisartan binding

Angiotensin II type 1 (AT(1)) receptor belongs to the super-family of G-protein-coupled receptors, and antagonists of the AT(1) receptor are effectively used in the treatment of hypertension. To understand the molecular interactions of these antagonists, such as losartan and telmisartan, with the AT(1) receptor, a homology model of the human AT(1) (hAT(1)) receptor with all connecting loops was constructed from the 2.6 A resolution crystal structure (PDB i.d., 1L9H) of bovine rhodopsin. The initial model generated by MODELLER was subjected to a stepwise ligand-supported model refinement. This protocol involved initial docking of non-peptide AT(1) antagonists in the putative binding site, followed by several rounds of iterative energy minimizations and molecular dynamics simulations. The final model was validated based on its correlation with several structure-activity relationships and site-directed mutagenesis data. The final model was also found to be in agreement with a previously reported AT(1) antagonist pharmacophore model. Docking studies were performed for a series of non-peptide AT(1) receptor antagonists in the active site of the final hAT(1) receptor model. The docking was able to identify key molecular interactions for all the AT(1) antagonists studied. Reasonable correlation was observed between the interaction energy values and the corresponding binding affinities of these ligands, providing further validation for the model. In addition, an extensive unrestrained molecular dynamics simulation showed that the docking-derived bound pose of telmisartan is energetically stable. Knowledge gained from the present studies can be used in structure-based drug design for developing novel ligands for the AT(1) receptor.     

Mass spectrometry insights into a tandem ubiquitin‐binding domain hybrid engineered for the selective recognition of unanchored polyubiquitin

Unanchored polyubiquitin chains are emerging as important regulators of cellular physiology with diverse roles paralleling those of substrate‐conjugated polyubiquitin. However tools able to discriminate unanchored polyubiquitin chains of different isopeptide linkages have not been reported. We describe the design of a linker‐optimized ubiquitin‐binding domain hybrid (t‐UBD) containing two UBDs, a ZnF‐UBP domain in tandem with a linkage‐selective UBA domain, which exploits avidity effects to afford selective recognition of unanchored Lys48‐linked polyubiquitin chains. Utilizing native MS to quantitatively probe binding affinities we confirm cooperative binding of the UBDs within the synthetic protein, and desired binding specificity for Lys48‐linked ubiquitin dimers. Furthermore, MS/MS analyses indicate that the t‐UBD, when applied as an affinity enrichment reagent, can be used to favor the purification of endogenous unanchored Lys48‐linked polyubiquitin chains from mammalian cell extracts. Our study indicates that strategies for the rational design and engineering of polyubiquitin chain‐selective binding in nonbiological polymers are possible, paving the way for the generation of reagents to probe unanchored polyubiquitin chains of different linkages and more broadly the ‘ubiquitome’. All MS data have been deposited in the ProteomeXchange with identifier PXD004059 ( http://proteomecentral.proteomexchange.org/dataset/PXD004059 ).     

1.15 A crystal structure of the X. tropicalis Spred1 EVH1 domain suggests a fourth distinct peptide-binding mechanism within the EVH1 family

The recently described Spred protein family has been implicated in the modulation of receptor tyrosine kinase signalling. We report the crystal structure of the Enabled/vasodilator-stimulated phosphoprotein homology-1 (EVH1) domain from Xenopus tropicalis Spred1, solved to 1.15 A resolution. This structure confirms that the Spred EVH1 adopts the pleckstrin-homology fold, with a similar secondary structure to Enabled. A translation of one of the peptide-binding groove beta-strands narrows this groove, whilst one end of the groove shows structural flexibility. We propose that Spred1 will bind peptides that are less proline-rich than other EVH1 domains, with conformational changes indicating an induced fit.     

Structural insights into the binding of SARS-CoV-2, SARS-CoV,and hCoV-NL63 spike receptor-binding domain to horse ACE2

1. Download : Download high-res image (171KB)
2. Download : Download full-size image

     

Identification of a novel human MAST4 gene, a new member of human microtubule associated serine/threonine kinase family

Human protein kinases make up a large superfamily of homologous proteins, which are related by virtue of their kinase domains (also known as catalytic domains). Here, we report the cloning and characterization of a novel human MAST4 (microtubule associated serine/threonine kinase family member 4) gene, which locates on human chromosome 5q13. The MAST4 cDNA is 7587 base pairs in length and encodes a putative protein of 2435 amino acids which contains a serine/threonine kinase domain and a PDZ domain. MAST4 protein has 64, 63, 59, and 39% identical amino acid residues with MAST1, MAST2, MAST3, and MASTL, respectively. RT-PCR analysis revealed a relatively high expression level of MAST4 in most normal human tissues, with the exception of in testis, small intestine, colon, and peripheral blood leukocyte. Published in Russian in Molekulyarnaya Biologiya, 2006, Vol. 40, No. 5, pp. 808–815. The text was submitted by the authors in English. The nucleotide sequences reported in this paper have been submitted to GenBank under accession number: AY830839. These two authors contributed equally to this paper.