Ligand Lipophilicity Determines Molecular Mechanisms of Nanoparticle Adsorption to Lipid Bilayers

The interactions of ligand-functionalized nanoparticles with the cell membrane affect cellular uptake, cytotoxicity, and related behaviors, but relating these interactions to ligand properties remains challenging. In this work, we perform coarse-grained molecular dynamics simulations to study how the adsorption of ligand-functionalized cationic gold nanoparticles (NPs) to a single-component lipid bilayer (as a model cell membrane) is influenced by ligand end group lipophilicity. A set of 2 nm diameter NPs, each coated with a monolayer of organic ligands that differ only in their end groups, was simulated to mimic NPs recently studied experimentally. Metadynamics calculations were performed to determine key features of the free energy landscape for adsorption as a function of the distance of the NP from the bilayer and the number of NP-lipid contacts. These simulations revealed that NP adsorption is thermodynamically favorable for all NPs due to the extraction of lipids from the bilayer and into the NP monolayer. To resolve ligand-dependent differences in adsorption behavior, string method calculations were performed to compute minimum free energy pathways for adsorption. These calculations revealed a surprising nonmonotonic dependence of the free energy barrier for adsorption on ligand end group lipophilicity. Large free energy barriers are predicted for the least lipophilic end groups because favorable NP-lipid contacts are initiated only through the unfavorable protrusion of lipid tail groups out of the bilayer. The smallest free energy barriers are predicted for end groups of intermediate lipophilicity which promote NP-lipid contacts by intercalating within the bilayer. Unexpectedly, large free energy barriers are also predicted for the most lipophilic end groups which remain sequestered within the ligand monolayer rather than intercalating within the bilayer. These trends are broadly in agreement with past experimental measurements and reveal how subtle variations in ligand lipophilicity dictate adsorption mechanisms and associated kinetics by influencing the interplay of lipid-ligand interactions.

Effects of solvent conditions on the self-assembly of heterotrimeric collagen-like peptide (CLP) triple helices: a coarse-grained simulation study

We perform coarse-grained (CG) molecular dynamics (MD) simulations to investigate the self-assembly of collagen-like peptide (CLP) triple helices into fibrillar structures and percolated networks as a function of solvent quality. The focus of this study is on CLP triple helices whose strands are different lengths (i.e., heterotrimers), leading to dangling 'sticky ends'. These 'sticky ends' are segments of the CLP strands that have unbonded hydrogen-bonding donor/acceptor sites that drive heterotrimeric CLP triple helices to physically associate with one another, leading to assembly into higher-order structures. We use a validated CG model for CLP in implicit solvent and capture varying solvent quality through changing strength of attraction between CG beads representing the amino acids in the CLP strands. Our CG MD simulations show that, at lower CLP concentrations, CLP heterotrimers assemble into fibrils and, at higher CLP concentrations, into percolated networks. At higher concentrations, decreasing solvent quality causes (i) the formation of heterogeneous network structures with a lower degree of branching at network junctions and (ii) increases in the diameter of network strands and pore sizes. We also observe a nonmonotonic effect of solvent quality on distances between network junctions due to the balance between heterotrimer end-end associations driven by hydrogen bonding and side-side associations driven by worsening solvent quality. Below the percolation threshold, we observe that decreasing solvent quality leads to the formation of fibrils composed of multiple aligned CLP triple helices, while the number of 'sticky ends' governs the spatial extent (radius of gyration) of the assembled fibrils.

Combining Molecular Dynamics Simulations and Biophysical Characterization to Investigate Protein-Specific Excipient Effects on Reteplase during Freeze Drying

We performed molecular dynamics simulations of Reteplase in the presence of different excipients to study the stabilizing mechanisms and to identify the role of excipients during freeze drying. To simulate the freeze-drying process, we divided the process into five distinct steps: (i) protein-excipient formulations at room temperature, (ii) the ice-growth process, (iii)-(iv) the partially solvated and fully dried formulations, and (v) the reconstitution. Furthermore, coarse-grained (CG) simulations were employed to explore the protein-aggregation process in the presence of arginine. By using a coarse-grained representation, we could observe the collective behavior and interactions between protein molecules during the aggregation process. The CG simulations revealed that the presence of arginine prevented intermolecular interactions of the catalytic domain of Reteplase, thus reducing the aggregation propensity. This suggests that arginine played a stabilizing role by interacting with protein-specific regions. From the freeze-drying simulations, we could identify several protein-specific events: (i) collapse of the domain structure, (ii) recovery of the drying-induced damages during reconstitution, and (iii) stabilization of the local aggregation-prone region via direct interactions with excipients. Complementary to the simulations, we employed nanoDSF, size-exclusion chromatography, and CD spectroscopy to investigate the effect of the freeze-drying process on the protein structure and stability.

Implicit-Solvent Coarse-Grained Simulations of Linear-Dendritic Block Copolymer Micelles

The design of nanoassemblies can be conveniently achieved by tuning the strength of the hydrophobic interactions of block copolymers in selective solvents. These block copolymer micelles form supramolecular aggregates, which have attracted great attention in the area of drug delivery and imaging in biomedicine due to their easy-to-tune properties and straightforward large-scale production. In the present work, we have investigated the micellization process of linear-dendritic block copolymers in order to elucidate the effect of branching on the micellar properties. We focus on block copolymers formed by linear hydrophobic blocks attached to either dendritic neutral or charged hydrophilic blocks. We have implemented a simple protocol for determining the equilibrium micellar size, which permits the study of linear-dendritic block copolymers in a wide range of block morphologies in an efficient and parallelizable manner. We have explored the impact of different topological and charge properties of the hydrophilic blocks on the equilibrium micellar properties and compared them to predictions from self-consistent field theory and scaling theory. We have found that, at higher degrees of branching in the corona and for short polymer chains, excluded volume interactions strongly influence the micellar aggregation as well as their effective charge.

Differential interactions of resting, activated, and desensitized states of the α7 nicotinic acetylcholine receptor with lipidic modulators

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that modulates neuronal excitability, largely by allowing Ca2+ permeation. Agonist binding promotes transition from a resting state to an activated state, and then rapidly to a desensitized state. Recently, cryogenic electron microscopy (cryo-EM) structures of the human α7 receptor in nanodiscs were reported in multiple conformations. These were selectively stabilized by inhibitory, activating, or potentiating compounds. However, the functional annotation of these structures and their differential interactions with unresolved lipids and ligands remain incomplete. Here, we characterized their ion permeation, membrane interactions, and ligand binding using computational electrophysiology, free-energy calculations, and coarse-grained molecular dynamics. In contrast to nonconductive structures in apparent resting and desensitized states, the structure determined in the presence of the potentiator PNU-120596 was consistent with an activated state permeable to Ca2+. Transition to this state was associated with compression and rearrangement of the membrane, particularly in the vicinity of the peripheral MX helix. An intersubunit transmembrane site was implicated in selective binding of either PNU-120596 in the activated state or cholesterol in the desensitized state. This substantiates functional assignment of all three lipid-embedded α7-receptor structures with ion-permeation simulations. It also proposes testable models of their state-dependent interactions with lipophilic ligands, including a mechanism for allosteric modulation at the transmembrane subunit interface.

Hybrid Nanocrystals of Small Molecules and Chemically Disordered Polymers

Organic crystals formed by small molecules can be highly functional but are often brittle or insoluble structures with limited possibilities for use or processing from a liquid phase. A possible solution is the nanoscale integration of polymers into organic crystals without sacrificing long-range order and therefore function. This enables the organic crystals to benefit from the advantageous mechanical and chemical properties of the polymeric component. We report here on a strategy in which small molecules cocrystallize with side chains of chemically disordered polymers to create hybrid nanostructures containing a highly ordered lattice. Synchrotron X-ray scattering, absorption spectroscopy, and coarse-grained molecular dynamics simulations reveal that the polymer backbones form an "exo-crystalline" layer of disordered chains that wrap around the nanostructures, becoming a handle for interesting properties. The morphology of this "hybrid bonding polymer" nanostructure is dictated by the competition between the polymers' entropy and the enthalpy of the lattice allowing for control over the aspect ratio of the nanocrystal by changing the degree of polymer integration. We observed that nanostructures with an exo-crystalline layer of polymer exhibit enhanced fracture strength, self-healing capacity, and dispersion in water, which benefits their use as light-harvesting assemblies in photocatalysis. Guided by computation, future work could further explore these hybrid nanostructures as components for functional materials.

Combining Coarse-Grained Simulations and Single Molecule Analysis Reveals a Three-State Folding Model of the Guanidine-II Riboswitch

Riboswitch RNAs regulate gene expression by conformational changes induced by environmental conditions and specific ligand binding. The guanidine-II riboswitch is proposed to bind the small molecule guanidinium and to subsequently form a kissing loop interaction between the P1 and P2 hairpins. While an interaction was shown for isolated hairpins in crystallization and electron paramagnetic resonance experiments, an intrastrand kissing loop formation has not been demonstrated. Here, we report the first evidence of this interaction in cis in a ligand and Mg²⁺ dependent manner. Using single-molecule FRET spectroscopy and detailed structural information from coarse-grained simulations, we observe and characterize three interconvertible states representing an open and kissing loop conformation as well as a novel Mg²⁺ dependent state for the guanidine-II riboswitch from E. coli. The results further substantiate the proposed switching mechanism and provide detailed insight into the regulation mechanism for the guanidine-II riboswitch class. Combining single molecule experiments and coarse-grained simulations therefore provides a promising perspective in resolving the conformational changes induced by environmental conditions and to yield molecular insights into RNA regulation.

Impact of collagen-like peptide (CLP) heterotrimeric triple helix design on helical thermal stability and hierarchical assembly: a coarse-grained molecular dynamics simulation study

Collagen-like peptides (CLP) are multifunctional materials garnering a lot of recent interest from the biomaterials community due to their hierarchical assembly and tunable physicochemical properties. In this work, we present a computational study that links the design of CLP heterotrimers to the thermal stability of the triple helix and their self-assembly into fibrillar aggregates and percolated networks. Unlike homotrimeric helices, the CLP heterotrimeric triple helices in this study are made of CLP strands of different chain lengths that result in 'sticky' ends with available hydrogen bonding groups. These 'sticky' ends at one end or both ends of the CLP heterotrimer then facilitate inter-helix hydrogen bonding leading to self-assembly into fibrils (clusters) and percolated networks. We consider the cases of three sticky end lengths - two, four, and six repeat units - present entirely on one end or split between two ends of the CLP heterotrimer. We observe in CLP heterotrimer melting curves generated using coarse grained Langevin dynamics simulations at low CLP concentration that increasing sticky end length results in lower melting temperatures for both one and two sticky ended CLP designs. At higher CLP concentrations, we observe non-monotonic trends in cluster sizes with increasing sticky end length with one sticky end but not for two sticky ends with the same number of available hydrogen bonding groups as the one sticky end; this nonmonotonicity stems from the formation of turn structures stabilized by hydrogen bonds at the single, sticky end for sticky end lengths greater than four repeat units. With increasing CLP concentration, heterotrimers also form percolated networks with increasing sticky end length with a minimum sticky end length of four repeat units required to observe percolation. Overall, this work informs the design of thermoresponsive, peptide-based biomaterials with desired morphologies using strand length and dispersity as a handle for tuning thermal stability and formation of supramolecular structures.

Characterization of the Interfacial Orientation and Molecular Conformation in a Glass-Forming Organic Semiconductor

The ability to control structure in molecular glasses has enabled them to play a key role in modern technology; in particular, they are ubiquitous in organic light-emitting diodes. While the interplay between bulk structure and optoelectronic properties has been extensively investigated, few studies have examined molecular orientation near buried interfaces despite its critical role in emergent functionality. Direct, quantitative measurements of buried molecular orientation are inherently challenging, and many methods are insensitive to orientation in amorphous soft matter or lack the necessary spatial resolution. To overcome these challenges, we use polarized resonant soft X-ray reflectivity (p-RSoXR) to measure nanometer-resolved, molecular orientation depth profiles of vapor-deposited thin films of an organic semiconductor Tris(4-carbazoyl-9-ylphenyl)amine (TCTA). Our depth profiling approach characterizes the vertical distribution of molecular orientation and reveals that molecules near the inorganic substrate and free surface have a different, nearly isotropic orientation compared to those of the anisotropic bulk. Comparison of p-RSoXR results with near-edge X-ray absorption fine structure spectroscopy and optical spectroscopies reveals that TCTA molecules away from the interfaces are predominantly planar, which may contribute to their attractive charge transport qualities. Buried interfaces are further investigated in a TCTA bilayer (each layer deposited under separate conditions resulting in different orientations) in which we find a narrow interface between orientationally distinct layers extending across ≈1 nm. Coupling this result with molecular dynamics simulations provides additional insight into the formation of interfacial structure. This study characterizes the local molecular orientation at various types of buried interfaces in vapor-deposited glasses and provides a foundation for future studies to develop critical structure-function relationships.

Unified Nanotechnology Format: One Way to Store Them All

The domains of DNA and RNA nanotechnology are steadily gaining in popularity while proving their value with various successful results, including biosensing robots and drug delivery cages. Nowadays, the nanotechnology design pipeline usually relies on computer-based design (CAD) approaches to design and simulate the desired structure before the wet lab assembly. To aid with these tasks, various software tools exist and are often used in conjunction. However, their interoperability is hindered by a lack of a common file format that is fully descriptive of the many design paradigms. Therefore, in this paper, we propose a Unified Nanotechnology Format (UNF) designed specifically for the biomimetic nanotechnology field. UNF allows storage of both design and simulation data in a single file, including free-form and lattice-based DNA structures. By defining a logical and versatile format, we hope it will become a widely accepted and used file format for the nucleic acid nanotechnology community, facilitating the future work of researchers and software developers. Together with the format description and publicly available documentation, we provide a set of converters from existing file formats to simplify the transition. Finally, we present several use cases visualizing example structures stored in UNF, showcasing the various types of data UNF can handle.

Channels with Helical Modulation Display Stereospecific Sensitivity for Chiral Superstructures

By means of coarse-grained molecular dynamics simulations, we explore chiral sensitivity of confining spaces modelled as helical channels to chiral superstructures represented by polymer knots. The simulations show that helical channels exhibit stereosensitivity to chiral knots localized on linear chains by effect of external pulling force and also to knots embedded on circular chains. The magnitude of the stereoselective effect is stronger for torus knots, the effect is weaker in the case of twist knots, and amphichiral knots do exhibit no chiral effects. The magnitude of the effect can be tuned by the so-far investigated radius of the helix, the pitch of the helix and the strength of the pulling force. The model is aimed to simulate and address a range of practical situations that may occur in experimental settings such as designing of nanotechnological devices for the detection of topological state of molecules, preparation of new gels with tailor made stereoselective properties, or diffusion of knotted DNA in biological conditions.

Computational Approaches to Explore Bacterial Toxin Entry into the Host Cell

Many bacteria secrete toxic protein complexes that modify and disrupt essential processes in the infected cell that can lead to cell death. To conduct their action, these toxins often need to cross the cell membrane and reach a specific substrate inside the cell. The investigation of these protein complexes is essential not only for understanding their biological functions but also for the rational design of targeted drug delivery vehicles that must navigate across the cell membrane to deliver their therapeutic payload. Despite the immense advances in experimental techniques, the investigations of the toxin entry mechanism have remained challenging. Computer simulations are robust complementary tools that allow for the exploration of biological processes in exceptional detail. In this review, we first highlight the strength of computational methods, with a special focus on all-atom molecular dynamics, coarse-grained, and mesoscopic models, for exploring different stages of the toxin protein entry mechanism. We then summarize recent developments that are significantly advancing our understanding, notably of the glycolipid–lectin (GL-Lect) endocytosis of bacterial Shiga and cholera toxins. The methods discussed here are also applicable to the design of membrane-penetrating nanoparticles and the study of the phenomenon of protein phase separation at the surface of the membrane. Finally, we discuss other likely routes for future development.

Synergistic Effect of Chemical Penetration Enhancers on Lidocaine Permeability Revealed by Coarse-Grained Molecular Dynamics Simulations

The search for new formulations for transdermal drug delivery (TDD) is an important field in medicine and cosmetology. Molecules with specific physicochemical properties which can increase the permeability of active ingredients across the stratum corneum (SC) are called chemical penetration enhancers (CPEs), and it was shown that some CPEs can act synergistically. In this study, we performed coarse-grained (CG) molecular dynamics (MD) simulations of the lidocaine delivery facilitated by two CPEs-linoleic acid (LA) and ethanol-through the SC model membrane containing cholesterol, N-Stearoylsphingosine (DCPE), and behenic acid. In our simulations, we probed the effects of individual CPEs as well as their combination on various properties of the SC membrane and the lidocaine penetration across it. We demonstrated that the addition of both CPEs decreases the membrane thickness and the order parameters of the DPCE hydrocarbon chains. Moreover, LA also enhances diffusion of the SC membrane components, especially cholesterol. The estimated potential of mean force (PMF) profiles for the lidocaine translocation across SC in the presence/absence of two individual CPEs and their combination demonstrated that while ethanol lowers the free energy barrier for lidocaine to enter SC, LA decreases the depth of the free energy minima for lidocaine inside SC. These two effects supposedly result in synergistic penetration enhancement of drugs. Altogether, the present simulations provide a detailed molecular picture of CPEs' action and their synergistic effect on the penetration of small molecular weight therapeutics that can be beneficial for the design of novel drug and cosmetics formulations.

The effect of polydispersity, shape fluctuations and curvature on small unilamellar vesicle small-angle X-ray scattering curves

It is a challenge to distinguish the effect of shape fluctuations and size polydispersity on experimental small-angle X-ray scattering curves of small unilamellar vesicles. Here it is shown that both effects have distinguishable spectral patterns, and an efficient simulation tool is presented for simulating and analysing experimental data. The importance of curvature-induced electron-density profile asymmetry for estimating the vesicle size from SAXS scattering curves is also demonstrated.

Small unilamellar vesicles (20–100 nm diameter) are model systems for strongly curved lipid membranes, in particular for cell organelles. Routinely, small-angle X-ray scattering (SAXS) is employed to study their size and electron-density profile (EDP). Current SAXS analysis of small unilamellar vesicles (SUVs) often employs a factorization into the structure factor (vesicle shape) and the form factor (lipid bilayer electron-density profile) and invokes additional idealizations: (i) an effective polydispersity distribution of vesicle radii, (ii) a spherical vesicle shape and (iii) an approximate account of membrane asymmetry, a feature particularly relevant for strongly curved membranes. These idealizations do not account for thermal shape fluctuations and also break down for strong salt- or protein-induced deformations, as well as vesicle adhesion and fusion, which complicate the analysis of the lipid bilayer structure. Presented here are simulations of SAXS curves of SUVs with experimentally relevant size, shape and EDPs of the curved bilayer, inferred from coarse-grained simulations and elasticity considerations, to quantify the effects of size polydispersity, thermal fluctuations of the SUV shape and membrane asymmetry. It is observed that the factorization approximation of the scattering intensity holds even for small vesicle radii (∼30 nm). However, the simulations show that, for very small vesicles, a curvature-induced asymmetry arises in the EDP, with sizeable effects on the SAXS curve. It is also demonstrated that thermal fluctuations in shape and the size polydispersity have distinguishable signatures in the SAXS intensity. Polydispersity gives rise to low-q features, whereas thermal fluctuations predominantly affect the scattering at larger q, related to membrane bending rigidity. Finally, it is shown that simulation of fluctuating vesicle ensembles can be used for analysis of experimental SAXS curves.

Combining simulations and experiments for the molecular engineering of multifunctional collagen mimetic peptide-based materials

Assembling peptides allow the creation of structurally complex materials, where amino acid selection influences resulting properties. We present a synergistic approach of experiments and simulations for examining the influence of natural and non-natural amino acid substitutions via incorporation of charged residues and a reactive handle on the thermal stability and assembly of multifunctional collagen mimetic peptides (CMPs). Experimentally, we observed inclusion of charged residues significantly decreased the melting temperature of CMP triple helices with further destabilization upon inclusion of the reactive handle. Atomistic simulations of a single CMP triple helix in explicit water showed increased residue-level and helical structural fluctuations caused by the inclusion of the reactive handle; however, these atomistic simulations cannot be used to predict changes in CMP melting transition. Coarse-grained (CG) simulations of CMPs at experimentally relevant solution conditions, showed, qualitatively, the same trends as experiments in CMP melting transition temperature with CMP design. These simulations show that when charged residues are included electrostatic repulsions significantly destabilize the CMP triple helix and that an additional inclusion of a reactive handle does not significantly change the melting transition. Based on findings from both experiments and simulations, the sequence design was refined for increased CMP triple helix thermal stability, and the reactive handle was utilized for the incorporation of the assembled CMPs within covalently crosslinked hydrogels. Overall, a unique approach was established for predicting stability of CMP triple helices for various sequences prior to synthesis, providing molecular insights for sequence design towards the creation of bulk nanostructured soft biomaterials.

Entropic Competition between Supercoiled and Torsionally Relaxed Chromatin Fibers Drives Loop Extrusion through Pseudo-Topologically Bound Cohesin

Simple Summary

Chromatin dynamics and chromatin structure are a two-way relationship governed by polymer physics and active biological processes. Thanks to the research in the field of computational biology and modeling, computer simulations became indispensable in studying these complex relationships. It is now generally accepted that looped structures occurring in the intermediate range of ordering of chromatin are formed by a loop extrusion mechanism involving specialized proteins (structural maintenance complexes or SMCs). Although the motor activity of SMCs has been speculated for a long time, the motor activity of cohesin was discovered only recently (Davidson 2019). While evidence of the cohesin’s motor activity is missing, other mechanisms that could efficiently drive the loop extrusion without motor activity of SMCs have been discovered by computer simulations. These mechanisms account for transcriptionally driven loop extrusion or entropically driven loop extrusion by osmotic pressure. In our previous model, we have shown that the cohesin in handcuffed conformation can be pushed mechanically by emerging plectoneme formed during transcription, exerting pressure on the joint section of handcuffs. In the current work, we use coarse-grained molecular simulation to further explore the extrusion driven by supercoiling while employing much lower levels of supercoiling. Moreover, recent works favor non-topological binding of cohesin on fibers, which would solve a range of topological problems while bypassing other molecular machinery sitting on DNA. We show by means of computer simulations that supercoiling can drive loop extrusion without taking advantage of mechanic push on the joint section of cohesin handcuffs. As such, the work addresses current problems in molecular biology and employs advanced methods and original solutions in the study.

Abstract

We propose a model for cohesin-mediated loop extrusion, where the loop extrusion is driven entropically by the energy difference between supercoiled and torsionally relaxed chromatin fibers. Different levels of negative supercoiling are controlled by varying imposed friction between the cohesin ring and the chromatin fiber. The speed of generation of negative supercoiling by RNA polymerase associated with TOP1 is kept constant and corresponds to 10 rotations per second. The model was tested by coarse-grained molecular simulations for a wide range of frictions between 2 to 200 folds of that of generic fiber and the surrounding medium. The higher friction allowed for the accumulation of higher levels of supercoiling, while the resulting extrusion rate also increased. The obtained extrusion rates for the given range of investigated frictions were between 1 and 10 kbps, but also a saturation of the rate at high frictions was observed. The calculated contact maps indicate a qualitative improvement obtained at lower levels of supercoiling. The fits of mathematical equations qualitatively reproduce the loop sizes and levels of supercoiling obtained from simulations and support the proposed mechanism of entropically driven extrusion. The cohesin ring is bound on the fibers pseudo-topologically, and the model suggests that the topological binding is not necessary.

Residue-Level Contact Reveals Modular Domain Interactions of PICK1 Are Driven by Both Electrostatic and Hydrophobic Forces

PICK1 is a multi-domain scaffolding protein that is uniquely comprised of both a PDZ domain and a BAR domain. While previous experiments have shown that the PDZ domain and the linker positively regulate the BAR domain and the C-terminus negatively regulates the BAR domain, the details of internal regulation mechanisms are unknown. Molecular dynamics (MD) simulations have been proven to be a useful tool in revealing the intramolecular interactions at atomic-level resolution. PICK1 performs its biological functions in a dimeric form which is extremely computationally demanding to simulate with an all-atom force field. Here, we use coarse-grained MD simulations to expose the key residues and driving forces in the internal regulations of PICK1. While the PDZ and BAR domains do not form a stable complex, our simulations show the PDZ domain preferentially interacting with the concave surface of the BAR domain over other BAR domain regions. Furthermore, our simulations show that the short helix in the linker region can form interactions with the PDZ domain. Our results reveal that the surface of the βB-βC loop, βC strand, and αA-βD loop of the PDZ domain can form a group of hydrophobic interactions surrounding the linker helix. These interactions are driven by hydrophobic forces. In contrast, our simulations reveal a very dynamic C-terminus that most often resides on the convex surface of the BAR domain rather than the previously suspected concave surface. These interactions are driven by a combination of electrostatic and hydrophobic interactions.

Amyloid Aggregation under the Lens of Liquid-Liquid Phase Separation

Increasing experiments suggest that amyloid peptides can undergo liquid-liquid phase separation (LLPS) before the formation of amyloid fibrils. However, the exact role of LLPS in amyloid aggregation at the molecular level remains elusive. Here, we investigated the LLPS and amyloid fibrillization of a coarse-grained peptide, capable of capturing fundamental properties of amyloid aggregation over a wide range of concentrations in molecular dynamics simulations. On the basis of the Flory-Huggins theory of polymer solutions, we determined the binodal and spinodal concentrations of LLPS in the low-concentration regime, ϕ_BL and ϕ_SL, respectively. Only at concentrations above ϕ_BL, peptides formed metastable or stable oligomers corresponding to the high-density liquid phase (HDLP) in LLPS, out of which the nucleated conformational conversion to fibril seeds occurred. Below ϕ_SL, the HDLP was metastable and transient, and the subsequent fibrillization process followed the traditional nucleation and elongation mechanisms. Only above ϕ_SL, the HDLP became stable, and the initial fibril nucleation and growth were governed by the high local peptide concentrations. The predicted saturation of amyloid aggregation half-times with increasing peptide concentration to a constant, instead of the traditional power-law scaling to zero, was confirmed by simulations and by a thioflavin-T kinetic assay and the transmission electron microscopy of islet amyloid polypeptide (IAPP) aggregation. Our study provides a unified picture of amyloid aggregation for a wide range of concentrations within the framework of LLPS, which may help us better understand the etiology of amyloid diseases, where the amyloid protein concentration can vary by ∼9 orders of magnitude depending on the organ location and facilitate the engineering of novel amyloid-based functional materials.

Coarse-Grained Simulations Suggest the Epsin N-Terminal Homology Domain Can Sense Membrane Curvature without Its Terminal Amphipathic Helix

Nanoscale membrane curvature is a common feature in cell biology required for functions such as endocytosis, exocytosis and cell migration. These processes require the cytoskeleton to exert forces on the membrane to deform it. Cytosolic proteins contain specific motifs which bind to the membrane, connecting it to the internal cytoskeletal machinery. These motifs often bind charged phosphatidylinositol phosphate lipids present in the cell membrane which play significant roles in signaling. These lipids are important for membrane deforming processes, such as endocytosis, but much remains unknown about their role in the sensing of membrane nanocurvature by protein domains. Using coarse-grained molecular dynamics simulations, we investigated the interaction of a model curvature active protein domain, the epsin N-terminal homology domain (ENTH), with curved lipid membranes. The combination of anionic lipids (phosphatidylinositol 4,5-bisphosphate and phosphatidylserine) within the membrane, protein backbone flexibility, and structural changes within the domain were found to affect the domain's ability to sense, bind, and localize with nanoscale precision at curved membrane regions. The findings suggest that the ENTH domain can sense membrane curvature without the presence of its terminal amphipathic α helix via another structural region we have denoted as H3, re-emphasizing the critical relationship between nanoscale membrane curvature and protein function.

Physical mechanisms of amyloid nucleation on fluid membranes

Biological membranes can dramatically accelerate the aggregation of normally soluble protein molecules into amyloid fibrils and alter the fibril morphologies, yet the molecular mechanisms through which this accelerated nucleation takes place are not yet understood. Here, we develop a coarse-grained model to systematically explore the effect that the structural properties of the lipid membrane and the nature of protein-membrane interactions have on the nucleation rates of amyloid fibrils. We identify two physically distinct nucleation pathways-protein-rich and lipid-rich-and quantify how the membrane fluidity and protein-membrane affinity control the relative importance of those molecular pathways. We find that the membrane's susceptibility to reshaping and being incorporated into the fibrillar aggregates is a key determinant of its ability to promote protein aggregation. We then characterize the rates and the free-energy profile associated with this heterogeneous nucleation process, in which the surface itself participates in the aggregate structure. Finally, we compare quantitatively our data to experiments on membrane-catalyzed amyloid aggregation of α-synuclein, a protein implicated in Parkinson's disease that predominately nucleates on membranes. More generally, our results provide a framework for understanding macromolecular aggregation on lipid membranes in a broad biological and biotechnological context.

Clustering and Fibril Formation during GNNQQNY Aggregation: A Molecular Dynamics Study

The precise kinetic pathways of peptide clustering and fibril formation are not fully understood. Here we study the initial clustering kinetics and transient cluster morphologies during aggregation of the heptapeptide fragment GNNQQNY from the yeast prion protein Sup35. We use a mid-resolution coarse-grained molecular dynamics model of Bereau and Deserno to explore the aggregation pathways from the initial random distribution of free monomers to the formation of large clusters. By increasing the system size to 72 peptides we could follow directly the molecular events leading to the formation of stable fibril-like structures. To quantify those structures we developed a new cluster helicity parameter. We found that the formation of fibril-like structures is a cooperative processes that requires a critical number of monomers, M⋆≈25, in a cluster. The terminal tyrosine residue is the structural determinant in the formation of helical fibril-like structures. This work supports and quantifies the two-step aggregation model where the initially formed amorphous clusters grow and, when they are large enough, rearrange into mature twisted structures. However, in addition to the nucleated fibrillation, growing aggregates undergo further internal reorganization, which leads to more compact structures of large aggregates.

Molecular Dynamics Simulations on Interindividual Variability of Intestinal Fluids: Impact on Drug Solubilization

Efficient delivery of oral drugs is dependent on their solubility in human intestinal fluid, a complex and dynamic fluid that contains colloidal structures composed of small molecules. These structures solubilize poorly water-soluble compounds, increasing their apparent solubility, and possibly their bioavailability. In this study, we conducted coarse-grained molecular dynamics simulations with data from duodenal fluid samples previously acquired from five healthy volunteers. In these simulations, we observed the self-assembly of mixed micelles of bile salts, phospholipids, and free fatty acids. The micelles were ellipsoids with a size range of 4–7 nm. Next, we investigated micelle affinities of three model drugs. The affinities in our simulation showed the same trend as literature values for the solubility enhancement of drugs in human intestinal fluids. This type of simulations is useful for studies of events and interactions taking place in the small intestinal fluid.

Using Coarse-Grained Simulations to Characterize the Mechanisms of Protein-Protein Association

The formation of functionally versatile protein complexes underlies almost every biological process. The estimation of how fast these complexes can be formed has broad implications for unravelling the mechanism of biomolecular recognition. This kinetic property is traditionally quantified by association rates, which can be measured through various experimental techniques. To complement these time-consuming and labor-intensive approaches, we developed a coarse-grained simulation approach to study the physical processes of protein–protein association. We systematically calibrated our simulation method against a large-scale benchmark set. By combining a physics-based force field with a statistically-derived potential in the simulation, we found that the association rates of more than 80% of protein complexes can be correctly predicted within one order of magnitude relative to their experimental measurements. We further showed that a mixture of force fields derived from complementary sources was able to describe the process of protein–protein association with mechanistic details. For instance, we show that association of a protein complex contains multiple steps in which proteins continuously search their local binding orientations and form non-native-like intermediates through repeated dissociation and re-association. Moreover, with an ensemble of loosely bound encounter complexes observed around their native conformation, we suggest that the transition states of protein–protein association could be highly diverse on the structural level. Our study also supports the idea in which the association of a protein complex is driven by a “funnel-like” energy landscape. In summary, these results shed light on our understanding of how protein–protein recognition is kinetically modulated, and our coarse-grained simulation approach can serve as a useful addition to the existing experimental approaches that measure protein–protein association rates.

In Silico Study of the Mechanism of Binding of the N-Terminal Region of α Synuclein to Synaptic-Like Membranes

The membrane binding by α-synuclein (αS), a presynaptic protein whose aggregation is strongly linked with Parkinson’s disease, influences its biological behavior under functional and pathological conditions. This interaction requires a conformational transition from a disordered-unbound to a partially helical membrane-bound state of the protein. In the present study, we used enhanced coarse-grained MD simulations to characterize the sequence and conformational determinants of the binding to synaptic-like vesicles by the N-terminal region of αS. This region is the membrane anchor and is of crucial importance for the properties of the physiological monomeric state of αS as well as for its aberrant aggregates. These results identify the key factors that play a role in the binding of αS with synaptic lipid bilayers in both membrane-tethered and membrane-locked conformational states.

Over What Length Scale Does an Inorganic Substrate Perturb the Structure of a Glassy Organic Semiconductor?

While the bulk structure of vapor-deposited glasses has been extensively studied, structure at buried interfaces has received little attention, despite being important for organic electronic applications. To learn about glass structure at buried interfaces, we study the structure of vapor-deposited glasses of the organic semiconductor DSA-Ph (1,4-di-[4-(N,N-diphenyl)amino]styrylbenzene) as a function of film thickness; the structure is probed with grazing incidence X-ray scattering. We deposit on silicon and gold substrates and span a film thickness range of 10-600 nm. Our experiments demonstrate that interfacial molecular packing in vapor-deposited glasses of DSA-Ph is more disordered compared to the bulk. At a deposition temperature near room temperature, we estimate ∼8 nm near the substrate can have modified molecular packing. Molecular dynamics simulations of a coarse-grained representation of DSA-Ph reveal a similar length scale. In both the simulations and the experiments, deposition temperature controls glass structure beyond this interfacial layer of a few nanometers.

Tubulin tails and their modifications regulate protein diffusion on microtubules

Microtubules (MTs) are essential components of the eukaryotic cytoskeleton that serve as "highways" for intracellular trafficking. In addition to the well-known active transport of cargo by motor proteins, many MT-binding proteins seem to adopt diffusional motility as a transportation mechanism. However, because of the limited spatial resolution of current experimental techniques, the detailed mechanism of protein diffusion has not been elucidated. In particular, the precise role of tubulin tails and tail modifications in the diffusion process is unclear. Here, using coarse-grained molecular dynamics simulations validated against atomistic simulations, we explore the molecular mechanism of protein diffusion along MTs. We found that electrostatic interactions play a central role in protein diffusion; the disordered tubulin tails enhance affinity but slow down diffusion, and diffusion occurs in discrete steps. While diffusion along wild-type MT is performed in steps of dimeric tubulin, the removal of the tails results in a step of monomeric tubulin. We found that the energy barrier for diffusion is larger when diffusion on MTs is mediated primarily by the MT tails rather than the MT body. In addition, globular proteins (EB1 and PRC1) diffuse more slowly than an intrinsically disordered protein (Tau) on MTs. Finally, we found that polyglutamylation and polyglycylation of tubulin tails lead to slower protein diffusion along MTs, although polyglycylation leads to faster diffusion across MT protofilaments. Taken together, our results explain experimentally observed data and shed light on the roles played by disordered tubulin tails and tail modifications in the molecular mechanism of protein diffusion along MTs.

Kinetic Control of Amyloidogenesis Calls for Unconventional Drugs To Fight Alzheimer's Disease

Computer simulations had predicted that amyloid fibrillogenesis is governed by free energy barriers and kinetic traps (kinetic control), rather than the free energy of the final aggregates. The simulations suggested that the diversity in fibril morphologies can originate from variations in the number of protofilaments which has been confirmed by recent cryo-electron microscopy studies of amyloid-β fibrils derived from brain tissue of Alzheimer's patients. The kinetic control of fibril formation and polymorphism imply that chemical substances with new mechanisms of action are needed to fight Alzheimer's disease.

Structural insights into Acyl-coenzyme A binding domain containing 3 (ACBD3) protein hijacking by picornaviruses

Many picornaviruses hijack the Golgi resident Acyl-coenzyme A binding domain containing 3 (ACBD3) protein in order to recruit the phosphatidylinositol 4-kinase B (PI4KB) to viral replication organelles (ROs). PI4KB, once recruited and activated by ACBD3 protein, produces the lipid phosphatidylinositol 4-phosphate (PI4P), which is a key step in the biogenesis of viral ROs. To do so, picornaviruses use their small nonstructural protein 3A that binds the Golgi dynamics domain of the ACBD3 protein. Here, we present the analysis of the highly flexible ACBD3 proteins and the viral 3A protein in solution using small-angle X-ray scattering and computer simulations. Our analysis revealed that both the ACBD3 protein and the 3A:ACBD3 protein complex have an extended and flexible conformation in solution.

Manipulation of Confined Polyelectrolyte Conformations through Dielectric Mismatch

We demonstrate that a highly charged polyelectrolyte confined in a spherical cavity undergoes reversible transformations between amorphous conformations and a four-fold symmetry morphology as a function of dielectric mismatch between the media inside and outside the cavity. Surface polarization due to dielectric mismatch exhibits an extra "confinement" effect, which is most pronounced within a certain range of the cavity radius and the electrostatic strength between the monomers and counterions and multivalent counterions. For cavities with a charged surface, surface polarization leads to an increased amount of counterions adsorbed in the outer side, further compressing the confined polyelectrolyte into a four-fold symmetry morphology. The equilibrium conformation of the chain is dependent upon several key factors including the relative permittivities of the media inside and outside the cavity, multivalent counterion concentration, cavity radius relative to the chain length, and interface charge density. Our findings offer insights into the effects of dielectric mismatch in packaging and delivery of polyelectrolytes across media with different relative permittivities. Moreover, the reversible transformation of the polyelectrolyte conformations in response to environmental permittivity allows for potential applications in biosensing and medical monitoring.

Domain topology, stability, and translation speed determine mechanical force generation on the ribosome

The concomitant folding of a nascent protein domain with its synthesis can generate mechanical forces that act on the ribosome and alter translation speed. Such changes in speed can affect the structure and function of the newly synthesized protein as well as cellular phenotype. The domain properties that govern force generation have yet to be identified and understood, and the influence of translation speed is unknown because all reported measurements have been carried out on arrested ribosomes. Here, using coarse-grained molecular simulations and statistical mechanical modeling of protein synthesis, we demonstrate that force generation is determined by a domain's stability and topology, as well as translation speed. The statistical mechanical models we create predict how force profiles depend on these properties. These results indicate that force measurements on arrested ribosomes will not always accurately reflect what happens in a cell, especially for slow-folding domains, and suggest the possibility that certain domain properties may be enriched or depleted across the structural proteome of organisms through evolutionary selection pressures to modulate protein synthesis speed and posttranslational protein behavior.

Structural Model of the mIgM B-Cell Receptor Transmembrane Domain From Self-Association Molecular Dynamics Simulations

Antigen binding to B-cell antigen receptors (BCRs) followed by signaling initiates the humoral immune response. The signaling is intimately coupled to nanoclustering of BCRs and their sorting to specific membrane domains, a process that is ruled by interactions between the BCR transmembrane domain and lipids. While the structure of the extracellular domains of BCRs has been resolved, little is known about the configuration of the constituting four immunoglobulin domains spanning the membrane. Here, we modeled the structure of the transmembrane (TM) domain of the IgM B-cell receptor using self-assembly coarse-grained molecular dynamics simulations. The obtained quaternary structure was validated against available experimental data and atomistic simulations. The IgM-BCR-TM domain configuration shows a 1:1 stoichiometry between the homodimeric membrane-bound domain of IgM (mIgM) and a Ig-α/Ig-β heterodimer. The mIgM homodimer is based on an asymmetric association of two mIgM domains. We show that a specific site of the Ig-α/Ig-β heterodimer is responsible for the association of IgM-BCRs with lipid rafts. Our results further suggest that this site is blocked in small-sized IgM-BCR clusters. The BCR TM structure provides a molecular basis for the previously suggested dissociation activation model of B-cell receptors. Self-assembly molecular dynamics simulations at the coarse-grained scale here proved as a versatile tool in the study of receptor complexes.

Organizing membrane-curving proteins: the emerging dynamical picture

Lipid membranes play key roles in cells, such as in trafficking, division, infection, remodeling of organelles, among others. The key step in all these processes is creating membrane curvature, typically under the control of many anchored, adhered or included proteins. However, it has become clear that the membrane itself can mediate the interactions among proteins to produce highly ordered assemblies. Computer simulations are ideally suited to investigate protein organization and the dynamics of membrane remodeling at near-micron scales, something that is extremely challenging to tackle experimentally. We review recent computational efforts in modeling protein-caused membrane deformation mechanisms, specifically focusing on coarse-grained simulations. We highlight work that exposed the membrane-mediated ordering of proteins into lines, meshwork, spirals and other assemblies, in what seems to be a very generic mechanism driven by a combination of short and long-ranged forces. Modulating the mechanical properties of membranes is an underexplored signaling mechanism in various processes deserving of more attention in the near future.

Monovalent ions modulate the flux through multiple folding pathways of an RNA pseudoknot

The assembly mechanism of RNA, vital to describing its functions, depends on both the sequence and the metal ion concentration. How the latter influences the folding trajectories remains an important unsolved problem. Here, we examine the folding pathways of an RNA pseudoknot (PK) with key functional roles in transcription and translation, using a combination of experiments and simulations. We demonstrate that the PK, consisting of two hairpins with differing stabilities, folds by parallel pathways. Surprisingly, the flux between them is modulated by monovalent salt concentration. Our work shows that the order of assembly of PKs is determined by the relative stability of the hairpins, implying that the folding landscape can be controlled by sequence and ion concentration.

The functions of RNA pseudoknots (PKs), which are minimal tertiary structural motifs and an integral part of several ribozymes and ribonucleoprotein complexes, are determined by their structure, stability, and dynamics. Therefore, it is important to elucidate the general principles governing their thermodynamics/folding mechanisms. Here, we combine laser temperature-jump experiments and coarse-grained simulations to determine the folding/unfolding pathways of VPK, a variant of the mouse mammary tumor virus (MMTV) PK involved in ribosomal frameshifting. Fluorescent nucleotide analogs (2-aminopurine and pyrrolocytidine) placed at different stem/loop positions in the PK serve as local probes allowing us to monitor the order of assembly of VPK that has two constituent hairpins with different intrinsic stabilities. We show that at 50 mM KCl, the dominant folding pathway populates only the more stable hairpin intermediate; as the salt concentration is increased, a parallel folding pathway emerges involving the less stable hairpin as an alternate intermediate. Notably, the flux between the pathways is modulated by the ionic strength. Our findings support the principle that the order of PK structure formation is determined by the relative stabilities of the hairpins, which can be altered by sequence variations or salt concentrations. The experimental results of salt effects on the partitioning between the two folding pathways are in remarkable agreement with simulations that were performed with no adjustable parameters. Our study not only unambiguously demonstrates that VPK folds by parallel pathways but also showcases the power of combining experiments and simulations for a more enriched description of RNA self-assembly.

Viral nanomechanics with a virtual atomic force microscope

One of the most important components of a virus is the protein shell or capsid that encloses its genetic material. The main role of the capsid is to protect the viral genome against external aggressions, facilitating its safe and efficient encapsulation and delivery. As a consequence, viral capsids have developed astonishing mechanical properties that are crucial for viral function. These remarkable properties have started to be unveiled in single-virus nanoindentation experiments, and are opening the door to the use of viral-derived artificial nanocages for promising bio- and nano-technological applications. However, the interpretation of nanoindentation experiments is often difficult, requiring the support of theoretical and simulation analysis. Here we present a 'Virtual AFM' (VAFM), a Brownian Dynamics simulation of a coarse-grained model of virus aimed to mimic the standard setup of atomic force microscopy (AFM) nanoindentation experiments. Despite the heavy level of coarse-graining, these simulations provide valuable information which is not accessible in experiments. Rather than focusing on a specific virus, the VAFM will be used to analyze how the mechanical response and breaking of viruses depend on different parameters controlling the effective interactions between capsid's structural units. In particular, we will discuss the influence of adsorption, the tip radius, and the rigidity and shape of the shell on its mechanical response.

Chemical Kinetics for Bridging Molecular Mechanisms and Macroscopic Measurements of Amyloid Fibril Formation

Understanding how normally soluble peptides and proteins aggregate to form amyloid fibrils is central to many areas of modern biomolecular science, ranging from the development of functional biomaterials to the design of rational therapeutic strategies against increasingly prevalent medical conditions such as Alzheimer's and Parkinson's diseases. As such, there is a great need to develop models to mechanistically describe how amyloid fibrils are formed from precursor peptides and proteins. Here we review and discuss how ideas and concepts from chemical reaction kinetics can help to achieve this objective. In particular, we show how a combination of theory, experiments, and computer simulations, based on chemical kinetics, provides a general formalism for uncovering, at the molecular level, the mechanistic steps that underlie the phenomenon of amyloid fibril formation.

Optical Voltage Sensing Using DNA Origami

We explore the potential of DNA nanotechnology for developing novel optical voltage sensing nanodevices that convert a local change of electric potential into optical signals. As a proof-of-concept of the sensing mechanism, we assembled voltage responsive DNA origami structures labeled with a single pair of FRET dyes. The DNA structures were reversibly immobilized on a nanocapillary tip and underwent controlled structural changes upon application of an electric field. The applied field was monitored through a change in FRET efficiency. By exchanging the position of a single dye, we could tune the voltage sensitivity of our DNA origami structure, demonstrating the flexibility and versatility of our approach. The experimental studies were complemented by coarse-grained simulations that characterized voltage-dependent elastic deformation of the DNA nanostructures and the associated change in the distance between the FRET pair. Our work opens a novel pathway for determining the mechanical properties of DNA origami structures and highlights potential applications of dynamic DNA nanostructures as voltage sensors.

Modeling the effects of lattice defects on microtubule breaking and healing

Microtubule reorganization often results from the loss of polymer induced through breakage or active destruction by energy-using enzymes. Pre-existing defects in the microtubule lattice likely lower structural integrity and aid filament destruction. Using large-scale molecular simulations, we model diverse microtubule fragments under forces generated at specific positions to locally crush the filament. We show that lattices with 2% defects are crushed and severed by forces three times smaller than defect-free ones. We validate our results with direct comparisons of microtubule kinking angles during severing. We find a high statistical correlation between the angle distributions from experiments and simulations indicating that they sample the same population of structures. Our simulations also indicate that the mechanical environment of the filament affects breaking: local mechanical support inhibits healing after severing, especially in the case of filaments with defects. These results recall reports of microtubule healing after flow-induced bending and corroborate prior experimental studies that show severing is more likely at locations where microtubules crossover in networks. Our results shed new light on mechanisms underlying the ability of microtubules to be destroyed and healed in the cell, either by external forces or by severing enzymes wedging dimers apart. © 2016 Wiley Periodicals, Inc.

Kinematics of the lever arm swing in myosin VI

Myosin VI (MVI) is the only known member of the myosin superfamily that, upon dimerization, walks processively toward the pointed end of the actin filament. The leading head of the dimer directs the trailing head forward with a power stroke, a conformational change of the motor domain exaggerated by the lever arm. Using a unique coarse-grained model for the power stroke of a single MVI, we provide the molecular basis for its motility. We show that the power stroke occurs in two major steps. First, the motor domain attains the poststroke conformation without directing the lever arm forward; and second, the lever arm reaches the poststroke orientation by undergoing a rotational diffusion. From the analysis of the trajectories, we discover that the potential that directs the rotating lever arm toward the poststroke conformation is almost flat, implying that the lever arm rotation is mostly uncoupled from the motor domain. Because a backward load comparable to the largest interhead tension in a MVI dimer prevents the rotation of the lever arm, our model suggests that the leading-head lever arm of a MVI dimer is uncoupled, in accord with the inference drawn from polarized total internal reflection fluorescence (polTIRF) experiments. Without any adjustable parameter, our simulations lead to quantitative agreement with polTIRF experiments, which validates the structural insights. Finally, in addition to making testable predictions, we also discuss the implications of our model in explaining the broad step-size distribution of the MVI stepping pattern.

Carbon Nanotubes Mediate Fusion of Lipid Vesicles

The fusion of lipid membranes is opposed by high energetic barriers. In living organisms, complex protein machineries carry out this biologically essential process. Here we show that membrane-spanning carbon nanotubes (CNTs) can trigger spontaneous fusion of small lipid vesicles. In coarse-grained molecular dynamics simulations, we find that a CNT bridging between two vesicles locally perturbs their lipid structure. Their outer leaflets merge as the CNT pulls lipids out of the membranes, creating an hourglass-shaped fusion intermediate with still intact inner leaflets. As the CNT moves away from the symmetry axis connecting the vesicle centers, the inner leaflets merge, forming a pore that completes fusion. The distinct mechanism of CNT-mediated membrane fusion may be transferable, providing guidance in the development of fusion agents, e.g., for the targeted delivery of drugs or nucleic acids.

Assembly of Influenza Hemagglutinin Fusion Peptides in a Phospholipid Bilayer by Coarse-grained Computer Simulations

Membrane fusion is critical to eukaryotic cellular function and crucial to the entry of enveloped viruses such as influenza and human immunodeficiency virus. Influenza viral entry in the host cell is mediated by a 20–23 amino acid long sequence, called the fusion peptide (FP). Recently, possible structures for the fusion peptide (ranging from an inverted V shaped α-helical structure to an α-helical hairpin, or to a complete α-helix) and their implication in the membrane fusion initiation have been proposed. Despite the large number of studies devoted to the structure of the FP, the mechanism of action of this peptide remains unclear with several mechanisms having been suggested, including the induction of local disorder, promoting membrane curvature, and/or altering local membrane composition. In recent years, several research groups have employed atomistic and/or coarse-grained molecular dynamics (MD) simulations to investigate the matter. In all previous works, the behavior of a single FP monomer was studied, while in this manuscript, we use a simplified model of a tripeptide (TP) monomer of FP (TFP) instead of a single FP monomer because each Influenza Hemagglutinin contains three FP molecules in the biological system. In this manuscript we report findings targeted at understanding the fusogenic properties and the collective behavior of these trimers of FP peptides on a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine model membrane. Here we show how the TFP monomers self-assemble into differently sized oligomers in the presence of the membrane. We measure the perturbation to the structure of the phospholipid membrane caused by the presence of these TFP oligomers. Our work (i) shows how self-assembly of TFP in the presence of the membrane induces non negligible deformation to the membrane and (ii) could be a useful starting point to stimulate discussion and further work targeted to fusion pore formation.

Designed protein reveals structural determinants of extreme kinetic stability

The design of stable, functional proteins is difficult. Improved design requires a deeper knowledge of the molecular basis for design outcomes and properties. We previously used a bioinformatics and energy function method to design a symmetric superfold protein composed of repeating structural elements with multivalent carbohydrate-binding function, called ThreeFoil. This and similar methods have produced a notably high yield of stable proteins. Using a battery of experimental and computational analyses we show that despite its small size and lack of disulfide bonds, ThreeFoil has remarkably high kinetic stability and its folding is specifically chaperoned by carbohydrate binding. It is also extremely stable against thermal and chemical denaturation and proteolytic degradation. We demonstrate that the kinetic stability can be predicted and modeled using absolute contact order (ACO) and long-range order (LRO), as well as coarse-grained simulations; the stability arises from a topology that includes many long-range contacts which create a large and highly cooperative energy barrier for unfolding and folding. Extensive data from proteomic screens and other experiments reveal that a high ACO/LRO is a general feature of proteins with strong resistances to denaturation and degradation. These results provide tractable approaches for predicting resistance and designing proteins with sufficient topological complexity and long-range interactions to accommodate destabilizing functional features as well as withstand chemical and proteolytic challenge.

Multiple stiffening effects of nanoscale knobs on human red blood cells infected with Plasmodium falciparum malaria parasite

Our coarse-grained molecular dynamics (CGMD) simulations show that the deposition of nanoscale knobs, rather than spectrin network remodeling, is the primary cause of the dramatically increased stiffness of the Plasmodium falciparum (Pf)-infected red blood cell (RBC) membranes. Our analyses further reveal that the knobs stiffen the RBC membrane in a unique manner by simultaneously harnessing composite strengthening, strain hardening, and knob density-dependent vertical coupling effects. In addition to providing a fundamental understanding of the stiffening mechanism of Pf-infected RBCs, our simulation results suggest potential targets for antimalarial therapies.

During its asexual development within the red blood cell (RBC), Plasmodium falciparum (Pf), the most virulent human malaria parasite, exports proteins that modify the host RBC membrane. The attendant increase in cell stiffness and cytoadherence leads to sequestration of infected RBCs in microvasculature, which enables the parasite to evade the spleen, and leads to organ dysfunction in severe cases of malaria. Despite progress in understanding malaria pathogenesis, the molecular mechanisms responsible for the dramatic loss of deformability of Pf-infected RBCs have remained elusive. By recourse to a coarse-grained (CG) model that captures the molecular structures of Pf-infected RBC membrane, here we show that nanoscale surface protrusions, known as “knobs,” introduce multiple stiffening mechanisms through composite strengthening, strain hardening, and knob density-dependent vertical coupling. On one hand, the knobs act as structural strengtheners for the spectrin network; on the other, the presence of knobs results in strain inhomogeneity in the spectrin network with elevated shear strain in the knob-free regions, which, given its strain-hardening property, effectively stiffens the network. From the trophozoite to the schizont stage that ensues within 24–48 h of parasite invasion into the RBC, the rise in the knob density results in the increased number of vertical constraints between the spectrin network and the lipid bilayer, which further stiffens the membrane. The shear moduli of Pf-infected RBCs predicted by the CG model at different stages of parasite maturation are in agreement with experimental results. In addition to providing a fundamental understanding of the stiffening mechanisms of Pf-infected RBCs, our simulation results suggest potential targets for antimalarial therapies.

Crucial role of nonspecific interactions in amyloid nucleation

Protein oligomers have been implicated as toxic agents in a wide range of amyloid-related diseases. However, it has remained unsolved whether the oligomers are a necessary step in the formation of amyloid fibrils or just a dangerous byproduct. Analogously, it has not been resolved if the amyloid nucleation process is a classical one-step nucleation process or a two-step process involving prenucleation clusters. We use coarse-grained computer simulations to study the effect of nonspecific attractions between peptides on the primary nucleation process underlying amyloid fibrillization. We find that, for peptides that do not attract, the classical one-step nucleation mechanism is possible but only at nonphysiologically high peptide concentrations. At low peptide concentrations, which mimic the physiologically relevant regime, attractive interpeptide interactions are essential for fibril formation. Nucleation then inevitably takes place through a two-step mechanism involving prefibrillar oligomers. We show that oligomers not only help peptides meet each other but also, create an environment that facilitates the conversion of monomers into the β-sheet-rich form characteristic of fibrils. Nucleation typically does not proceed through the most prevalent oligomers but through an oligomer size that is only observed in rare fluctuations, which is why such aggregates might be hard to capture experimentally. Finally, we find that the nucleation of amyloid fibrils cannot be described by classical nucleation theory: in the two-step mechanism, the critical nucleus size increases with increases in both concentration and interpeptide interactions, which is in direct contrast with predictions from classical nucleation theory.

Local Lipid Reorganization by a Transmembrane Protein Domain

Membrane proteins interact with their lipid bilayer environment via both a transmembrane helix and juxtamembrane regions. The effect of juxtamembrane regions and membrane lipid composition on these interactions has been explored by multiscale molecular dynamics simulations. The consequences of anionic lipids within the inner leaflet of a membrane were studied in combination with membrane spanning protein models differing in their juxtamembrane domains. The simulations reveal sensitivity of the protein-lipid interactions to membrane lipid composition and charged amino acid side chains. Basic residues on the intracellular side of the protein facilitated interactions with anionic lipids. Protein systems without basic residues do not show selectivity for anionic compared with zwitterionic lipids. This reveals the sensitivity to the composition of both the membrane and the protein system when studying membrane-embedded proteins. The results presented here illustrate how even a simple transmembrane domain is able to induce lipid reorganization in a mixed asymmetric bilayer.

Coarse-grained molecular dynamics simulation of DNA translocation in chemically modified nanopores

Solid-state nanopores provide a direct means to detect and analyze DNA and proteins. In a typical setup, the DNA molecules travel through a nanopore under electrophoretic voltage bias. The nanopore is sandwiched between two chambers that are filled with ionic solution. A major challenge in using solid-state nanopores for DNA sequencing and gene detection is to improve their selectivity and detection sensitivity. To achieve these goals, one solution is to functionalize the nanopores by chemically modifying the pore walls with silanes or nucleic acids. However, little is known about molecular interactions in functionalized nanopores. This paper presents DNA translocation dynamics and the mechanism of DNA sequencing in a functionalized nanopore through a coarse-grained molecular dynamics model. The DNA nucleotide is coarse-grained into two interaction sites: one site corresponds to the base group and the other encompasses the phosphate and sugar groups. The water molecules are included in the model implicitly through Langevin dynamics. The coarse-grained model immensely improves the computational efficiency while still capturing the essential translocation dynamics. The model characterizes important physical properties of functionalized nanopores such as the effective pore diameter and effect of biasing voltage on the DNA translocation dynamics. The model reveals a nonlinear relationship between translocation speed of DNA and applied voltage. Moreover, DNA translocation in nanopores functionalized with hairpin-loop (HPL) DNA and single-stranded DNA (ss-DNA) shows significant differences: a target DNA is found to translocate through a ss-DNA coated nanopore 9 times faster than through an HPL coated one at a bias of 100 mV, putatively from lower stiffness of ss-DNA than that for HPL. The DNA translocation speed is also largely influenced by interaction potential between the DNA and surface-tethered molecules. The results reveal that such selective translocation, distinctly different translocation dynamics of target DNA molecules largely stem from the flexibility and orientation of the surface-tethered molecules. These findings can significantly impact the rational design of DNA transport experiments leading to rapid molecule-level diagnostics.

Probing the origin of tubulin rigidity with molecular simulations

Tubulin heterodimers are the building blocks of microtubules, a major component of the cytoskeleton, whose mechanical properties are fundamental for the life of the cell. We uncover the microscopic origins of the mechanical response in microtubules by probing features of the energy landscape of the tubulin monomers and tubulin heterodimer. To elucidate the structures of the unfolding pathways and reveal the multiple unfolding routes, we performed simulations of a self-organized polymer (SOP) model of tubulin. The SOP representation, which is a coarse-grained description of chains, allows us to perform force-induced simulations at loading rates and time scales that closely match those used in single-molecule experiments. We show that the forced unfolding of each monomer involves a bifurcation in the pathways to the stretched state. After the unfolding of the C-term domain, the unraveling continues either from the N-term domain or from the middle domain, depending on the monomer and the pathway. In contrast to the unfolding complexity of the monomers, the dimer unfolds according to only one route corresponding to the unraveling of the C-term domain and part of the middle domain of beta-tubulin. We find that this surprising behavior is due to the viscoelastic properties of the interface between the monomers. We map precise features of the complex energy landscape of tubulin by surveying the structures of the various metastable intermediates, which, in the dimer case, are characterized only by changes in the beta-tubulin monomer.