Computational analyses of the catalytic and heparin-binding sites and their interactions with glycosaminoglycans in glycoside hydrolase family 79 endo-d-glucuronidase (heparanase)

Mammalian heparanase is an endo-β-glucuronidase associated with cell invasion in cancer metastasis, angiogenesis and inflammation. Heparanase cleaves heparan sulfate proteoglycans in the extracellular matrix and basement membrane, releasing heparin/heparan sulfate oligosaccharides of appreciable size. This in turn causes the release of growth factors, which accelerate tumor growth and metastasis. Heparanase has two glycosaminoglycan-binding domains; however, no three-dimensional structure information is available for human heparanase that can provide insights into how the two domains interact to degrade heparin fragments. We have constructed a new homology model of heparanase that takes into account the most recent structural and bioinformatics data available. Heparin analogs and glycosaminoglycan mimetics were computationally docked into the active site with energetically stable ring conformations and their interaction energies were compared. The resulting docked structures were used to propose a model for substrates and conformer selectivity based on the dimensions of the active site. The docking of substrates and inhibitors indicates the existence of a large binding site extending at least two saccharide units beyond the cleavage site (toward the nonreducing end) and at least three saccharides toward the reducing end (toward heparin-binding site 2). The docking of substrates suggests that heparanase recognizes the N-sulfated and O-sulfated glucosamines at subsite +1 and glucuronic acid at the cleavage site, whereas in the absence of 6-O-sulfation in glucosamine, glucuronic acid is docked at subsite +2. These findings will help us to focus on the rational design of heparanase-inhibiting molecules for anticancer drug development by targeting the two heparin/heparan sulfate recognition domains.

The structure of glycosaminoglycans and their interactions with proteins

Glycosaminoglycans (GAGs) are important complex carbohydrates that participate in many biological processes through the regulation of their various protein partners. Biochemical, structural biology and molecular modelling approaches have assisted in understanding the molecular basis of such interactions, creating an opportunity to capitalize on the large structural diversity of GAGs in the discovery of new drugs. The complexity of GAG–protein interactions is in part due to the conformational flexibility and underlying sulphation patterns of GAGs, the role of metal ions and the effect of pH on the affinity of binding. Current understanding of the structure of GAGs and their interactions with proteins is here reviewed: the basic structures and functions of GAGs and their proteoglycans, their clinical significance, the three-dimensional features of GAGs, their interactions with proteins and the molecular modelling of heparin binding sites and GAG–protein interactions. This review focuses on some key aspects of GAG structure–function relationships using classical examples that illustrate the specificity of GAG–protein interactions, such as growth factors, anti-thrombin, cytokines and cell adhesion molecules. New approaches to the development of GAG mimetics as possible new glycotherapeutics are also briefly covered.

Computational methods for the prediction of the structure and interactions of coiled-coil peptides

The past several years have seen significant advances in the development of computational methods for the prediction of the structure and interactions of coiled-coil peptides. These methods are generally based on pairwise correlations of amino acids, helical propensity, thermal melts and the energetics of sidechain interactions, as well as statistical patterns based on Hidden Markov Model (HMM) and Support Vector Machine (SVM) techniques. These methods are complemented by a number of public databases that contain sequences, motifs, domains and other details of coiled-coil structures identified by various algorithms. Some of these computational methods have been developed to make predictions of coiled-coil structure on the basis of sequence information; however, structural predictions of the oligomerisation state of these peptides still remains largely an open question due to the dynamic behaviour of these molecules. This review focuses on existing in silico methods for the prediction of coiled-coil peptides of functional importance using sequence and/or three-dimensional structural data.

Heparin/heparan sulphate-based drugs

Glycosaminoglycans (GAGs) are an untapped source of novel chemical entities and, therefore, offer exciting new opportunities for the development of novel drug molecules because of their unique physical and biological properties. Advances in the functional understanding of GAG–protein interactions are enabling the development of GAG mimetics for use as anti-angiogenic, anti-metastatic, anti-inflammatory, anticoagulant and anti-thrombotic agents. Many anti-thrombotic molecules, such as Fondaparinux and Idraparinux, have been successful in clinical trials, and a new generation of heparin mimetic oligosaccharides and small molecules are currently in different stages of clinical development. In particular, the recent increased activity in the development of new mimetics by altering the composition of sulphated GAGs is very encouraging. This article reviews structurally defined heparin-mimetic oligosaccharides and small molecules currently in development or clinical trials...

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Prediction of heparin binding sites in bone morphogenetic proteins (BMPs)

Heparin is a glycosaminoglycan known to bind bone morphogenetic proteins (BMPs) and the growth and differentiation factors (GDFs) and has strong and variable effects on BMP osteogenic activity. In this paper we report our predictions of the likely heparin binding sites for BMP-2 and 14. The N-terminal sequences upstream of TGF-β-type cysteine-knot domains in BMP-2, 7 and 14 contain the basic residues arginine and lysine, which are key components of the heparin/HS-binding sites, with these residues being highly non-conserved. Importantly, evolutionary conserved surfaces on the beta sheets are required for interactions with receptors and antagonists. Furthermore, BMP-2 has electropositive surfaces on two sides compared to BMP-7 and BMP-14. Molecular docking simulations suggest the presence of high and low affinity binding sites in dimeric BMP-2. Histidines were found to play a role in the interactions of BMP-2 with heparin; however, a pKa analysis suggests that histidines are likely not protonated. This is indicative that interactions of BMP-2 with heparin do not require acidic pH. Taken together, non-conserved amino acid residues in the N-terminus and residues protruding from the beta sheet (not overlapping with the receptor binding sites and the dimeric interface) and not C-terminal are found to be important for heparin–BMP interactions.

Characterization of sequence and structural features of the Candida krusei Enolase

The incidence of human infections by the fungal pathogen Candida species has been increasing in recent years. Enolase is an essential protein in fungal metabolism. Sequence data is available for human and a number of medically important fungal species. An understanding of the structural and functional features of fungal enolases may provide the structural basis for their use as a target for the development of new anti-fungal drugs. We have obtained the sequence of the enolase of Candida krusei (C. krusei), as it is a significant medically important fungal pathogen. We have then used multiple sequence alignments with various enolase isoforms in order to identify C. krusei specific amino acid residues. The phylogenetic tree of enolases shows that the C. krusei enolase assembles on the tree with the fungal genes. Importantly, C. krusei lacks four amino acids in the active site compared to human enolase, as revealed by multiple sequence alignments. These differences in the substrate binding site may be exploited for the design of new anti-fungal drugs to selectively block this enzyme. The lack of the important amino acids in the active site also indicates that C. krusei enolase might have evolved as a member of a mechanistically diverse enolase superfamily catalying somewhat different reactions.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.

A comparative structural bioinformatics analysis of the insulin receptor family ectodomain based on phylogenetic information

The insulin receptor (IR), the insulin-like growth factor 1 receptor (IGF1R) and the insulin receptor-related receptor (IRR) are covalently-linked homodimers made up of several structural domains. The molecular mechanism of ligand binding to the ectodomain of these receptors and the resulting activation of their tyrosine kinase domain is still not well understood. We have carried out an amino acid residue conservation analysis in order to reconstruct the phylogeny of the IR Family. We have confirmed the location of ligand binding site 1 of the IGF1R and IR. Importantly, we have also predicted the likely location of the insulin binding site 2 on the surface of the fibronectin type III domains of the IR. An evolutionary conserved surface on the second leucine-rich domain that may interact with the ligand could not be detected. We suggest a possible mechanical trigger of the activation of the IR that involves a slight ‘twist’ rotation of the last two fibronectin type III domains in order to face the likely location of insulin. Finally, a strong selective pressure was found amongst the IRR orthologous sequences, suggesting that this orphan receptor has a yet unknown physiological role which may be conserved from amphibians to mammals.

Molecular modeling of Bt Cry1Ac (DI–DII)–ASAL (Allium sativum lectin)–fusion protein and its interaction with aminopeptidase N (APN) receptor of Manduca sexta

Genetic engineering of Bacillus thuringiensis (Bt) Cry proteins has resulted in the synthesis of various novel toxin proteins with enhanced insecticidal activity and specificity towards different insect pests. In this study, a fusion protein consisting of the DI–DII domains of Cry1Ac and garlic lectin (ASAL) has been designed in silico by replacing the DIII domain of Cry1Ac with ASAL. The binding interface between the DI–DII domains of Cry1Ac and lectin has been identified using protein–protein docking studies. Free energy of binding calculations and interaction profiles between the Cry1Ac and lectin domains confirmed the stability of fusion protein. A total of 18 hydrogen bonds was observed in the DI–DII–lectin fusion protein compared to 11 hydrogen bonds in the Cry1Ac (DI–DII–DIII) protein. Molecular mechanics/Poisson–Boltzmann (generalized-Born) surface area [MM/PB (GB) SA] methods were used for predicting free energy of interactions of the fusion proteins. Protein–protein docking studies based on the number of hydrogen bonds, hydrophobic interactions, aromatic–aromatic, aromatic–sulphur, cation–pi interactions and binding energy of Cry1Ac/fusion proteins with the aminopeptidase N (APN) of Manduca sexta rationalised the higher binding affinity of the fusion protein with the APN receptor compared to that of the Cry1Ac–APN complex, as predicted by ZDOCK, Rosetta and ClusPro analysis. The molecular binding interface between the fusion protein and the APN receptor is well packed, analogously to that of the Cry1Ac–APN complex. These findings offer scope for the design and development of customized fusion molecules for improved pest management in crop plants.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.

The global burden of injury: Incidence, mortality, disability-adjusted life years and time trends from the Global Burden of Disease study 2013

Background The Global Burden of Diseases (GBD), Injuries, and Risk Factors study used the disability-adjusted life year (DALY) to quantify the burden of diseases, injuries, and risk factors. This paper provides an overview of injury estimates from the 2013 update of GBD, with detailed information on incidence, mortality, DALYs and rates of change from 1990 to 2013 for 26 causes of injury, globally, by region and by country. Methods Injury mortality was estimated using the extensive GBD mortality database, corrections for ill-defined cause of death and the cause of death ensemble modelling tool. Morbidity estimation was based on inpatient and outpatient data sets, 26 cause-of-injury and 47 nature-of-injury categories, and seven follow-up studies with patient-reported long-term outcome measures. Results In 2013, 973 million (uncertainty interval (UI) 942 to 993) people sustained injuries that warranted some type of healthcare and 4.8 million (UI 4.5 to 5.1) people died from injuries. Between 1990 and 2013 the global age-standardised injury DALY rate decreased by 31% (UI 26% to 35%). The rate of decline in DALY rates was significant for 22 cause-of-injury categories, including all the major injuries. Conclusions Injuries continue to be an important cause of morbidity and mortality in the developed and developing world. The decline in rates for almost all injuries is so prominent that it warrants a general statement that the world is becoming a safer place to live in. However, the patterns vary widely by cause, age, sex, region and time and there are still large improvements that need to be made.

Mapping the risk of soil-transmitted helminthic infections in the Philippines

Background In order to increase the efficient allocation of soil-transmitted helminth (STH) disease control resources in the Philippines, we aimed to describe for the first time the spatial variation in the prevalence of A. lumbricoides, T. trichiura and hookworm across the country, quantify the association between the physical environment and spatial variation of STH infection and develop predictive risk maps for each infection. Methodology/Principal Findings Data on STH infection from 35,573 individuals across the country were geolocated at the barangay level and included in the analysis. The analysis was stratified geographically in two major regions: 1) Luzon and the Visayas and 2) Mindanao. Bayesian geostatistical models of STH prevalence were developed, including age and sex of individuals and environmental variables (rainfall, land surface temperature and distance to inland water bodies) as predictors, and diagnostic uncertainty was incorporated. The role of environmental variables was different between regions of the Philippines. This analysis revealed that while A. lumbricoides and T. trichiura infections were widespread and highly endemic, hookworm infections were more circumscribed to smaller foci in the Visayas and Mindanao. Conclusions/Significance This analysis revealed significant spatial variation in STH infection prevalence within provinces of the Philippines. This suggests that a spatially targeted approach to STH interventions, including mass drug administration, is warranted. When financially possible, additional STH surveys should be prioritized to high-risk areas identified by our study in Luzon.