The role of extrahelical peptides in stabilization of collagen fibrils

A quantitative model for fibril assembly of type I collagen was extended to include the explicit effect of extrahelical peptides. The collagen molecule was simulated by rod-like sequences to which short, rigid tails were connected by "nondimensional" flexible joints. Three collagen structures were studied: (1) intact collagen, simulated by a rod of axial ratio 200 (The axial ratio x was taken as a segment length) with two tails of length x = 1 and x = 2, respectively, appended to each end; (2) pepsin-digested collagen, simulated by one rigid segment of length 200 and one tail of length 1; and (3) pronase-digested collagen, by a single rigid segment of length x = 200. Phase equilibria of such structures were calculated, using a lattice theory of Matheson and Flory, and the relation of the polymer-solvent interaction parameter chi to the equilibrium solubility was determined. The chi for each collagen species was then related to temperature (T) and ionic strength (I), based on the approximation that local (per segment) stabilization of collagen fibrils was due to hydrophobic and electrostatic forces only. Solubility vs temperature curves for all three collagen species were computed and compared to published experimental data. From the chi factors for each species, the composite chi was resolved into components representing energetic contributions of the extrahelical peptides relative to the helix, which were interpreted in terms of hydrophobic or electrostatic interactions stabilizing the collagen fibril.     

The role of hydrophobic bonding in collagen fibril formation: a quantitative model

A model has been developed for approximating the free energy of collagen fibril formation (ΔF_f) and the equilibrium solubility of collagen under physiological conditions. The model utilizes an expression of Flory for rodlike polymers, with the modification that the “pure” anisotropic phase is defined as a collagen fibril containing about 0.3 g water/g collagen. The model also assumes that χ₁, the polymer–solvent interaction term, is entirely due to hydrophobic effects. χ₁ is estimated from hydrophobic bond energies of amino acid side chains, using the results of Némethy and Scheraga. The temperature dependence of χ₁ is utilized to calculate equilibrium solubilities and ΔF_f as a function of temperature.     

Early assembly pathways of type I collagen.

A method was developed for computing the free energy (delta Fi) of aggregates of type I collagen. The method was based on a treatment of Matheson and Flory describing phase equilibria of rigid rod polymers. It included a polymer-solvent interaction term that depended on near neighbor transfer energies. Extrahelical portions of the molecule were assigned local interaction energies differing from that assigned to the helix. Free energies of reaction for successive steps along assembly pathways (delta Fi-i+1) were computed. When allowance was made for specific pairing between extrahelical and helical domains, the so-called D-staggered (D = 670 A) alignment of molecules was preferred, as opposed to a nonstaggered, or nematic, alignment. Based on delta Fi-i+1 alone, it appeared that 1D-staggered oligomers arise first in assembly, followed later by addition of molecules in 4D alignment. Neither 4D dimers nor 4D-8D trimers were predicted to be major intermediates in assembly. This result is contrary to previous hypotheses. When energies of activation were included in the analysis, the prediction was less certain, and specific circumstances were identified in which 4D dimers and 4D-8D trimers were the earliest aggregated species in assembly.     

Interaction of collagen molecules from the aspect of fibril formation: acid-soluble, alkali-treated, and MMP1-digested fragments of type I collagen.

Collagen type I extracted with acid or digested with pepsin forms fibrils under physiological conditions, but this ability is lost when the collagen is treated with alkaline solution or digested with matrix metalloproteinase 1 (MMP1). When acid-soluble collagen was incubated with alkali-treated collagen, the fibril formation of acid-soluble collagen was inhibited. At 37 degrees C, at which alkali-treated collagen is denatured, the lag time was prolonged but the growth rate of fibrils was not affected. At 30 degrees C, at which the triple helical conformation of alkali-treated collagen is retained, the lag time was prolonged and the growth rate reduced. Heat-denatured alkali-treated collagen and MMP1-digested fragments have no inhibitory effect on the fibril formation of acid-soluble collagen. This means that the triple helical conformation and the molecular length are important factors in the interaction of collagen molecules and that alkali-treated collagen acts as a competitive inhibitor for fibril formation of collagen. We found that alkali-treated collagen and MMP1-digested fragments form fibrils that lack the D periodic banding pattern and twisted morphology under acidic conditions at the appropriate ionic strength. We also calculated the relative strengths of hydrophobic and electrostatic interactions between collagen molecules. When the hydrophobic interaction between linear collagen molecules was considered, we found a pattern of periodic maximization of the interactive force including the D period. On the other hand, the electrostatic interaction did not show the periodic pattern, but the overall interaction score affected fibril formation.     

Membrane binding of the colicin E1 channel: activity requires an electrostatic interaction of intermediate magnitude.

In vitro channel activity of the C-terminal colicin E1 channel polypeptide under conditions of variable electrostatic interaction with synthetic lipid membranes showed distinct maxima with respect to pH and membrane surface potential. The membrane binding energy was determined from fluorescence quenching of the intrinsic tryptophans of the channel polypeptide by liposomes containing N-trinitrophenyl-phosphatidylethanolamine. Maximum in vitro colicin channel activity correlated with an intermediate magnitude of the electrostatic interaction. For conditions associated with maximum activity (40% anionic lipid, I = 0.12 M, pH 4.0), the free energy of binding was delta G approximately -9 kcal/mol, with nonelectrostatic and electrostatic components, delta Gnel approximately -5 kcal/mol and delta Gel approximately -4 kcal/mol, and an effective binding charge of +7 at pH 4.0. Binding of the channel polypeptide to negative membranes at pH 8 is minimal, whereas initial binding at pH 4 followed by a shift to pH 8 causes only 3-10% reversal of binding, implying that it is kinetically trapped, probably by a hydrophobic interaction. It was inferred that membrane binding and insertion involves an initial electrostatic interaction responsible for concentration and binding to the membrane surface. This is followed by insertion into the bilayer driven by hydrophobic forces, which are countered in the case of excessive electrostatic binding.     

Information contained in the amino acid sequence of theα1(I)-chain of collagen and its consequences upon the formation of the triple helix,of fibrils and crosslinks

The molecule of type I collagen from skin consists of two alpha1(I)-chains and one alpha2-chain. The sequence of the entire alpha1-chain comprising 1052 residues is summarily presented and discussed. Apart from the 279 residues of alpha1(I)-CB8 whose sequence has been established for rat skin collagen, all sequences have been determined for calf skin collagen. In order to facilitate sequence analysis, the alpha1-chain was cleaved into defined fragments by cyanogen bromide or hydroxylamine or limited collagenase digestion. Most of the sequence was established by automated stepwise Edman degradation. The alpha1-chain contains two basically different types of sequences: the triple helical region of 1011 amino acid residues in which every third position is occupied by glycine and the N- and C-terminal regions not displaying this type of regularity. Both of these non-triple helical regions carry oxidizable lysine or hydroxylysine residues as functional sites for the intermolecular crosslink formation. Implications of the amino acid sequence for the stability of the triple helix and the fibril as well as for formation of crosslinks are discussed. Evaluation of the sequence in connection with electron microscopical investigations yielded the parameters of the axial arrangement of the molecules within the fibrils. Axial stagger of the molecules by a distance D = 670 angstrom = 233 amino acid residues results in maximal interaction of polar sequence regions of adjacent molecules and similarly of regions of hydrophobic residues. Ordered aggregation of molecules into fibrils is, therefore, regulated by electrostatic and electrophobic forces. Possible loci of intermolecular crosslinks between the alpha1-chains of adjacent molecules may be deduced from the dimensions of the axial aggregation of molecules.     

Role of the amino-terminal extrahelical region of type I collagen in directing the 4D overlap in fibrillogenesis

The amino-terminal telopeptide of the collagen α1(I) chain has a highly conserved sequence. This sequence was analyzed by the Chou-Fasman criteria, and a folded β-sheet conformation, including a β-turn, was predicted. This folded “hairpin” region favors both ionic and hydrophobic intermolecular interactions with α1(I) chain residues 930–938 on a neighboring, end-overlapped molecule. An end-overlap interaction of this nature could direct the initial step in fibril formation. The predicted structure also places the potential crosslink-forming lysyl residue, 9^N, in a unique site at the β-turn end of the telopeptide.     

On the stability of nucleic acid structures in solution: enthalpy-entropy compensations, internal rotations and reversibility.

The thermodynamics of self-association (stacking) of free bases and nucleotides, intramolecular stacking in dinucleotides, nearest-neighbour base pair stacking interactions in duplex DNA and RNA, and the formation of hairpin loops illustrate enthalpy/entropy compensations. Large stacking exothermicities are associated with large negative entropy changes that ensure that delta G is small, permitting readily reversible associations in solution. We rationalise enthalpy/entropy compensations with reference to residual motions and torsional vibrations which make a larger entropic contribution to binding when - delta H approximately kT (thermal energy at room temperature), than when - delta H >> kT. We present a factorisation of experimental free energies for helix formation in terms of approximate contributions from the restriction of rotations, hydrophobic interactions, electrostatic interactions due to base stacking, and contributions from hydrogen bonding, and estimate the adverse free energy cost per rotor (mainly entropy) of ordering the phosphate backbone as between 1.9 and 5.4 kJ mol-1 [averaged over 12 rotors per base pair for A-U on A-U stacking (lower limit), and G-C on C-G stacking (upper limit)]. The largest cost is associated with the most exothermic stacking interactions, while the range of values is consistent with earlier conclusions from data on the fusion of hydrocarbon chains (lower value), and with entropy changes in covalent isomerisations of small molecules involving severe restrictions (upper value).     

Collagen fibril formation in the presence of sodium dodecyl sulphate.

Sodium dodecyl sulphate (SDS) was used to weaken both the electrostatic and the hydrophobic interactions during collagen fibrillogenesis in vitro. The rate and extent of fibril formation as well as fibril morphology were affected by SDS concentration. Both the formation of large fibrils at 0.3 mM-SDS and the complete cessation of fibril formation at 0.5 mM-SDS were considered to be the result of SDS-induced conformational changes in the non-helical telopeptides. A possible mechanism of SDS interaction with the N-terminal and the distal region of the C-terminal telopeptides is offered.     

Calculations of antibody-antigen interactions: microscopic and semi-microscopic evaluation of the free energies of binding of phosphorylcholine analogs to McPC603.

The study of antibody-antigen interactions should greatly benefit from the development of quantitative models for the evaluation of binding free energies in proteins. The present work addresses this challenge by considering the test case of the binding free energies of phosphorylcholine analogs to the murine myeloma protein McPC603. This includes the evaluation of the differential binding energy as well as the absolute binding energies and their corresponding electrostatic contributions. Four different approaches are examined: the Protein Dipoles Langevin Dipoles (PDLD) method, the semi-microscopic PDLD (PDLD/S) method, a free energy perturbation (FEP) method based on an adiabatic charging procedure and a linear response approximation that accelerates the FEP calculation. The PDLD electrostatic calculations are augmented by estimates of the relevant hydrophobic and steric contributions. The determination of the hydrophobic energy involves an approach which considers the modification of the effective surface area of the solute by local field effects. The steric contributions are analyzed in terms of the corresponding reorganization energies. This treatment, which considers the protein as a harmonic system, views the steric forces as the restoring forces for the electrostatic interactions. The FEP method is found to give unreliable results with regular cut-off radii and starts to give quantitative results only in very expensive treatment with very large cut-off radii. The PDLD and PDLD/S methods are much faster than the FEP approach and give reasonable results for both the relative and absolute binding energies. The speed and simplicity of the PDLD/S method make it an effective strategy for interactive docking studies and indeed such an option is incorporated in the program MOLARIS. A component analysis of the different energy contributions of the FEP treatment and a similar PDLD analysis indicate that electrostatic effects provide the largest contribution to the differential binding energy, while the hydrophobic and steric contributions are much smaller. This finding lends further support to the idea that electrostatic interactions play a major role in determining the antigen specificity of McPC603.     

Semi-empirical conformational analysis of propranolol interacting with dipalmitoylphosphatidylcholine

A semi-empirical conformational analysis is used to compute the conformation of (+)-propranolol inserted in dipalmitoylphosphatidylcholine. In a first step, the minimal conformational energy of the isolated molecule at the hydrocarbon-water interface is calculated as the sum of the contributions resulting from the Van der Waals, the torsional, the electrostatic and the transfer energies. Five pairs of conformers of minimal energy are determined. They are compared to data available from other experimental approaches. In a second step, they are assembled with dipalmitoylphosphatidylcholine at the interface. Although propranolol is considered in its protonated form, the electrostatic interaction with dipalmitoylphosphatidylcholine is negligible as compared to the Van der Waals interaction. The area occupied per propranolol molecule is between 0.53 and 0.64 nm2/molecule. In the most probable modes of insertion of propranolol into the lipid layer, the naphthyl moiety of the compound interacts with the lipid acyl chains. The protonated amino group is located in the vicinity of the phosphate residue possibly causing an electrostatic interaction.     

Estimation of the binding force of the collagen molecule-decorin core protein complex in collagen fibril

Decorin belongs to the small leucine proteoglycans family and is considered to play an important role in extracellular matrix organization. Experimental studies suggest that decorin is required for the assembly of collagen fibrils, as well as for the development of proper tissue mechanical properties. In tendons, decorins tie adjoining collagen fibrils together and probably guarantee the mechanical coupling of fibrils. The decorin molecule consists of one core protein and one glycosaminoglycan chain covalently linked to a serine residue of the core protein. Several studies have indicated that each core protein binds to the surface of collagen fibrils every 67 nm, by interacting non-covalently to one collagen molecule of the fibril surface, while the decorin glycosaminoglycans extend from the core protein to connect to another decorin core protein laying on adjacent fibril surface. The present paper investigates the complex composed of one decorin core protein and one collagen molecule in order to obtain their binding force. For this purpose, molecular models of collagen molecules type I and decorin core protein were developed and their interaction energies were evaluated by means of the molecular mechanics approach. Results show that the complex is characterized by a maximum binding force of about 12.4 x 10(3) nN and a binding stiffness of 8.33 x 10(-8) N/nm; the attained binding force is greater than the glycosaminoglycan chain's ultimate strength, thus indicating that overloads are likely to damage the collagen fibre's mechanical integrity by disrupting the glycosaminoglycan chains rather than by causing decorin core protein detachment from the collagen fibril.     

Decorin Core Protein (Decoron) Shape Complements Collagen Fibril Surface Structure and Mediates Its Binding

Decorin is the archetypal small leucine rich repeat proteoglycan of the vertebrate extracellular matrix (ECM). With its glycosaminoglycuronan chain, it is responsible for stabilizing inter-fibrillar organization. Type I collagen is the predominant member of the fibrillar collagen family, fulfilling both organizational and structural roles in animal ECMs. In this study, interactions between decoron (the decorin core protein) and binding sites in the d and e₁ bands of the type I collagen fibril were investigated through molecular modeling of their respective X-ray diffraction structures. Previously, it was proposed that a model-based, highly curved concave decoron interacts with a single collagen molecule, which would form extensive van der Waals contacts and give rise to strong non-specific binding. However, the large well-ordered aggregate that is the collagen fibril places significant restraints on modes of ligand binding and necessitates multi-collagen molecular contacts. We present here a relatively high-resolution model of the decoron-fibril collagen complex. We find that the respective crystal structures complement each other well, although it is the monomeric form of decoron that shows the most appropriate shape complementarity with the fibril surface and favorable calculated energies of interaction. One molecule of decoron interacts with four to six collagen molecules, and the binding specificity relies on a large number of hydrogen bonds and electrostatic interactions, primarily with the collagen motifs KXGDRGE and AKGDRGE (d and e₁ bands). This work helps us to understand collagen-decorin interactions and the molecular architecture of the fibrillar ECM in health and disease.     

A combined bioinformatic approach oriented to the analysis and design of peptides with high affinity to MHC class I molecules

We report on a new method to compute the antigenic degree of peptides from available experimental data on peptide binding affinity to class I MHC molecules. The methodology is a combination of two strategies at different levels of information. The first, at the primary structure level, consists in expressing the peptides binding activity as a profile of amino acid contributions, amino acid similarity being accounted for by their characteristic physicochemical properties and their position within the sequence. The higher level of the strategy is based on a meticulous analysis of the contact interface of the peptides with the cleft constituting the receptor region of a particular class I MHC molecule. Interaction interfaces are inferred by docking the peptide onto the receptor groove of the MHC molecule; evaluation of the affinity of the peptide to the receptor is then performed by analysis of the electrostatic and hydrophobic energies on points of the interaction interface.The result is a robust system for analysis of peptide affinity to class I MHC molecules since while the first analysis dictates the composition of active sequences at the amino acid level, the second translates this information to the atomic level, where the molecular interaction can be analyzed in terms of the intrinsic interatomic forces and energies. Evaluation results for the methodology are encouraging since high affinity peptides are reflected by high scores at both levels of information, and are proportionally lower for peptides of medium and lower affinity for which interaction surfaces show relatively lower electrostatic complementarity and hydrophobic correlation than for the former.     

Empirical calculation of the relative free energies of peptide binding to the molecular chaperone DnaK

We describe a methodology to calculate the relative free energies of protein-peptide complex formation. The interaction energy was decomposed into nonpolar, electrostatic and entropic contributions. A free energy-surface area relationship served to calculate the nonpolar free energy term. The electrostatic free energy was calculated with the finite difference Poisson-Boltzmann method and the entropic contribution was estimated from the loss in the conformational entropy of the peptide side chains. We applied this methodology to a series of DnaK*peptide complexes. On the basis of the single known crystal structure of the peptide-binding domain of DnaK with a bound heptapeptide, we modeled ten other DnaK*heptapeptide complexes with experimentally measured K(d) values from 0.06 microM to 11 microM, using molecular dynamics to refine the structures of the complexes. Molecular dynamic trajectories, after equilibration, were used for calculating the energies with greater accuracy. The calculated relative binding free energies were compared with the experimentally determined free energies. Linear scaling of the calculated terms was applied to fit them to the experimental values. The calculated binding free energies were between -7.1 kcal/mol and - 9.4 kcal/mol with a correlation coefficient of 0.86. The calculated nonpolar contributions are mainly due to the central hydrophobic binding pocket of DnaK for three amino acid residues. Negative electrostatic fields generated by the protein increase the binding affinity for basic residues flanking the hydrophobic core of the peptide ligand. Analysis of the individual energy contributions indicated that the nonpolar contributions are predominant compared to the other energy terms even for peptides with low affinity and that inclusion of the change in conformational entropy of the peptide side chains does not improve the discriminative power of the calculation. The method seems to be useful for predicting relative binding energies of peptide ligands of DnaK and might be applicable to other protein-peptide systems, particularly if only the structure of one protein-ligand complex is available.     

Probing the initial stage of aggregation of the Abeta(10-35)-protein: assessing the propensity for peptide dimerization

Characterization of the early stages of peptide aggregation is of fundamental importance in elucidating the mechanism of the formation of deposits associated with amyloid disease. The initial step in the pathway of aggregation of the Abeta-protein, whose monomeric NMR structure is known, was studied through the simulation of the structure and stability of the peptide dimer in aqueous solution. A protocol based on shape complementarity was used to generate an assortment of possible dimer structures. The structures generated based on shape complementarity were evaluated using rapidly computed estimates of the desolvation and electrostatic interaction energies to identify a putative stable dimer structure. The potential of mean force associated with the dimerization of the peptides in aqueous solution was computed for both the hydrophobic and the electrostatic driven forces using umbrella sampling and classical molecular dynamics simulation at constant temperature and pressure with explicit solvent and periodic boundary conditions. The comparison of the two free energy profiles suggests that the structure of the peptide dimer is determined by the favorable desolvation of the hydrophobic residues at the interface. Molecular dynamics trajectories originating from two putative dimer structures indicate that the peptide dimer is stabilized primarily through hydrophobic interactions, while the conformations of the peptide monomers undergo substantial structural reorganization in the dimerization process. The finding that the phi-dimer may constitute the ensemble of stable Abeta(10-35) dimer has important implications for fibril formation. In particular, the expulsion of water molecules at the interface might be a key event, just as in the oligomerization of Abeta(16-22) fragments. We conjecture that events prior to the nucleation process themselves might involve crossing free energy barriers which depend on the peptide-peptide and peptide-water interactions. Consistent with existing experimental studies, the peptides within the ensemble of aggregated states show no signs of formation of secondary structure.     

Solution nonideality related to solute molecular characteristics of amino acids.

By measuring the freezing-point depression for dilute, aqueous solutions of all water-soluble amino acids, we test the hypothesis that nonideality in aqueous solutions is due to solute-induced water structuring near hydrophobic surfaces and solute-induced water destructuring in the dipolar electric fields generated by the solute. Nonideality is expressed with a single solute/solvent interaction parameter I, calculated from experimental measure of delta T. A related parameter, I(n), gives a method of directly relating solute characteristics to solute-induced water structuring or destructuring. I(n)-values correlate directly with hydrophobic surface area and inversely with dipolar strength. By comparing the nonideality of amino acids with progressively larger hydrophobic side chains, structuring is shown to increase with hydrophobic surface area at a rate of one perturbed water molecule per 8.8 square angstroms, implying monolayer coverage. Destructuring is attributed to dielectric realignment as described by the Debye-Hückel theory, but with a constant separation of charges in the amino-carboxyl dipole. By using dimers and trimers of glycine and alanine, this destructuring is shown to increase with increasing dipole strength using increased separation of fixed dipolar charges. The capacity to predict nonideal solution behavior on the basis of amino acid characteristics will permit prediction of free energy of transfer to water, which may help predict the energetics of folding and unfolding of proteins based on the characteristics of constituent amino acids.     

Predicting the rate enhancement of protein complex formation from the electrostatic energy of interaction

The rate of association of proteins is dictated by diffusion, but can be enhanced by favorable electrostatic forces. Here the relationship between the electrostatic energy of interaction, and the kinetics of protein-complex formation was analyzed for the protein pairs of: hirudin-thrombin, acetylcholinesterase-fasciculin and barnase-barstar, and for a panel of point mutants of these proteins. Electrostatic energies of interaction were calculated as the difference between the electrostatic energy of the complex and the sum of the energies of the two individual proteins, using the computer simulation package DelPhi. Calculated electrostatic energies of interaction were compared to experimentally determined rates of association. One kcal/mol of Coulombic interaction energy increased the rate of association by a factor of 2.8, independent of the protein-complex or mutant analyzed. Electrostatic energies of interaction were also determined from the salt dependence of the association rate constant, using the same basic equation as for the theoretical calculation. A Br?nsted analysis of the electrostatic energies of interactions plotted versus experimentally determined ln(rate)s of association shows a linear relation between the two, with a beta value close to 1. This is interpreted as the energy of the transition state varies according to the electrostatic interaction energy, fitting a two state model for the association reaction. Calculating electrostatic rate enhancement from the electrostatic interaction energy can be used as a powerful tool to design protein complexes with altered rates of association and affinities.     

Unhydroxylated triple helical collagen I produced in transgenic plants provides new clues on the role of hydroxyproline in collagen folding and fibril formation

Human unhydroxylated homotrimeric triple-helical collagen I produced in transgenic plants was used as an experimental model to provide insights into the role of hydroxyproline in molecular folding and fibril formation. By using chemically cross-linked molecules, we show here that the absence of hydroxyproline residues does not prevent correct folding of the recombinant collagen although it markedly slows down the propagation rate compared with bovine fully hydroxylated homotrimeric collagen I. Relatively slow cis-trans-isomerization in the absence of hydroxyproline likely represents the rate-limiting factor in the propagation of the unhydroxylated collagen helix. Because of the lack of hydroxylation, recombinant collagen molecules showed increased flexibility as well as a reduced melting temperature compared with native homotrimers and heterotrimers, whereas the distribution of charged amino acids was unchanged. However, unlike with bovine collagen I, the recombinant collagen did not self-assemble into banded fibrils in physiological ionic strength buffer at 20 degrees C. Striated fibrils were only obtained with low ionic strength buffer. We propose that, under physiological ionic strength conditions, the hydroxyl groups in the native molecule retain water more efficiently thus favoring correct fibril formation. The importance of hydroxyproline in collagen self-assembly suggested by others from the crystal structures of collagen model peptides is thus confirmed experimentally on the entire collagen molecule.