Synthetic NAC 71-82 Peptides Designed to Produce Fibrils with Different Protofilament Interface Contacts

Alpha-synucleinopathies are featured by fibrillar inclusions in brain cells. Although α-synuclein fibrils display structural diversity, the origin of this diversity is not fully understood. We used molecular dynamics simulations to design synthetic peptides, based on the NAC 71-82 amino acid fragment of α-synuclein, that govern protofilament contacts and generation of twisted fibrillar polymorphs. Four peptides with structures based on either single or double fragments and capped or non-capped ends were selected for further analysis. We determined the fibrillar yield and the structures from these peptides found in the solution after fibrillisation using protein concentration determination assay and circular dichroism spectroscopy. In addition, we characterised secondary structures formed by individual fibrillar complexes using laser-tweezers Raman spectroscopy. Results suggest less mature fibrils, based on the lower relative β-sheet content for double- than single-fragment peptide fibrils. We confirmed this structural difference by TEM analysis which revealed, in addition to short protofibrils, more elongated, twisted and rod-like fibril structures in non-capped and capped double-fragment peptide systems, respectively. Finally, time-correlated single-photon counting demonstrated a difference in the Thioflavin T fluorescence lifetime profiles upon fibril binding. It could be proposed that this difference originated from morphological differences in the fibril samples. Altogether, these results highlight the potential of using peptide models for the generation of fibrils that share morphological features relevant for disease, e.g., twisted and rod-like polymorphs.     

Role of water in protein aggregation and amyloid polymorphism

A variety of neurodegenerative diseases are associated with amyloid plaques, which begin as soluble protein oligomers but develop into amyloid fibrils. Our incomplete understanding of this process underscores the need to decipher the principles governing protein aggregation. Mechanisms of in vivo amyloid formation involve a number of coconspirators and complex interactions with membranes. Nevertheless, understanding the biophysical basis of simpler in vitro amyloid formation is considered important for discovering ligands that preferentially bind regions harboring amyloidogenic tendencies. The determination of the fibril structure of many peptides has set the stage for probing the dynamics of oligomer formation and amyloid growth through computer simulations. Most experimental and simulation studies, however, have been interpreted largely from the perspective of proteins: the role of solvent has been relatively overlooked in oligomer formation and assembly to protofilaments and amyloid fibrils. In this Account, we provide a perspective on how interactions with water affect folding landscapes of amyloid beta (Aβ) monomers, oligomer formation in the Aβ16-22 fragment, and protofilament formation in a peptide from yeast prion Sup35. Explicit molecular dynamics simulations illustrate how water controls the self-assembly of higher order structures, providing a structural basis for understanding the kinetics of oligomer and fibril growth. Simulations show that monomers of Aβ peptides sample a number of compact conformations. The formation of aggregation-prone structures (N*) with a salt bridge, strikingly similar to the structure in the fibril, requires overcoming a high desolvation barrier. In general, sequences for which N* structures are not significantly populated are unlikely to aggregate. Oligomers and fibrils generally form in two steps. First, water is expelled from the region between peptides rich in hydrophobic residues (for example, Aβ16-22), resulting in disordered oligomers. Then the peptides align along a preferred axis to form ordered structures with anti-parallel β-strand arrangement. The rate-limiting step in the ordered assembly is the rearrangement of the peptides within a confining volume. The mechanism of protofilament formation in a polar peptide fragment from the yeast prion, in which the two sheets are packed against each other and create a dry interface, illustrates that water dramatically slows self-assembly. As the sheets approach each other, two perfectly ordered one-dimensional water wires form. They are stabilized by hydrogen bonds to the amide groups of the polar side chains, resulting in the formation of long-lived metastable structures. Release of trapped water from the pore creates a helically twisted protofilament with a dry interface. Similarly, the driving force for addition of a solvated monomer to a preformed fibril is water release; the entropy gain and favorable interpeptide hydrogen bond formation compensate for entropy loss in the peptides. We conclude by offering evidence that a two-step model, similar to that postulated for protein crystallization, must also hold for higher order amyloid structure formation starting from N*. Distinct water-laden polymorphic structures result from multiple N* structures. Water plays multifarious roles in all of these protein aggregations. In predominantly hydrophobic sequences, water accelerates fibril formation. In contrast, water-stabilized metastable intermediates dramatically slow fibril growth rates in hydrophilic sequences.     

Phosphorylation as a tool to modulate aggregation propensity and to predict fibril architecture

Despite the importance of post-translational modifications in controlling the solubility and conformational properties of proteins and peptides, precisely how the aggregation propensity of different peptide sequences is modulated by chemical modification remains unclear. Here we have investigated the effect of phosphorylation on the aggregation propensity of a 13-residue synthetic peptide incorporating one or more phosphate groups at seven different sites at various pH values. Fibril formation was shown to be inhibited when a single phosphate group was introduced at all seven locations in the peptide sequence at pH 7.5, when the phosphate group is fully charged. By contrast, when the same peptides were analysed at pH 1.1, when the phosphate is fully protonated, fibrils from all seven peptide sequences form rapidly. At intermediate pH values (pH 3.6) when the phosphate group is mono-anionic, the aggregation propensity of the peptides was found to be highly dependent on the position of the phosphate group in the peptide sequence. Using this information, combined with molecular dynamics (MD) simulations of the peptide sequence, we provide evidence consistent with the peptide forming amyloid fibrils with a class 7 architecture. The results highlight the potential utility of phosphorylation as a method of reversibly controlling the aggregation kinetics of peptide sequences both during and after synthesis. Moreover, by exploiting the ability of the phosphate group to adopt different charge states as a function of pH, and combining experimental insights with atomistic information calculated from MD simulations as pH is varied, we show how the resulting information can be used to predict fibril structures consistent with both datasets, and use these to rationalise their sensitivity of fibrillation kinetics both to the location of the phosphate group and its charge state.     

Aggregation Condition–Structure Relationship of Mouse Prion Protein Fibrils

Prion diseases are associated with conformational conversion of cellular prion protein into a misfolded pathogenic form, which resembles many properties of amyloid fibrils. The same prion protein sequence can misfold into different conformations, which are responsible for variations in prion disease phenotypes (prion strains). In this work, we use atomic force microscopy, FTIR spectroscopy and magic-angle spinning NMR to devise structural models of mouse prion protein fibrils prepared in three different denaturing conditions. We find that the fibril core region as well as the structure of its N- and C-terminal parts is almost identical between the three fibrils. In contrast, the central part differs in length of β-strands and the arrangement of charged residues. We propose that the denaturant ionic strength plays a major role in determining the structure of fibrils obtained in a particular condition by stabilizing fibril core interior-facing glutamic acid residues.     

Self-folding and aggregation of amyloid nanofibrils

Amyloids are highly organized protein filaments, rich in β-sheet secondary structures that self-assemble to form dense plaques in brain tissues affected by severe neurodegenerative disorders (e.g. Alzheimer's Disease). Identified as natural functional materials in bacteria, in addition to their remarkable mechanical properties, amyloids have also been proposed as a platform for novel biomaterials in nanotechnology applications including nanowires, liquid crystals, scaffolds and thin films. Despite recent progress in understanding amyloid structure and behavior, the latent self-assembly mechanism and the underlying adhesion forces that drive the aggregation process remain poorly understood. On the basis of previous full atomistic simulations, here we report a simple coarse-grain model to analyze the competition between adhesive forces and elastic deformation of amyloid fibrils. We use simple model system to investigate self-assembly mechanisms of fibrils, focused on the formation of self-folded nanorackets and nanorings, and thereby address a critical issue in linking the biochemical (Angstrom) to micrometre scales relevant for larger-scale states of functional amyloid materials. We investigate the effect of varying the interfibril adhesion energy on the structure and stability of self-folded nanorackets and nanorings and demonstrate that these aggregated amyloid fibrils are stable in such states even when the fibril-fibril interaction is relatively weak, given that the constituting amyloid fibril length exceeds a critical fibril length-scale of several hundred nanometres. We further present a simple approach to directly determine the interfibril adhesion strength from geometric measures. In addition to providing insight into the physics of aggregation of amyloid fibrils our model enables the analysis of large-scale amyloid plaques and presents a new method for the estimation and engineering of the adhesive forces responsible of the self-assembly process of amyloid nanostructures, filling a gap that previously existed between full atomistic simulations of primarily ultra-short fibrils and much larger micrometre-scale amyloid aggregates. Via direct simulation of large-scale amyloid aggregates consisting of hundreds of fibrils we demonstrate that the fibril length has a profound impact on their structure and mechanical properties, where the critical fibril length-scale derived from our analysis of self-folded nanorackets and nanorings defines the structure of amyloid aggregates. A multi-scale modeling approach as used here, bridging the scales from Angstroms to micrometres, opens a wide range of possible nanotechnology applications by presenting a holistic framework that balances mechanical properties of individual fibrils, hierarchical self-assembly, and the adhesive forces determining their stability to facilitate the design of de novo amyloid materials.     

Effects of Intramolecular Distance between Amyloidogenic Domains on Amyloid Aggregation

Peptide/protein aggregation is implicated in many amyloid diseases. Some amyloidogenic peptides/proteins, such as those implicated in Alzheimer’s and Parkinson’s diseases, contain multiple amyloidogenic domains connected by “linker” sequences displaying high propensities to form turn structures. Recent studies have demonstrated the importance of physicochemical properties of each amino acid contained in the polypeptide sequences in amyloid aggregation. However, effects on aggregation related to the intramolecular distance between amyloidogenic domains, which may be determined by a linker length, have yet to be examined. In the study presented here, we created peptides containing two copies of KFFE, a simple four-residue amyloidogenic domain, connected by GS-rich linker sequences with different lengths yet similar physicochemical properties. Our experimental results indicate that aggregation occurred most rapidly when KFFE domains were connected by a linker of an intermediate length. Our experimental findings were consistent with estimated entropic contribution of a linker length toward formation of (partially) structured intermediates on the aggregation pathway. Moreover, inclusion of a relatively short linker was found to inhibit formation of aggregates with mature fibril morphology. When the results are assimilated, our study demonstrates that intramolecular distance between amyloidogenic domains is an important yet overlooked factor affecting amyloid aggregation.     

Searching for the Best Transthyretin Aggregation Protocol to Study Amyloid Fibril Disruption

Several degenerative amyloid diseases, with no fully effective treatment, affect millions of people worldwide. These pathologies—amyloidoses—are known to be associated with the formation of ordered protein aggregates and highly stable and insoluble amyloid fibrils, which are deposited in multiple tissues and organs. The disruption of preformed amyloid aggregates and fibrils is one possible therapeutic strategy against amyloidosis; however, only a few compounds have been identified as possible fibril disruptors in vivo to date. To properly identify chemical compounds as potential fibril disruptors, a reliable, fast, and economic screening protocol must be developed. For this purpose, three amyloid fibril formation protocols using transthyretin (TTR), a plasma protein involved in several amyloidoses, were studied using thioflavin-T fluorescence assays, circular dichroism (CD), turbidity, dynamic light scattering (DLS), and transmission electron microscopy (TEM), in order to characterize and select the most appropriate fibril formation protocol. Saturation transfer difference nuclear magnetic resonance spectroscopy (STD NMR) was successfully used to study the interaction of doxycycline, a known amyloid fibril disruptor, with preformed wild-type TTR (TTRwt) aggregates and fibrils. DLS and TEM were also used to characterize the effect of doxycycline on TTRwt amyloid species disaggregation. A comparison of the TTR amyloid morphology formed in different experimental conditions is also presented.     

Temperature-Dependent Structural Variability of Prion Protein Amyloid Fibrils

Prion protein aggregation into amyloid fibrils is associated with the onset and progression of prion diseases—a group of neurodegenerative amyloidoses. The process of such aggregate formation is still not fully understood, especially regarding their polymorphism, an event where the same type of protein forms multiple, conformationally and morphologically distinct structures. Considering that such structural variations can greatly complicate the search for potential antiamyloid compounds, either by having specific propagation properties or stability, it is important to better understand this aggregation event. We have recently reported the ability of prion protein fibrils to obtain at least two distinct conformations under identical conditions, which raised the question if this occurrence is tied to only certain environmental conditions. In this work, we examined a large sample size of prion protein aggregation reactions under a range of temperatures and analyzed the resulting fibril dye-binding, secondary structure and morphological properties. We show that all temperature conditions lead to the formation of more than one fibril type and that this variability may depend on the state of the initial prion protein molecules.     

Effect of 1-Ethyl-3-methylimidazolium Tetrafluoroborate and Acetate Ionic Liquids on Stability and Amyloid Aggregation of Lysozyme

Amyloid fibrils draw attention as potential novel biomaterials due to their high stability, strength, elasticity or resistance against degradation. Therefore, the controlled and fast fibrillization process is of great interest, which raises the demand for effective tools capable of regulating amyloid fibrillization. Ionic liquids (ILs) were identified as effective modulators of amyloid aggregation. The present work is focused on the study of the effect of 1-ethyl-3-methyl imidazolium-based ILs with kosmotropic anion acetate (EMIM-ac) and chaotropic cation tetrafluoroborate (EMIM-BF₄) on the kinetics of lysozyme amyloid aggregation and morphology of formed fibrils using fluorescence and CD spectroscopy, differential scanning calorimetry, AFM with statistical image analysis and docking calculations. We have found that both ILs decrease the thermal stability of lysozyme and significantly accelerate amyloid fibrillization in a dose-dependent manner at concentrations of 0.5%, 1% and 5% (v/v) in conditions and time-frames when no fibrils are formed in ILs-free solvent. The effect of EMIM-BF₄ is more prominent than EMIM-ac due to the different specific interactions of the anionic part with the protein surface. Although both ILs induced formation of amyloid fibrils with typical needle-like morphology, a higher variability of fibril morphology consisting of a different number of intertwining protofilaments was identified for EMIM-BF_4.     

Molecular Simulations of Amyloid Structures,Toxicity, and Inhibition

The misfolding and aggregation of proteins and peptides into amyloid fibrils are believed to be responsible for the dysfunction and death of neuron cells in many neurodegenerative diseases. Resolving the atomic structures of amyloid peptides at different aggregation stages by molecular simulations has opened new ways to probe the molecular mechanisms of amyloid aggregation, toxicity, and inhibition, as well as to validate computational data with available experimental ones. In this review article, we summarize some recent and important findings on: 1) a number of atomic structures of amyloid oligomers with typical β-sheet-rich conformations, related to amyloid aggregation; 2) different amyloid peptide-induced membrane-disruption mechanisms, related to amyloid toxicity; and 3) rational design of different amyloid inhibitors capable of preventing amyloid aggregation and toxicity, related to amyloid inhibition. All these findings will provide some mechanistic implications for molecular mechanisms of amyloid aggregation, toxicity, and inhibition, which are fundamentally and practically important for the treatment of amyloid diseases.     

Amyloid architecture: complementary assembly of heterogeneous combinations of three or four peptides into amyloid fibrils

The amyloid fibril is a misfolded and undesirable state for proteins that has been proposed to be a causative agent for a variety of fatal diseases known as amyloid diseases, such as Alzheimer's and prion diseases. However, the fibril has a highly ordered tertiary structure in which numerous beta-strand polypeptide chains align in a regular pattern. Thus, this kind of fibril has the potential to be engineered into proteinaceous materials. Amyloid fibrils of misfolded proteins primarily comprise a single polypeptide species, that is, the self-assembly is homogeneous. We here found that three or four designed peptides can assemble heterogeneously and cooperatively into amyloid fibrils, a process accompanied by a drastic secondary structural transition from alpha helix to beta sheet. Heterogeneous assembly into fibrils is accomplished by complementary electrostatic interactions between three or four peptide species, each of which is not able to self-assemble homogeneously. These findings will lead to a novel way to study the molecular details of amyloid formation and also to design beta-sheet peptidyl materials.     

Influence of Cortisol on the Fibril Formation Kinetics of Aβ42 Peptide: A Multi-Technical Approach

Amyloid-β peptide (Aβ) aggregates are known to be correlated with pathological neurodegenerative diseases. The fibril formation process of such peptides in solution is influenced by several factors, such as the ionic strength of the buffer, concentration, pH, and presence of other molecules, just to mention a few. In this paper, we report a detailed analysis of in vitro Aβ42 fibril formation in the presence of cortisol at different relative concentrations. The thioflavin T fluorescence assay allowed us to monitor the fibril formation kinetics, while a morphological characterization of the aggregates was obtained by atomic force microscopy. Moreover, infrared absorption spectroscopy was exploited to investigate the secondary structure changes along the fibril formation path. Molecular dynamics calculations allowed us to understand the intermolecular interactions with cortisol. The combined results demonstrated the influence of cortisol on the fibril formation process: indeed, at cortisol-Aβ42 concentration ratio (ρ) close to 0.1 a faster organization of Aβ42 fragments into fibrils is promoted, while for ρ = 1 the formation of fibrils is completely inhibited.     

In Silico Screening of Novel α1-GABAA Receptor PAMs towards Schizophrenia Based on Combined Modeling Studies of Imidazo [1,2-a]-Pyridines

The ionotropic GABA_A receptor (GABA_AR) has been proven to be an important target of atypical antipsychotics. A novel series of imidazo [1,2-a]-pyridine derivatives, as selective positive allosteric modulators (PAMs) of α1-containing GABA_ARs with potent antipsychotic activities, have been reported recently. To better clarify the pharmacological essentiality of these PAMs and explore novel antipsychotics hits, three-dimensional quantitative structure–activity relationships (3D-QSAR), molecular docking, pharmacophore modeling, and molecular dynamics (MD) were performed on 33 imidazo [1,2-a]-pyridines. The constructed 3D-QSAR models exhibited good predictive abilities. The dockings results and MD simulations demonstrated that hydrogen bonds, π–π stackings, and hydrophobic interactions play essential roles in the binding of these novel PAMs in the GABA_AR binding pocket. Four hit compounds (DS01–04) were then screened out by the combination of the constructed models and computations, including the pharmacophore model, Topomer Search, molecular dockings, ADME/T predictions, and MD simulations. The compounds DS03 and DS04, with higher docking scores and better predicted activities, were also found to be relatively stable in the binding pocket by MD simulations. These results might provide a significant theoretical direction or information for the rational design and development of novel α1-GABA_AR PAMs with antipsychotic activities.     

The Effect of Octapeptide Repeats on Prion Folding and Misfolding

Misfolding of prion protein (PrP) into amyloid aggregates is the central feature of prion diseases. PrP has an amyloidogenic C-terminal domain with three α-helices and a flexible tail in the N-terminal domain in which multiple octapeptide repeats are present in most mammals. The role of the octapeptides in prion diseases has previously been underestimated because the octapeptides are not located in the amyloidogenic domain. Correlation between the number of octapeptide repeats and age of onset suggests the critical role of octapeptide repeats in prion diseases. In this study, we have investigated four PrP variants without any octapeptides and with 1, 5 and 8 octapeptide repeats. From the comparison of the protein structure and the thermal stability of these proteins, as well as the characterization of amyloids converted from these PrP variants, we found that octapeptide repeats affect both folding and misfolding of PrP creating amyloid fibrils with distinct structures. Deletion of octapeptides forms fewer twisted fibrils and weakens the cytotoxicity. Insertion of octapeptides enhances the formation of typical silk-like fibrils but it does not increase the cytotoxicity. There might be some threshold effect and increasing the number of peptides beyond a certain limit has no further effect on the cell viability, though the reasons are unclear at this stage. Overall, the results of this study elucidate the molecular mechanism of octapeptides at the onset of prion diseases.     

Heredity and difference of multiple-scale microstructures in PAN-based carbon fibers and their precursor fibers

Multiple-scale microstructures, including skin-core structure, fibril structure, lamellar structure, crystal/amorphous structure, were found co-present in the fibers during the whole production process of polyacrylonitrile (PAN)-based carbon fibers. The structural heredity and difference among them were systematically investigated for the first time by scanning electron microscope, optical microscope, transmission electron microscope, and X-ray diffraction. The relations between the four kinds of structures and their formation mechanisms were analyzed. The skin-core structure is contributed to inhomogeneous distribution of composition and chemical structure along the radial direction in oxidative stabilized fibers, which is proved different in PAN precursor fibers and carbon fibers. Fibrils are successively separated from low-temperature oxidative stabilized fibers by ultrasonic etching in dimethylsulfoxide solution. The separation of individual fibril becomes harder and even impossible in those fibers prepared at temperatures higher than 245°C. This suggests a stronger bonding force between fibrils in high-temperature oxidative stabilized fibers and carbon fibers. The lamellar structures within fibrils are observed in all of these fibers but with thicker lamella width with increasing temperature. They are unlikely due to the alternatively alignment of crystal regions and amorphous regions as reported by many previous literature, because the oxidative stabilized fibers are amorphous but have lamellar structures. The (002) diffraction arc gives the evidence that the lamellar structure in carbon fibers is not strictly perpendicular to the fiber axis, but have an angle of about 45° with it. © 2012 Wiley Periodicals, Inc. J Appl Polym Sci, 2012     

Methylcellulose fibrils: a mini review

This review covers the recently discovered methylcellulose (MC) fibrils and their impact on MC thermogelation and overall gel strength. The thermogelation of MC aqueous solutions has been studied for nearly 100 years, but it has only recently been attributed to fibril formation at elevated temperatures. The assembly of a fibrillar network explains the general features of MC gels: the unusually high modulus, the turbidity, and the nonlinear rheological behavior at high strains. Recent scattering, atomic force microscopy (AFM) and cryogenic transmission electron microscopy (cryo‐TEM) experimental efforts have led to critical discoveries about the fibril structure. MC fibrils are 60% water by volume when in solution. The length of the fibrils, especially for shorter chains, is closely correlated with the MC chain contour length. On the other hand, the fibril radius does not depend on the temperature, molecular weight or concentration of MC chains but does seem to be correlated with chain stiffness. The discovery of MC fibrils represents a significant departure from the entanglement of polymer chains and localized phase separation based physical model of MC gelation, opens new questions about the mechanism of fibril formation and presents new pathways for stimulus responsive material design. © 2019 Society of Chemical Industry     

Structural insights into peptide self-assembly using photo-induced crosslinking experiments and discontinuous molecular dynamics

Determining the structure of the (oligomeric) intermediates that form during the self-assembly of amyloidogenic peptides is challenging because of their heterogeneous and dynamic nature. Thus, there is need for methodology to analyze the underlying molecular structure of these transient species. In this work, a combination of fluorescence quenching, photo-induced crosslinking (PIC) and molecular dynamics simulation was used to study the assembly of a synthetic amyloid-forming peptide, Aβ_16-22. A PIC amino acid containing a trifluormethyldiazirine (TFMD) group—Fmoc(TFMD)Phe—was incorporated into the sequence (Aβ*_16–22). Electrospray ionization ion-mobility spectrometry mass-spectrometry (ESI-IMS-MS) analysis of the PIC products confirmed that Aβ*_16–22 forms assemblies with the monomers arranged as anti-parallel, in-register β-strands at all time points during the aggregation assay. The assembly process was also monitored separately using fluorescence quenching to profile the fibril assembly reaction. The molecular picture resulting from discontinuous molecule dynamics simulations showed that Aβ_16-22 assembles through a single-step nucleation into a β-sheet fibril in agreement with these experimental observations. This study provides detailed structural insights into the Aβ_16-22 self-assembly processes, paving the way to explore the self-assembly mechanism of larger, more complex peptides, including those whose aggregation is responsible for human disease.     

Mapping the Binding Interface of PET Tracer Molecules and Alzheimer Disease Aβ Fibrils by Using MAS Solid-State NMR Spectroscopy

Positron emission tomography (PET) tracer molecules like thioflavin T specifically recognize amyloid deposition in brain tissue by selective binding to hydrophobic or aromatic surface grooves on the β-sheet surface along the fibril axis. The molecular basis of this interaction is, however, not well understood. We have employed magic angle spinning (MAS) solid-state NMR spectroscopy to characterize Aβ-PET tracer complexes at atomic resolution. We established a titration protocol by using bovine serum albumin as a carrier to transfer hydrophobic small molecules to Aβ(1-40) fibrillar aggregates. The same Aβ(1-40) amyloid fibril sample was employed in subsequent titrations to minimize systematic errors that potentially arise from sample preparation. In the experiments, the small molecules ¹³C-methylated Pittsburgh compound B (PiB) as well as a novel Aβ tracer based on a diarylbithiazole (DABTA) scaffold were employed. Classical ¹³C-detected as well as proton-detected spectra of protonated and perdeuterated samples with back-substituted protons, respectively, were acquired and analyzed. After titration of the tracers, chemical-shift perturbations were observed in the loop region involving residues Gly25-Lys28 and Ile32-Gly33, thus suggesting that the PET tracer molecules interact with the loop region connecting β-sheets β1 and β2 in Aβ fibrils. We found that titration of the PiB derivatives suppressed fibril polymorphism and stabilized the amyloid fibril structure.     

Contrasting Local and Macroscopic Effects of Collagen Hydroxylation

Collagen is heavily hydroxylated. Experiments show that proline hydroxylation is important to triple helix (monomer) stability, fibril assembly, and interaction of fibrils with other molecules. Nevertheless, experiments also show that even without hydroxylation, type I collagen does assemble into its native D-banded fibrillar structure. This raises two questions. Firstly, even though hydroxylation removal marginally affects macroscopic structure, how does such an extensive chemical change, which is expected to substantially reduce hydrogen bonding capacity, affect local structure? Secondly, how does such a chemical perturbation, which is expected to substantially decrease electrostatic attraction between monomers, affect collagen’s mechanical properties? To address these issues, we conduct a benchmarked molecular dynamics study of rat type I fibrils in the presence and absence of hydroxylation. Our simulations reproduce the experimental observation that hydroxylation removal has a minimal effect on collagen’s D-band length. We also find that the gap-overlap ratio, monomer width and monomer length are minimally affected. Surprisingly, we find that de-hydroxylation also has a minor effect on the fibril’s Young’s modulus, and elastic stress build up is also accompanied by tightening of triple-helix windings. In terms of local structure, de-hydroxylation does result in a substantial drop (23%) in inter-monomer hydrogen bonding. However, at the same time, the local structures and inter-monomer hydrogen bonding networks of non-hydroxylated amino acids are also affected. It seems that it is this intrinsic plasticity in inter-monomer interactions that preclude fibrils from undergoing any large changes in macroscopic properties. Nevertheless, changes in local structure can be expected to directly impact collagen’s interaction with extra-cellular matrix proteins. In general, this study highlights a key challenge in tissue engineering and medicine related to mapping collagen chemistry to macroscopic properties but suggests a path forward to address it using molecular dynamics simulations.     

Methods and Applications of In Silico Aptamer Design and Modeling

Aptamers are nucleic acid analogues of antibodies with high affinity to different targets, such as cells, viruses, proteins, inorganic materials, and coenzymes. Empirical approaches allow the design of in vitro aptamers that bind particularly to a target molecule with high affinity and selectivity. Theoretical methods allow significant expansion of the possibilities of aptamer design. In this study, we review theoretical and joint theoretical-experimental studies dedicated to aptamer design and modeling. We consider aptamers with different targets, such as proteins, antibiotics, organophosphates, nucleobases, amino acids, and drugs. During nucleic acid modeling and in silico design, a full set of in silico methods can be applied, such as docking, molecular dynamics (MD), and statistical analysis. The typical modeling workflow starts with structure prediction. Then, docking of target and aptamer is performed. Next, MD simulations are performed, which allows for an evaluation of the stability of aptamer/ligand complexes and determination of the binding energies with higher accuracy. Then, aptamer/ligand interactions are analyzed, and mutations of studied aptamers made. Subsequently, the whole procedure of molecular modeling can be reiterated. Thus, the interactions between aptamers and their ligands are complex and difficult to understand using only experimental approaches. Docking and MD are irreplaceable when aptamers are studied in silico.