The effect of point mutations on structure and mechanical properties of collagen-like fibril: A molecular dynamics study

Designing collagens to shed light on the multi-scale structure-function mapping of matrix disorders

Collagens are the most abundant structural proteins in the extracellular matrix of animals and play crucial roles in maintaining the structural integrity and mechanical properties of tissues and organs while mediating important biological processes. Fibrillar collagens have a unique triple helix structure with a characteristic repeating sequence of (Gly-X-Y)n. Variations within the repetitive sequence can cause misfolding of the triple helix, resulting in heritable connective tissue disorders. The most common variations are single-point missense mutations that lead to the substitution of a glycine residue with a bulkier amino acid (Gly → X). In this review, we will first discuss the importance of collagen's triple helix structure and how single Gly substitutions can impact its folding, structure, secretion, assembly into higher-order structures, and biological functions. We will review the role of "designer collagens," i.e., synthetic collagen-mimetic peptides and recombinant bacterial collagen as model systems to include Gly → X substitutions observed in collagen disorders and investigate their impact on structure and function utilizing in vitro studies. Lastly, we will explore how computational modeling of collagen peptides, especially molecular and steered molecular dynamics, has been instrumental in probing the effects of Gly substitutions on structure, receptor binding, and mechanical stability across multiple length scales.

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.

Investigating the impact of missense mutations in hCES1 by in silico structure-based approaches

Genetic variations in drug-metabolizing enzymes have been reported to influence pharmacokinetics, drug dosage, and other aspects that affect therapeutic outcomes. Most particularly, non-synonymous single-nucleotide polymorphisms (nsSNPs) resulting in amino acid changes disrupt potential functional sites responsible for protein activity, structure, or stability, which can account for individual susceptibility to disease and drug response. Investigating the impact of nsSNPs at a protein's structural level is a key step in understanding the relationship between genetic variants and the resulting phenotypic changes. For this purpose, in silico structure-based approaches have proven their relevance in providing an atomic-level description of the underlying mechanisms. The present review focuses on nsSNPs in human carboxylesterase 1 (hCES1), an enzyme involved in drug metabolism. We highlight how prioritization of functional nsSNPs through computational prediction techniques in combination with structure-based approaches, namely molecular docking and molecular dynamics simulations, is a powerful tool in providing insight into the underlying molecular mechanisms of nsSNPs phenotypic effects at microscopic level. Examples of in silico studies of carboxylesterases (CESs) are discussed, ranging from exploring the effect of mutations on enzyme activity to predicting the metabolism of new hCES1 substrates as well as to guiding rational design of CES-selective inhibitors.

A systematic framework for molecular dynamics simulations of protein post-translational modifications

By directly affecting structure, dynamics and interaction networks of their targets, post-translational modifications (PTMs) of proteins play a key role in different cellular processes ranging from enzymatic activation to regulation of signal transduction to cell-cycle control. Despite the great importance of understanding how PTMs affect proteins at the atomistic level, a systematic framework for treating post-translationally modified amino acids by molecular dynamics (MD) simulations, a premier high-resolution computational biology tool, has never been developed. Here, we report and validate force field parameters (GROMOS 45a3 and 54a7) required to run and analyze MD simulations of more than 250 different types of enzymatic and non-enzymatic PTMs. The newly developed GROMOS 54a7 parameters in particular exhibit near chemical accuracy in matching experimentally measured hydration free energies (RMSE = 4.2 kJ/mol over the validation set). Using this tool, we quantitatively show that the majority of PTMs greatly alter the hydrophobicity and other physico-chemical properties of target amino acids, with the extent of change in many cases being comparable to the complete range spanned by native amino acids.

Post-translational modifications, i.e. chemical changes of protein amino acids, play a key role in different cellular processes, ranging from enzymatic activation to transcription and translation regulation to disease development and aging. However, our understanding of their effects on protein structure, dynamics and interaction networks at the atomistic level is still largely incomplete. In particular, molecular dynamics simulations, despite their power to provide a high-resolution insight into biomolecular function and underlying mechanisms, have been limited to unmodified, native proteins due to a surprising deficiency of suitable tools and systematically developed parameters for treating modified proteins. To fill this gap, we develop and validate force field parameters, an essential part of the molecular dynamics method, for more than 250 different types of enzymatic and non-enzymatic post-translational modifications. Additionally, using this tool, we quantitatively show that microscopic properties of target amino acids, such as hydrophobicity, are greatly affected by the majority of modifications. The parameters presented in this study greatly expand the range of applicability of computational methods, and in particular molecular dynamics simulations, to a large set of new systems with utmost biological and biomedical importance.

Vienna-PTM web server: a toolkit for MD simulations of protein post-translational modifications

Post-translational modifications (PTMs) play a key role in numerous cellular processes by directly affecting structure, dynamics and interaction networks of target proteins. Despite their importance, our understanding of protein PTMs at the atomistic level is still largely incomplete. Molecular dynamics (MD) simulations, which provide high-resolution insight into biomolecular function and underlying mechanisms, are in principle ideally suited to tackle this problem. However, because of the challenges associated with the development of novel MD parameters and a general lack of suitable computational tools for incorporating PTMs in target protein structures, MD simulations of post-translationally modified proteins have historically lagged significantly behind the studies of unmodified proteins. Here, we present Vienna-PTM web server (http://vienna-ptm.univie.ac.at), a platform for automated introduction of PTMs of choice to protein 3D structures (PDB files) in a user-friendly visual environment. With 256 different enzymatic and non-enzymatic PTMs available, the server performs geometrically realistic introduction of modifications at sites of interests, as well as subsequent energy minimization. Finally, the server makes available force field parameters and input files needed to run MD simulations of modified proteins within the framework of the widely used GROMOS 54A7 and 45A3 force fields and GROMACS simulation package.