Zachary A. Levine

Regulation and aggregation of intrinsically disordered peptides

Significance The microtubule-regulating protein tau is a prototypical intrinsically disordered protein (IDP) that plays an important physiological role in the human body; however, aggregates of tau are a pathological hallmark of Alzheimer’s disease. Here we demonstrate through simulations and experiments with an aggregating tau fragment that cosolvent interactions can significantly affect the balance between hydrogen bonds and salt bridge formation in IDPs, subsequently determining their preferred conformations. These subtle perturbations can dramatically shift IDPs from compact ensembles to extended ones, thereby influencing aggregate formation. These results lend considerable insight into the biophysics of the regulation and aggregation of IDPs. Intrinsically disordered proteins (IDPs) are a unique class of proteins that have no stable native structure, a feature that allows them to adopt a wide variety of extended and compact conformations that facilitate a large number of vital physiological functions. One of the most well-known IDPs is the microtubule-associated tau protein, which regulates microtubule growth in the nervous system. However, dysfunctions in tau can lead to tau oligomerization, fibril formation, and neurodegenerative disease, including Alzheimer’s disease. Using a combination of simulations and experiments, we explore the role of osmolytes in regulating the conformation and aggregation propensities of the R2/wt peptide, a fragment of tau containing the aggregating paired helical filament (PHF6*). We show that the osmolytes urea and trimethylamine N-oxide (TMAO) shift the population of IDP monomer structures, but that no new conformational ensembles emerge. Although urea halts aggregation, TMAO promotes the formation of compact oligomers (including helical oligomers) through a newly proposed mechanism of redistribution of water around the perimeter of the peptide. We put forth a “superposition of ensembles” hypothesis to rationalize the mechanism by which IDP structure and aggregation is regulated in the cell.

Electric Field-Driven Water Dipoles: Nanoscale Architecture of Electroporation

Electroporation is the formation of permeabilizing structures in the cell membrane under the influence of an externally imposed electric field. The resulting increased permeability of the membrane enables a wide range of biological applications, including the delivery of normally excluded substances into cells. While electroporation is used extensively in biology, biotechnology, and medicine, its molecular mechanism is not well understood. This lack of knowledge limits the ability to control and fine-tune the process. In this article we propose a novel molecular mechanism for the electroporation of a lipid bilayer based on energetics analysis. Using molecular dynamics simulations we demonstrate that pore formation is driven by the reorganization of the interfacial water molecules. Our energetics analysis and comparisons of simulations with and without the lipid bilayer show that the process of poration is driven by field-induced reorganization of water dipoles at the water-lipid or water-vacuum interfaces into more energetically favorable configurations, with their molecular dipoles oriented in the external field. Although the contributing role of water in electroporation has been noted previously, here we propose that interfacial water molecules are the main players in the process, its initiators and drivers. The role of the lipid layer, to a first-order approximation, is then reduced to a relatively passive barrier. This new view of electroporation simplifies the study of the problem, and opens up new opportunities in both theoretical modeling of the process and experimental research to better control or to use it in new, innovative ways.

Determination of biomembrane bending moduli in fully atomistic simulations.

The bilayer bending modulus (Kc) is one of the most important physical constants characterizing lipid membranes, but precisely measuring it is a challenge, both experimentally and computationally. Experimental measurements on chemically identical bilayers often differ depending upon the techniques employed, and robust simulation results have previously been limited to coarse-grained models (at varying levels of resolution). This Communication demonstrates the extraction of Kc from fully atomistic molecular dynamics simulations for three different single-component lipid bilayers (DPPC, DOPC, and DOPE). The results agree quantitatively with experiments that measure thermal shape fluctuations in giant unilamellar vesicles. Lipid tilt, twist, and compression moduli are also reported.

Surface force measurements and simulations of mussel-derived peptide adhesives on wet organic surfaces

Significance The need for bio-inspired wet adhesives has significantly increased in the past few decades (e.g., for dental and medical transplants, coronary artery coatings, cell encapsulants, etc.). However, the molecular basis behind catechol-facilitated adhesion to organic surfaces remains unclear, thus hindering synthesis and optimization of novel underwater adhesives. The present combined experimental and theoretical study reconciles bioadhesion measurements of novel catechol-containing peptides to self-assembled monolayers (SAMs) with all-atom molecular dynamics simulations, yielding a comprehensive framework that explicitly identifies the basis for underwater adhesion. Simulations and surface forces apparatus measurements agree with one another, and both approaches show strong peptide adhesion to hydrophobic SAMs, and weak peptide adhesion to hydrophilic SAMs, providing a starting point for the development of next-generation underwater glues. Translating sticky biological molecules—such as mussel foot proteins (MFPs)—into synthetic, cost-effective underwater adhesives with adjustable nano- and macroscale characteristics requires an intimate understanding of the glue’s molecular interactions. To help facilitate the next generation of aqueous adhesives, we performed a combination of surface forces apparatus (SFA) measurements and replica-exchange molecular dynamics (REMD) simulations on a synthetic, easy to prepare, Dopa-containing peptide (MFP-3s peptide), which adheres to organic surfaces just as effectively as its wild-type protein analog. Experiments and simulations both show significant differences in peptide adsorption on CH3-terminated (hydrophobic) and OH-terminated (hydrophilic) self-assembled monolayers (SAMs), where adsorption is strongest on hydrophobic SAMs because of orientationally specific interactions with Dopa. Additional umbrella-sampling simulations yield free-energy profiles that quantitatively agree with SFA measurements and are used to extract the adhesive properties of individual amino acids within the context of MFP-3s peptide adhesion, revealing a delicate balance between van der Waals, hydrophobic, and electrostatic forces.

Simulations of disordered proteins and systems with conformational heterogeneity

Intrinsically disordered proteins (IDPs) and protein regions can facilitate a wide variety of complex physiological processes such as binding, signaling, and formation of membraneless organelles. They can however also play pathological roles by aggregating into cytotoxic oligomers and fibrils. Characterizing the structure and function of disordered proteins is an onerous task, primarily because these proteins adopt transient structures, which are difficult to capture in experiments. Simulations have emerged as a powerful tool for interpreting and augmenting experimental measurements of IDPs. In this review we focus on computer simulations of disordered protein structures, functions, assemblies, and emerging questions that, taken together, give an overview of the field as it exists today.

Quantitative Limits on Small Molecule Transport via the Electropermeome — Measuring and Modeling Single Nanosecond Perturbations

The detailed molecular mechanisms underlying the permeabilization of cell membranes by pulsed electric fields (electroporation) remain obscure despite decades of investigative effort. To advance beyond descriptive schematics to the development of robust, predictive models, empirical parameters in existing models must be replaced with physics- and biology-based terms anchored in experimental observations. We report here absolute values for the uptake of YO-PRO-1, a small-molecule fluorescent indicator of membrane integrity, into cells after a single electric pulse lasting only 6 ns. We correlate these measured values, based on fluorescence microphotometry of hundreds of individual cells, with a diffusion-based geometric analysis of pore-mediated transport and with molecular simulations of transport across electropores in a phospholipid bilayer. The results challenge the “drift and diffusion through a pore” model that dominates conventional explanatory schemes for the electroporative transfer of small molecules into cells and point to the necessity for a more complex model.

To What Extent Does Surface Hydrophobicity Dictate Peptide Folding and Stability near Surfaces

Protein-surface interactions are ubiquitous in both the cellular setting and in modern bioengineering devices, but how such interactions impact protein stability is not well understood. We investigate the folding of the GB1 hairpin peptide in the presence of self-assembled monolayers and graphite like surfaces using replica exchange molecular dynamics simulations. By varying surface hydrophobicity, and decoupling direct protein-surface interactions from water-mediated interactions, we show that surface wettability plays a surprisingly minor role in dictating protein stability. For both the β-hairpin GB1 and the helical miniprotein TrpCage, adsorption and stability is largely dictated by the nature of the direct chemical interactions between the protein and the surface. Independent of the surface hydrophobicity profile, strong protein-surface interactions destabilize the folded structure while weak interactions stabilize it.

Significant Performance Enhancement of Polymer Resins by Bioinspired Dynamic Bonding

Marine mussels use catechol-rich interfacial mussel foot proteins (mfps) as primers that attach to mineral surfaces via hydrogen, metal coordination, electrostatic, ionic, or hydrophobic bonds, creating a secondary surface that promotes bonding to the bulk mfps. Inspired by this biological adhesive primer, it is shown that a ≈1 nm thick catecholic single-molecule priming layer increases the adhesion strength of crosslinked polymethacrylate resin on mineral surfaces by up to an order of magnitude when compared with conventional primers such as noncatecholic silane- and phosphate-based grafts. Molecular dynamics simulations confirm that catechol groups anchor to a variety of mineral surfaces and shed light on the binding mode of each molecule. Here, a ≈50% toughness enhancement is achieved in a stiff load-bearing polymer network, demonstrating the utility of mussel-inspired bonding for processing a wide range of polymeric interfaces, including structural, load-bearing materials.

Studying the Early Stages of Protein Aggregation Using Replica Exchange Molecular Dynamics Simulations.

The simulation of protein aggregation poses several computational challenges due to the disparate time and lengths scales that are involved. This chapter focuses on the use of atomistically detailed simulations to probe the initial steps of aggregation, with an emphasis on the Tau peptide as a model system, run under a replica exchange molecular dynamics protocol.