Insights into the structure and dynamics of lysyl oxidase propeptide, a flexible protein with numerous partners

Lysyl oxidase (LOX) catalyzes the oxidative deamination of lysine and hydroxylysine residues in collagens and elastin, which is the first step of the cross-linking of these extracellular matrix proteins. It is secreted as a proenzyme activated by bone morphogenetic protein-1, which releases the LOX catalytic domain and its bioactive N-terminal propeptide. We characterized the recombinant human propeptide by circular dichroism, dynamic light scattering, and small-angle X-ray scattering (SAXS), and showed that it is elongated, monomeric, disordered and flexible (Dmax: 11.7 nm, Rg: 3.7 nm). We generated 3D models of the propeptide by coarse-grained molecular dynamics simulations restrained by SAXS data, which were used for docking experiments. Furthermore, we have identified 17 new binding partners of the propeptide by label-free assays. They include four glycosaminoglycans (hyaluronan, chondroitin, dermatan and heparan sulfate), collagen I, cross-linking and proteolytic enzymes (lysyl oxidase-like 2, transglutaminase-2, matrix metalloproteinase-2), a proteoglycan (fibromodulin), one growth factor (Epidermal Growth Factor, EGF), and one membrane protein (tumor endothelial marker-8). This suggests new roles for the propeptide in EGF signaling pathway

Mechanism of fiber assembly: Treatment of Aβ peptide aggregation with a coarse-grained united-residue force field

The growth mechanism of β-amyloid (Aβ) peptide fibrils was studied by a physics-based coarse-grained united-residue model and molecular dynamics (MD) simulations. To identify the mechanism of monomer addition to an Aβ1-40 fibril, we placed an unstructured monomer at a distance of 20 Å from a fibril template and allowed it to interact freely with the latter. The monomer was not biased towards fibril conformation by either the force field or the MD algorithm. With the use of a coarse-grained model with replica-exchange molecular dynamics, a longer timescale was accessible, making it possible to observe how the monomers probe different binding modes during their search for the fibril conformation. Although different assembly pathways were seen, they all follow a dock-lock mechanism with two distinct locking stages, consistent with experimental data on fibril elongation. Whereas these experiments have not been able to characterize the conformations populating the different stages, we have been able to describe these different stages explicitly by following free monomers as they dock onto a fibril template and to adopt the fibril conformation (i.e., we describe fibril elongation step by step at the molecular level). During the first stage of the assembly ( docking ), the monomer tries different conformations. After docking, the monomer is locked into the fibril through two different locking stages. In the first stage, the monomer forms hydrogen bonds with the fibril template along one of the strands in a two-stranded β-hairpin; in the second stage, hydrogen bonds are formed along the second strand, locking the monomer into the fibril structure. The data reveal a free-energy barrier separating the two locking stages. The importance of hydrophobic interactions and hydrogen bonds in the stability of the Aβ fibril structure was examined by carrying out additional canonical MD simulations of oligomers with different numbers of chains (4-16 chains), with the fibril structure as the initial conformation. The data confirm that the structures are stabilized largely by hydrophobic interactions and show that intermolecular hydrogen bonds are highly stable and contribute to the stability of the oligomers as well. © 2010 Elsevier Ltd

Improvements and new functionalities of UNRES server for coarse-grained modeling of protein structure, dynamics, and interactions

In this paper we report the improvements and extensions of the UNRES server (https://unres-server.chem.ug.edu.pl) for physics-based simulations with the coarse-grained UNRES model of polypeptide chains. The improvements include the replacement of the old code with the recently optimized one and adding the recent scale-consistent variant of the UNRES force field, which performs better in the modeling of proteins with the β and the α+β structures. The scope of applications of the package was extended to data-assisted simulations with restraints from nuclear magnetic resonance (NMR) and chemical crosslink mass-spectroscopy (XL-MS) measurements. NMR restraints can be input in the NMR Exchange Format (NEF), which has become a standard. Ambiguous NMR restraints are handled without expert intervention owing to a specially designed penalty function. The server can be used to run smaller jobs directly or to prepare input data to run larger production jobs by using standalone installations of UNRES

Protein-Ligand Interaction Energy-Based Entropy Calculations: Fundamental Challenges For Flexible Systems

Entropy calculations represent one of the most challenging steps in obtaining the binding free energy in biomolecular systems. A novel computationally effective approach (IE) was recently proposed to calculate the entropy based on the computation of protein-ligand interaction energy directly from molecular dynamics (MD) simulations. We present a study focused on the application of this method to flexible molecular systems and compare its performance with well-established normal mode (NM) and quasiharmonic (QH) entropy calculation approaches. Our results raise substantial concerns on the general applicability of IE in terms of reproducibility, reasonable absolute values of the entropy and agreement with NM and QM approaches. IE shows significant variation in the computed entropy values depending on the MD frames chosen for calculations. These deviations render reproducibility of IE calculations to be far from sufficient. We conclude that IE is recommended to be used after substantial modifications with respect to its sampling methodology

Ab initio theory of helix-coil phase transition

In this paper we suggest a theoretical method based on the statistical mechanics for treating the alpha-helix-random coil transition in alanine polypeptides. We consider this process as a first-order phase transition and develop a theory which is free of model parameters and is based solely on fundamental physical principles. It describes essential thermodynamical properties of the system such as heat capacity, the phase transition temperature and others from the analysis of the polypeptide potential energy surface calculated as a function of two dihedral angles, responsible for the polypeptide twisting. The suggested theory is general and with some modification can be applied for the description of phase transitions in other complex molecular systems (e.g. proteins, DNA, nanotubes, atomic clusters, fullerenes).Comment: 24 pages, 3 figure

High resolution protein folding with a transferable potential

A generalized computational method for folding proteins with a fully transferable potential and geometrically realistic all-atom model is presented and tested on seven different helix bundle proteins. The protocol, which includes graph-theoretical analysis of the ensemble of resulting folded conformations, was systematically applied and consistently produced structure predictions of approximately 3 Angstroms without any knowledge of the native state. To measure and understand the significance of the results, extensive control simulations were conducted. Graph theoretic analysis provides a means for systematically identifying the native fold and provides physical insight, conceptually linking the results to modern theoretical views of protein folding. In addition to presenting a method for prediction of structure and folding mechanism, our model suggests that a accurate all-atom amino acid representation coupled with a physically reasonable atomic interaction potential (that does not require optimization to the test set) and hydrogen bonding are essential features for a realistic protein model.Comment: submitted to PNAS 2005-03-1

Assessment of chemical-crosslink-assisted protein structure modeling in CASP13

International audienceWith the advance of experimental procedures obtaining chemical crosslinking information is becoming a fast and routine practice. Information on crosslinks can greatly enhance the accuracy of protein structure modeling. Here, we review the current state of the art in modeling protein structures with the assistance of experimentally determined chemical crosslinks within the framework of the 13th meeting of Critical Assessment of Structure Prediction approaches. This largest‐to‐date blind assessment reveals benefits of using data assistance in difficult to model protein structure prediction cases. However, in a broader context, it also suggests that with the unprecedented advance in accuracy to predict contacts in recent years, experimental crosslinks will be useful only if their specificity and accuracy further improved and they are better integrated into computational workflows

An analysis and evaluation of the WeFold collaborative for protein structure prediction and its pipelines in CASP11 and CASP12

Every two years groups worldwide participate in the Critical Assessment of Protein Structure Prediction (CASP) experiment to blindly test the strengths and weaknesses of their computational methods. CASP has significantly advanced the field but many hurdles still remain, which may require new ideas and collaborations. In 2012 a web-based effort called WeFold, was initiated to promote collaboration within the CASP community and attract researchers from other fields to contribute new ideas to CASP. Members of the WeFold coopetition (cooperation and competition) participated in CASP as individual teams, but also shared components of their methods to create hybrid pipelines and actively contributed to this effort. We assert that the scale and diversity of integrative prediction pipelines could not have been achieved by any individual lab or even by any collaboration among a few partners. The models contributed by the participating groups and generated by the pipelines are publicly available at the WeFold website providing a wealth of data that remains to be tapped. Here, we analyze the results of the 2014 and 2016 pipelines showing improvements according to the CASP assessment as well as areas that require further adjustments and research

Topology of Type II REases revisited; structural classes and the common conserved core

Type II restriction endonucleases (REases) are deoxyribonucleases that cleave DNA sequences with remarkable specificity. Type II REases are highly divergent in sequence as well as in topology, i.e. the connectivity of secondary structure elements. A widely held assumption is that a structural core of five β-strands flanked by two α-helices is common to these enzymes. We introduce a systematic procedure to enumerate secondary structure elements in an unambiguous and reproducible way, and use it to analyze the currently available X-ray structures of Type II REases. Based on this analysis, we propose an alternative definition of the core, which we term the αβα-core. The αβα-core includes the most frequently observed secondary structure elements and is not a sandwich, as it consists of a five-strand β-sheet and two α-helices on the same face of the β-sheet. We use the αβα-core connectivity as a basis for grouping the Type II REases into distinct structural classes. In these new structural classes, the connectivity correlates with the angles between the secondary structure elements and with the cleavage patterns of the REases. We show that there exists a substructure of the αβα-core, namely a common conserved core, ccc, defined here as one α-helix and four β-strands common to all Type II REase of known structure