A phosphorylation-induced turn defines the Alzheimer's disease AT8 antibody epitope on the tau protein

Post mortem biochemical staging of Alzheimer’s disease is currently based on immunochemical analysis of brain slices with the AT8 antibody. The epitope of AT8 is described around the pSer202/pThr205 region of the hyperphosphorylated form of the neuronal protein tau. In this study, NMR spectroscopy was used to precisely map the AT8 epitope on phosphorylated tau, and derive its defining structural features by a combination of NMR analyses and molecular dynamics. A particular turn conformation is stabilized by a hydrogen bond of the phosphorylated Thr205 residue to the amide proton of Gly207, and is further stabilized by the two Arg residues opposing the pSer202/pThr205.

Free energy calculations of glycosaminoglycan - Protein interactions

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5-6 of PECAM-1 and a region of high affinity in domains 2-3, consistent with experimental data and ligand-protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules.

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.

Identification of IDUA and WNT16 phosphorylation-related non-synonymous polymorphisms for bone mineral density in meta-analyses of genome-wide association studies

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single-nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome-wide association (GWA) study, we conducted a three-stage meta-analysis targeting phosphorylation-related SNPs (phosSNPs) for femoral neck (FN)-bone mineral density (BMD), total hip (HIP)-BMD, and lumbar spine (LS)-BMD phenotypes. In stage 1, 9593 phosSNPs were meta-analyzed in 11,140 individuals of various ancestries. Genome-wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10–6 (0.05/9593) and 1.00 × 10–4, respectively. In stage 2, nine stage 1–discovered phosSNPs (based on α = 1.00 × 10–4) were in silico meta-analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10–3, 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN-BMD (p = 8.36 × 10–10, p = 5.26 × 10–10, and p = 3.01 × 10–10, respectively) and HIP-BMD (p = 3.26 × 10–6, p = 1.97 × 10–6, and p = 1.63 × 10–12, respectively). Although in vitro studies demonstrated no differences in expressions of wild-type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD-associated non-synonymous variants. © 2015 American Society for Bone and Mineral Research.

Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies

MicroRNAs (miRNAs) are critical post-transcriptional regulators. Based on a previous genome-wide association (GWA) scan, we conducted a polymorphism in microRNAs' Target Sites (poly-miRTS)-centric multistage meta-analysis for lumbar spine (LS)-, total hip (HIP)-, and femoral neck (FN)-bone mineral density (BMD). In stage I, 41,102 poly-miRTSs were meta-analyzed in 7 cohorts with a genome-wide significance (GWS) α=0.05/41,102=1.22×10-6. By applying α=5×10-5 (suggestive significance), 11 poly-miRTSs were selected, with FGFRL1 rs4647940 and PRR5 rs3213550 as top signals for FN-BMD (P-value=7.67×10-6 and 1.58×10-5) in gender-combined sample. In stage II in silico replication (two cohorts), FGFRL1 rs4647940 was the only signal marginally replicated for FN-BMD (P-value=5.08×10-3) at α=0.10/11=9.09×10-3. PRR5 rs3213550 was also selected based on biological significance. In stage III de novo genotyping replication (two cohorts), FGFRL1 rs4647940 was the only signal significantly replicated for FN-BMD (P-value=7.55×10-6) at α=0.05/2=0.025 in gender-combined sample. Aggregating three stages, FGFRL1 rs4647940 was the single stage I-discovered and stages II- and III-replicated signal attaining GWS for FN-BMD (P-value=8.87×10-12). Dual-luciferase reporter assays demonstrated that FGFRL1 3' untranslated region harboring rs4647940 appears to be hsa-miR-140-5p's target site. In a zebrafish microinjection experiment, dre-miR-140-5p is shown to exert a dramatic impact on craniofacial skeleton formation. Taken together, we provided functional evidence for a novel FGFRL1 poly-miRTS rs4647940 in a previously known 4p16.3 locus, and experimental and clinical genetics studies have shown both FGFRL1 and hsa-miR-140-5p are important for bone formation. © The Author 2015. Published by Oxford University Press. All rights reserved.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10-8) level: 14q24.2 (rs227425, P-value 3.98 × 10-13, SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10-9, CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n 5 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis. © The Author 2013. Published by Oxford University Press.

Evaluating viewpoint entropy for ribbon representation of protein structure

While many measures of viewpoint goodness have been proposed in computer graphics, none have been evaluated for ribbon representations of protein secondary structure. To fill this gap, we conducted a user study on Amazon’s Mechanical Turk platform, collecting human viewpoint preferences from 65 participants for 4 representative su- perfamilies of protein domains. In particular, we evaluated viewpoint entropy, which was previously shown to be a good predictor for human viewpoint preference of other, mostly non-abstract objects. In a second study, we asked 7 molecular biology experts to find the best viewpoint of the same protein domains and compared their choices with viewpoint entropy. Our results show that viewpoint entropy overall is a significant predictor of human viewpoint preference for ribbon representations of protein secondary structure. However, the accuracy is highly dependent on the complexity of the structure: while most participants agree on good viewpoints for small, non-globular structures with few secondary structure elements, viewpoint preference varies considerably for complex structures. Finally, experts tend to choose viewpoints of both low and high viewpoint entropy to emphasize different aspects of the respective structure.

Nano-confined synthesis of fullerene mesoporous carbon (C60-FMC) with bimodal pores: XRD, TEM, structural properties, NMR, and protein immobilization

Nanoconfined synthesized crystalline fullerene mesoporous carbon (C60-FMC) with bimodal pore architectures of 4.95 nm and 10-15 nm pore sizes characterized by XRD, TEM, nitrogen adsorption/ desorption isotherm and solid-state NMR, and the material was used for protein immobilization. The solid-state 13C NMR spectrum of C60-FMC along with XRD, BET and TEM confirms the formation of fullerene mesoporous carbon structure C60-FMC. The immobilization of albumin (from bovine serum, BSA) protein biomolecule in a buffer solution at pH 4.7 was used to determine the adsorption properties of the C60-FMC material and its structural changes investigated by FT-IR. We demonstrated that the C60-FMC with high surface area and pore volumes have excellent adsorption capacity towards BSA protein molecule. Protein adsorption experiments clearly showed that the C60-FMC with bimodal pore architectures (4.95 nm and 10-15 nm) are suitable material to be used for protein adsorption

Neck length and mean arterial pressure in the sauropod dinosaurs

How blood was able to reach the heads of the long-necked sauropod dinosaurs has long been a matter of debate and several hypotheses have been presented. For example, it has been proposed that sauropods had exceptionally large hearts, multiple ‘normal’ sized hearts spaced at regular intervals up the neck or held their necks horizontal, or that the siphon effect was in operation. By means of an experimental model, we demonstrate that the siphon principle is able to explain how blood was able to adequately perfuse the sauropod brain. The return venous circulation may have been protected from complete collapse by a structure akin to the vertebral venous plexus. We derive an equation relating neck height and mean arterial pressure, which indicates that with a mean arterial pressure similar to that of the giraffe, the maximum safe vertical distance between heart and head would have been about 12 m. A hypothesis is presented that the maximum neck length in the fossil record is due to the siphon height limit. The equation indicates that to migrate over high ground, sauropods would have had to either significantly increase their mean arterial pressure or keep their necks below a certain height dependent on altitude.

Structural analysis of the roles of influenza A virus membrane-associated proteins in assembly and morphology

The assembly of influenza A virus at the plasma membrane of infected cells leads to release of enveloped virions that are typically round in tissue culture-adapted strains but filamentous in strains isolated from patients. The viral proteins hemagglutinin (HA), neuraminidase (NA), matrix protein 1 (M1), and M2 ion channel all contribute to virus assembly. When expressed individually or in combination in cells, they can all, under certain conditions, mediate release of membrane-enveloped particles, but their relative roles in virus assembly, release, and morphology remain unclear. To investigate these roles, we produced membrane-enveloped particles by plasmid-derived expression of combinations of HA, NA, and M proteins (M1 and M2) or by infection with influenza A virus. We monitored particle release, particle morphology, and plasma membrane morphology by using biochemical methods, electron microscopy, electron tomography, and cryo-electron tomography. Our data suggest that HA, NA, or HANA (HA plus NA) expression leads to particle release through nonspecific induction of membrane curvature. In contrast, coexpression with the M proteins clusters the glycoproteins into filamentous membrane protrusions, which can be released as particles by formation of a constricted neck at the base. HA and NA are preferentially distributed to differently curved membranes within these particles. Both the budding intermediates and the released particles are morphologically similar to those produced during infection with influenza A virus. Together, our data provide new insights into influenza virus assembly and show that the M segment together with either of the glycoproteins is the minimal requirement to assemble and release membrane-enveloped particles that are truly virus-like.