HyperBeta: characterizing the structural dynamics of proteins and self-assembling peptides.

Self-assembling processes are ubiquitous phenomena that drive the organization and the hierarchical formation of complex molecular systems. The investigation of assembling dynamics, emerging from the interactions among biomolecules like amino-acids and polypeptides, is fundamental to determine how a mixture of simple objects can yield a complex structure at the nano-scale level. In this paper we present HyperBeta, a novel open-source software that exploits an innovative algorithm based on hyper-graphs to efficiently identify and graphically represent the dynamics of [Formula: see text]-sheets formation. Differently from the existing tools, HyperBeta directly manipulates data generated by means of coarse-grained molecular dynamics simulation tools (GROMACS), performed using the MARTINI force field. Coarse-grained molecular structures are visualized using HyperBeta 's proprietary real-time high-quality 3D engine, which provides a plethora of analysis tools and statistical information, controlled by means of an intuitive event-based graphical user interface. The high-quality renderer relies on a variety of visual cues to improve the readability and interpretability of distance and depth relationships between peptides. We show that HyperBeta is able to track the [Formula: see text]-sheets formation in coarse-grained molecular dynamics simulations, and provides a completely new and efficient mean for the investigation of the kinetics of these nano-structures. HyperBeta will therefore facilitate biotechnological and medical research where these structural elements play a crucial role, such as the development of novel high-performance biomaterials in tissue engineering, or a better comprehension of the molecular mechanisms at the basis of complex pathologies like Alzheimer's disease.

Elucidating Self-Assembling Peptide Aggregation via Morphoscanner: A New Tool for Protein-Peptide Structural Characterization.

Self-assembling and molecular folding are ubiquitous in Nature: they drive the organization of systems ranging from living creatures to DNA molecules. Elucidating the complex dynamics underlying these phenomena is of crucial importance. However, a tool for the analysis of the various phenomena involved in protein/peptide aggregation is still missing. Here, an innovative software is developed and validated for the identification and visualization of b-structuring and b-sheet formation in both simulated systems and crystal structures of proteins and peptides. The novel software suite, dubbed Morphoscanner, is designed to identify and intuitively represent b-structuring and b-sheet formation during molecular dynamics trajectories, paying attention to temporary strand-strand alignment, suboligomer formation and evolution of local order. Self-assembling peptides (SAPs) constitute a promising class of biomaterials and an interesting model to study the spontaneous assembly of molecular systems in vitro. With the help of coarse-grained molecular dynamics the self-assembling of diverse SAPs is simulated into molten aggregates. When applied to these systems, Morphoscanner highlights different b-structuring schemes and kinetics related to SAP sequences. It is demonstrated that Morphoscanner is a novel versatile tool designed to probe the aggregation dynamics of self-assembling systems, adaptable to the analysis of differently coarsened simulations of a variety of biomolecules.

Nanomaterials design and tests for neural tissue engineering.

Nanostructured scaffolds recently showed great promise in tissue engineering: nanomaterials can be tailored at the molecular level and scaffold morphology may more closely resemble features of extracellular matrix components in terms of porosity, framing and biofunctionalities. As a consequence, both biomechanical properties of scaffold microenvironments and biomaterial-protein interactions can be tuned, allowing for improved transplanted cell engraftment and better controlled diffusion of drugs. Easier said than done, a nanotech-based regenerative approach encompasses different fields of know-how, ranging from in silico simulations, nanomaterial synthesis and characterization at the nano-, micro- and mesoscales to random library screening methods (e.g. phage display), in vitro cellular-based experiments and validation in animal models of the target injury. All of these steps of the "assembly line" of nanostructured scaffolds are tightly interconnected both in their standard analysis techniques and in their most recent breakthroughs: indeed their efforts have to jointly provide the deepest possible analyses of the diverse facets of the challenging field of neural tissue engineering. The purpose of this review is therefore to provide a critical overview of the recent advances in and drawbacks and potential of each mentioned field, contributing to the realization of effective nanotech-based therapies for the regeneration of peripheral nerve transections, spinal cord injuries and brain traumatic injuries. Far from being the ultimate overview of such a number of topics, the reader will acknowledge the intrinsic complexity of the goal of nanotech tissue engineering for a conscious approach to the development of a regenerative therapy and, by deciphering the thread connecting all steps of the research, will gain the necessary view of its tremendous potential if each piece of stone is correctly placed to work synergically in this impressive mosaic.

Conformational behaviour determines the low-relaxivity state of a conditional MRI contrast agent.

The conformational behaviour in aqueous solution of the EgadMe complex, a conditional gadolinium-based contrast agent sensitive to beta-galactosidase enzymatic activity, is investigated by means of ab initio calculations and classical molecular dynamics simulations. Furthermore, force field parameterization of gadolinium-ligand interactions is performed, and its reliability is tested on the bench mark [Gd(DOTA)](-) system by MD simulations. Both computational methods highlight the presence in EgadMe of two main conformational isomers. The lowest energy conformation is a "close" form, corresponding to a state of low-relaxivity (MRI "inactive"), in which the ninth coordination site of the gadolinium ion is occupied by one oxygen atom of the galactopyranose residue. The second isomer, which is 2.9 (at ab initio level) and 4.2 (at MD level) kcal mol(-1) above the global minimum, presents an "open" form, corresponding to a state of high-relaxivity (MRI "active") in which one water molecule coordinates the ion. These results are consistent with experimental findings reported for EgadMe, and show that competition at the ninth coordination site of gadolinium ion, between the intra (the galactopyranose residue) and inter (water molecules) molecular interactions, affects the relaxivity of this system.

Spontaneous beta-helical fold in prion protein: the case of PrP(82-146).

The presence of amyloid is a hallmark of Gerstmann-Sträussler-Scheinker (GSS) disease, which is a prion disease caused by germ line mutations in the PRNP gene. The major component of amyloid is a fragment spanning residues from 81-82 to 144-153, part of the minimal sequence thought to play a crucial role in the conversion reaction and to sustain prion replication. We present here a molecular dynamics study on the 82-146 peptide from the human prion protein. The aim is to identify its aggregation-prone folds. The 82-146 prion sequence corresponds to a naturally occurring prion peptide able to form fibrils rich in parallel beta-sheets. A spontaneous right-handed beta-helical arrangement with 13 residues per turn can be observed in the 103-135 segment of the 82-146 peptide. The observed fold is in accordance with the evidence of a parallel beta-sheet organization in amyloid and with experiments on 82-146 discussed in the literature. To elucidate the conformational properties that trigger this peptide's aggregation propensity, the conformational behavior of peptides of different length (106-126 and 113-120 prion segments) was also investigated. Simulation analysis has led to some interesting considerations on sequence specific flexibility and the effects of growth. Comparing peptides of different length allows the localization of the origin of the beta-helix conformational propensity in the 106-126 segment, though longer sequences appear necessary for a clear beta-helical arrangement. Structural features of the observed 82-146 beta-helical fold are compatible with the "dock and lock" mechanism proposed to interpret peptide aggregation kinetics.

The anti-fibrillogenic activity of tetracyclines on PrP 106-126: a 3D-QSAR study.

There is evidence that Tetracyclines are potentially useful drugs to treat prion disease, the fatal neurodegenerative disease in which cellular prion proteins change in conformation to become a disease-specific species (PrP(Sc)). Based on an in vitro anti-fibrillogenesis test, and using the peptide PrP106-126 in the presence of tetracycline and 14 derivatives, we carried out a three-dimensional quantitative structure-activity relationship (3D-QSAR) study to investigate the stereoelectronic features required for anti-fibrillogenic activity. A preliminary variable reduction technique was used to search for grid points where statistical indexes of interaction potential distributions present local maximum (or minimum) values. Variable selection genetic algorithms were then used to search for the best 3D-QSAR models. A 6-variable model showed the best predictability of the anti-fibrillogenic activity that highlighted the best tetracycline substitution patterns: hydroxyl group presence in positions 5 and 6, electrodonor substituents on the aromatic D-ring, alkylamine substituent at the amidic group in position 2 and non-epi configuration of the NMe2 group.

Computational approaches to shed light on molecular mechanisms in biological processes.

Computational approaches based on Molecular Dynamics simulations, Quantum Mechanical methods and 3D Quantitative Structure-Activity Relationships were employed by computational chemistry groups at the University of Milano-Bicocca to study biological processes at the molecular level. The paper reports the methodologies adopted and the results obtained on Aryl hydrocarbon Receptor and homologous PAS proteins mechanisms, the properties of prion protein peptides, the reaction pathway of hydrogenase and peroxidase enzymes and the defibrillogenic activity of tetracyclines.

Conformational polymorphism of the PrP106-126 peptide in different environments: a molecular dynamics study.

Extensive molecular dynamic simulations (approximately 240 ns) have been used to investigate the conformational behavior of PrP106-126 prion peptide in four different environments (water, dimethyl sulfoxide, hexane, and trifluoroethanol) and under both neutral and acidic conditions. The conformational polymorphism of PrP106-126 in solution observed in the simulations supports the role of this fragment in the structural transition of the native to the abnormal form of prion protein in response to changes in the local environmental conditions. The peptide in solution is primarily unstructured. The simulations show an increased presence of helical structure in an apolar solvent, in agreement with the results from circular dichroism spectroscopy. In water solution, beta-sheet elements were observed between residues 108-112 and either residues 115-121 or 121-126. An alpha-beta transition was observed under neutral conditions. In DMSO, the peptide adopted an extended conformation, in agreement with nuclear magnetic resonance experiments.