Structural role of compensatory amino acid replacements in the α-synuclein protein

A subset of familial Parkinson's disease (PD) cases is associated with the presence of disease-causing point mutations in human α-synuclein [huAS(wt)], including A53T. Surprisingly, the human neurotoxic amino acid 53T is present in non-primate, wild-type sequences of α-synucleins, including that expressed by mice [mAS(wt)]. Because huAS(A53T) causes neurodegeneration when expressed in rodents, the amino acid changes between the wild-type human protein [huAS(wt)] and mAS(wt) might act as intramolecular suppressors of A53T toxicity in the mouse protein, restoring its physiological structure and function. The lack of structural information for mAS(wt) in aqueous solution has prompted us to conduct a comparative molecular dynamics study of huAS(wt), huAS(A53T), and mAS(wt) in water at 300 K. The calculations are based on an ensemble of nuclear magnetic resonance-derived huAS(wt) structures. huAS(A53T) turns out to be more flexible and less compact than huAS(wt). Its central (NAC) region, involved in fibril formation by the protein, is more solvent-exposed than that of the wild-type protein, in agreement with nuclear magnetic resonance data. The compactness of mAS(wt) is similar to that of the human protein. In addition, its NAC region is less solvent-exposed and more rigid than that of huAS(A53T). All of these features may be caused by an increase in the level of intramolecular interactions on passing from huAS(A53T) to mAS(wt). We conclude that the presence of "compensatory replacements" in the mouse protein causes a significant change in the protein relative to huAS(A53T), restoring features not too dissimilar to those of the human protein.

A soft-core Gay-Berne model for the simulation of liquid crystals by Hamiltonian replica exchange

The Gay-Berne (GB) potential has proved highly successful in the simulation of liquid crystal phases, although it is fairly demanding in terms of resources for simulations of large (e.g., N>10(5)) systems, as increasingly required in applications. Here, we introduce a soft-core GB model, which exhibits both liquid crystal phase behavior and rapid equilibration. We show that the Hamiltonian replica exchange method, coupled with the newly introduced soft-core GB model, can effectively speed up the equilibration of a GB liquid crystal phase by frequent exchange of configurations between replicas, while still recovering the mesogenic properties of the standard GB potential.

Theoretical characterization of the structural and hole transport dynamics in liquid-crystalline phthalocyanine stacks

We present a joint molecular dynamics (MD)/kinetic Monte Carlo (KMC) study aimed at the atomistic description of charge transport in stacks of liquid-crystalline tetraalkoxy-substituted, metal-free phthalocyanines. The molecular dynamics simulations reproduce the major structural features of the mesophases, in particular, a phase transition around 340 K between the rectangular and hexagonal phases. Charge transport simulations based on a Monte Carlo algorithm show an increase by 2 orders of magnitude in the hole mobility when accounting for the rotational and translational dynamics. The results point to the formation of dynamical structural defects along the columns.

Computer simulations of biaxial nematics

Biaxial nematic (N(b)) liquid crystals are a fascinating condensed matter phase that has baffled, for more than thirty years, scientists engaged in the challenge of demonstrating its actual existence, and which has only recently been experimentally found. During this period computer simulations of model N(b) have played an important role, both in providing the basic physical properties to be expected from these systems, and in giving clues about the molecular features essential for the thermodynamic stability of N(b) phases. However, simulation studies are expected to be even more crucial in the future for unravelling the structural features of biaxial mesogens at the molecular level, and for helping in the design and optimization of devices towards the technological deployment of N(b) materials. This review article gives an overview of the simulation work performed so far, and relying on the recent experimental findings, focuses on the still unanswered questions which will determine the future challenges in the field.

Controlling surface defect valence in colloids

We perform large-scale Monte Carlo simulations of orientational ordering in nematic shells and study the type and position of topological defects when an external electric field (homogeneous or quadrupolar) is applied. The field-induced variation of the defect number (and strength) can be used to change the valence of colloidal particles coated with a nematic layer.

Unveiling the role of histidine and tyrosine residues on the conformation of the avian prion hexarepeat domain

The prion protein (PrPC) is a glycoprotein that in mammals, differently from avians, can lead to prion diseases, by misfolding into a beta-sheet-rich pathogenic isoform (PrPSc). Mammal and avian proteins show different N-terminal tandem repeats: PHGGGWGQ and PHNPGY, both containing histidine, whereas tyrosine is included only in the primary sequence of the avian protein. Here, by means of potentiometric, circular dichroism (CD), and molecular dynamics (MD) studies at different pH values, we have investigated the conformation of the avian tetrahexarepeat (PHNPGY)4 (TetraHexaPY) with both N- and C-termini blocked by acetylation and amidation, respectively. We have found, also with the help of a recently proposed protein chirality indicator (Pietropaolo, A.; Muccioli, L.; Berardi, R.; Zannoni, C. Proteins 2008, 70, 667-677), a conformational dependence on the protonation states of histidine and tyrosine residues: the turn formation is pH driven, and at physiological pH a pivotal role is played by the tyrosine OH groups which give rise to a very compact bent structure of backbone upon forming a hydrogen-bond network.

Biaxial liquid-crystal elastomers: a lattice model

We present a simple coarse-grained lattice model for monodomain biaxial liquid-crystal elastomers and perform large-scale Monte Carlo simulations in the proposed model system. Orientational ordering--uniaxial or biaxial--reflects in sample deformations on cooling the system. The simulation output is used to predict calorimetry data and deuterium magnetic resonance spectra.

A chirality index for investigating protein secondary structures and their time evolution

We propose a methodology for the description of the secondary structure of proteins, based on assigning a chirality parameter to short aminoacid sequences according to their arrangement in space at a certain time. We validated the method on ideal and crystalline structures, showing that it can assign secondary structures and that this assignment is robust with respect to random conformational perturbations. From the values of the index and its pattern along a sequence it is possible to recognize many structural motifs of a protein, and in particular poly-L-proline II left-handed helices, often not detected by secondary structure assignment algorithms. Assigning an instantaneous chirality index to the fragments also allows the dynamics to be studied. With this purpose, molecular dynamics simulations were carried out in water for selected hemoglobin (110 ns) and immunoglobulin antigen fragments (50 ns), showing the capability of the chiral index in identifying the stable secondary structure elements, as well as in following their time evolution and conformational changes during the trajectory.

Core charge distribution and self assembly of columnar phases: the case of triphenylenes and azatriphenylenes

BACKGROUND

The relation between the structure of discotic molecules and columnar properties, a crucial point for the realization of new advanced materials, is still largely unknown. A paradigmatic case is that hexa-alkyl-thio substituted triphenylenes present mesogenic behavior while the corresponding azatriphenylenes, similar in shape and chemical structure, but with a different core charge distribution, do not form any liquid crystalline mesophase. This study is aimed at investigating, with the help of computer simulations techniques, the effects on phase behaviour of changes of the charge distribution in the discotic core.

RESULTS

We described the shape and the pair, dispersive and electrostatic, interactions of hexa alkyl triphenylenes by uniaxial Gay-Berne discs with embedded point charges. Gay-Berne parameters were deduced by fitting the dispersive energies obtained from an atomistic molecular dynamics simulation of a small sample of hexa-octyl-thio triphenylene molecules in columnar phase, while a genetic algorithm was used to get a minimal set of point charges that properly reproduces the ab anitio electrostatic potential. We performed Monte Carlo simulations of three molecular models: the pure Gay-Berne disc, used as a reference, the Gay-Berne disc with hexa-thio triphenylene point charges, the Gay-Berne disc with hexa-thio azatriphenylene point charges. The phase diagram of the pure model evidences a rich polymorphism, with isotropic, columnar and crystalline phases at low pressure, and the appearance of nematic phase at higher pressure.

CONCLUSION

We found that the intermolecular electrostatic potential among the cores is fundamental in sta-bilizing/destabilizing columnar phases; in particular the triphenylene charge distribution stabilizes the columnar structure, while the azatriphenylene distribution suppresses its formation in favor of the nematic phase. We believe the present model could be successfully employed as the basis for coarse-grained level simulations of a wider class of triphenylene derivatives.

Expected and Unexpected Behavior of the Orientational Order and Dynamics Induced by Azobenzene Solutes in a Nematic

We have explored the changes in the phase stability, orientational order, and dynamics of the nematic 4-cyano-4'-n-pentylbiphenyl (5CB) doped with either the trans or the cis form of different p-azobenzene derivatives using the ESR spin-probe technique. In particular, we have studied the effects induced by each of the seven nonmesogenic 4-R-phenylazobenzenes (R = H, F, Br, CH3, CF3, On-Bu, Ot-Bu) at 1% and 7% mole fraction on the order parameter <P2> and on the shift of the nematic-isotropic transition temperature (TNI), as reported by a nitroxide spin probe, and we have tried to relate them to the solute shape and charge distribution. In all the cases the presence of the azo-derivative causes a depression of T(NI), more pronounced for the cis isomers. The dependence of <P2> on the reduced temperature T* = T/T(NI) remains the same as that of pure 5CB in all trans-doped samples at 1% and 7% and decreases only slightly in the cis at 1%. However, we observe different and in some cases large variations (up to 25%) in <P2> for the cis at 7%, showing solute effects that go beyond the shift in T(NI). Surprisingly enough, even at the highest concentration, the probe dynamics appears to be essentially independent of the nature, the configuration, and the concentration of the different solutes and very similar to that observed in the pure 5CB.