Joshua T. Berryman

Deciphering drift time measurements from travelling wave ion mobility spectrometry-mass spectrometry studies.

Detailed knowledge of the tertiary and quaternary structure of proteins and protein complexes is of immense importance in understanding their functionality. Similarly, variations in the conformational states of proteins form the underlying mechanisms behind many biomolecular processes, numerous of which are disease-related. Thus, the availability of reliable and accurate biophysical techniques that can provide detailed information concerning these issues is of paramount importance. Ion mobility spectrometry (IMS) coupled to mass spectrometry (MS) offers a unique opportunity to separate multi-component biomolecular entities and to measure the molecular mass and collision cross-section of individual components in a single, rapid (≤ 2 min) experiment, providing 3D-architectural information directly. Here we report a method of calibrating a commercially available electrospray ionisation (ESI)-travelling wave ion mobility spectrometry (TWIMS)–mass spectrometer using known cross-sectional areas determined for a range of biomolecules by conventional IMS-MS. Using this method of calibration, we have analysed a range of proteins of differing mass and 3D architecture in their native conformations by ESI-TWIMS-MS and found that the cross-sectional areas measured in this way compare extremely favourably with cross-sectional areas calculated using an in-house computing method based on Protein Data Bank NMR-derived co-ordinates. This not only provides a high degree of confidence in the calibration method, but also suggests that the gas phase ESI-TWIMS-MS measurements relate well to solution-based measurements derived from other biophysical techniques. In order to determine which instrumental parameters affect the ESI-TWIMS-MS cross-sectional area calibration, a systematic study of the parameters used to optimise TWIMS drift time separations has been carried out, observing the effect each parameter has on drift times and IMS resolution. Finally, the ESI-TWIMS-MS cross-sectional area calibration has been applied to the analysis of the amyloidogenic protein β2-microglobulin and measurements for three co-populated conformational families, present under denaturing conditions, have been made: the folded, partially unfolded and unfolded states.

Thermodynamic Description of Polymorphism in Q- and N-Rich Peptide Aggregates Revealed by Atomistic Simulation

Amyloid fibrils are long, helically symmetric protein aggregates that can display substantial variation (polymorphism), including alterations in twist and structure at the beta-strand and protofilament levels, even when grown under the same experimental conditions. The structural and thermodynamic origins of this behavior are not yet understood. We performed molecular-dynamics simulations to determine the thermodynamic properties of different polymorphs of the peptide GNNQQNY, modeling fibrils containing different numbers of protofilaments based on the structure of amyloid-like cross-beta crystals of this peptide. We also modeled fibrils with new orientations of the side chains, as well as a de novo designed structure based on antiparallel beta-strands. The simulations show that these polymorphs are approximately isoenergetic under a range of conditions. Structural analysis reveals a dynamic reorganization of electrostatics and hydrogen bonding in the main and side chains of the Gln and Asn residues that characterize this peptide sequence. Q/N-rich stretches are found in several amyloidogenic proteins and peptides, including the yeast prions Sup35-N and Ure2p, as well as in the human poly-Q disease proteins, including the ataxins and huntingtin. Based on our results, we propose that these residues imbue a unique structural plasticity to the amyloid fibrils that they comprise, rationalizing the ability of proteins enriched in these amino acids to form prion strains with heritable and different phenotypic traits.

Systematic Examination of Polymorphism in Amyloid Fibrils by Molecular-Dynamics Simulation

Amyloid fibrils often exhibit polymorphism. Polymorphs are formed when proteins or peptides with identical sequences self-assemble into fibrils containing substantially different arrangements of the β-strands. We used atomistic molecular-dynamics simulation to examine the thermodynamic stability of a amyloid fibrils in different polymorphic forms by performing a systematic investigation of sequence and symmetry space for a series of peptides with a range of physicochemical properties. We show that the stability of fibrils depends on both sequence and the symmetry because these factors determine the availability of favorable interactions between the peptide strands within a sheet and in intersheet packing. By performing a detailed analysis of these interactions as a function of symmetry, we obtained a series of simple design rules that can be used to determine which polymorphs of a given sequence are most likely to form thermodynamically stable fibrils. These rules can potentially be employed to design peptide sequences that aggregate into a preferred polymorphic form for nanotechnological purposes.

ILQINS Hexapeptide, Identified in Lysozyme Left-Handed Helical Ribbons and Nanotubes, Forms Right-Handed Helical Ribbons and Crystals

Amyloid fibrils are implicated in over 20 neurodegenerative diseases. The mechanisms of fibril structuring and formation are not only of medical and biological importance but are also relevant for material science and nanotechnologies due to the unique structural and physical properties of amyloids. We previously found that hen egg white lysozyme, homologous to the disease-related human lysozyme, can form left-handed giant ribbons, closing into nanotubes. By using matrix-assisted laser desorption ionization mass spectrometry analysis, we here identify a key component of such structures: the ILQINS hexapeptide. By combining atomic force microscopy and circular dichorism, we find that this fragment, synthesized by solid-phase peptide synthesis, also forms fibrillar structures in water at pH 2. However, all fibrillar structures formed possess an unexpected right-handed twist, a rare chirality within the corpus of amyloid experimental observations. We confirm by small- and wide-angle X-ray scattering and molecular dynamics simulations that these fibrils are composed of conventional left-handed β-sheets, but that packing stresses between adjacent sheets create this twist of unusual handedness. We also show that the right-handed fibrils represent a metastable state toward β-sheet-based microcrystals formation.

Sampling rare events in nonequilibrium and nonstationary systems

Although many computational methods for rare event sampling exist, this type of calculation is not usually practical for general nonequilibrium conditions, with macroscopically irreversible dynamics and away from both stationary and metastable states. A novel method for calculating the time-series of the probability of a rare event is presented which is designed for these conditions. The method is validated for the cases of the Glauber-Ising model under time-varying shear flow, the Kawasaki-Ising model after a quench into the region between nucleation dominated and spinodal decomposition dominated phase change dynamics, and the parallel open asymmetric exclusion process. The method requires a subdivision of the phase space of the system: it is benchmarked and found to scale well for increasingly fine subdivisions, meaning that it can be applied without detailed foreknowledge of the physically important reaction pathways.

Crystal nucleation mechanism in melts of short polymer chains under quiescent conditions and under shear flow

We present a molecular dynamics simulation study of crystal nucleation from undercooled melts of n-alkanes, and we identify the molecular mechanism of homogeneous crystal nucleation under quiescent conditions and under shear flow. We compare results for n-eicosane (C20) and n-pentacontahectane (C150), i.e., one system below the entanglement length and one above, at 20%-30% undercooling. Under quiescent conditions, we observe that entanglement does not have an effect on the nucleation mechanism. For both chain lengths, the chains first align and then straighten locally, then the local density increases and finally positional ordering sets in. At low shear rates the nucleation mechanism is the same as under quiescent conditions, while at high shear rates the chains align and straighten at the same time. We report on the effects of shear rate and temperature on the nucleation rates and estimate the critical shear rates, beyond which the nucleation rates increase with the shear rate. In agreement with previous experimental observation and theoretical work, we find that the critical shear rate corresponds to a Weissenberg number of order 1. Finally, we show that the viscosity of the system is not affected by the crystalline nuclei.

Competition between crystal and fibril formation in molecular mutations of amyloidogenic peptides

Amyloidogenic model peptides are invaluable for investigating assembly mechanisms in disease related amyloids and in protein folding. During aggregation, such peptides can undergo bifurcation leading to fibrils or crystals, however the mechanisms of fibril-to-crystal conversion are unclear. We navigate herein the energy landscape of amyloidogenic peptides by studying a homologous series of hexapeptides found in animal, human and disease related proteins. We observe fibril-to-crystal conversion occurring within single aggregates via untwisting of twisted ribbon fibrils possessing saddle-like curvature and cross-sectional aspect ratios approaching unity. Changing sequence, pH or concentration shifts the growth towards larger aspect ratio species assembling into stable helical ribbons possessing mean-curvature. By comparing atomistic calculations of desolvation energies for association of peptides we parameterise a kinetic model, providing a physical explanation of fibril-to-crystal interconversion. These results shed light on the self-assembly of amyloidogenic peptides, suggesting amyloid crystals, not fibrils, represent the ground state of the protein folding energy landscape.Aggregation of amyloidogenic peptides into fibrils and crystals has incidence in several amyloid-related diseases. Here, the authors investigate the origins of the fibril-to-crystal conversion in amyloidogenic hexapeptides pointing to the amyloid crystals as the ground state in the protein folding energy landscape.

Free Energies by Thermodynamic Integration Relative to an Exact Solution, Used to Find the Handedness-Switching Salt Concentration for DNA

Sets of free energy differences are useful for finding the equilibria of chemical reactions, while absolute free energies have little physical meaning. However finding the relative free energy between two macrostates by subtraction of their absolute free energies is a valuable strategy in certain important cases. We present calculations of absolute free energies of biomolecules, using a combination of the well-known Einstein molecule method (for treating the solute) with a conceptually related method of recent genesis for computing free energies of liquids (to treat the solvent and counterions). The approach is based on thermodynamic integration from a detailed atomistic model to one which is simplified but analytically solvable, thereby giving the absolute free energy as that of the tractable model plus a correction term found numerically. An example calculation giving the free energy with respect to salt concentration for the B- and Z-isomers of all-atom duplex DNA in explicit solvent and counterions is presented. The coexistence salt concentration is found with unprecedented accuracy.