Mapping intermolecular interactions and active site conformations: from human MMP-1 crystal structure to molecular dynamics free energy calculations

The zinc-dependent Matrix Metalloproteinases (MMPs) found within the extracellular matrix (ECM) of vertebrates are linked to pathological processes such as arthritis, skin ulceration and cancer. Although a general backbone proteolytic mechanism is understood, crystallographic data continue to suggest an active site that is too narrow to encompass the respective substrate. We present a fully parameterised Molecular Dynamics (MD) study of the structural properties of an MMP-1-collagen crystallographic structure (Protein Data Bank – 4AUO), followed by an exploration of the free energy surface of a collagen polypeptide chain entering the active site, using a combined meta-dynamics and umbrella sampling (MDUS) approach. We conclude that the interactions between MMP-1 and the collagen substrate are in good agreement with a number of experimental studies. As such, our unrestrained MD simulations and our MDUS results, which indicate an energetic barrier for a local uncoiling and insertion event, can inform future investigations of the collagen-peptide non-bonded association steps with the active site prior to proteolytic mechanisms. The elucidation of such free energy barriers provides a better understanding of the role of the enzyme in the ECM and is important in the design of future MMP inhibitors

An equine tendon model for studying intra-tendinous shear in tendons that have more than one muscle contribution

Human Achilles tendon is composed of three smaller sub-tendons and exhibits non-uniform internal displacements, which decline with age and after injury, suggesting a potential role in the development of tendinopathies. Studying internal sliding behaviour is therefore important but difficult in human Achilles tendon. Here we propose the equine deep digital flexor tendon (DDFT) and its accessory ligament (AL) as a model to understand the sliding mechanism. The AL-DDFT has a comparable sub-bundle structure, is subjected to high and frequent asymmetric loads and is a natural site of injury similar to human Achilles tendons. Equine AL-DDFT were collected and underwent whole tendon level (n=7) and fascicle level (n=7) quasi-static mechanical testing. Whole tendon level testing was performed by sequentially loading through the proximal AL and subsequently through the proximal DDFT and recording regional strain in the free structures and joined DDFT and AL. Fascicle level testing was performed with focus on the inter-sub-bundle matrix between the two structures at the junction. Our results demonstrate a significant difference in the regional strain between the joined DDFT and AL and a greater transmission of force from the AL to the DDFT than vice versa. These results can be partially explained by the mechanical properties and geometry of the two structures and by differences in the properties of the interfascicular matrices. In conclusion, this tendon model successfully demonstrates that high displacement discrepancy occurs between the two structures and can be used as an easy-access model for study intra-tendinous shear mechanics at the sub-tendon level. STATEMENT OF SIGNIFICANCE: : Our study provides a naturally occurring and easily accessible equine model to study the complex behaviour of sub-tendons within the human Achilles tendon, which is likely to play a critical role in the pathogenesis of tendon disease. Our results demonstrate that the difference in material stiffness between the equine AL and DDFT stems largely from differences in the inter-fascicular matrix and furthermore that differences in strain are maintained in distal parts of the tightly joined structure. Furthermore, our results suggest that distribution of load between sub-structures is highly dependent on the morphological relationship between them; a finding that has important implications for understanding Achilles tendon mechanical behaviour, injury mechanisms and rehabilitation

Detection of age-related changes in tendon molecular composition by Raman spectroscopy—potential for rapid, non-invasive assessment of susceptibility to injury

The lack of clinical detection tools at the molecular level hinders our progression in preventing age-related tendon pathologies. Raman spectroscopy can rapidly and non-invasively detect tissue molecular compositions and has great potential for in vivo applications. In biological tissues, a highly fluorescent background masks the Raman spectral features and is usually removed during data processing, but including this background could help age differentiation since fluorescence level in tendons increases with age. Therefore, we conducted a stepwise analysis of fluorescence and Raman combined spectra for better understanding of the chemical differences between young and old tendons. Spectra were collected from random locations of vacuum-dried young and old equine tendon samples (superficial digital flexor tendon (SDFT) and deep digital flexor tendon (DDFT), total n = 15) under identical instrumental settings. The fluorescence-Raman spectra showed an increase in old tendons as expected. Normalising the fluorescence-Raman spectra further indicated a potential change in intra-tendinous fluorophores as tendon ages. After fluorescence removal, the pure Raman spectra demonstrated between-group differences in CH2 bending (1450 cm-1) and various ring-structure and carbohydrate-associated bands (1000-1100 cm-1), possibly relating to a decline in cellular numbers and an accumulation of advanced glycation end products in old tendons. These results demonstrated that Raman spectroscopy can successfully detect age-related tendon molecular differences

Lysine-Arginine Advanced Glycation End-Product Cross-links and the Effect on Collagen Structure: A Molecular Dynamics Study

The accumulation of advanced glycation end-products is a fundamental process that is central to age-related decline in musculoskeletal tissues and locomotor system function and other collagen-rich tissues. However, although computational studies of advanced glycation end-product cross-links could be immensely valuable, this area remains largely unexplored given the limited availability of structural parameters for the derivation of force fields for Molecular Dynamics simulations. In this article, we present the bonded force constants, atomic partial charges and geometry of the arginine-lysine cross-links DOGDIC, GODIC, and MODIC. We have performed in vacuo Molecular Dynamics simulations to validate their implementation against quantum mechanical frequency calculations. A DOGDIC advanced glycation endproduct cross-link was then inserted into a model collagen fibril to explore structural changes of collagen and dynamics in interstitial water. Unlike our previous studies of glucosepane, our findings suggest that intra-collagen DOGDIC cross-links furthers intra-collagen peptide hydrogen-bonding and does not promote the diffusion of water through the collagen triple helices

Intra-molecular lysine-arginine derived advanced glycation end-product cross-linking in Type I collagen: A molecular dynamics simulation study

Covalently cross-linked advanced glycation end products (AGE) are among the major post-translational modifications to proteins as a result of non-enzymatic glycation. The formation of AGEs has been shown to have adverse effects on the properties of the collagenous tissue; they are even linked to a number of age related disorders. Little is known about the sites at which these AGEs form or why certain sites within the collagen are energetically more favourable than others. In this study we have used a proven fully atomistic molecular dynamics approach to identify six sites where the formation of the intra-molecular 3-deoxyglucosone-derived imidazolium cross-link (DOGDIC) is energetically favourable. We have also conducted a comparison of these positions with those of the more abundant glucosepane cross-link, to determine any site preference. We show that when we consider both lysine and arginine AGEs, they exhibit a prevalence to form within the gap region of the collagen fibril

Individual variation in Achilles tendon morphology and geometry changes susceptibility to injury

The unique structure of the Achilles tendon, combining three smaller sub-tendons, enhances movement efficiency by allowing individual control from connected muscles. This requires compliant interfaces between sub-tendons, but compliance decreases with age and may account for increased injury frequency. Current understanding of sub-tendon sliding and its role in the whole Achilles tendon function is limited. Here we show changing the degree of sliding greatly affects the tendon mechanical behaviour. Our in vitro testing discovered distinct sub-tendon mechanical properties in keeping with their mechanical demands. In silico study based on measured properties, subject-specific tendon geometry, and modified sliding capacity demonstrated age-related displacement reduction similar to our in vivo ultrasonography measurements. Peak stress magnitude and distribution within the whole Achilles tendon are affected by individual tendon geometries, the sliding capacity between sub-tendons, and different muscle loading conditions. These results suggest clinical possibilities to identify patients at risk and design personalised rehabilitation protocols

Effect on the mechanical properties of type I collagen of intra-molecular lysine-arginine derived advanced glycation end-product cross-linking

Non-enzymatic advanced glycation end product (AGE) cross-linking of collagen molecules has been hypothesised to result in significant changes to the mechanical properties of the connective tissues within the body, potentially resulting in a number of age related diseases. We have investigated the effect of two of these cross-links, glucosepane and DOGDIC, on the tensile and lateral moduli of the collagen molecule through the use of a steered molecular dynamics approach, using previously identified preferential formation sites for intra-molecular cross-links. Our results show that the presence of intra-molecular AGE cross-links increases the tensile and lateral Young’s moduli in the low strain domain by between 3.0–8.5% and 2.9–60.3% respectively, with little effect exhibited at higher strains

ForceGen: atomic covalent bond value derivation for Gromacs

A large number of crystallographic protein structures include ligands, small molecules and post-translational modifications. Atomic bond force values for computational atomistic models of post-translational or non-standard amino acids, metal binding active sites, small molecules and drug molecules are not readily available in most simulation software packages. We present ForceGen, a Java tool that extracts the bond stretch and bond angle force values and equilibrium values from the Hessian of a Gaussian vibrational frequency analysis. The parameters are compatible with force fields derived using the second order tensor of the Hessian. The output is formatted with the Gromacs topology in mind. This study further demonstrates the use of ForceGen over the quantum mechanically derived structures of a small organic solvent, a naturally occurring protein crosslink derived from two amino acids following post-translational modification and the amino acid ligands of a zinc ion. We then derive Laplacian bond orders to understand how the resulting force values relate back to the quantum mechanical model. The parameterisation of the organic solvent, toluene, was verified using Molecular Mechanics simulations. The structural data from the simulation compared well with the quantum mechanical structure and the system density compared well with experimental values

Relative orientation of collagen molecules within a fibril: a homology model for homo sapiens type I collagen

Type I collagen is an essential extracellular protein that plays an important structural role in tissues that require high tensile strength. However, owing to the molecule’s size, to date no experimental structural data are available for the Homo sapiens species. Therefore, there is a real need to develop a reliable homology model and a method to study the packing of the collagen molecules within the fibril. Through the use of the homology model and implementation of a novel simulation technique, we have ascertained the orientations of the collagen molecules within a fibril, which is currently below the resolution limit of experimental techniques. The longitudinal orientation of collagen molecules within a fibril has a significant effect on the mechanical and biological properties of the fibril, owing to the different amino acid side chains available at the interface between the molecules

Computational study of glucosepane–water and hydrogen bond formation: an electron topology and orbital analysis

The collagen protein provides tensile strength to the extracellular matrix in addition to localising cells, proteins and protein cofactors. Collagen is susceptible to a build up of glycation modifications as a result of an exceptionally long half-life. Glucosepane is a collagen cross-linking advanced glycation end product; the structural and mechanical effects of glucosepane are still the subjects of much debate. With the prospect of an ageing population, the management and treatment of age-related diseases is becoming a pressing concern. One area of interest is the isolation of hydrated glucosepane, which has yet to be reported at an atomistic level. This study presents a series of glucosepane–water complexes within an implicit aqueous environment. Electronic structure calculations were performed using density functional theory and a high level basis set. Hydrogen bonds between glucosepane and explicit water were identified by monitoring changes to covalent bonds, calculating levels of electron donation from Natural Bonding Orbital analysis and the detection of bond critical points. Hydrogen bond strength was calculated using second-order perturbation calculations. The combined results suggest that glucosepane is very hydrophilic, with the imidazole feature being energetically more attractive to water than either hydroxyl group, although all hydrogen bonds, regardless of bond strength, were electrostatic in nature. Our results are in growing support of an earlier hypothesis that cross-links may result in an increase in interstitial water retention, which would permit the collagen fibril to swell, thereby potentially affecting the tensile and compression properties and biological function of connective tissues