Medical, Pharma, Engineering, Science, Technology and Business

**Jiapu Zhang ^{*}**

Faulty of Science, School of Science, Informatics Technology and Engineering, Centre of Informatics and Applied Optimisation, The Federation University of Australia, MT Helen Campus, Mt Helen, Ballarat, Victoria 3353, Australia

- *Corresponding Author:
- Jiapu Zhang

Faulty of Science, School of Science, Informatics Technology and Engineering

Centre of Informatics and Applied Optimisation, The Federation University of Australia

MT Helen Campus, Mt Helen, Ballarat, Victoria 3353, Australia

**Tel:**61423487360, 61353276335

**E-mail:**[email protected], [email protected]

**Received Date:** May 22, 2014; **Accepted Date:** June 02, 2014; **Published Date:** June 05, 2014

**Citation:** Zhang J (2014) The Sensor Network in Molecular Structures of PrP(113-120) AGAAAAGA Amyloid Fibrils. J Comput Sci Syst Biol 7: 137-143. doi: 10.4172/jcsb.1000148

**Copyright:** © 2014 Zhang J. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Computer Science & Systems Biology

The problem of locating sensors in telecommunication networks is a distance geometry problem (DGP). In such a case, the positions of some sensors are known (which are called anchors) and some of the distances between sensors (which can or cannot be anchors) are known. The DGP is to locate the positions of all the sensors. Molecular DGP (MDGP) looks sensors as atoms and their telecommunication network as a molecule for the determination of its three-dimensional (3D) structure. This Chapter defines some sensor networks for determining molecular structures of PrP(113-120) AGAAAAGA amyloid fibrils, which are unstable, noncrystalline, insoluble and hard to be determined in NMR or X-ray experimental laboratories. The amyloid fibril structure is the common structure associated with some 20 neurodegenerative amyloid diseases (including Parkinson’s, Alzheimer’s, Huntington’s, and Prions’), and other diseases such as Type II diabetes, etc. The sensor networks established in this Chapter will benefit the study of 3D molecular structures of all these diseases and will be useful in the research areas such as structural materials, computer-aided or structure-based drug design, and the computational theory of molecular dynamics, and quantum mechanics/molecular mechanics.

Sensor networks; Distance Geometry Problem (DGP); Molecular Distance Geometry Problem (MDGP); Amyloid fibril MDGP models; Neurodegenerative diseases

Neurodegenerative diseases including Parkinson’s, Alzheimer’s, Huntington’s, and Prions’ are amyloid fibril diseases [1], where such as mad cow disease in cattle, scrapie in sheep and kuru and Creutzfeldt- Jakob disease (CJD) in humans are prion diseases. Creutzfeldt-Jakob disease (CJD), Gerstmann-Straussler-Scheinker (GSS) syndrome, Fatal familial insomnia (FFI) and Kuru are the main prion diseases affecting humans. The atomic structures of all these amyloid fibrils revealed steric zippers, with strong van der Waals (vdWs) interactions between β-sheets and hydrogen bonds (HBs) to maintain the β-strands [2]. As we all know, prion proteins have two regions: the N-terminal unstructured region and the C-terminal structured region. The prion AGAAAAGA palindrome hydrophobic region PrP^{C}(113-120) just falls within the N-terminal unstructured region PrP^{C}(1-123). According to Brown [3], AGAAAAGA is important for fibril formation and is an inhibitor of PrPSc neurotoxicity. However, to date the structural data on AGAAAAGA is very little, and traditional X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy experimental methods cannot be easily used to get the structural data of AGAAAAGA. Thus, the advanced computational methods introducing novel mathematical formulations and physical concepts into MDGPs of the molecular medicine field have some advantages under such a circumstance.

The atomic structures of all amyloid fibrils revealed steric zippers, with strong vdW interactions between β-sheets and HBs to maintain the β-strands [2]. The “amyloid fibril” problem can be looked as a MDGP [4], which arises in the interpretation of NMR (nuclear magnetic resonance) data and in the determination of protein structure [as an example to understand MDGP, the problem of locating sensors in telecommunication networks is a DGP. In such a case, the positions of some sensors are known (which are called anchors) and some of the distances between sensors (which can be anchors or not) are known: the DGP is to locate the positions of all the sensors. Here we look sensors as atoms and their telecommunication network as a molecule]. The three dimensional structure of a molecule with n atoms can be described by specifying the 3-Dimensional coordinate positions x_{1}, x_{2}, …, x_{n} in R^{3} of all its atoms. Given bond lengths d_{ij} between a subset S of the atom pairs, the determination of the molecular structure is

(P) to find subject to in S,

where denotes a norm in a real vector space and it is calculated as the Euclidean distance 2-norm in this article. (P) can be reformulated as a mathematical global optimization problem

in the terms of finding the global minimum of the function P(X), where in S are positive weights, in R^{3n} [5] and usually S has many fewer than n2/2 elements due to the error in the theoretical or experimental data [4,6]. There may even not exist any solution to satisfy the distance constraints in the orginal problem (P), for example when data for atoms i, j, k in S violate the triangle inequality; in this case, we may add a perturbation term - ε^{T}X to P(X) :

where

In some cases, instead exact values in S can be found, we can only specify lower and upper bounds on the distances: ; in such cases we may penalize all the unsatisfied constraints into the objective function of (Pε) by adding

into Pε (X) [4,6], where we may let d_{ij} be the interatomic distance (less than 6 angstroms) for the pair in successive residues of a protein and set

[4]. If we look at the prion AGAAAAGA molecular modeling problem as a MDGP with two anchors and two sensors, we can easily construct the prion AGAAAAGA amyloid fibril models. In fact we may let the coordinates of these two anchors being variable. But, these two anchors belong to one body of a chain of a molecule, and the two sensors belong to another body of another chain of the molecule. This is a simple two-body problem model of theoretical physics, i.e. Einstein’s absolute relative theory. Hence, we may look the coordinates of two anchors being fixed.

This article will use do some homology modellings for the PrPSc unstructured region (113-120), AGAAAAGA. The atomic-resolution amyloid fibril structure of this peptide is a steric zipper, with strong vdW contacts and HBs to maintain the molecular structure. The minimization problem of vdW and HB contacts can be looked as a MDGP. In the section of Methods, we present the templates from the Protein Data Bank (http://www.rcsb.org, PDB), which will be used to construct the amyloid fibril models. How to build the models will also be presented in the Methods’ section. Next, in the section of Results, we will show the optimized amyloid fibril models. At last, we give some comments in the section of Conclusions.

The constructions will be based on the most recently released experimental molecular structures in the PDB bank:

• 3NHC.pdb,

• 3MD4.pdb, 3MD5.pdb,

• 3NVE.pdb, 3NVF.pdb, 3NVG.pdb, 3NVH.pdb,

where

• 3NHC.pdb is a crystal structure of human prion segment (127- 132) GYMLGS having Met at residue 129 [7]. This 6-residue segment of PrP centered on residue 129 is a “steric zipper”, pairs of interacting β-sheets. This structure of the “homozygous steric zipper” reveals direct intermolecular interactions between Met in one sheet and the identical residue in the mating sheet. This structure, plus a structure-based model of the heterozygous Met steric zipper, suggests an explanation for the previously observed effects of this locus on prion disease susceptibility and progression;

• 3MD4.pdb and 3MD5.pdb are the crystal cross-beta spine structures of prion 127-132 peptides; and the following four X-ray crystal atomic structures of the same segment 137-143 from human, mouse, and hamster PrP, which is critical for forming amyloid and confers species specificity in PrP seeding experiments:

• 3NVE - MMHFGN segment 138-143 from Syrian Hamster prion [8];

• 3NVF - IIHFGS hexapeptide (residues 138-143) from human prion protein [8];

• 3NVG - MIHFGN hexapeptide (residues 137-142) from mouse prion protein [8];

• 3NVH - MIHFGND (residues 137-143) from mouse prion protein [8].

The models for prion 113-120 AGAAAAGA amyloid fibrils were built by the mutations in the use of SPDBV.4.01.PC (which make all the vdW contacts between the two β-sheets are very far), any Optimization Solver (which will remove the bad vdW/HB contacts) and Optimization program of Amber 11 (which furthermore refines the models and removes all the bad contacts and relax the models into a perfect way) [9]. The models after the mutations of SPDBV.4.01.PC are illuminated in the following **Figure 1** (3NHC.pdb), **Figure 2** (3MD4.pdb), **Figure 3** (3MD5.pdb), **Figure 4** (3NVE.pdb), **Figure 5** (3NVF.pdb), **Figure 6** (3NVG.pdb), and **Figure 7** (3NVH.pdb), respectively.

The optimization problems (Pε) for each Figures are listed as follows. For **Figure 1**, fixing the coordinates of A.ALA3.CB and B.ALA4. CB (two anchors) ((6.014, 5.917, 0.065), (5.658, 1.630, -0.797)) and letting the coordinates of G.ALA4.CB and H.ALA3.CB (two sensors) be variables, there is a Pε-GOP:

For **Figure 2**, fixing the coordinate of B.ALA130.CB (the anchor) ((11.796, 7.063, 8.213)) and letting the coordinates of A.ALA130.CB (a sensor) or A.ALA132.CB (a sensor) be variables, we may get a simple MDGP with 3 variables:

For **Figure 3**, fixing the coordinates of B.ALA130.CB and B.ALA128. CB (two anchors) ((11.796, 7.063, 8.213), (5.849, 6.703, 6.531)) and letting the coordinates of A.ALA130.CB (a sensor) and A.ALA128. CB (another sensor) be variables, we may get a simple MDGP with 6 variables and its dual with 3 variable:

For **Figure 4**, fixing the coordinates of A.ALA4.CB and B.ALA4. CB (two anchors) ((5.150, -1.295, 5.372), (-0.586, -0.444, 12.766)) for 3NVE-AGAAAAG model, fixing the coordinates of A.ALA4. CB, A.ALA2.CB (two anchors) ((5.150, -1.295, 5.372), (4.624, -0.327, 12.132)) for 3NVE-AGAAAAG, 3NVE-GAAAAGA models, letting the coordinates of G.ALA3.CB of 3NVE-models (one sensor) be variables, we may get a simple MDGP with 3 variables and its dual with 2 variables for 3nve-AGAAAAG model:

And for 3NVE-GAAAAG, AAAAGA Models, similarly we may get a simple MDGP with 3 variables and its dual with 2 variables:

For **Figure 5**, We know that for 3NVF-AGAAAA Model 1 at least the VDw interaction between A.GLY2.CA-H.GLY2.CA, A.ALA4.CB-H. GLY2.CA should be maintained, for 3NVF-GAAAAG Model 2 at least three vdW interactions between A.ALA4.CBH-ALA2.CB, A.ALA2. CB-H.ALA2.CB, A.ALA2.CB-H.ALA4.CB should be maintained, and for 3nvf-Model 3 at least three vdw interactions between A.ALA2. CB-H.ALA2.CB, A.ALA2.CB-H.ALA4.CB, A.ALA4.CB-H.ALA2. CB should be maintained. Fixing the coordinates of A.GLY2.CA and A.ALA4.CB (two anchors) ((-10.919, -3.862, -1.487), (6.357, 1.461, -1.905)) for 3NVF-AGAAAA Model 1, fixing the coordinates of A.ALA2.CB and A.ALA4.CB (two anchors) ((11.959, -2.844, -1.977), (6.357, 1.461, -1.905)) for 3NVF-GAAAAG, AAAAGA Models 2-3, letting the coordinates of H.GLY2.CA of 3NVF-AGAAAA Model 1 (two sensors) and the coordinates of H.ALA2.CB and H.ALA4.CB of 3NVF-GAAAAG, AAAAGA Models 2-3 (two sensors) be variables, we may get a simple MDGP with 3/6 variables and its dual with 2/3 variables for 3NVF-AGAAAA Model 1:

For 3NVF-GAAAAG, AAAAGA Models 2-3, similarly we may get a simple MDGP with 6 variables and its dual with 3 variables:

For **Figure 6**, We know that for 3NVG-AGAAAA Model 1 at least the three vdW interaction between A.GLY2.CA-H.GLY2.CA, A.GLY2. CA-H.ALA4.CB, A.ALA4.CB-H.GLY2.CA should be maintained, for 3NVG-GAAAAG Model 2 at least the three vdw interactions between A.ALA2.CB-H.ALA2.CB, A.ALA2.CB-H.ALA4.CB, A.ALA4.CB-H. ALA2.CB should be maintained, and for 3NVG-AAAAGA Model 3 at least the three vdw interactions between A.ALA2.CB-H.ALA2. CB, A.ALA2.CB-H.ALA4.CB, A.ALA4.CB-H.ALA2.CB should be maintained. Fixing the coordinates of A.GLY2.CA and A.ALA4.CB (two anchors) ((-11.159, -2.241, 4.126), (-5.865, -2.618, 8.696)) for 3NVGAGAAAA Model 1, fixing the coordinates of A.ALA2.CB and A.ALA4. CB (two anchors) ((-12.040, -2.675, 5.307), (-5.865, -2.618, 8.696)) for 3NVG-AGAAAA, GAAAAG Models 1-2, letting the coordinates of H.GLY2.CA and H.ALA4.CB of 3NVG-AGAAAA Model 1 (two sensors) and the coordinates of H.ALA2.CB and H.ALA4.CB of 3NVGGAAAAG, AAAAGA Models 2-3 (two sensors) be variables, we may get a simple MDGP with 6 variables and its dual with 3 variables for 3NVG-AGAAAA Model 1:

For 3nvg-Models 2~3, similarly we may get a simple MDGP with 6 variables and its dual with 3 variables:

For **Figure 7**, we may know that for 3NVH-AGAAAAG, GAAAAGA Models 1-2 at least one vdW interaction between A.ALA4.CB-H.ALA4. CB should be maintained. Fixing the coordinates of A.ALA4.CB (the anchor) ((1.731, -1.514, -7.980)) and letting the coordinate of H.ALA4. CB (one sensor) be variables, we may get a simple MDGP with 3 variables and its dual with 1 variable:

Thus, we got the above ten MDGPs.

Solving the ten MDGPs in the above section by any Optimization Solver (which will remove the bad vdW/HB contacts) [9-16] and then refining all the models by the Optimization program of Amber 11 [9,17]. At last we got the optimized prion 113-120 GAAAAGA amyloid fibril models, which are illuminated in **Figures 8-14**.

Sensor Network Problems (SNPs) can be looked as DGPs. MDGPs are SNPs, which look sensors as atoms and their telecommunication network as a molecule for the determination of its three-dimensional (3D) structure. Recently, the 2013 Workshop on Distance Geometry and Applications (2013DGA, http://dga2013.icomp.ufam.edu.br ) united the DGPs and SNPs. Inspired by the reports of 2013DGA, this article presented the amyloid fibril structures of all the neurodegenerative diseases as the SNPs, and then optimized and refined the models of prion 113-120 AGAAAAGA amyloid fibril models. These models can be acted as a reference for the laboratory researches on this region and helpful for the treatment of neurodegenerative prion diseases.

In conclusion, this paper presented molecular structure of amyloid fibril of all the neurodegenerative diseases as a Sensor Network Problem and further presented the solutions of its corresponding Distance Geometry Problem using optimization method. As the presented models to act as reference for future laboratory research, in the result section, the graphs presented not necessarily captured all the details of the optimized model result and is hard to use as reference directly without running the optimization again; thus we will deposit the optimized models as Support Materials into a standard online repository to share it (also be able to find from website https://sites. google.com/site/jiapuzhang/), so that the models would serve the reference purpose better.

The author appreciates the anonymous referees for the numerous insightful comments, which have improved this paper greatly to readers.

20 prion AGAAAGA amyloid fibril models (20 pdb files : optimized_models.zip)

3MD4_AAAAGA_model

3MD4_AGAAAA_model

3MD4_GAAAAG_model

3MD5_AAAAGA_model

3MD5_AGAAAA_model

3MD5_GAAAAG_model

3NHC_AAAAGA_model

3NHC_AGAAAA_model

3NHC_GAAAAG_model

3NVE_AAAAGA_model

3NVE_AGAAAA_model

3NVE_GAAAAG_model

3NVF_AAAAGA_model

3NVF_AGAAAA_model

3NVF_GAAAAG_model

3NVG_AAAAGA_model

3NVG_AGAAAA_model

3NVG_GAAAAG_model

3NVH_AGAAAAG_model

3NVH_GAAAAGA_model

- Graether SP, Slupsky CM, Sykes BD (2003) Freezing of a fish antifreeze protein results in amyloid fibril formation. Biophys J 84: 552-557.
- Sawaya MR, Sambashivan S, Nelson R, Ivanova MI, Sievers SA, et al. (2007) Atomic structures of amyloid cross-beta spines reveal varied steric zippers. Nature 447: 453-457.
- Brown DR (2000) Prion protein peptides: optimal toxicity and peptide blockade of toxicity. Mol Cell Neurosci 15: 66-78.
- Grosso A, Locatelli M, Schoen F (2009) Solving molecular distance geometry problems by global optimization algorithms. ComputOptimAppl 43: 23-37.
- More JJ, Wu ZJ (1997) Global continuation for distance geometry problems. SIAM J Optim 7: 814-836.
- Zou ZH, Bird RH, Schnabel RB (1997) A stochastic/perturbation global optimization algorithm for distance geometry problems. J Glob Optim 11: 91-105.
- Apostol MI, Sawaya MR, Cascio D, Eisenberg D (2010) Crystallographic studies of prion protein (PrP) segments suggest how structural changes encoded by polymorphism at residue 129 modulate susceptibility to human prion disease. J BiolChem 285: 29671-29675.
- Apostol MI, Wiltzius JJ, Sawaya MR, Cascio D, Eisenberg D (2011) Atomic structures suggest determinants of transmission barriers in mammalian prion disease. Biochemistry 50: 2456-2463.
- Case DA, Darden TA, Cheatham TE, Simmerling III CL, Wang J, et al. (2010)Assisted Model Building with Energy Refinement. AMBER 11, University of California, San Francisco.
- Zhang JP, Yating H, Yiju W, Changyu W, Xiangsun Z (2012) The LBFGS quasi-Newtonian method for molecular modeling prion AGAAAAGA amyloid fibril. Natural Science 4: 1097-1108.
- Zhang J, Gao DY, Yearwood J (2011) A novel canonical dual computational approach for prion AGAAAAGA amyloid fibril molecular modeling. J TheorBiol 284: 149-157.
- Zhang J, Sun J, Wu C (2011) Optimal atomic-resolution structures of prion AGAAAAGA amyloid fibrils. J TheorBiol 279: 17-28.
- Zhang J (2011) Optimal molecular structures of prion AGAAAAGA amyloid fibrils formatted by simulated annealing. J Mol Model 17: 173-179.
- Zhang JP(2011) An effective simulated annealing refined replica exchange Markov chain Monte Carlo algorithm for the infectious disease model of H1N1 influenza pandemic. World Journal of Modelling and Simulation 7: 29-39.
- Zhang JP(2012) Atomic-resolution structures of prion AGAAAAGA amyloid fibrils. Amyloids: Composition, Functions and Pathology 10: 177-186.
- Zhang JP (2012) Computational studies on the prion AGAAAAGA amyloid fibrils molecular structures. Recent Advances in Crystallography 12: 297-312.
- Zhang J, Zhang Y (2013) Molecular dynamics studies on 3D structures of the hydrophobic region PrP(109-136). ActaBiochimBiophys Sin (Shanghai) 45: 509-519.

Select your language of interest to view the total content in your interested language

- Advanced DNA Sequencing
- Algorithm
- Animal and Tissue Engineering
- Applications of Bioinformatics
- Artificial Intelligence Studies
- Artificial intelligence
- Artificial neural networks
- Big data
- Bioinformatics Algorithms
- Bioinformatics Databases
- Bioinformatics Modeling
- Bioinformatics Tools
- Biology Engineering
- Biostatistics: Current Trends
- Cancer Proteomics
- Chemistry of Biology
- Clinical Proteomics
- Cloud Computation
- Cluster analysis
- Comparative genomics
- Comparative proteomics
- Computational Chemistry
- Computational Sciences
- Computational drug design
- Computer Science
- Computer-aided design (CAD)
- Current Proteomics
- Data Mining Current Research
- Data algorithms
- Data mining applications in genomics
- Data mining applications in proteomics
- Data mining in drug discovery
- Data mining tools
- Data modelling and intellegence
- Data warehousing
- Ethics in Synthetic Biology
- Evolution of social network
- Evolutionary Optimisation
- Evolutionary algorithm
- Evolutionary algorithm in datamining
- Evolutionary computation
- Evolutionary science
- Experimental Physics
- Findings on Machine Learning
- Gene Synthesis
- Genome annotation
- Genomic data mining
- Genomic data warehousing
- Handover
- Human Proteome Project Applications
- Hybrid soft computing
- Industrial Biotechnology
- Knowledge modelling
- Machine Learninng
- Mapping of genomes
- Mass Spectrometry in Proteomics
- Mathematical Modeling
- Mathematics for Computer Science
- Meta genomics
- Microarray Proteomics
- Models of Science
- Molecular and Cellular Proteomics
- Multi Objective Programming
- Neural Network
- Ontology Engineering
- P4 medicine
- Physics Models
- Protein Sequence Analysis
- Proteogenomics
- Proteome Profiling
- Proteomic Analysis
- Proteomic Biomarkers
- Proteomics Clinical Applications
- Proteomics Research
- Proteomics Science
- Proteomics data warehousing
- Python for Bioinformatics
- Quantitative Proteomics
- Robotics Research
- Scientific Computing
- Simulation Computer Science
- Soft Computing
- Statistical data mining
- Studies on Computational Biology
- Swarm Robotics
- Swarm intelligence
- Synthetic Biology
- Synthetic Biology medicine
- Synthetic Biotechnology
- Synthetic Genomics
- Synthetic biology drugs
- Systems Biology
- Technologies in Computer Science
- Theoretical Chemistry
- Theoretical Computer Science
- Theoretical Issues in Ergonomics Science
- Theoretical Methods
- Theoretical and Applied Science

- International Conference on
**Proteomics**and**Bioinformatics**

May 16-17, 2018 Singapore City, Singapore - International Conference on
**Metabolomics**and**Systems Biology**

June 11-12, 2018 London, UK - International Conference on
**Computational Biology**and**Bioinformatics**

Sep 05-06 2018 Tokyo, Japan - International Conference on Advancements in
**Bioinformatics**and**Drug Discovery**

November 26-27, 2018 Dublin, Ireland

- Total views:
**11683** - [From(publication date):

July-2014 - Mar 22, 2018] - Breakdown by view type
- HTML page views :
**7910** - PDF downloads :
**3773**

Peer Reviewed Journals

International Conferences
2018-19