Computation meets experiment: identification of highly efficient fibrillating peptides.

Self-assembling peptides are of huge interest for biological, medical and nanotechnological applications. The enormous chemical variety that is available from the 20 amino acids offers potentially unlimited peptide sequences, but it is currently an issue to predict their supramolecular behavior in a reliable and cheap way. Herein we report a computational method to screen and forecast the aqueous self-assembly propensity of amyloidogenic pentapeptides. This method was found also as an interesting tool to predict peptide crystallinity, which may be of interest for the development of peptide based drugs.

In Silico Analysis of Nanoplastics' and ß-amyloid Fibrils' Interactions.

Plastic pollution has become a global environmental threat, which leads to an increasing concern over the consequences of plastic exposition on global health. Plastic nanoparticles have been shown to influence the folding of proteins and influence the formation of aberrant amyloid proteins, therefore potentially triggering the development of systemic and local amyloidosis. This work aims to study the interaction between nanoplastics and ß-amyloid fibrils to better understand the potential role of nanoplastics in the outbreak of neurodegenerative disorders. Using microsecond-long coarse-grained molecular dynamics simulations, we investigated the interactions between neutral and charged nanoparticles made of the most common plastic materials (i.e., polyethylene, polypropylene, and polystyrene) and ß-amyloid fibrils. We observe that the occurrence of contacts, region of amyloid fibril involved, and specific amino acids mediating the interaction depend on the type and charge of the nanoparticles.

Emergence of Elastic Properties in a Minimalist Resilin-Derived Heptapeptide upon Bromination.

Bromination is herein exploited to promote the emergence of elastic behavior in a short peptide-SDSYGAP-derived from resilin, a rubber-like protein exerting its role in the jumping and flight systems of insects. Elastic and resilient hydrogels are obtained, which also show self-healing behavior, thanks to the promoted non-covalent interactions that limit deformations and contribute to the structural recovery of the peptide-based hydrogel. In particular, halogen bonds may stabilize the ß-sheet organization working as non-covalent cross-links between nearby peptide strands. Importantly, the unmodified peptide (i.e., wild type) does not show such properties. Thus, SDSY(3,5-Br)GAP is a novel minimalist peptide elastomer.

Micro- and Nanoplastics' Effects on Protein Folding and Amyloidosis.

A significant portion of the world's plastic is not properly disposed of and, through various processes, is degraded into microscopic particles termed micro- and nanoplastics. Marine and terrestrial faunae, including humans, inevitably get in contact and may inhale and ingest these microscopic plastics which can deposit throughout the body, potentially altering cellular and molecular functions in the nervous and other systems. For instance, at the cellular level, studies in animal models have shown that plastic particles can cross the blood-brain barrier and interact with neurons, and thus affect cognition. At the molecular level, plastics may specifically influence the folding of proteins, induce the formation of aberrant amyloid proteins, and therefore potentially trigger the development of systemic and local amyloidosis. In this review, we discuss the general issue of plastic micro- and nanoparticle generation, with a focus on their effects on protein folding, misfolding, and their possible clinical implications.

Rational backbone redesign of a fructosyl peptide oxidase to widen its active site access tunnel.

Fructosyl peptide oxidases (FPOXs) are enzymes currently used in enzymatic assays to measure the concentration of glycated hemoglobin and albumin in blood samples, which serve as biomarkers of diabetes. However, since FPOX are unable to work directly on glycated proteins, current enzymatic assays are based on a preliminary proteolytic digestion of the target proteins. Herein, to improve the speed and costs of the enzymatic assays for diabetes testing, we applied a rational design approach to engineer a novel enzyme with a wider access tunnel to the catalytic site, using a combination of Rosetta design and molecular dynamics simulations. Our final design, L3_35A, shows a significantly wider and shorter access tunnel, resulting from the deletion of five-amino acids lining the gate structures and from a total of 35 point mutations relative to the wild-type (WT) enzyme. Indeed, upon experimental testing, our engineered enzyme shows good structural stability and maintains significant activity relative to the WT.

The Anti-Amyloidogenic Action of Doxycycline: A Molecular Dynamics Study on the Interaction with Aß42.

The pathological aggregation of amyloidogenic proteins is a hallmark of many neurological diseases, including Alzheimer's disease and prion diseases. We have shown both in vitro and in vivo that doxycycline can inhibit the aggregation of Aß42 amyloid fibrils and disassemble mature amyloid fibrils. However, the molecular mechanisms of the drug's anti-amyloidogenic property are not understood. In this study, a series of molecular dynamics simulations were performed to explain the molecular mechanism of the destabilization of Aß42 fibrils by doxycycline and to compare the action of doxycycline with those of iododoxorubicin (a toxic structural homolog of tetracyclines), curcumin (known to have anti-amyloidogenic activity) and gentamicin (an antibiotic with no experimental evidence of anti-amyloidogenic properties). We found that doxycycline tightly binds the exposed hydrophobic amino acids of the Aß42 amyloid fibrils, partly leading to destabilization of the fibrillar structure. Clarifying the molecular determinants of doxycycline binding to Aß42 may help devise further strategies for structure-based drug design for Alzheimer's disease.

Nanostructure and stability of calcitonin amyloids.

Calcitonin is a 32-amino acid thyroid hormone that can form amyloid fibrils. The structural basis of the fibril formation and stabilization is still debated and poorly understood. The reason is that NMR data strongly suggest antiparallel ß-sheet calcitonin assembly, whereas modeling studies on the short DFNKF peptide (corresponding to the sequence from Asp15 to Phe19 of human calcitonin and reported as the minimal amyloidogenic module) show that it assembles with parallel ß-sheets. In this work, we first predict the structure of human calcitonin through two complementary molecular dynamics (MD) methods, finding that human calcitonin forms an α-helix. We use extensive MD simulations to compare previously proposed calcitonin fibril structures. We find that two conformations, the parallel arrangement and one of the possible antiparallel structures (with Asp15 and Phe19 aligned), are highly stable and ordered. Nonetheless, fibrils with parallel molecules show bulky loops formed by residues 1 to 7 located on the same side, which could limit or prevent the formation of larger amyloids. We investigate fibrils formed by the DFNKF peptide by simulating different arrangements of this amyloidogenic core sequence. We show that DFNKF fibrils are highly stable when assembled in parallel ß-sheets, whereas they quickly unfold in antiparallel conformation. Our results indicate that the DFNKF peptide represents only partially the full-length calcitonin behavior. Contrary to the full-length polypeptide, in fact, the DFNKF sequence is not stable in antiparallel conformation, suggesting that the residue flanking the amyloidogenic peptide contributes to the stabilization of the experimentally observed antiparallel ß-sheet packing.

Crystal structure of the deglycating enzyme Amadoriase I in its free form and substrate-bound complex.

Amadoriases, also known as fructosyl amine oxidases (FAOX), are enzymes that catalyze the de-glycosylation of fructosyl amino acids. As such, they are excellent candidates for the development of enzyme-based diagnostic and therapeutic tools against age- and diabetes-induced protein glycation. However, mostly because of the lack of a complete structural characterization of the different members of the family, the molecular bases of their substrate specificity have yet to be fully understood. The high resolution crystal structures of the free and the substrate-bound form of Amadoriase I shown herein allow for the identification of key structural features that account for the diverse substrate specificity shown by this class of enzymes. This is of particular importance in the context of the rather limited and partially incomplete structural information that has so far been available in the literature on the members of the FAOX family. Moreover, using molecular dynamics simulations, we describe the tunnel conformation and the free energy profile experienced by the ligand in going from bulk water to the catalytic cavity, showing the presence of four gating helices/loops, followed by an "L-shaped" narrow cavity. In summary, the tridimensional architecture of Amadoriase I presented herein provides a reference structural framework for the design of novel enzymes for diabetes monitoring and protein deglycation. Proteins 2016; 84:744-758. © 2016 Wiley Periodicals, Inc.

Role of intrafibrillar collagen mineralization in defining the compressive properties of nascent bone.

Bone is the sole biological material found in the human body that is able to sustain compressive loads. However, although the structure of bone is well-known (it is a natural composite of collagen protein and hydroxyapatite mineral with a complex hierarchical organization), the details about the mechanisms that govern deformation at the molecular scale under compressive loading are still not completely understood. To investigate the molecular origins of bone's unique compressive properties, we perform full atomistic simulations of the three-dimensional molecular structure of a mineralized collagen fibril, focusing on the role of intrafibrillar mineral densities in dictating the mechanical performance under compressive loading. We find that as the mineral density increases, the compressive modulus of the mineralized collagen increases monotonically and well beyond that of pure collagen fibrils. These findings reveal the mechanism by which bone is able to achieve superior load bearing characteristics beyond its individual constituents. Moreover, we find that intrafibrillar mineralization leads to compressive moduli that are one order of magnitude lower than the macroscale modulus of bone, indicating that extrafibrillar mineralization is mandatory for providing the load bearing properties of bone, consistent with recent experimental observations.

Molecular insights into nanoplastics-peptides binding and their interactions with the lipid membrane.

Micro- and nanoplastics have become a significant concern, due to their ubiquitous presence in the environment. These particles can be internalized by the human body through ingestion, inhalation, or dermal contact, and then they can interact with environmental or biological molecules, such as proteins, resulting in the formation of the protein corona. However, information on the role of protein corona in the human body is still missing. Coarse-grain models of the nanoplastics and pentapeptides were created and simulated at the microscale to study the role of protein corona. Additionally, a lipid bilayer coarse-grain model was reproduced to investigate the behavior of the coronated nanoplastics in proximity of a lipid bilayer. Hydrophobic and aromatic amino acids have a high tendency to create stable bonds with all nanoplastics. Moreover, polystyrene and polypropylene establish bonds with polar and charged amino acids. When the coronated nanoplastics are close to a lipid bilayer, different behaviors can be observed. Polyethylene creates a single polymeric chain, while polypropylene tends to break down into its single chains. Polystyrene can both separate into its individual chains and remain aggregated. The protein corona plays an important role when interacting with the nanoplastics and the lipid membrane. More studies are needed to validate the results and to enhance the complexity of the systems.

Unbiased in silico design of pH-sensitive tetrapeptides.

We used coarse-grain molecular dynamics simulations to screen all possible histidine-bearing tetrapeptide sequences, finding novel peptide sequences with pH-tunable assembly properties. These tetrapeptides could be used for various biological applications, such as triggered delivery of bioactive molecules.

Unsupervised deep learning for molecular dynamics simulations: a novel analysis of protein-ligand interactions in SARS-CoV-2 Mpro.

Molecular dynamics (MD) simulations, which are central to drug discovery, offer detailed insights into protein-ligand interactions. However, analyzing large MD datasets remains a challenge. Current machine-learning solutions are predominantly supervised and have data labelling and standardisation issues. In this study, we adopted an unsupervised deep-learning framework, previously benchmarked for rigid proteins, to study the more flexible SARS-CoV-2 main protease (Mpro). We ran MD simulations of Mpro with various ligands and refined the data by focusing on binding-site residues and time frames in stable protein conformations. The optimal descriptor chosen was the distance between the residues and the center of the binding pocket. Using this approach, a local dynamic ensemble was generated and fed into our neural network to compute Wasserstein distances across system pairs, revealing ligand-induced conformational differences in Mpro. Dimensionality reduction yielded an embedding map that correlated ligand-induced dynamics and binding affinity. Notably, the high-affinity compounds showed pronounced effects on the protein's conformations. We also identified the key residues that contributed to these differences. Our findings emphasize the potential of combining unsupervised deep learning with MD simulations to extract valuable information and accelerate drug discovery.

Tailoring FPOX enzymes for enhanced stability and expanded substrate recognition.

Fructosyl peptide oxidases (FPOX) are deglycating enzymes that find application as key enzymatic components in diabetes monitoring devices. Indeed, their use with blood samples can provide a measurement of the concentration of glycated hemoglobin and glycated albumin, two well-known diabetes markers. However, the FPOX currently employed in enzymatic assays cannot directly detect whole glycated proteins, making it necessary to perform a preliminary proteolytic treatment of the target protein to generate small glycated peptides that can act as viable substrates for the enzyme. This is a costly and time consuming step. In this work, we used an in silico protein engineering approach to enhance the overall thermal stability of the enzyme and to improve its catalytic activity toward large substrates. The final design shows a marked improvement in thermal stability relative to the wild type enzyme, a distinct widening of its access tunnel and significant enzymatic activity towards a range of glycated substrates.

Thickness of hydroxyapatite nanocrystal controls mechanical properties of the collagen-hydroxyapatite interface.

Collagen-hydroxyapatite interfaces compose an important building block of bone structures. While it is known that the nanoscale structure of this elementary building block can affect the mechanical properties of bone, a systematic understanding of the effect of the geometry on the mechanical properties of this interface between protein and mineral is lacking. Here we study the effect of geometry, different crystal surfaces, and hydration on the mechanical properties of collagen-hydroxyapatite interfaces from an atomistic perspective, and discuss underlying deformation mechanisms. We find that the presence of hydroxyapatite significantly enhances the tensile modulus and strength compared with a tropocollagen molecule alone. The stiffening effect is strongly dependent on the thickness of the mineral crystal until a plateau is reached at 2 nm crystal thickness. We observe no significant differences due to the mineral surface (Ca surface vs OH surface) or due to the presence of water. Our result shows that the hydroxyapatite crystal with its thickness confined to the nanometer size efficiently increases the tensile modulus and strength of the collagen-hydroxyapatite composite, agreeing well with experimental observations that consistently show the existence of extremely thin mineral flakes in various types of bones. We also show that the collagen-hydroxyapatite interface can be modeled with an elastic network model which, based on the results of atomistic simulations, provides a good estimate of the surface energy and other mechanical features.

Hierarchical structure and nanomechanics of collagen microfibrils from the atomistic scale up.

Collagen constitutes one-third of the human proteome, providing mechanical stability, elasticity, and strength to organisms and is the prime construction material in biology. Collagen is also the dominating material in the extracellular matrix and its stiffness controls cell differentiation, growth, and pathology. However, the origin of the unique mechanical properties of collagenous tissues, and in particular its stiffness, extensibility, and nonlinear mechanical response at large deformation, remains unknown. By using X-ray diffraction data of a collagen fibril (Orgel, J. P. R. O. et al. Proc. Natl. Acad. Sci. 2006, 103, 9001) here we present an experimentally validated model of the nanomechanics of a collagen microfibril that incorporates the full biochemical details of the amino acid sequence of constituting molecules and the nanoscale molecular arrangement. We demonstrate by direct mechanical testing that hydrated (wet) collagen microfibrils feature a Young's modulus of ≈300 MPa at small, and ≈1.2 GPa at larger deformation in excess of 10% strain, which is in excellent agreement with experimental data. We find that dehydrated (dry) collagen microfibrils show a significantly increased Young's modulus of ≈1.8-2.25 GPa, which is in agreement with experimental measurements and owing to tighter molecular packing. Our results show that the unique mechanical properties of collagen microfibrils arise due to their hierarchical structure at the nanoscale, where key deformation mechanisms are straightening of twisted triple-helical molecules at small strains, followed by axial stretching and eventual molecular uncoiling. The establishment of a model of hierarchical deformation mechanisms explains the striking difference of the elastic modulus of collagen fibrils compared with single molecules, which is found in the range of 4.8 ± 2 GPa, or ≈10-20 times greater. We find that collagen molecules alone are not capable of providing the broad range of mechanical functionality required for physiological function of collagenous tissues. Rather, the existence of an array of deformation mechanisms, derived from the hierarchical makeup of the material, is critical to the material's ability to confer key mechanical properties, specifically large extensibility, strain hardening, and toughness, despite the limitation that collagenous materials are constructed from only few distinct amino acids. The atomistic model of collagen microfibril mechanics now enables the bottom-up elucidation of structure-property relationships in a broader class of collagen materials (e.g., tendon, bone, cornea), including studies of genetic disease where the incorporation of biochemical details is essential. The availability of a molecular-based model of collagen tissues may eventually result in novel nanomedicine approaches to develop treatments for a broad class of collagen diseases and the design of de novo biomaterials for regenerative medicine.

In Silico Engineering of Enzyme Access Tunnels.

Enzyme engineering is a tailoring process that allows the modification of naturally occurring enzymes to provide them with improved catalytic efficiency, stability, or specificity. By introducing partial modifications to their sequence and to their structural features, enzyme engineering can transform natural enzymes into more efficient, specific and resistant biocatalysts and render them suitable for virtually countless industrial processes. Current enzyme engineering methods mostly target the active site of the enzyme, where the catalytic reaction takes place. Nonetheless, the tunnel that often connects the surface of an enzyme with its buried active site plays a key role in the activity of the enzyme as it acts as a gatekeeper and regulates the access of the substrate to the catalytic pocket. Hence, there is an increasing interest in targeting the sequence and the structure of substrate entrance tunnels in order to fine-tune enzymatic activity, regulate substrate specificity, or control reaction promiscuity.In this chapter, we describe the use of a rational in silico design and screening method to engineer the access tunnel of a fructosyl peptide oxidase with the aim to facilitate access to its catalytic site and to expand its substrate range. Our goal is to engineer this class of enzymes in order to utilize them for the direct detection of glycated proteins in diabetes monitoring devices. The design strategy involves remodeling of the backbone structure of the enzyme , a feature that is not possible with conventional enzyme engineering techniques such as single-point mutagenesis and that is highly unlikely to occur using a directed evolution approach.The proposed strategy, which results in a significant reduction in cost and time for the experimental production and characterization of candidate enzyme variants, represents a promising approach to the expedited identification of novel and improved enzymes. Rational enzyme design aims to provide in silico strategies for the fast, accurate, and inexpensive development of biocatalysts that can meet the needs of multiple industrial sectors, thus ultimately promoting the use of green chemistry and improving the efficiency of chemical processes.

Composite Peptide-Agarose Hydrogels for Robust and High-Sensitivity 3D Immunoassays.

Canonical immunoassays rely on highly sensitive and specific capturing of circulating biomarkers by interacting biomolecular baits. In this frame, bioprobe immobilization in spatially discrete three-dimensional (3D) spots onto analytical surfaces by hydrogel encapsulation was shown to provide relevant advantages over conventional two-dimensional (2D) platforms. Yet, the broad application of 3D systems is still hampered by hurdles in matching their straightforward fabrication with optimal functional properties. Herein, we report on a composite hydrogel obtained by combining a self-assembling peptide (namely, Q3 peptide) with low-temperature gelling agarose that is proved to have simple and robust application in the fabrication of microdroplet arrays, overcoming hurdles and limitations commonly associated with 3D hydrogel assays. We demonstrate the real-case scenario feasibility of our 3D system in the profiling of Covid-19 patients' serum IgG immunoreactivity, which showed remarkably improved signal-to-noise ratio over canonical assays in the 2D format and exquisite specificity. Overall, the new two-component hydrogel widens the perspectives of hydrogel-based arrays and represents a step forward towards their routine use in analytical practices.

In silico prediction of the in vitro behavior of polymeric gene delivery vectors.

Non-viral gene delivery vectors have increasingly come under the spotlight, but their performaces are still far from being satisfactory. Therefore, there is an urgent need for forecasting tools and screening methods to enable the development of ever more effective transfectants. Here, coarse-grained (CG) models of gold standard transfectant poly(ethylene imine)s (PEIs) have been profitably used to investigate and highlight the effect of experimentally-relevant parameters, namely molecular weight (2 vs. 10 kDa) and topologies (linear vs. branched), protonation state, and ammine-to-phosphate ratios (N/Ps), on the complexation and the gene silencing efficiency of siRNA molecules. The results from the in vitro screening of cationic polymers and conditions were used to validate the in silico platform that we developed, such that the hits which came out of the CG models were of high practical relevance. We show that our in silico platform enables to foresee the most suitable conditions for the complexation of relevant siRNA-polycation assemblies, thereby providing a reliable predictive tool to test bench transfectants in silico, and foster the design and development of gene delivery vectors.

Molecular and mesoscale mechanisms of osteogenesis imperfecta disease in collagen fibrils.

Osteogenesis imperfecta (OI) is a genetic disorder in collagen characterized by mechanically weakened tendon, fragile bones, skeletal deformities, and in severe cases, prenatal death. Although many studies have attempted to associate specific mutation types with phenotypic severity, the molecular and mesoscale mechanisms by which a single point mutation influences the mechanical behavior of tissues at multiple length scales remain unknown. We show by a hierarchy of full atomistic and mesoscale simulation that OI mutations severely compromise the mechanical properties of collagenous tissues at multiple scales, from single molecules to collagen fibrils. Mutations that lead to the most severe OI phenotype correlate with the strongest effects, leading to weakened intermolecular adhesion, increased intermolecular spacing, reduced stiffness, as well as a reduced failure strength of collagen fibrils. We find that these molecular-level changes lead to an alteration of the stress distribution in mutated collagen fibrils, causing the formation of stress concentrations that induce material failure via intermolecular slip. We believe that our findings provide insight into the microscopic mechanisms of this disease and lead to explanations of characteristic OI tissue features such as reduced mechanical strength and a lower cross-link density. Our study explains how single point mutations can control the breakdown of tissue at much larger length scales, a question of great relevance for a broad class of genetic diseases.

Alport syndrome mutations in type IV tropocollagen alter molecular structure and nanomechanical properties.

Alport Syndrome is a genetic disease characterized by breakdown of the glomerular basement membrane (GBM) around blood vessels in the kidney, leading to kidney failure in most patients. It is the second most inherited kidney disease in the US, and many other symptoms are associated with the disease, including hearing loss and ocular lesions. Here we probe the molecular level structure-property relationships of this disease using a bottom-up computational materiomics approach implemented through large-scale molecular dynamics simulation. Since the GBM is under constant mechanical loading due to blood flow, changes in mechanical properties due to amino acid mutations may be critical in the symptomatic GBM breakdown seen in Alport Syndrome patients. Through full-atomistic simulations in explicit solvent, the effects of single-residue glycine substitution mutations of varying clinical severity are studied in short segments of type IV tropocollagen molecules. The segments with physiological amino acid sequences are equilibrated and then subjected to tensile loading. Major changes are observed at the single molecule level of the mutated sequence, including a bent shape of the structures after equilibration (with the kink located at the mutation site) and a significant alteration of the molecules' stress-strain responses and stiffnesses. These results suggest that localized structural changes at amino acid level induce severe alterations of the molecular properties. Our study opens a new approach in pursuing a bottom-up multi-scale analysis of this disease.