Multiscale Molecular Dynamics Simulations of an Active Self-Assembling Material.

Inspired by the adaptability observed in biological materials, self-assembly processes have attracted significant interest for their potential to yield novel materials with unique properties. However, experimental methods have often fallen short in capturing the molecular details of the assembly process. In this study, we employ a multiscale molecular dynamics simulation approach, complemented by NMR quantification, to investigate the mechanism of self-assembly in a redox-fueled bioinspired system. Contrary to conventional assumptions, we have uncovered a significant role played by the monomer precursor in the assembly process, with its presence varying with concentration and the extent of conversion of the monomer to the dimer. Experimental confirmation through NMR quantification underscores the concentration-dependent incorporation of monomers into the fibrous structures. Furthermore, our simulations also shed light on the diverse intermolecular interactions, including T-shaped and parallel π stacking, as well as hydrogen bonds, in stabilizing the aggregates. Overall, the open conformation of the dimer is preferred within these aggregates. However, inside the aggregates, the distribution of conformations shifts slightly to the closed conformation compared to on the surface. These findings contribute to the growing field of bioinspired materials science by providing valuable mechanistic and structural insights to guide the design and development of self-assembling materials with biomimetic functionalities.

A ß-barrel-like tetramer formed by a ß-hairpin derived from Aß.

ß-Hairpins formed by the ß-amyloid peptide Aß are building blocks of Aß oligomers. Three different alignments of ß-hairpins have been observed in the structures of Aß oligomers or fibrils. Differences in ß-hairpin alignment likely contribute to the heterogeneity of Aß oligomers and thus impede their study at high-resolution. Here, we designed, synthesized, and studied a series of ß-hairpin peptides derived from Aß12-40 in one of these three alignments and investigated their solution-phase assembly and folding. These assays reveal the formation of tetramers and octamers that are stabilized by intermolecular hydrogen bonding interactions between Aß residues 12-14 and 38-40 as part of an extended ß-hairpin conformation. X-ray crystallographic studies of one peptide from this series reveal the formation of ß-barrel-like tetramers and octamers that are stabilized by edge-to-edge hydrogen bonding and hydrophobic packing. Dye-leakage and caspase 3/7 activation assays using tetramer and octamer forming peptides from this series reveal membrane-damaging and apoptotic properties. A molecular dynamics simulation of the ß-barrel-like tetramer embedded in a lipid bilayer shows membrane disruption and water permeation. The tetramers and octamers described herein provide additional models of how Aß may assemble into oligomers and supports the hypothesis that ß-hairpin alignment and topology may contribute directly to oligomer heterogeneity.

Native American ataxia medicines rescue ataxia-linked mutant potassium channel activity via binding to the voltage sensing domain.

There are currently no drugs known to rescue the function of Kv1.1 voltage-gated potassium channels carrying loss-of-function sequence variants underlying the inherited movement disorder, Episodic Ataxia 1 (EA1). The Kwakwaka'wakw First Nations of the Pacific Northwest Coast used Fucus gardneri (bladderwrack kelp), Physocarpus capitatus (Pacific ninebark) and Urtica dioica (common nettle) to treat locomotor ataxia. Here, we show that extracts of these plants enhance wild-type Kv1.1 current, especially at subthreshold potentials. Screening of their constituents revealed that gallic acid and tannic acid similarly augment wild-type Kv1.1 current, with submicromolar potency. Crucially, the extracts and their constituents also enhance activity of Kv1.1 channels containing EA1-linked sequence variants. Molecular dynamics simulations reveal that gallic acid augments Kv1.1 activity via a small-molecule binding site in the extracellular S1-S2 linker. Thus, traditional Native American ataxia treatments utilize a molecular mechanistic foundation that can inform small-molecule approaches to therapeutically correcting EA1 and potentially other Kv1.1-linked channelopathies.

Insights into the solubility of γ D-crystallin from multiscale atomistic simulations.

The molecular basis underlying the rich phase behavior of globular proteins remains poorly understood. We use atomistic multiscale molecular simulations to model the solution-state conformational dynamics and interprotein interactions of γ D-crystallin and its P23T-R36S mutant, which drastically limits the protein solubility, at both infinite dilution and at a concentration where the mutant fluid phase and crystalline phase coexist. We find that while the mutant conserves the protein fold, changes to the surface exposure of residues in the neighborhood of residue-36 enhance protein-protein interactions and develop specific protein-protein contacts found in the protein crystal lattice.

Network Hamiltonian Models for Unstructured Protein Aggregates, with Application to γD-Crystallin.

Network Hamiltonian models (NHMs) are a framework for topological coarse-graining of protein-protein interactions, in which each node corresponds to a protein, and edges are drawn between nodes representing proteins that are noncovalently bound. Here, this framework is applied to aggregates of γD-crystallin, a structural protein of the eye lens implicated in cataract disease. The NHMs in this study are generated from atomistic simulations of equilibrium distributions of wild-type and the cataract-causing variant W42R in solution, performed by Wong, E. K.; Prytkova, V.; Freites, J. A.; Butts, C. T.; Tobias, D. J. Molecular Mechanism of Aggregation of the Cataract-Related γD-Crystallin W42R Variant from Multiscale Atomistic Simulations. Biochemistry2019, 58 (35), 3691-3699. Network models are shown to successfully reproduce the aggregate size and structure observed in the atomistic simulation, and provide information about the transient protein-protein interactions therein. The system size is scaled from the original 375 monomers to a system of 10000 monomers, revealing a lowering of the upper tail of the aggregate size distribution of the W42R variant. Extrapolation to higher and lower concentrations is also performed. These results provide an example of the utility of NHMs for coarse-grained simulation of protein systems, as well as their ability to scale to large system sizes and high concentrations, reducing computational costs while retaining topological information about the system.

Ca2+ -dependent interactions between lipids and the tumor-targeting peptide pHLIP.

Cancerous tissues undergo extensive changes to their cellular environments that differentiate them from healthy tissues. These changes include changes in extracellular pH and Ca2+ concentrations, and the exposure of phosphatidylserine (PS) to the extracellular environment, which can modulate the interaction of peptides and proteins with the plasma membrane. Deciphering the molecular mechanisms of such interactions is critical for advancing the knowledge-based design of cancer-targeting molecular tools, such as pH-low insertion peptide (pHLIP). Here, we explore the effects of PS, Ca2+ , and peptide protonation state on the interactions of pHLIP with lipid membranes. Cellular studies demonstrate that exposed PS on the plasma membrane promotes pHLIP targeting. The magnitude of this effect is dependent on extracellular Ca2+ concentration, indicating that divalent cations play an important role in pHLIP targeting in vivo. The targeting mechanism is further explored with a combination of fluorescence and circular dichroism experiments in model membranes and microsecond-timescale all-atom molecular dynamics simulations. Our results demonstrate that Ca2+ is engaged in coupling peptide-lipid interactions in the unprotonated transmembrane conformation of pHLIP. The simulations reveal that while the pH-induced insertion leads to a strong depletion of PS around pHLIP, the Ca2+ -induced insertion has the opposite effect. Thus, extracellular levels of Ca2+ are crucial to linking cellular changes in membrane lipid composition with the selective targeting and insertion of pHLIP. The characterized Ca2+ -dependent coupling between pHLIP sidechains and PS provides atomistic insights into the general mechanism for lipid-coupled regulation of protein-membrane insertion by divalent cations.

Effects of Cardiolipin on the Conformational Dynamics of Membrane-Anchored Bcl-xL.

The anti-apoptotic protein Bcl-xL regulates apoptosis by preventing the permeation of the mitochondrial outer membrane by pro-apoptotic pore-forming proteins, which release apoptotic factors into the cytosol that ultimately lead to cell death. Two different membrane-integrated Bcl-xL constructs have been identified: a membrane-anchored and a membrane-inserted conformation. Here, we use molecular dynamics simulations to study the effect of the mitochondrial specific lipid cardiolipin and the protein protonation state on the conformational dynamics of membrane-anchored Bcl-xL. The analysis reveals that the protonation state of the protein and cardiolipin content of the membrane modulate the orientation of the soluble head region (helices α1 through α7) and hence the exposure of its BH3-binding groove, which is required for its interaction with pro-apoptotic proteins.

Anomalous Diffusion of Peripheral Membrane Signaling Proteins from All-Atom Molecular Dynamics Simulations.

Peripheral membrane proteins bind transiently to membrane surfaces as part of many signaling pathways. The bound proteins perform two-dimensional (2-D) diffusion on the membrane surface during the recruitment function. To better understand the interplay between the 2-D diffusion of these protein domains and their membrane binding modes, we performed multimicrosecond all-atom molecular dynamics simulations of two regulatory domains, a C2 domain and a pleckstrin homology (PH) domain, in their experimentally determined bound configuration to a lipid bilayer. The protein bound configurations are preserved throughout the simulation trajectories. Both protein domains exhibit anomalous diffusion with distinct features in their dynamics that reflect their different modes of binding. An analysis of their diffusive behavior reveals common features with the diffusion of lipid molecules in lipid bilayers, suggesting that the 2-D motion of protein domains bound to the membrane surface is modulated by the viscoelastic nature of the lipid bilayer.

Thermodynamics and Mechanism of the Membrane Permeation of Hv1 Channel Blockers.

The voltage-gated proton channel Hv1 mediates efflux of protons from the cell. Hv1 integrally contributes to various physiological processes including pH homeostasis and the respiratory burst of phagocytes. Inhibition of Hv1 may provide therapeutic avenues for the treatment of inflammatory diseases, breast cancer, and ischemic brain damage. In this work, we investigate two prototypical Hv1 inhibitors, 2-guanidinobenzimidazole (2GBI), and 5-chloro-2-guanidinobenzimidazole (GBIC), from an experimentally screened class of guanidine derivatives. Both compounds block proton conduction by binding the same site located on the intracellular side of the channel. However, when added to the extracellular medium, the compounds strongly differ in their ability to inhibit proton conduction, suggesting substantial differences in membrane permeability. Here, we compute the potential of mean force for each compound to permeate through the membrane using atomistic molecular dynamics simulations with the adaptive biasing force method. Our results rationalize the putative distinction between these two blockers with respect to their abilities to permeate the cellular membrane.

Insights on small molecule binding to the Hv1 proton channel from free energy calculations with molecular dynamics simulations.

Hv1 is a voltage-gated proton channel whose main function is to facilitate extrusion of protons from the cell. The development of effective channel blockers for Hv1 can lead to new therapeutics for the treatment of maladies related to Hv1 dysfunction. Although the mechanism of proton permeation in Hv1 remains to be elucidated, a series of small molecules have been discovered to inhibit Hv1. Here, we computed relative binding free energies of a prototypical Hv1 blocker on a model of human Hv1 in an open state. We used alchemical free energy perturbation techniques based on atomistic molecular dynamics simulations. The results support our proposed open state model and shed light on the preferred tautomeric state of the channel blocker. This work lays the groundwork for future studies on adapting the blocker molecule for more effective inhibition of Hv1.

Voltage-dependent structural models of the human Hv1 proton channel from long-timescale molecular dynamics simulations.

The voltage-gated Hv1 proton channel is a ubiquitous membrane protein that has roles in a variety of cellular processes, including proton extrusion, pH regulation, production of reactive oxygen species, proliferation of cancer cells, and increased brain damage during ischemic stroke. A crystal structure of an Hv1 construct in a putative closed state has been reported, and structural models for the channel open state have been proposed, but a complete characterization of the Hv1 conformational dynamics under an applied membrane potential has been elusive. We report structural models of the Hv1 voltage-sensing domain (VSD), both in a hyperpolarized state and a depolarized state resulting from voltage-dependent conformational changes during a 10-µs-timescale atomistic molecular dynamics simulation in an explicit membrane environment. In response to a depolarizing membrane potential, the S4 helix undergoes an outward displacement, leading to changes in the VSD internal salt-bridge network, resulting in a reshaping of the permeation pathway and a significant increase in hydrogen bond connectivity throughout the channel. The total gating charge displacement associated with this transition is consistent with experimental estimates. Molecular docking calculations confirm the proposed mechanism for the inhibitory action of 2-guanidinobenzimidazole (2GBI) derived from electrophysiological measurements and mutagenesis. The depolarized structural model is also consistent with the formation of a metal bridge between residues located in the core of the VSD. Taken together, our results suggest that these structural models are representative of the closed and open states of the Hv1 channel.

Human αB-crystallin discriminates between aggregation-prone and function-preserving variants of a client protein.

BACKGROUND: The eye lens crystallins are highly soluble proteins that are required to last the lifespan of an organism due to low protein turnover in the lens. Crystallin aggregation leads to formation of light-scattering aggregates known as cataract. The G18V mutation of human γS-crystallin (γS-G18V), which is associated with childhood-onset cataract, causes structural changes throughout the N-terminal domain and increases aggregation propensity. The holdase chaperone protein αB-crystallin does not interact with wild-type γS-crystallin, but does bind its G18V variant. The specific molecular determinants of αB-crystallin binding to client proteins is incompletely charcterized. Here, a new variant of γS, γS-G18A, was created to test the limits of αB-crystallin selectivity. METHODS: Molecular dynamics simulations were used to investigate the structure and dynamics of γS-G18A. The overall fold of γS-G18A was assessed by circular dichroism (CD) spectroscopy and intrinsic tryptophan fluorescence. Its thermal unfolding temperature and aggregation propensity were characterized by CD and DLS, respectively. Solution-state NMR was used to characterize interactions between αB-crystallin and γS-G18A. RESULTS: γS-G18A exhibits minimal structural changes, but has compromised thermal stability relative to γS-WT. The placement of alanine, rather than valine, at this highly conserved glycine position produces minor changes in hydrophobic surface exposure. However, human αB-crystallin does not bind the G18A variant, in contrast to previous observations for γS-G18V, which aggregates at physiological temperature. CONCLUSIONS: αB-crystallin is capable of distinguishing between aggregation-prone and function-preserving variants, and recognizing the transient unfolding or minor conformers that lead to aggregation in the disease-related variant. GENERAL SIGNIFICANCE: Human αB-crystallin distinguishes between highly similar variants of a structural crystallin, binding the cataract-related γS-G18V variant, but not the function-preserving γS-G18A variant, which is monomeric at physiological temperature.

Molecular Mechanism of Aggregation of the Cataract-Related γD-Crystallin W42R Variant from Multiscale Atomistic Simulations.

The mechanisms leading to aggregation of the crystallin proteins of the eye lens remain largely unknown. We use atomistic multiscale molecular simulations to model the solution-state conformational dynamics of γD-crystallin and its cataract-related W42R variant at both infinite dilution and physiologically relevant concentrations. We find that the W42R variant assumes a distinct conformation in solution that leaves the Greek key domains of the native fold largely unaltered but lacks the hydrophobic interdomain interface that is key to the stability of wild-type γD-crystallin. At physiologically relevant concentrations, exposed hydrophobic regions in this alternative conformation become primary sites for enhanced interprotein interactions leading to large-scale aggregation.

Voltage-Dependent Profile Structures of a Kv-Channel via Time-Resolved Neutron Interferometry.

Available experimental techniques cannot determine high-resolution three-dimensional structures of membrane proteins under a transmembrane voltage. Hence, the mechanism by which voltage-gated cation channels couple conformational changes within the four voltage sensor domains, in response to either depolarizing or polarizing transmembrane voltages, to opening or closing of the pore domain's ion channel remains unresolved. Single-membrane specimens, composed of a phospholipid bilayer containing a vectorially oriented voltage-gated K+ channel protein at high in-plane density tethered to the surface of an inorganic multilayer substrate, were developed to allow the application of transmembrane voltages in an electrochemical cell. Time-resolved neutron reflectivity experiments, enhanced by interferometry enabled by the multilayer substrate, were employed to provide directly the low-resolution profile structures of the membrane containing the vectorially oriented voltage-gated K+ channel for the activated, open and deactivated, closed states of the channel under depolarizing and hyperpolarizing transmembrane voltages applied cyclically. The profile structures of these single membranes were dominated by the voltage-gated K+ channel protein because of the high in-plane density. Importantly, the use of neutrons allowed the determination of the voltage-dependent changes in both the profile structure of the membrane and the distribution of water within the profile structure. These two key experimental results were then compared to those predicted by three computational modeling approaches for the activated, open and deactivated, closed states of three different voltage-gated K+ channels in hydrated phospholipid bilayer membrane environments. Of the three modeling approaches investigated, only one state-of-the-art molecular dynamics simulation that directly predicted the response of a voltage-gated K+ channel within a phospholipid bilayer membrane to applied transmembrane voltages by utilizing very long trajectories was found to be in agreement with the two key experimental results provided by the time-resolved neutron interferometry experiments.

Experimental and Simulation Studies of Aquaporin 0 Water Permeability and Regulation.

We begin with the history of aquaporin zero (AQP0), the most prevalent membrane protein in the eye lens, from the early days when AQP0 was a protein of unknown function known as Major Intrinsic Protein 26. We progress through its joining the aquaporin family as a water channel in its own right and discuss how regulation of its water permeability by pH and calcium came to be discovered experimentally and linked to lens homeostasis and development. We review the development of molecular dynamics (MD) simulations of lipid bilayers and membrane proteins, including aquaporins, with an emphasis on simulation studies that have elucidated the mechanisms of water conduction, selectivity, and proton exclusion by aquaporins in general. We also review experimental and theoretical progress toward understanding why mammalian AQP0 has a lower water permeability than other aquaporins and the evolution of our present understanding of how its water permeability is regulated by pH and calcium. Finally, we discuss how MD simulations have elucidated the nature of lipid interactions with AQP0.

Cooperativity and allostery in aquaporin 0 regulation by Ca2.

Aquaporin 0 (AQP0) is essential for eye lens homeostasis as is regulation of its water permeability by Ca2+, which occurs through interactions with calmodulin (CaM), but the underlying molecular mechanisms are not well understood. Here, we use molecular dynamics (MD) simulations on the microsecond timescale under an osmotic gradient to explicitly model water permeation through the AQP0 channel. To identify any structural features that are specific to water permeation through AQP0, we also performed simulations of aquaporin 1 (AQP1) and a pure mixed lipid bilayer under the same conditions. The relative single-channel water osmotic permeability coefficients (pf) calculated from all of our simulations are in reasonable agreement with experiment. Our simulations allowed us to characterize the dynamics of the key structural elements that modulate the diffusion of water single-files through the AQP0 and AQP1 pores. We find that CaM binding influences the collective dynamics of the whole AQP0 tetramer, promoting the closing of both the extracellular and intracellular gates by inducing cooperativity between neighboring subunits.

Role of Conformational Flexibility in Monte Carlo Simulations of Many-Protein Systems.

Efficient computational modeling of biological systems characterized by high concentrations of macromolecules often relies on rigid-body Brownian Dynamics or Monte Carlo (MC) simulations. However, the accuracy of rigid-body models is limited by the fixed conformation of the simulated biomolecules. Multi-conformation Monte Carlo (mcMC) simulations of protein solutions incorporate conformational flexibility via a conformational swap trial move within a predetermined library of discrete protein structures, thereby alleviating artifacts arising from the use of a single protein conformation. Here, we investigate the impact of the number of distinct protein structures in the conformational library and the extent of conformational sampling used in its generation on structural observables computed from simulations of hen egg white lysozyme (HEWL), human γD-Crystallin, and bovine γB-Crystallin solutions. We find that the importance of specific protocols for the construction of the protein structure library is strongly dependent on the nature of the simulated system.

Refining Protein Penetration into the Lipid Bilayer Using Fluorescence Quenching and Molecular Dynamics Simulations: The Case of Diphtheria Toxin Translocation Domain.

Dynamic disorder of the lipid bilayer presents a challenge for establishing structure-function relationships in membranous systems. The resulting structural heterogeneity is especially evident for peripheral and spontaneously inserting membrane proteins, which are not constrained by the well-defined transmembrane topology and exert their action in the context of intimate interaction with lipids. Here, we propose a concerted approach combining depth-dependent fluorescence quenching with Molecular Dynamics simulation to decipher dynamic interactions of membrane proteins with the lipid bilayers. We apply this approach to characterize membrane-mediated action of the diphtheria toxin translocation domain. First, we use a combination of the steady-state and time-resolved fluorescence spectroscopy to characterize bilayer penetration of the NBD probe selectively attached to different sites of the protein into membranes containing lipid-attached nitroxyl quenching groups. The constructed quenching profiles are analyzed with the Distribution Analysis methodology allowing for accurate determination of transverse distribution of the probe. The results obtained for 12 NBD-labeled single-Cys mutants are consistent with the so-called Open-Channel topology model. The experimentally determined quenching profiles for labeling sites corresponding to L350, N373, and P378 were used as initial constraints for positioning TH8-9 hairpin into the lipid bilayer for Molecular Dynamics simulation. Finally, we used alchemical free energy calculations to characterize protonation of E362 in soluble translocation domain and membrane-inserted conformation of its TH8-9 fragment. Our results indicate that membrane partitioning of the neutral E362 is more favorable energetically (by ~ 6 kcal/mol), but causes stronger perturbation of the bilayer, than the charged E362.

Specific cation effects at aqueous solution-vapor interfaces: Surfactant-like behavior of Li+ revealed by experiments and simulations.

It is now well established by numerous experimental and computational studies that the adsorption propensities of inorganic anions conform to the Hofmeister series. The adsorption propensities of inorganic cations, such as the alkali metal cations, have received relatively little attention. Here we use a combination of liquid-jet X-ray photoelectron experiments and molecular dynamics simulations to investigate the behavior of K+ and Li+ ions near the interfaces of their aqueous solutions with halide ions. Both the experiments and the simulations show that Li+ adsorbs to the aqueous solution-vapor interface, while K+ does not. Thus, we provide experimental validation of the "surfactant-like" behavior of Li+ predicted by previous simulation studies. Furthermore, we use our simulations to trace the difference in the adsorption of K+ and Li+ ions to a difference in the resilience of their hydration shells.

YidC Insertase of Escherichia coli: Water Accessibility and Membrane Shaping.

The YidC/Oxa1/Alb3 family of membrane proteins function to insert proteins into membranes in bacteria, mitochondria, and chloroplasts. Recent X-ray structures of YidC from Bacillus halodurans and Escherichia coli revealed a hydrophilic groove that is accessible from the lipid bilayer and the cytoplasm. Here, we explore the water accessibility within the conserved core region of the E. coli YidC using in vivo cysteine alkylation scanning and molecular dynamics (MD) simulations of YidC in POPE/POPG membranes. As expected from the structure, YidC possesses an aqueous membrane cavity localized to the membrane inner leaflet. Both the scanning data and the MD simulations show that the lipid-exposed transmembrane helices 3, 4, and 5 are short, leading to membrane thinning around YidC. Close examination of the MD data reveals previously unrecognized structural features that are likely important for protein stability and function.