Analyses of protein cores reveal fundamental differences between solution and crystal structures.

There have been several studies suggesting that protein structures solved by NMR spectroscopy and X-ray crystallography show significant differences. To understand the origin of these differences, we assembled a database of high-quality protein structures solved by both methods. We also find significant differences between NMR and crystal structures-in the root-mean-square deviations of the C α atomic positions, identities of core amino acids, backbone, and side-chain dihedral angles, and packing fraction of core residues. In contrast to prior studies, we identify the physical basis for these differences by modeling protein cores as jammed packings of amino acid-shaped particles. We find that we can tune the jammed packing fraction by varying the degree of thermalization used to generate the packings. For an athermal protocol, we find that the average jammed packing fraction is identical to that observed in the cores of protein structures solved by X-ray crystallography. In contrast, highly thermalized packing-generation protocols yield jammed packing fractions that are even higher than those observed in NMR structures. These results indicate that thermalized systems can pack more densely than athermal systems, which suggests a physical basis for the structural differences between protein structures solved by NMR and X-ray crystallography.

The Mitochondrial Peptide Humanin Targets but Does Not Denature Amyloid Oligomers in Type II Diabetes.

Mitochondrially derived peptides (MDPs) such as humanin (HN) have shown a remarkable ability to modulate neurological amyloids and apoptosis-associated proteins in cells and animal models. Recently, we found that humanin-like peptides also inhibit amyloid formation outside of neural environments in islet amyloid polypeptide (IAPP) fibrils and plaques, which are hallmarks of Type II diabetes. However, the biochemical basis for regulating amyloids through endogenous MDPs remains elusive. One hypothesis is that MDPs stabilize intermediate amyloid oligomers and discourage the formation of insoluble fibrils. To test this hypothesis, we carried out simulations and experiments to extract the dominant interactions between the S14G-HN mutant (HNG) and a diverse set of IAPP structures. Replica-exchange molecular dynamics suggests that MDPs cap the growth of amyloid oligomers. Simulations also indicate that HNG-IAPP heterodimers are 10 times more stable than IAPP homodimers, which explains the substoichiometric ability of HNG to inhibit amyloid growth. Despite this strong attraction, HNG does not denature IAPP. Instead, HNG binds IAPP near the disordered NFGAIL motif, wedging itself between amyloidogenic fragments. Shielding of NFGAIL-flanking fragments reduces the formation of parallel IAPP ß-sheets and subsequent nucleation of mature amyloid fibrils. From ThT spectroscopy and electron microscopy, we found that HNG does not deconstruct mature IAPP fibrils and oligomers, consistent with the simulations and our proposed hypothesis. Taken together, this work provides new mechanistic insight into how endogenous MDPs regulate pathological amyloid growth at the molecular level and in highly substoichiometric quantities, which can be exploited through peptidomimetics in diabetes or Alzheimer's disease.

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

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.

Regulation and aggregation of intrinsically disordered peptides.

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.

Molecularly Smooth Self-Assembled Monolayer for High-Mobility Organic Field-Effect Transistors.

Despite the need for molecularly smooth self-assembled monolayers (SAMs) on silicon dioxide surfaces (the most common dielectric surface), current techniques are limited to nonideal silane grafting. Here, we show unique bioinspired zwitterionic molecules forming a molecularly smooth and uniformly thin SAM in "water" in <1 min on various dielectric surfaces, which enables a dip-coating process that is essential for organic electronics to become reality. This monomolecular layer leads to high mobility of organic field-effect transistors (OFETs) based on various organic semiconductors and source/drain electrodes. A combination of experimental and computational techniques confirms strong adsorption (Wad > 20 mJ m-2), uniform thickness (∼0.5 or ∼1 nm) and orientation (all catechol head groups facing the oxide surface) of the "monomolecular" layers. This robust (strong adsorption), rapid, and green SAM represents a promising advancement toward the next generation of nanofabrication compared to the current nonuniform and inconsistent polysiloxane-based SAM involving toxic chemicals, long processing time (>10 h), or heat (>80 °C).

Phospholipid and Hydrocarbon Interactions with a Charged Electrode Interface.

Using a combination of molecular dynamics simulations and experiments we examined the interactions of alkanes and phospholipids at charged interfaces in order to understand how interfacial charge densities affect the association of these two representative molecules with electrodes. Consistent with theory and experiment, these model systems reveal interfacial associations mediated through a combination of Coulombic and van der Waals forces. van der Waals forces, in particular, mediate rapid binding of decane to neutral electrodes. No decane binding was observed at high surface charge densities because of interfacial water polarization, which screens hydrophobic attractions. The positively charged choline moiety of the phospholipid palmitoyloleoylphosphatidylcholine (POPC) is primarily responsible for POPC attraction by a moderately negatively charged electrode. The hydrocarbon tails of POPC interact with the hydrophobic electrode interface similarly to decane. Previously reported electrochemical results confirm these findings by demonstrating bipolar displacement currents from PC vesicles adhering to moderately negatively charged interfaces, originating from the choline interactions observed in simulations. At more negatively charged interfaces, choline-to-surface binding was stronger. In both simulations and experiments the maximal interaction of anionic PS occurs with a positively charged interface, provided that the electrostatic forces outweigh local Lennard-Jones interactions. Direct comparisons between the binding affinities measured in experiments and those obtained in simulations reveal previously unobserved atomic interactions that facilitate lipid vesicle adhesion to charged interfaces. Moreover, the implementation of a charged interface in molecular dynamics simulations provides an alternative method for the generation of large electric fields across phospholipid bilayers, especially for systems with periodic boundary conditions, and may be useful for simulations of membrane electropermeabilization.

Picosecond and Terahertz Perturbation of Interfacial Water and Electropermeabilization of Biological Membranes.

Non-thermal probing and stimulation with subnanosecond electric pulses and terahertz electromagnetic radiation may lead to new, minimally invasive diagnostic and therapeutic procedures and to methods for remote monitoring and analysis of biological systems, including plants, animals, and humans. To effectively engineer these still-emerging tools, we need an understanding of the biophysical mechanisms underlying the responses that have been reported to these novel stimuli. We show here that subnanosecond (≤500 ps) electric pulses induce action potentials in neurons and cause calcium transients in neuroblastoma-glioma hybrid cells, and we report complementary molecular dynamics simulations of phospholipid bilayers in electric fields in which membrane permeabilization occurs in less than 1 ns. Water dipoles in the interior of these model membranes respond in less than 1 ps to permeabilizing electric potentials by aligning in the direction of the field, and they re-orient at terahertz frequencies to field reversals. The mechanism for subnanosecond lipid electropore formation is similar to that observed on longer time scales-energy-minimizing intrusions of interfacial water into the membrane interior and subsequent reorganization of the bilayer into hydrophilic, conductive structures.

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.

Water influx and cell swelling after nanosecond electropermeabilization.

Pulsed electric fields are used to permeabilize cell membranes in biotechnology and the clinic. Although molecular and continuum models provide compelling representations of the mechanisms underlying this phenomenon, a clear structural link between the biomolecular transformations displayed in molecular dynamics (MD) simulations and the micro- and macroscale cellular responses observed in the laboratory has not been established. In this paper, plasma membrane electropermeabilization is characterized by exposing Jurkat T lymphoblasts to pulsed electric fields less than 10ns long (including single pulse exposures), and by monitoring the resulting osmotically driven cell swelling as a function of pulse number and pulse repetition rate. In this way, we reduce the complexity of the experimental system and lay a foundation for gauging the correspondence between measured and simulated values for water and ion transport through electropermeabilized membranes. We find that a single 10MV/m pulse of 5ns duration produces measurable swelling of Jurkat T lymphoblasts in growth medium, and we estimate from the swelling kinetics the ion and water flux that follows the electropermeabilization of the membrane. From these observations we set boundaries on the net conductance of the permeabilized membrane, and we show how this is consistent with model predictions for the conductance and areal density of nanoelectropulse-induced lipid nanopores.

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.

Humanin variant P3S is associated with longevity in APOE4 carriers and resists APOE4-induced brain pathology.

The APOE4 allele is recognized as a significant genetic risk factor to Alzheimer's disease (AD) and influences longevity. Nonetheless, some APOE4 carriers exhibit resistance to AD even in advanced age. Humanin, a mitochondrial-derived peptide comprising 24 amino acids, has variants linked to cognitive resilience and longevity. Our research uncovered a unique humanin variant, P3S, specifically enriched in centenarians with the APOE4 allele. Through in silico analyses and subsequent experimental validation, we demonstrated a strong affinity between humanin P3S and APOE4. Utilizing an APOE4-centric mouse model of amyloidosis (APP/PS1/APOE4), we observed that humanin P3S significantly attenuated brain amyloid-beta accumulation compared to the wild-type humanin. Transcriptomic assessments of mice treated with humanin P3S highlighted its potential mechanism involving the enhancement of amyloid beta phagocytosis. Additionally, in vitro studies corroborated humanin P3S's efficacy in promoting amyloid-beta clearance. Notably, in the temporal cortex of APOE4 carriers, humanin expression is correlated with genes associated with phagocytosis. Our findings suggest a role of the rare humanin variant P3S, especially prevalent among individuals of Ashkenazi descent, in mitigating amyloid beta pathology and facilitating phagocytosis in APOE4-linked amyloidosis, underscoring its significance in longevity and cognitive health among APOE4 carriers.

Nanoscale, electric field-driven water bridges in vacuum gaps and lipid bilayers.

Formation of a water bridge across the lipid bilayer is the first stage of pore formation in molecular dynamic (MD) simulations of electroporation, suggesting that the intrusion of individual water molecules into the membrane interior is the initiation event in a sequence that leads to the formation of a conductive membrane pore. To delineate more clearly the role of water in membrane permeabilization, we conducted extensive MD simulations of water bridge formation, stabilization, and collapse in palmitoyloleoylphosphatidylcholine bilayers and in water-vacuum-water systems, in which two groups of water molecules are separated by a 2.8 nm vacuum gap, a simple analog of a phospholipid bilayer. Certain features, such as the exponential decrease in water bridge initiation time with increased external electric field, are similar in both systems. Other features, such as the relationship between water bridge lifetime and the diameter of the water bridge, are quite different between the two systems. Data such as these contribute to a better and more quantitative understanding of the relative roles of water and lipid in membrane electropore creation and annihilation, facilitating a mechanism-driven development of electroporation protocols. These methods can be extended to more complex, heterogeneous systems that include membrane proteins and intracellular and extracellular membrane attachments, leading to more accurate models of living cells in electric fields.

Calcium and phosphatidylserine inhibit lipid electropore formation and reduce pore lifetime.

Molecular dynamics simulations of electroporation of homogeneous phospholipid bilayers show that the pore creation time is strongly dependent on the magnitude of the applied electric field. Here, we investigated whether heterogeneous bilayers containing phospholipids with zwitterionic and anionic headgroups exhibit a similar dependence. To facilitate this analysis we divide the life cycle of an electropore into several stages, marking the sequence of steps for pore creation and pore annihilation (restoration of the bilayer after removal of the electric field). We also report simulations of calcium binding isotherms and the effects of calcium ions on the electroporation of heterogeneous lipid bilayers. Calcium binding simulations are consistent with experimental data using a 1:2 Langmuir binding isotherm. We find that calcium ions and phosphatidylserine increase pore creation time and decrease pore annihilation time. For all systems tested, pore creation time was inversely proportional to the bilayer internal electric field.

Defining the Neuropathological Aggresome across in Silico, in Vitro, and ex Vivo Experiments.

The loss of proteostasis over the life course is associated with a wide range of debilitating degenerative diseases and is a central hallmark of human aging. When left unchecked, proteins that are intrinsically disordered can pathologically aggregate into highly ordered fibrils, plaques, and tangles (termed amyloids), which are associated with countless disorders such as Alzheimer's disease, Parkinson's disease, type II diabetes, cancer, and even certain viral infections. However, despite significant advances in protein folding and solution biophysics techniques, determining the molecular cause of these conditions in humans has remained elusive. This has been due, in part, to recent discoveries showing that soluble protein oligomers, not insoluble fibrils or plaques, drive the majority of pathological processes. This has subsequently led researchers to focus instead on heterogeneous and often promiscuous protein oligomers. Unfortunately, significant gaps remain in how to prepare, model, experimentally corroborate, and extract amyloid oligomers relevant to human disease in a systematic manner. This Review will report on each of these techniques and their successes and shortcomings in an attempt to standardize comparisons between protein oligomers across disciplines, especially in the context of neurodegeneration. By standardizing multiple techniques and identifying their common overlap, a clearer picture of the soluble neuropathological aggresome can be constructed and used as a baseline for studying human disease and aging.

Molecular Context of Dopa Influences Adhesion of Mussel-Inspired Peptides.

Improving adhesives for wet surfaces is an ongoing challenge. While the adhesive proteins of marine mussels have inspired many synthetic wet adhesives, the mechanisms of mussel adhesion are still not fully understood. Using surface forces apparatus (SFA) measurements and replica-exchange and umbrella-sampling molecular dynamics simulations, we probed the relationships between the sequence, structure, and adhesion of mussel-inspired peptides. Experimental and computational results reveal that peptides derived from mussel foot protein 3 slow (mfp-3s) containing 3,4-dihydroxyphenylalanine (Dopa), a post-translationally modified variant of tyrosine commonly found in mussel foot proteins, form adhesive monolayers on mica. In contrast, peptides with tyrosine adsorb as weakly adhesive clusters. We further considered simulations of mfp-3s derivatives on a range of hydrophobic and hydrophilic organic and inorganic surfaces (including silica, self-assembled monolayers, and a lipid bilayer) and demonstrated that the chemical character of the target surface and proximity of cationic and hydrophobic residues to Dopa affect peptide adsorption and adhesion. Collectively, our results suggest that conversion of tyrosine to Dopa in hydrophobic, sparsely charged peptides influences peptide self-association and ultimately dictates their adhesive performance.

Age-dependent ataxia and neurodegeneration caused by an αII spectrin mutation with impaired regulation of its calpain sensitivity.

The neuronal membrane-associated periodic spectrin skeleton (MPS) contributes to neuronal development, remodeling, and organization. Post-translational modifications impinge on spectrin, the major component of the MPS, but their role remains poorly understood. One modification targeting spectrin is cleavage by calpains, a family of calcium-activated proteases. Spectrin cleavage is regulated by activated calpain, but also by the calcium-dependent binding of calmodulin (CaM) to spectrin. The physiologic significance of this balance between calpain activation and substrate-level regulation of spectrin cleavage is unknown. We report a strain of C57BL/6J mice harboring a single αII spectrin point mutation (Sptan1 c.3293G > A:p.R1098Q) with reduced CaM affinity and intrinsically enhanced sensitivity to calpain proteolysis. Homozygotes are embryonic lethal. Newborn heterozygotes of either gender appear normal, but soon develop a progressive ataxia characterized biochemically by accelerated calpain-mediated spectrin cleavage and morphologically by disruption of axonal and dendritic integrity and global neurodegeneration. Molecular modeling predicts unconstrained exposure of the mutant spectrin's calpain-cleavage site. These results reveal the critical importance of substrate-level regulation of spectrin cleavage for the maintenance of neuronal integrity. Given that excessive activation of calpain proteases is a common feature of neurodegenerative disease and traumatic encephalopathy, we propose that damage to the spectrin MPS may contribute to the neuropathology of many disorders.

α-CGRP disrupts amylin fibrillization and regulates insulin secretion: implications on diabetes and migraine.

Despite being relatively benign and not an indicative signature of toxicity, fibril formation and fibrillar structures continue to be key factors in assessing the structure-function relationship in protein aggregation diseases. The inability to capture molecular cross-talk among key players at the tissue level before fibril formation greatly accounts for the missing link toward the development of an efficacious therapeutic intervention for Type II diabetes mellitus (T2DM). We show that human α-calcitonin gene-related peptide (α-CGRP) remodeled amylin fibrillization. Furthermore, while CGRP and/or amylin monomers reduce the secretion of both mouse Ins1 and Ins2 proteins, CGRP oligomers have a reverse effect on Ins1. Genetically reduced Ins2, the orthologous version of human insulin, has been shown to enhance insulin sensitivity and extend the life-span in old female mice. Beyond the mechanistic insights, our data suggest that CGRP regulates insulin secretion and lowers the risk of T2DM. Our result rationalizes how migraine might be protective against T2DM. We envision the new paradigm of CGRP : amylin interactions as a pivotal aspect for T2DM diagnostics and therapeutics. Maintaining a low level of amylin while increasing the level of CGRP could become a viable approach toward T2DM prevention and treatment.

Life cycle of an electropore: field-dependent and field-independent steps in pore creation and annihilation.

Electropermeabilization, an electric field-induced modification of the barrier functions of the cell membrane, is widely used in laboratories and increasingly in the clinic; but the mechanisms and physical structures associated with the electromanipulation of membrane permeability have not been definitively characterized. Indirect experimental observations of electrical conductance and small molecule transport as well as molecular dynamics simulations have led to models in which hydrophilic pores form in phospholipid bilayers with increased probability in the presence of an electric field. Presently available methods do not permit the direct, nanoscale examination of electroporated membranes that would confirm the existence of these structures. To facilitate the reconciliation of poration models with the observed properties of electropermeabilized lipid bilayers and cell membranes, we propose a scheme for characterizing the stages of electropore formation and resealing. This electropore life cycle, based on molecular dynamics simulations of phospholipid bilayers, defines a sequence of discrete steps in the electric field-driven restructuring of the membrane that leads to the formation of a head group-lined, aqueous pore and then, after the field is removed, to the dismantling of the pore and reassembly of the intact bilayer. Utilizing this scheme we can systematically analyze the interactions between the electric field and the bilayer components involved in pore initiation, construction and resealing. We find that the pore creation time depends strongly on the electric field gradient across the membrane interface and that the pore annihilation time is at least weakly dependent on the magnitude of the pore-initiating electric field and, in general, much longer than the pore creation time.

Using physical features of protein core packing to distinguish real proteins from decoys.

The ability to consistently distinguish real protein structures from computationally generated model decoys is not yet a solved problem. One route to distinguish real protein structures from decoys is to delineate the important physical features that specify a real protein. For example, it has long been appreciated that the hydrophobic cores of proteins contribute significantly to their stability. We used two sources to obtain datasets of decoys to compare with real protein structures: submissions to the biennial Critical Assessment of protein Structure Prediction competition, in which researchers attempt to predict the structure of a protein only knowing its amino acid sequence, and also decoys generated by 3DRobot, which have user-specified global root-mean-squared deviations from experimentally determined structures. Our analysis revealed that both sets of decoys possess cores that do not recapitulate the key features that define real protein cores. In particular, the model structures appear more densely packed (because of energetically unfavorable atomic overlaps), contain too few residues in the core, and have improper distributions of hydrophobic residues throughout the structure. Based on these observations, we developed a feed-forward neural network, which incorporates key physical features of protein cores, to predict how well a computational model recapitulates the real protein structure without knowledge of the structure of the target sequence. By identifying the important features of protein structure, our method is able to rank decoy structures with similar accuracy to that obtained by state-of-the-art methods that incorporate many additional features. The small number of physical features makes our model interpretable, emphasizing the importance of protein packing and hydrophobicity in protein structure prediction.

Targeting the Intrinsically Disordered Proteome Using Small-Molecule Ligands.

Intrinsically disordered proteins (IDPs) and regions (IDRs) make up a significant part of the proteome and facilitate a wide range of physiological and pathological functions that are only beginning to be understood. As such, they are highly attractive targets for drug development and bioengineering. However, their inability to adopt well-defined structures provides significant obstacles for developing ligands that regulate their behaviors. In this chapter, we review how the conformational flexibility of IDPs and their propensity to phase separate make them tractable targets for small-molecule manipulation. We also describe both theoretical and experimental approaches to characterize disordered proteins, including novel thermodynamic and single-molecule techniques that help identify complimentary partners of IDPs and their ability to shift protein ensembles toward preferred conformations.