Skip to main content
  • Research article
  • Open access
  • Published:

Bioinformatics in crosslinking chemistry of collagen with selective cross linkers

Abstract

Background

Identifying the molecular interactions using bioinformatics tools before venturing into wet lab studies saves the energy and time considerably. The present study summarizes, molecular interactions and binding energy calculations made for major structural protein, collagen of Type I and Type III with the chosen cross-linkers, namely, coenzyme Q10, dopaquinone, embelin, embelin complex-1 & 2, idebenone, 5-O-methyl embelin, potassium embelate and vilangin.

Results

Molecular descriptive analyses suggest, dopaquinone, embelin, idebenone, 5-O-methyl embelin, and potassium embelate display nil violations. And results of docking analyses revealed, best affinity for Type I (- 4.74 kcal/mol) and type III (-4.94 kcal/mol) collagen was with dopaquinone.

Conclusions

Among the selected cross-linkers, dopaquinone, embelin, potassium embelate and 5-O-methyl embelin were the suitable cross-linkers for both Type I and Type III collagen and stabilizes the collagen at the expected level.

Background

Collagen is the most abundant fibrillar protein in multicellular animals and the protein of importance for current scenario on development of tissue engineering materials for therapeutic applications. It has unique protein motif containing three super-coiled polyproline II helices with high percentage of imino acids. Thus, it contains repeating units of Gly-Xaa-Yaa, where, proline and 4-hydroxyproline occupy Xaa and Yaa positions respectively. Almost all prolines in the Yaa position of vertebrate collagen post-translationally modified to hydroxyproline in the presence of enzyme (prolyl-4-hydroxylase). Ever since the report of triple helical structure of collagen by Ramachandran and Kartha [1] nearly fifty years ago, there has been a considerable interest in the study of molecular properties and structure of the collagen. Molecular modeling of collagen has a long fascinating history, as Miller and Scheraga [2] reported computational model of collagen, especially evaluate the effects of specific side chains on conformation. Later Chen et al. [3] performed molecular docking to form a microfibril template based on the smith model [4]. Piez and Truzs [5] constructed three-dimensional energy-minimized model for calf-skin type-I collagen. Qi et al. [6] constructed the N-terminal te1opeptides model.

Cross linking of collagen is often necessary to improve the stability as well as resistance against enzymatic degradation. The existing cross linking agents such as formaldehyde, glutaraldehyde and epoxy compounds are all identified as cytotoxic, which limits the clinical applications and provokes scientific community to look for alternative natural cross linking/stabilizing agents. In addition, awareness on collagen modification and fibril coating in tanning has recently gained interest among Leather Chemist's. Ramasami [7] reported interaction of the collagen with small molecules like water, polyphenols and chromium (III) salt. Wu et al. [8] observed covalent interactions between glutaraldehyde and collagen. With regard to quinone tannage, Thomas and Kelly [9] reported, effect of pH on quinone fixation by collagen [9]. Stecker and Highberger [10] observed, nature of the buffer system plays an important role in collagen-quinone reactions. Theis [11] reported, collagen upon treatment with quinone, an appreciable increase in shrinkage temperature ultimately increased the structural strength of the collagen due to the bonding between polypeptide chains by means of the reaction of the reactive centers of quinone with amino groups in juxtaposition.

With regard to the bioinformatics on cross-linking chemistry of collagen, only very few reports are available collagen with polyphenols [12], and with alginic acid [13]. In order to select the suitable cross-linkers for the preparation of collagen based biomaterials, in the present study we made an attempt on bioinformatics in cross linking of selected cross linkers coenzyme Q10, dopaquinone, embelin, embelin complex-1 & 2, idebenone, 5-O-methyl embelin, potassium embelate (unless otherwise mentioned without potassium metal) and vilangin with collagen of both Type I and Type III. Since most of the selected cross-linkers are of from natural sources, the results of the study will find application in wet lab conditions.

Methods

Protein preparation

Restricting the number of repeating units is necessary in the modeling and simulation of a large macromolecule like collagen. In the present study, 24-mer collagen triple helix was constructed by Object Technology Framework (OTF) using the GENCOLLAGEN package [14]. The 24-residue long triple helix constructed corresponds to the residues 193 to 216 (2a 1 and 1a 2 chains) of the native type I collagen except residue 204 of the a 1 chain, where Ala of native collagen was replaced by lysine in order to study the interaction of ligands with the side chains of basic amino acids. In case of type III collagen, all the three are identically α 1chains (3 α 1chains).

Ligand preparation

Chemical structures of ligands such as embelin [CID no: 3218], vilangin [CID no: 417182], coenzyme Q10 [CID no: 5281915], idebenone [CID No: 3686], potassium embelate [CID no: 23677950], 5-O-methyl embelin [CID no: 171489] and dopaquinone [CID no: 439316] were retrieved from Pubchem compound database [15]. Unavailable three dimensional structures such as embelin complex-1 & 2 were generated using ACD/ChemSketch [16].

Docking setup

Docking was performed using Autodock 4. Autodock combines energy evaluation through precalculated grids of affinity potential employing various search algorithms to find the suitable binding position for a ligand on a given protein [17]. Kollman united atom charges and polar hydrogens were added to the protein PDB using Autodock tools [17]. All rotatable bonds in the ligands were kept free to allow for flexible docking. Grid size was set to 40 × 40 × 40 grid points (x, y and z), with spacing between grid points kept at 0.375 Å. The Lamarckian genetic algorithm was chosen to search for the best conformers. Standard docking protocol was applied. One hundred independent docking runs were carried out for each ligand was generated by using genetic algorithm searches.

Molecular descriptors calculation

Quantitative structure-activity relationships (QSARs) correlate the response with molecular properties of compounds under interest. Any compound to be considered as a lead must possess acceptable scores for all of the descriptors. Molinspiration [18] was used to calculate thirteen descriptors such as logP, polar surface area, molecular weight, number of atoms, number of O or N, number of OH or NH, number of rotatable bonds, volume, drug likeness (includes GPCR ligand, ion channel modulator, kinase inhibitor and nuclear receptor ligand) and number of violations to Lipinski's rule for all ligands taken for the analysis [19].

Results and discussion

Molecular Physicochemical and the Drug-Likeness are the two properties that are important for considering a compound to become a successful drug. Table 1 depicts the International Union of Pure and Applied Chemistry (IUPAC) name and Simplified Molecular Input Line Entry Specification (SMILES) of ligands. Two dimensional structure of dopaquinone, embelin, potassium embelate and 5-O-methyl embelin also represented in additional files 1, 2, 3 and 4. A chemical nomenclature is a set of rules followed to generate systematic names for chemical compounds. IUPAC nomenclature is a universal chemical nomenclature. It is developed and kept up to date under the auspices of the International Union of Pure and Applied Chemistry (IUPAC). The simplified molecular input line entry specification (SMILES) is a specification for unambiguously describing the structure of chemical molecules using short ASCII strings. SMILES strings can be imported by most molecule editors for conversion back into two-dimensional drawings or three-dimensional models of the molecules.

Table 1 Ligand Molecules - IUPAC Name and SMILES

The rule formulated by Christopher A. Lipinski et al. [20] considered as the thumb rule thumb rule to evaluate drug likeness, or determine if a chemical compound with a certain pharmacological or biological activity has properties that would make it a likely orally active drug in humans. The rule describes molecular properties important for a drug's pharmacokinetics in the human body, including their absorption, distribution, metabolism, and excretion ("ADME"). The rule is important for drug development where a pharmacologically active lead structure is optimized step-wise for increased activity and selectivity, as well as drug-like properties as described by Lipinski's rule.

LogP (Octanol-water partition coefficient) is used as important tool in both quantitative structure-activity relationship (QSAR) studies and rational drug design as a measure of molecular hydrophobicity. Hydrophobicity affects drug absorption, bioavailability, hydrophobic drug-receptor interactions, metabolism of molecules, as well as their toxicity. LogP has become a key parameter in studies of the environmental fate of chemicals. In the present study LogP value of dopaquinone was -2.684, which indicates more hydrophilic nature, whereas LogP value of coenzyme Q 10 was 10.509, indicates more lipophilic or hydrophobic nature. On other hand potassium embelate, idebenone, embelin and 5-O-methyl embelin demonstrated LogP value between 2.2-4.8 as shown in Table-2. LogP value less than 5 will be preferred for drug likeness property.

The preferred range of molecular weight for drug likeness property was 160-480 g/mol as reported by Tambunan and Wulandari [21]. The molecular weight of the selected cross-linkers was calculated as 195.17 g/mol (dopaquinone), 293.38 g/mol (potassium embelate), 294.39 g/mol (embelin), was 308.41 g/mol (5-O-methyl embelin) and 338.44 g/mol (idebenone). With regard to the preferred number of N, O (hydrogen bond acceptors) and OH & NH (hydrogen bond donors) 10 and or less than 10 and 5 and or less than 5 respectively compliance with the rule. From table 2, it has been observed that embelin complex-2 showed greater than 10 with respect to N and O, whereas the embelin complex-1 showed greater than 5 with respect to OH and NH. Further, the preferred number of rotatable bonds (rotb) is 15 and or less than 15, and we observed that it was greater than 15 for the cross linkers vilangin, embelin complex-1, embelin complex-2 and coenzyme Q 10 (Table 2). The preferred number of Violations is 0, and we observed nil violations for the dopaquinone, embelin, idebenone, 5-O-methyl embelin, and potassium embelate and suggest all the five chosen ligands satisfy well with thumb rule (Table-2).

Table 2 The descriptor analysis helped in the identification of the better ligand

With regard to affinity and binding energy calculations, each chosen ligands displayed different affinities with the collagen types (I & III). Lead '1' dopaquinone alone showed the best affinity with both the type-I & III collagen (-4.74 kcal/mol and -4.94 kcal/mol) followed by potassium embelate, embelin, and 5-O-methyl embelin for type I collagen (Table 3 & Figure 1). However, with type III collagen, followed by dopaquinone, embelin showed the second best affinity, which was followed by potassium embelate.

Table 3 The binding energy helped in the identification of the better ligand
Figure 1
figure 1

Docking results of the four lead molecules namely dopaquinone, embelin, potassium embelate and 5-O-methyl embelin. [Where a-e (docking results of the four lead molecules namely dopaquinone, embelin, potassium embelate and 5-O-methyl embelin) with type I collagen. Where f-i (docking results of the four lead molecules namely dopaquinone, embelin, potassium embelate and 5-O-methyl embelin) with type III collagen].

With regard to the interaction sites, bond sites and bond distance, calculated based on the bioinformatics tool for the cross-linking between collagen (Type I) and four ligands (four lead molecules) as shown in table 4, emphasizes, the first lead molecule 'dopaquinone' interacts with collagen (Type I) through non-polar aminoacids (Proline and Valine) and polar aminoacid (Lysine and Asparagine). Second lead molecule 'potassium embelate' interacts with collagen (Type I) through non-polar aminoacid (Valine). Third lead molecule 'embelin' interacts with collagen (Type I) through non-polar aminoacids (Proline, Glycine and Valine). Fourth lead molecule '5-O-methyl embelin' interacts with collagen (Type I) through non-polar aminoacids (Proline and Valine) and polar aminoacid (Tyrosine).

Table 4 Interaction sites, bond sites and bond distance, between collagen (Type I) and four ligands.

And with reference to type III collagen, the interaction sites, bond sites and bond distance, calculated based on the bioinformatics tool for the cross-linking of four ligands (four lead molecules) Table 5 emphasizes, the first lead molecule 'dopaquinone' interacts with collagen (Type III) through non-polar aminoacids (Proline and Alanine) and polar aminoacid (Lysine). Second lead molecule 'embelin' interacts with collagen (Type III) through non-polar aminoacid (Glycine) and polar aminoacid (Lysine). Third lead molecule 'potassium embelate' interacts with collagen (Type III) through non-polar aminoacids (Proline and Glycine). Fourth lead molecule '5-O-methyl embelin' interacts with collagen (Type III) through non-polar aminoacids (Proline) and polar aminoacid (Lysine).

Table 5 Interaction sites, bond sites and bond distance, between collagen (Type III) and four ligands

The reason behind the choice of quinone based compounds is described below. Meunier and Seyewetz [22] identified the remarkable stabilizing (tanning) potential of p-benzoquinone and after them it has been named as quinone tannage. Later, Thomas and Kelly [9] comprehensively investigated the stabilization of collagen by quinone. Wilson [23] reported benzoquinone tans/stabilizes well in alcoholic solution, and well correlated with many aspects of formaldehyde tannage. Theis [11] reported collagen treated with quinone has increased shrinkage temperature and the structural strength. Later Suparno [24] reported two salient features of quinone-tannage; viz., an increase in shrinkage temperature (Ts) (> 90°C) and high resistant to the proteolytic degradation. Further, Covington [25] also reported increase in shrinkage temperature upon treating collagen with quinone alone; however, he said the toxicity has to be considered before going for commercial use.

With regard to molecular interaction studies using bioinformatics tools, Vaidyanathan et al [26] studied the interactions of five ligands (2-Hydroxyethyl methacrylate (HEMA), Glutaraldehyde-HEMA adduct, Glyceryl dimethacrylate, Methacryloyloxyethyl maleate and Acryloyloxyethyl citraconate) with collagen and reported that steric and electrostatic complementarity interactions form the potential basis of binding between dentin adhesive ligands and type 1 collagen. Madhan et al. [27] studied stabilization of collagen by catechin, reported on the interactions of a catechin with 24-mer collagen triple helix through hydrogen bonding interaction. Mitra et al. [13] studied, thermal stabilization of collagen by alginic acid and reported on the interactions of alginic acid with 24-mer collagen triple helix through hydrogen bonding interaction with binding energy of -7.28 kcal/mol.

Conclusion

Present study provides the molecular interaction view of quinones with both type I collagen and type III collagen, as a first eye opener on quinone and collagen interactions using bioinformatics tool. Out of nine quinone studied, dopaquinone, potassium embelate, embelin and 5-O-methyl embelin showed better affinities with both type I collagen and type III collagen. Hence dopa quinone, embelin, potassium embelate and 5-O-methyl embelin could be developed as potential cross-linking/stabilization agent of collagen preparation and found application as wound dressing sheet in clinical applications.

References

  1. Ramachandran GN, Kartha G: Structure of collagen. Nature. 1955, 176: 593-595. 10.1038/176593a0.

    Article  PubMed  CAS  Google Scholar 

  2. Miller MH, Scheraga HA: Calculation of the structures of collagen models. J Polymer Sci Polym Symp. 1976, 54: 171-200.

    CAS  Google Scholar 

  3. Chen JM, Kung CE, Feairheller SH, Brown EM: An energetic evaluation of a "Smith" collagen microfibril model. J Protein Chem. 1991, 10: 535-552. 10.1007/BF01025482.

    Article  PubMed  CAS  Google Scholar 

  4. Smith JW: Molecular pattern in native collagen. Nature. 1968, 219: 157-158. 10.1038/219157a0.

    Article  PubMed  CAS  Google Scholar 

  5. Piez KA, Trus BL: Sequence regularities and packing of collagen molecules. J Mol Biol. 1978, 122: 419-432. 10.1016/0022-2836(78)90419-9.

    Article  PubMed  CAS  Google Scholar 

  6. Qi PX, Brown EM: Molecular modeling of N-terminal telopeptides of bovine type I collagen. J Am Leather Chem Assoc. 2002, 97: 235-242.

    CAS  Google Scholar 

  7. Ramasami T: Approach towards a unified theory for tanning: Wilson's dream. J Am Leather Chem Assoc. 2001, 96: 290-304.

    CAS  Google Scholar 

  8. Wu HC, Wang TW, Kang PL, Tsuang YH, Sun JS, Lin FH: Coculture of endothelial and smooth muscle cells on a collagen membrane in the development of a small-diameter vascular graft. Biomat. 2007, 28 (7): 1385-1392. 10.1016/j.biomaterials.2006.11.012.

    Article  CAS  Google Scholar 

  9. Thomas AW, Kelly MW: Further studies of quinone tannage. Indust Eng Chem. 1926, 18: 383-385. 10.1021/ie50196a014.

    Article  CAS  Google Scholar 

  10. Stecker HC, Highberger JH: Quinone tannage. I. The weight increase of collagen in buffered quinone systems and its relation to oxidation. J Am Leather Chem Assoc. 1942, 36: 226-

    Google Scholar 

  11. Theis ER: The collagen-quinone reaction: I. Fixation and thermolability as a function of pH values. J Biol Chem. 1945, 157: 23-34.

    CAS  Google Scholar 

  12. Madhan B, Dhathathreyan A, Subramanian V, Ramasami T: Investigations on geometrical features in induced ordering of collagen by small molecules. Proc Indian Acad Sci (Chem Sci). 2003, 115 (5 & 6): 751-766.

    Article  CAS  Google Scholar 

  13. Mitra T, Sailakshmi G, Gnanamani A, Thirupathi Kumara Raja S, Thiruselvi T, Mangala Gowri V, Naga Vignesh S, Ramesh G, Mandal AB: Preparation and characterization of a thermostable and biodegradable biopolymers using natural cross-linker. Int J Biol Macromol. 2011, 48: 276-285. 10.1016/j.ijbiomac.2010.11.011.

    Article  PubMed  CAS  Google Scholar 

  14. The GenCollagen Database. [http://www.cgl.ucsf.edu./cgi-bin/gencollagen.py]

  15. The PubChem Structure search. [http://pubchem.ncbi.nlm.nih.gov/search/search.cgi]

  16. ACD, ChemSketch Version 12: 2009, Advanced Chemistry Development, Inc., Toronto, ON, Canada, [http://www.acdlabs.com/download/chemsk.html]

  17. Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, Olson AJ: AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comput Chem. 2009, 30: 2785-27910. 10.1002/jcc.21256.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. The Molinspiration Database. [http://www.molinspiration.com]

  19. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ: Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev. 2001, 46: 3-26. 10.1016/S0169-409X(00)00129-0.

    Article  PubMed  CAS  Google Scholar 

  20. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ: Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev. 1997, 23: 3-25. 10.1016/S0169-409X(96)00423-1.

    Article  CAS  Google Scholar 

  21. Tambunan USF, Wulandari EK: Identification of a better Homo sapiens Class II HDAC inhibitor through binding energy calculations and descriptor analysis. BMC Bioinform. 2010, 11 (Suppl 7): S16-10.1186/1471-2105-11-S7-S16.

    Article  Google Scholar 

  22. Meunier L, Seyewetz A: A new method of tanning. Compt rend Acad. 1908, 146: 987-989.

    Google Scholar 

  23. Wilson JA: The chemistry of leather manufacture. 1929, New York: The Chemical catalog company, Inc, 2

    Google Scholar 

  24. Suparno O: Quinone tannage. 2008, [http://web.ipb.ac.id/~ono.suparno/html/quinone_tanning.html]

    Google Scholar 

  25. Covington AD: Modern tanning chemistry. Chem Soc Rev. 1997, 26: 111-126. 10.1039/cs9972600111.

    Article  CAS  Google Scholar 

  26. Vaidyanathan J, Vaidyanathan TK, Yadav P, Linaras CE: Collagen-ligand interaction in dentinal adhesion: computer visualization and analysis. Biomat. 2001, 22: 2911-2920. 10.1016/S0142-9612(01)00038-2.

    Article  CAS  Google Scholar 

  27. Madhan B, Subramanian V, Raghava Rao J, Nair BU, Ramasami T: Stabilization of collagen using plant polyphenol: Role of catechin. Int J Biol Macromol. 2005, 37: 47-53. 10.1016/j.ijbiomac.2005.08.005.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

One of the authors thanks the Council of Scientific and Industrial Research (CSIR), New Delhi, India for financial assistance in the form Senior Research Fellowship (SRF) is gratefully acknowledged.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to A Gnanamani.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RN, who designed and executed the experiment as part of Doctorate programme. GA & ABM Research supervisor and Director, who guided the experiment and interpretation. SS and RG who, helped to carry out the experiment and interpretation. All authors read and approved the final manuscript.

Electronic supplementary material

13104_2011_1173_MOESM1_ESM.CDX

Additional file 1:2D structure of dopaquinone. Figure represents the two dimensional structure of dopaquinone. (CDX 2 KB)

Additional file 2:2D structure of embelin. Figure represents the two dimensional structure of embelin. (CDX 2 KB)

13104_2011_1173_MOESM3_ESM.CDX

Additional file 3:2D structure of potassium embelate. Figure represents the two dimensional structure of potassium embelate. (CDX 2 KB)

13104_2011_1173_MOESM4_ESM.CDX

Additional file 4:2D structure of 5-O-methyl embelin. Figure represents the two dimensional structure of 5-O-methyl embelin. (CDX 2 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Swamy, R.N., Gnanamani, A., Shanmugasamy, S. et al. Bioinformatics in crosslinking chemistry of collagen with selective cross linkers. BMC Res Notes 4, 399 (2011). https://doi.org/10.1186/1756-0500-4-399

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1756-0500-4-399

Keywords