Discovering pathways through ribozyme fitness landscapes using information theoretic quantification of epistasis.

The identification of catalytic RNAs is typically achieved through primarily experimental means. However, only a small fraction of sequence space can be analyzed even with high-throughput techniques. Methods to extrapolate from a limited data set to predict additional ribozyme sequences, particularly in a human-interpretable fashion, could be useful both for designing new functional RNAs and for generating greater understanding about a ribozyme fitness landscape. Using information theory, we express the effects of epistasis (i.e., deviations from additivity) on a ribozyme. This representation was incorporated into a simple model of the epistatic fitness landscape, which identified potentially exploitable combinations of mutations. We used this model to theoretically predict mutants of high activity for a self-aminoacylating ribozyme, identifying potentially active triple and quadruple mutants beyond the experimental data set of single and double mutants. The predictions were validated experimentally, with nine out of nine sequences being accurately predicted to have high activity. This set of sequences included mutants that form a previously unknown evolutionary "bridge" between two ribozyme families that share a common motif. Individual steps in the method could be examined, understood, and guided by a human, combining interpretability and performance in a simple model to predict ribozyme sequences by extrapolation.

Thermodynamic Modeling and Experimental Data Reveal That Sugars Stabilize Proteins According to an Excluded Volume Mechanism.

We present a new thermodynamic model to investigate the relative effects of excluded volume and soft interaction contributions in determining whether a cosolute will either destabilize or stabilize a protein in solution. This model is unique in considering an atomistically detailed model of the protein and accounting for the preferential accumulation/exclusion of the osmolyte molecules from the protein surface. Importantly, we use molecular dynamics simulations and experiments to validate the model. The experimental approach presents a unique means of decoupling excluded volume and soft interaction contributions using a linear polymeric series of cosolutes with different numbers of glucose subunits, from 1 (glucose) to 8 (maltooctaose), as well as an 8-mer of glucose units in the closed form (γ-CD). By studying the stabilizing effect of cosolutes along this polymeric series using lysozyme as a model protein, we validate the thermodynamic model and show that sugars stabilize proteins according to an excluded volume mechanism.

Competition between ß-Sheet and Coacervate Domains Yields Diverse Morphologies in Mixtures of Oppositely Charged Homochiral Polypeptides.

Biomolecular assembly processes involving competition between specific intermolecular interactions and thermodynamic phase instability have been implicated in a number of pathological states and technological applications of biomaterials. As a model for such processes, aqueous mixtures of oppositely charged homochiral polypeptides such as poly-l-lysine and poly-l-glutamic acid have been reported to form either ß-sheet-rich solid-like precipitates or liquid-like coacervate droplets depending on competing hydrogen bonding interactions. Herein, we report studies of polypeptide mixtures that reveal unexpectedly diverse morphologies ranging from partially coalescing and aggregated droplets to bulk precipitates, as well as a previously unreported re-entrant liquid-liquid phase separation at high polypeptide concentration and ionic strength. Combining our experimental results with all-atom molecular dynamics simulations of folded polypeptide complexes reveals a concentration dependence of ß-sheet-rich secondary structure, whose relative composition correlates with the observed macroscale morphologies of the mixtures. These results elucidate a crucial balance of interactions that are important for controlling morphology during coacervation in these and potentially similar biologically relevant systems.

Mechanical Behavior of Octopus Egg Tethers Composed of Topologically Constrained, Tandemly Repeated EGF Domains.

Whether and how intramolecular crosslinks in polymeric materials contribute to mechanical properties is debated in both experimental and theoretical arenas. The tethering threads of Octopus bimaculoides egg cases provide a rare window to investigate this question in a biomaterial. The only detectable component of the load-bearing fibers in octopus threads is a 135 kDa protein, octovafibrin, comprising 29 tandem repeats of epidermal growth factor (EGF) each of which contains 3 intramolecular disulfide linkages. The N- and C-terminal C-type lectins mediate linear end-to-end octovafibrin self-assembly. Mechanical testing of threads shows that the regularly spaced disulfide linkages result in improved stiffness, toughness, and energy dissipation. In response to applied loads, molecular dynamics and X-ray scattering show that EGF-like domains deform by recruiting two hidden length ß-sheet structures nested between the disulfides. The results of this study further the understanding of intramolecular crosslinking in polymers and provide a foundation for the mechanical contributions of EGF domains to the extracellular matrix.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis.

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aß, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Tachykinin Neuropeptides and Amyloid ß (25-35) Assembly: Friend or Foe?

Amyloid ß (Aß) protein is responsible for Alzheimer's disease, and one of its important fragments, Aß(25-35), is found in the brain and has been shown to be neurotoxic. Tachykinin neuropeptides, including Neuromedin K (NK), Kassinin, and Substance P, have been reported to reduce Aß(25-35)'s toxicity in cells even though they share similar primary structures with Aß(25-35). Here, we seek to understand the molecular mechanisms of how these peptides interact with Aß(25-35) and to shed light on why some peptides with similar primary structures are toxic and others nontoxic. We use both experimental and computational methods, including ion mobility mass spectrometry and enhanced-sampling replica-exchange molecular dynamics simulations, to study the aggregation pathways of Aß(25-35), NK, Kassinin, Substance P, and mixtures of the latter three with Aß(25-35). NK and Substance P were observed to remove the higher-order oligomers (i.e., hexamers and dodecamers) of Aß(25-35), which are related to its toxicity, although Substance P did so more slowly. In contrast, Kassinin was found to promote the formation of these higher-order oligomers. This result conflicts with what is expected and is elaborated on in the text. We also observe that even though they have significant structural homology with Aß(25-35), NK, Kassinin, and Substance P do not form hexamers with a ß-sheet structure like Aß(25-35). The hexamer structure of Aß(25-35) has been identified as a cylindrin, and this structure has been strongly correlated to toxic species. The reasons why the three tachykinin peptides behave so differently when mixed with Aß(25-35) are discussed.

Molecular design of self-coacervation phenomena in block polyampholytes.

Coacervation is a common phenomenon in natural polymers and has been applied to synthetic materials systems for coatings, adhesives, and encapsulants. Single-component coacervates are formed when block polyampholytes exhibit self-coacervation, phase separating into a dense liquid coacervate phase rich in the polyampholyte coexisting with a dilute supernatant phase, a process implicated in the liquid-liquid phase separation of intrinsically disordered proteins. Using fully fluctuating field-theoretic simulations using complex Langevin sampling and complementary molecular-dynamics simulations, we develop molecular design principles to connect the sequenced charge pattern of a polyampholyte with its self-coacervation behavior in solution. In particular, the lengthscale of charged blocks and number of connections between oppositely charged blocks are shown to have a dramatic effect on the tendency to phase separate and on the accessible chain conformations. The field and particle-based simulation results are compared with analytical predictions from the random phase approximation (RPA) and postulated scaling relationships. The qualitative trends are mostly captured by the RPA, but the approximation fails catastrophically at low concentration.

Catalytic Cross Talk between Key Peptide Fragments That Couple Alzheimer's Disease with Amyotrophic Lateral Sclerosis.

Protein aggregation is a common feature in prominent neurodegenerative diseases, usually thought to be due to the assembly of a single peptide or protein. Recent studies have challenged this notion and suggested several proteins may be involved in promoting and amplifying disease. For example, the TDP-43 protein associated with Amyotrophic Lateral Sclerosis has been found in the brain along with Aß assemblies associated with Alzheimer's disease, and those patients that show the presence of TDP-43 are 10 times more likely to demonstrate cognitive impairment compared to TDP-43-negative Alzheimer's patients. Here we examine the interactions between the amyloidogenic core of TDP-43, TDP-43307-319, and a neurotoxic physiologically observed fragment of Aß, Aß25-35. Utilizing ion mobility mass spectrometry in concert with atomic force microscopy and molecular dynamics simulations, we investigate which oligomers are involved in seeding aggregation across these two different protein systems and gain insight into which structures initiate and result from these interactions. Studies were conducted by mixing Aß25-35 with the toxic wild type TDP-43307-319 peptide and with the nontoxic synthetic TDP-43307-319 mutant, G314V. Our findings identify a strong catalytic effect of TDP-43307-319 WT monomer in the acceleration of Aß25-35 aggregation to its toxic cylindrin and ß barrel forms. This observation is unprecedented in both its speed and specificity. Interestingly, the nontoxic G314V mutant of TDP-43307-319 and dimers or higher order oligomers of WT TDP-43307-319 do not promote aggregation of Aß25-35 but rather dissociate preformed toxic higher order oligomers of Aß25-35. Reasons for these very different behaviors are reported.

CORE-MD II: A fast, adaptive, and accurate enhanced sampling method.

In this paper, we present a fast and adaptive correlation guided enhanced sampling method (CORE-MD II). The CORE-MD II technique relies, in part, on partitioning of the entire pathway into short trajectories that we refer to as instances. The sampling within each instance is accelerated by adaptive path-dependent metadynamics simulations. The second part of this approach involves kinetic Monte Carlo (kMC) sampling between the different states that have been accessed during each instance. Through the combination of the partition of the total simulation into short non-equilibrium simulations and the kMC sampling, the CORE-MD II method is capable of sampling protein folding without any a priori definitions of reaction pathways and additional parameters. In the validation simulations, we applied the CORE-MD II on the dialanine peptide and the folding of two peptides: TrpCage and TrpZip2. In a comparison with long time equilibrium Molecular Dynamics (MD), 1 µs replica exchange MD (REMD), and CORE-MD I simulations, we find that the level of convergence of the CORE-MD II method is improved by a factor of 8.8, while the CORE-MD II method reaches acceleration factors of ∼120. In the CORE-MD II simulation of TrpZip2, we observe the formation of the native state in contrast to the REMD and the CORE-MD I simulations. The method is broadly applicable for MD simulations and is not restricted to simulations of protein folding or even biomolecules but also applicable to simulations of protein aggregation, protein signaling, or even materials science simulations.

Heightened Cold-Denaturation of Proteins at the Ice-Water Interface.

The process of freezing proteins is widely used in applications ranging from processing and storage of biopharmaceuticals to cryo-EM analysis of protein complexes. The formation of an ice-water interface is a critical destabilization factor for the protein, which can be offset by the use of cryo-protectants. Using molecular dynamics simulation, we demonstrate that the presence of the ice-water interface leads to a lowering of the free-energy barrier for unfolding, resulting in rapid unfolding of the protein. The unfolding process does not require direct adsorption of the protein to the surface, but is rather mediated by nearby liquid molecules that show an increased tendency for hydrating nonpolar groups. The observed enhancement in the cold denaturation process upon ice formation can be mitigated by addition of glucose, which acts as a cryoprotectant through preferential exclusion from side chains of the protein.

CORE-MD, a path correlated molecular dynamics simulation method.

We present an enhanced Molecular Dynamics (MD) simulation method, which is free from the requirement of a priori structural information of the system. The technique is capable of folding proteins with very low computational effort and requires only an energy parameter. The path correlated MD (CORE-MD) method uses the autocorrelation of the path integral over the reduced action and propagates the system along the history dependent path correlation. We validate the new technique in simulations of the conformational landscapes of dialanine and the TrpCage mini-peptide. We find that the novel method accelerates the sampling by three orders of magnitude and observe convergence of the conformational sampling in both cases. We conclude that the new method is broadly applicable for the enhanced sampling in MD simulations. The CORE-MD algorithm reaches a high accuracy compared with long time equilibrium MD simulations.

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.

Small ion effects on self-coacervation phenomena in block polyampholytes.

Self-coacervation is a phenomenon in which a solution of polyampholytes spontaneously phase separates into a dense liquid coacervate phase, rich in the polyampholyte, coexisting with a dilute supernatant phase. Such coacervation results in the formation of membraneless organelles in vivo and has further been applied industrially as synthetic encapsulants and coatings. It has been suggested that coacervation is primarily driven by the entropy gain from releasing counter-ions upon complexation. Using fully fluctuating field-theoretic simulations employing complex Langevin sampling and complementary molecular dynamics simulations, we have determined that the small ions contribute only weakly to the self-coacervation behavior of charge-symmetric block polyampholytes in solution. Salt partitioning between the supernatant and coacervate is also found to be negligible in the weak-binding regime at low electrostatic strengths. Asymmetries in charge distribution along the polyampholytes can cause net-charges that lead to "tadpole" configurations in dilute solution and the suppression of phase separation at low salt content. The field and particle-based simulation results are compared with analytical predictions from the random phase approximation (RPA) and postulated scaling relationships. The qualitative trends are mostly captured by the RPA, but the approximation fails at low concentration.

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.

Phospholipid Chain Interactions with Cholesterol Drive Domain Formation in Lipid Membranes.

Cholesterol is a key component of eukaryotic membranes, but its role in cellular biology in general and in lipid rafts in particular remains controversial. Model membranes are used extensively to determine the phase behavior of ternary mixtures of cholesterol, a saturated lipid, and an unsaturated lipid with liquid-ordered and liquid-disordered phase coexistence. Despite many different experiments that determine lipid-phase diagrams, we lack an understanding of the molecular-level driving forces for liquid phase coexistence in bilayers with cholesterol. Here, we use atomistic molecular dynamics computer simulations to address the driving forces for phase coexistence in ternary lipid mixtures. Domain formation is directly observed in a long-timescale simulation of a mixture of 1,2-distearoyl-sn-glycero-3-phosphocholine, unsaturated 1,2-dilinoleoyl-sn-glycero-3-phosphocholine, and cholesterol. Free-energy calculations for the exchange of the saturated and unsaturated lipids between the ordered and disordered phases give insight into the mixing behavior. We show that a large energetic contribution to domain formation is favorable enthalpic interactions of the saturated lipid in the ordered phase. This favorable energy for forming an ordered, cholesterol-rich phase is opposed by a large unfavorable entropy. Martini coarse-grained simulations capture the unfavorable free energy of mixing but do not reproduce the entropic contribution because of the reduced representation of the phospholipid tails. Phospholipid tails and their degree of unsaturation are key energetic contributors to lipid phase separation.

Trimethylamine N-oxide Counteracts Urea Denaturation by Inhibiting Protein-Urea Preferential Interaction.

Osmolytes are small organic molecules that can modulate the stability and function of cellular proteins by altering the chemical environment of the cell. Some of these osmolytes work in conjunction, via mechanisms that are poorly understood. An example is the naturally occurring protein-protective osmolyte trimethylamine N-oxide (TMAO) that stabilizes cellular proteins in marine organisms against the detrimental denaturing effects of another naturally occurring osmolyte, urea. From a computational standpoint, our understanding of this counteraction mechanism is hampered by the fact that existing force fields fail to capture the correct balance of TMAO and urea interactions in ternary solutions. Using molecular dynamics simulations and Kirkwood-Buff theory of solutions, we have developed an optimized force field that reproduces experimental Kirkwood-Buff integrals. We show through the study of two model systems, a 15-residue polyalanine chain and the R2-fragment (273GKVQIINKKLDL284) of the Tau protein, that TMAO can counteract the denaturing effects of urea by inhibiting protein-urea preferential interaction. The extent to which counteraction can occur is seen to depend heavily on the amino acid composition of the peptide.

Conformational investigation of the structure-activity relationship of GdFFD and its analogues on an achatin-like neuropeptide receptor of Aplysia californica involved in the feeding circuit.

Proteins and peptides in nature are almost exclusively made from l-amino acids, and this is even more absolute in the metazoan. With the advent of modern bioanalytical techniques, however, previously unappreciated roles for d-amino acids in biological processes have been revealed. Over 30 d-amino acid containing peptides (DAACPs) have been discovered in animals where at least one l-residue has been isomerized to the d-form via an enzyme-catalyzed process. In Aplysia californica, GdFFD and GdYFD (the lower-case letter "d" indicates a d-amino acid residue) modulate the feeding behavior by activating the Aplysia achatin-like neuropeptide receptor (apALNR). However, little is known about how the three-dimensional conformation of DAACPs influences activity at the receptor, and the role that d-residues play in these peptide conformations. Here, we use a combination of computational modeling, drift-tube ion-mobility mass spectrometry, and receptor activation assays to create a simple model that predicts bioactivities for a series of GdFFD analogs. Our results suggest that the active conformations of GdFFD and GdYFD are similar to their lowest energy conformations in solution. Our model helps connect the predicted structures of GdFFD analogs to their activities, and highlights a steric effect on peptide activity at position 1 on the GdFFD receptor apALNR. Overall, these methods allow us to understand ligand-receptor interactions in the absence of high-resolution structural data.