Strain-Promoted Cycloadditions in Lipid Bilayers Triggered by Liposome Fusion.

Due to the variety of roles served by the cell membrane, its composition and structure are complex, making it difficult to study. Bioorthogonal reactions, such as the strain promoted azide-alkyne cycloaddition (SPAAC), are powerful tools for exploring the function of biomolecules in their native environment but have been largely unexplored within the context of lipid bilayers. Here, we developed a new approach to study the SPAAC reaction in liposomal membranes using azide- and strained alkyne-functionalized Förster resonance energy transfer (FRET) dye pairs. This study represents the first characterization of the SPAAC reaction between diffusing molecules inside liposomal membranes. Potential applications of this work include in situ bioorthogonal labeling of membrane proteins, improved understanding of membrane dynamics and fluidity, and the generation of new probes for biosensing assays.

Engineering Flexible Metal-Phenolic Networks with Guest Responsiveness via Intermolecular Interactions.

Flexible metal-organic materials are of growing interest owing to their ability to undergo reversible structural transformations under external stimuli. Here, we report flexible metal-phenolic networks (MPNs) featuring stimuli-responsive behavior to diverse solute guests. The competitive coordination of metal ions to phenolic ligands of multiple coordination sites and solute guests (e.g., glucose) primarily determines the responsive behavior of the MPNs, as revealed experimentally and computationally. Glucose molecules can be embedded into the dynamic MPNs upon mixing, leading to the reconfiguration of the metal-organic networks and thus changes in their physicochemical properties for targeting applications. This study expands the library of stimuli-responsive flexible metal-organic materials and the understanding of intermolecular interactions between metal-organic materials and solute guests, which is essential for the rational design of responsive materials for various applications.

Block Length-Dependent Protein Fouling on Poly(2-oxazoline)-Based Polymersomes: Influence on Macrophage Association and Circulation Behavior.

Polymersomes are vesicular structures self-assembled from amphiphilic block copolymers and are considered an alternative to liposomes for applications in drug delivery, immunotherapy, biosensing, and as nanoreactors and artificial organelles. However, the limited availability of systematic stability, protein fouling (protein corona formation), and blood circulation studies hampers their clinical translation. Poly(2-oxazoline)s (POx) are valuable antifouling hydrophilic polymers that can replace the current gold-standard, poly(ethylene glycol) (PEG), yet investigations of POx functionality on nanoparticles are relatively sparse. Herein, a systematic study is reported of the structural, dynamic and antifouling properties of polymersomes made of poly(2-methyl-2-oxazoline)-block-poly(dimethylsiloxane)-block-poly(2-methyl-2-oxazoline) (PMOXA-b-PDMS-b-PMOXA). The study relates in vitro antifouling performance of the polymersomes to atomistic molecular dynamics simulations of polymersome membrane hydration behavior. These observations support the experimentally demonstrated benefit of maximizing the length of PMOXA (degree of polymerization (DP) > 6) while keeping PDMS at a minimal length that still provides sufficient membrane stability (DP > 19). In vitro macrophage association and in vivo blood circulation evaluation of polymersomes in zebrafish embryos corroborate these findings. They further suggest that single copolymer presentation on polymersomes is outperformed by blends of varied copolymer lengths. This study helps to rationalize design rules for stable and low-fouling polymersomes for future medical applications.

Membrane Permeating Macrocycles: Design Guidelines from Machine Learning.

The ability to predict cell-permeable candidate molecules has great potential to assist drug discovery projects. Large molecules that lie beyond the Rule of Five (bRo5) are increasingly important as drug candidates and tool molecules for chemical biology. However, such large molecules usually do not cross cell membranes and cannot access intracellular targets or be developed as orally bioavailable drugs. Here, we describe a random forest (RF) machine learning model for the prediction of passive membrane permeation rates developed using a set of over 1000 bRo5 macrocyclic compounds. The model is based on easily calculated chemical features/descriptors as independent variables. Our random forest (RF) model substantially outperforms a multiple linear regression model based on the same features and achieves better performance metrics than previously reported models using the same underlying data. These features include: (1) polar surface area in water, (2) the octanol-water partitioning coefficient, (3) the number of hydrogen-bond donors, (4) the sum of the topological distances between nitrogen atoms, (5) the sum of the topological distances between nitrogen and oxygen atoms, and (6) the multiple molecular path count of order 2. The last three features represent molecular flexibility, the ability of the molecule to adopt different conformations in the aqueous and membrane interior phases, and the molecular "chameleonicity." Guided by the model, we propose design guidelines for membrane-permeating macrocycles. It is anticipated that this model will be useful in guiding the design of large, bioactive molecules for medicinal chemistry and chemical biology applications.

Electromagnetic bioeffects: a multiscale molecular simulation perspective.

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

Site-Selective Coordination Assembly of Dynamic Metal-Phenolic Networks.

Coordination states of metal-organic materials are known to dictate their physicochemical properties and applications in various fields. However, understanding and controlling coordination sites in metal-organic systems is challenging. Herein, we report the synthesis of site-selective coordinated metal-phenolic networks (MPNs) using flavonoids as coordination modulators. The site-selective coordination was systematically investigated experimentally and computationally using ligands with one, two, and multiple different coordination sites. Tuning the multimodal Fe coordination with catechol, carbonyl, and hydroxyl groups within the MPNs enabled the facile engineering of diverse physicochemical properties including size, selective permeability (20-2000 kDa), and pH-dependent degradability. This study expands our understanding of metal-phenolic chemistry and provides new routes for the rational design of structurally tailorable coordination-based materials.

Fluorinated Metal-Organic Coatings with Selective Wettability.

Surface chemistry is a major factor that determines the wettability of materials, and devising broadly applicable coating strategies that afford tunable and selective surface properties required for next-generation materials remains a challenge. Herein, we report fluorinated metal-organic coatings that display water-wetting and oil-repelling characteristics, a wetting phenomenon different from responsive wetting induced by external stimuli. We demonstrate this selective wettability with a library of metal-organic coatings using catechol-based coordination and silanization (both fluorinated and fluorine-free), enabling sensing through interfacial reconfigurations in both gaseous and liquid environments, and establish a correlation between the coating wettability and polarity of the liquids. This selective wetting performance is substrate-independent, spontaneous, durable, and reversible and occurs over a range of polar and nonpolar liquids (60 studied). These results provide insight into advanced liquid-solid interactions and a pathway toward tuning interfacial affinities and realizing robust, selective superwettability according to the surrounding conditions.

Design of Lipid-Based Nanocarriers via Cation Modulation of Ethanol-Interdigitated Lipid Membranes.

Short-chain alcohols (i.e., ethanol) can induce membrane interdigitation in saturated-chain phosphatidylcholines (PCs). In this process, alcohol molecules intercalate between phosphate heads, increasing lateral separation and favoring hydrophobic interactions between opposing acyl chains, which interpenetrate forming an interdigitated phase. Unraveling mechanisms underlying the interactions between ethanol and model lipid membranes has implications for cell biology, biochemistry, and for the formulation of lipid-based nanocarriers. However, investigations of ethanol-lipid membrane systems have been carried out in deionized water, which limits their applicability. Here, using a combination of small- and wide-angle X-ray scattering, small-angle neutron scattering, and all-atom molecular dynamics simulations, we analyzed the effect of varying CaCl2 and NaCl concentrations on ethanol-induced interdigitation. We observed that while ethanol addition leads to the interdigitation of bulk phase 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers in the presence of CaCl2 and NaCl regardless of the salt concentration, the ethanol-induced interdigitation of vesicular DPPC depends on the choice of cation and its concentration. These findings unravel a key role for cations in the ethanol-induced interdigitation of lipid membranes in either bulk phase or vesicular form.

Cobalt-Directed Assembly of Antibodies onto Metal-Phenolic Networks for Enhanced Particle Targeting.

The orientation-specific immobilization of antibodies onto nanoparticles, to preserve antibody-antigen recognition, is a key challenge in developing targeted nanomedicines. Herein, we report the targeting ability of metal-phenolic network (MPN)-coated gold nanoparticles with surface-physisorbed antibodies against respective antigens. The MPN coatings were self-assembled from metal ions (FeIII, CoII, CuII, NiII, or ZnII) cross-linked with tannic acid. Upon physisorption of antibodies, all particle systems exhibited enhanced association with target antigens, with CoII systems demonstrating more than 2-fold greater association. These systems contained more metal atoms distributed in a way to specifically interact with antibodies, which were investigated by molecular dynamics simulations. A model antibody fragment crystallizable (Fc) region in solution with CoII-tannic acid complexes revealed that the solvent-exposed CoII can directly coordinate to the histidine-rich portion of the Fc region. This one-pot interaction suggests anchoring of the antibody Fc region to the MPN on nanoparticles, allowing for enhanced targeting.

Robust and Versatile Coatings Engineered via Simultaneous Covalent and Noncovalent Interactions.

Interfacial modular assembly has emerged as an adaptable strategy for engineering the surface properties of substrates in biomedicine, photonics, and catalysis. Herein, we report a versatile and robust coating (pBDT-TA), self-assembled from tannic acid (TA) and a self-polymerizing aromatic dithiol (i.e., benzene-1,4-dithiol, BDT), that can be engineered on diverse substrates with a precisely tuned thickness (5-40 nm) by varying the concentration of BDT used. The pBDT-TA coating is stabilized by covalent (disulfide) bonds and supramolecular (π-π) interactions, endowing the coating with high stability in various harsh aqueous environments across ionic strength, pH, temperature (e.g., 100 mM NaCl, HCl (pH 1) or NaOH (pH 13), and water at 100 °C), as well as surfactant solution (e.g., 100 mM Triton X-100) and biological buffer (e.g., Dulbecco's phosphate-buffered saline), as validated by experiments and simulations. Moreover, the reported pBDT-TA coating enables secondary reactions on the coating for engineering hybrid adlayers (e.g., ZIF-8 shells) via phenolic-mediated adhesion, and the facile integration of aromatic fluorescent dyes (e.g., rhodamine B) via π interactions without requiring elaborate synthetic processes.

Surface Dynamics and Ligand-Core Interactions of Quantum Sized Photoluminescent Gold Nanoclusters.

Quantum-sized metallic clusters protected by biological ligands represent a new class of luminescent materials; yet the understanding of structural information and photoluminescence origin of these ultrasmall clusters remains a challenge. Herein we systematically study the surface ligand dynamics and ligand-metal core interactions of peptide-protected gold nanoclusters (AuNCs) with combined experimental characterizations and theoretical molecular simulations. We show that the peptide sequence plays an important role in determining the surface peptide structuring, interfacial water dynamics and ligand-Au core interaction, which can be tailored by controlling peptide acetylation, constituent amino acid electron donating/withdrawing capacity, aromaticity/hydrophobicity and by adjusting environmental pH. Specifically, emission enhancement is achieved through increasing the electron density of surface ligands in proximity to the Au core, discouraging photoinduced quenching, and by reducing the amount of surface-bound water molecules. These findings provide key design principles for understanding the surface dynamics of peptide-protected nanoparticles and maximizing the photoluminescence of metallic clusters through the exploitation of biologically relevant ligand properties.

Intra- and Intersubunit Ion-Pair Interactions Determine the Ability of Apolipoprotein C-II Mutants To Form Hybrid Amyloid Fibrils.

The apolipoprotein family is structurally defined by amphipathic α-helical regions that interact with lipid surfaces. In the absence of lipid, human apolipoprotein (apo) C-II also forms well-defined amyloid fibrils with cross-ß structure. Formation of this ß-structure is accompanied by the burial of two charged residues, K30 and D69, that form an ion-pair within the amyloid fibril core. Molecular dynamics (MD) simulations indicate these buried residues form both intra- and intersubunit ion-pair interactions that stabilize the fibril. Mutations of the ion-pair (either K30D or D69K) reduce fibril stability and prevent fibril formation by K30D apoC-II under standard conditions. We investigated whether mixing K30D apoC-II with other mutants would overcome this loss of fibril forming ability. Co-incubation of equimolar mixtures of K30D apoC-II with wild-type, D69K, or double-mutant (K30D/D69K) apoC-II promoted the incorporation of K30D apoC-II into hybrid fibrils with increased stability. MD simulations showed an increase in the number of intersubunit ion-pair interactions accompanied the increased stability of the hybrid fibrils. These results demonstrate the important role of both intra- and intersubunit charge interactions in stabilizing apoC-II amyloid fibrils, a process that may be a key factor in determining the general ability of proteins to form amyloid fibrils.

Influence of Ionic Strength on the Deposition of Metal-Phenolic Networks.

Metal-phenolic networks (MPNs) are a versatile class of self-assembled materials that are able to form functional thin films on various substrates with potential applications in areas including drug delivery and catalysis. Different metal ions (e.g., FeIII, CuII) and phenols (e.g., tannic acid, gallic acid) have been investigated for MPN film assembly; however, a mechanistic understanding of the thermodynamics governing MPN formation remains largely unexplored. To date, MPNs have been deposited at low ionic strengths (<5 mM), resulting in films with typical thicknesses of ∼10 nm, and it is still unclear how a bulk complexation reaction results in homogeneous thin films when a substrate is present. Herein we explore the influence of ionic strength (0-2 M NaCl) on the conformation of MPN precursors in solution and how this determines the final thickness and morphology of MPN films. Specifically, the film thickness increases from 10 nm in 0 M NaCl to 12 nm in 0.5 M NaCl and 15 nm in 1 M NaCl, after which the films grow rougher rather than thicker. For example, the root-mean-square roughness values of the films are constant below 1 M NaCl at 1.5 nm; in contrast, the roughness is 3 nm at 1 M NaCl and increases to 5 nm at 2 M NaCl. Small-angle X-ray scattering and molecular dynamics simulations allow for comparisons to be made with chelated metals and polyelectrolyte thin films. For example, at a higher ionic strength (2 M NaCl), sodium ions shield the galloyl groups of tannic acid, allowing them to extend away from the FeIII center and interact with other MPN complexes in solution to form thicker and rougher films. As the properties of films determine their final performance and application, the ability to tune both thickness and roughness using salts may allow for new applications of MPNs.

Long-range dipolar order and dispersion forces in polar liquids.

Complex solvation phenomena, such as specific ion effects, occur in polar liquids. Interpretation of these effects in terms of structure and dispersion forces will lead to a greater understanding of solvation. Herein, using molecular dynamics, we probe the structure of polar liquids through specific dipolar pair correlation functions that contribute to the potential of mean force that is "felt" between thermally rotating dipole moments. It is shown that unique dipolar order exists at separations at least up to 20 Å for all liquids studied. When the structural order is compared with a dipolar dispersion force that arises from local co-operative enhancement of dipole moments, a strong agreement is found. Lifshitz theory of dispersion forces was compared with the structural order, where the theory is validated for all liquids that do not have significant local dipole correlations. For liquids that do have significant local dipole correlations, specifically liquid water, Lifshitz theory underestimates the dispersion force by a factor of 5-10, demonstrating that the force that leads to the increased structure in liquid water is missed by Lifshitz theory of van der Waals forces. We apply similar correlation functions to an ionic aqueous system, where long-range order between water's dipole moment and a single chloride ion is found to exist at 20 Å of separation, revealing a long-range perturbation of water's structure by an ion. Furthermore, we found that waters within the 1st, 2nd, and 3rd solvation shells of a chloride ion exhibit significantly enhanced dipolar interactions, particularly with waters at larger distances of separation. Our results provide a link between structures, dispersion forces, and specific ion effects, which may lead to a more robust understanding of solvation.

Solution Conditions Affect the Ability of the K30D Mutation To Prevent Amyloid Fibril Formation by Apolipoprotein C-II: Insights from Experiments and Theoretical Simulations.

Apolipoproteins form amphipathic helical structures that bind lipid surfaces. Paradoxically, lipid-free apolipoproteins display a strong propensity to form cross-ß structure and self-associate into disease-related amyloid fibrils. Studies of apolipoprotein C-II (apoC-II) amyloid fibrils suggest that a K30-D69 ion pair accounts for the dual abilities to form helix and cross-ß structure. Consistent with this is the observation that a K30D mutation prevents fibril formation under standard fibril forming conditions. However, we found that fibril formation by K30D apoC-II proceeded readily at low pH and a higher salt or protein concentration. Structural analysis demonstrated that K30D apoC-II fibrils at pH 7 have a structure similar to that of the wild-type fibrils but are less stable. Molecular dynamics simulations of the wild-type apoC-II fibril model at pH 7 and 3 showed that the loss of charge on D69 at pH 3 leads to greater separation between residues K30 and D69 within the fibril with a corresponding reduction in ß-strand content around residue 30. In contrast, in simulations of the K30D mutant model at pH 7 and 3, residues D30 and D69 moved closer at pH 3, accompanied by an increase in ß-strand content around residue 30. The simulations also demonstrated a strong dominance of inter- over intramolecular contacts between ionic residues of apoC-II and suggested a cooperative mechanism for forming favorable interactions between the individual strands under different conditions. These observations demonstrate the important role of the buried K30-D69 ion pair in the stability and solution properties of apoC-II amyloid fibrils.

Understanding and Designing the Gold-Bio Interface: Insights from Simulations.

Gold nanoparticles (AuNPs) are an integral part of many exciting and novel biomedical applications, sparking the urgent need for a thorough understanding of the physicochemical interactions occurring between these inorganic materials, their functional layers, and the biological species they interact with. Computational approaches are instrumental in providing the necessary molecular insight into the structural and dynamic behavior of the Au-bio interface with spatial and temporal resolutions not yet achievable in the laboratory, and are able to facilitate a rational approach to AuNP design for specific applications. A perspective of the current successes and challenges associated with the multiscale computational treatment of Au-bio interfacial systems, from electronic structure calculations to force field methods, is provided to illustrate the links between different approaches and their relationship to experiment and applications.

Gold Nanoparticles: Understanding and Designing the Gold-Bio Interface: Insights from Simulations (Small 18/2016).

Computational modeling is a virtual microscope providing molecular insight into the complex interactions occurring at nano-bio interfaces. The image shows the atomistic landscape of a peptide-adlayer on gold nanoparticles, providing the physicochemical information needed to tailor the design of efficient colorimetric biosensors. The successes and challenges of theoretical simulations in facilitating gold nanomaterial design for biomedicine are reviewed on page 2395 by I. Yarovsky and co-workers.

Surface heterogeneity: a friend or foe of protein adsorption - insights from theoretical simulations.

A lack in the detailed understanding of mechanisms through which proteins adsorb or are repelled at various solid/liquid interfaces limits the capacity to rationally design and produce more sophisticated surfaces with controlled protein adsorption in both biomedical and industrial settings. To date there are three main approaches to achieve anti biofouling efficacy, namely chemically adjusting the surface hydrophobicity and introducing various degrees of surface roughness, or a combination of both. More recently, surface nanostructuring has been shown to have an effect on protein adsorption. However, the current resolution of experimental techniques makes it difficult to investigate these three phase systems at the molecular level. In this molecular dynamics study we explore in all-atom detail the adsorption process of one of the most surface active proteins, EAS hydrophobin, known for its versatile ability to self-assemble on both hydrophobic and hydrophilic surfaces forming stable monolayers that facilitate further biofilm growth. We model the adsorption of this protein on organic ligand protected silica surfaces with varying degrees of chemical heterogeneity and roughness, including fully homogenous hydrophobic and hydrophilic surfaces for comparison. We present a detailed characterisation of the functionalised surface structure and dynamics for each of these systems, and the effect the ligands have on interfacial water, the adsorption process and conformational rearrangements of the protein. Results suggest that the ligand arrangement that produces the highest hydrophilic chain mobility and the lack of significant hydrophobic patches shows the most promising anti-fouling efficacy toward hydrophobin. However, the presence on the protein surface of a flexible loop with amphipathic character (the Cys3-Cys4 loop) is seen to facilitate EAS adsorption on all surfaces by enabling the protein to match the surface pattern.

Electromagnetic-field effects on structure and dynamics of amyloidogenic peptides.

Electromagnetic fields (EMFs) are ever-present, and so is the need to better understand their influence on human health and biological matter in general. The interaction between a molecular system and external EMF can alter the structure, and dynamical behaviour, and, hence, biological function of proteins with uncertain health consequences. This urges a detailed investigation of EMF-induced effects on basic protein biophysics. Here, we used all-atom non-equilibrium molecular dynamics simulations to understand and quantify the response mechanisms of the amyloidogenic apoC-II(60-70) peptides to non-ionising radiation by modelling their behaviour under external electromagnetic and electric fields of different strengths. Our simulations show high strength fields (>0.04 V/nm) cause structural changes in apoC-II(60-70) due to the peptide dipole alignment along the applied field direction, which disrupts the inherent ß-hairpin conformation known to be the intermediate state for fibril formation. The intermediate field-strength range (0.04-0.004 V/nm) causes a significant acceleration in peptide dynamics, which leads to the increased population of structures with fibril-inhibiting characteristics, such as the separated N- and C-termini and colocation of the aromatic residues at the same peptide face. In contrast, lower field strengths (<0.004 V/nm) promote the formation of the amyloid-prone hairpin structures relative to the ambient conditions. These findings suggest that intermediate-strength electromagnetic fields could be considered for designing alternative treatments of amyloid diseases, while the very high and low field strengths could be employed for engineering well-ordered fibrillar aggregates for non-medicinal applications.

Hydrogen/Deuterium Exchange and Molecular Dynamics Analysis of Amyloid Fibrils Formed by a D69K Charge-Pair Mutant of Human Apolipoprotein C-II.

Plasma apolipoproteins form amphipathic α helices in lipid environments but in the lipid-free state show a high propensity to form ß structure and self-associate into amyloid fibrils. The widespread occurrence of apolipoproteins in amyloid plaques suggests disease-related roles, specifically in atherosclerosis. To reconcile the dual abilities of apolipoproteins to form either α helices or cross-ß sheet structures, we examined fibrils formed by human apolipoprotein C-II (apoC-II). A structural model for apoC-II fibrils shows a cross-ß core with parallel ß strands, including a buried K30-D69 charge pair. We investigated the effect of abolishing this charge pair in mutant D69K apoC-II. Fluorescence studies indicated more rapid fibril formation and less solvent accessibility of tryptophan (W26) in D69K apoC-II fibrils than in wild-type (WT) fibrils. X-ray diffraction data of aligned D69K apoC-II fibrils yielded a typical cross-ß structure with increased ß sheet spacing compared to that of WT fibrils. Hydrogen/deuterium (H/D) exchange patterns were similar for D69K apoC-II fibrils compared to WT fibrils, albeit with an overall reduction in the level of slow H/D exchange, particularly around residues 29-32. Molecular dynamics simulations indicated reduced ß strand content for a model D69K apoC-II tetramer compared to the WT tetramer and confirmed an expansion of the cross-ß spacing that contributed to the formation of a stable charge pair between K69 and E27. The results highlight the importance of charge-pair interactions within the apoC-II fibril core, which together with numerous salt bridges in the flexible connecting loop play a major role in the ability of lipid-free apoC-II to form stable cross-ß fibrils.