Glycan-Driven Formation of Raft-Like Domains with Hierarchical Periodic Nanoarrays on Dendrimersome Synthetic Cells.

The accurate spatial segregation into distinct phases within cell membranes coordinates vital biochemical processes and functionalities in living organisms. One of nature's strategies to localize reactivity is the formation of dynamic raft domains. Most raft models rely on liquid-ordered L0 phases in a liquid-disordered Ld phase lacking correlation and remaining static, often necessitating external agents for phase separation. Here, we introduce a synthetic system of bicomponent glycodendrimersomes coassembled from Janus dendrimers and Janus glycodendrimers (JGDs), where lactose-lactose interactions exclusively drive lateral organization. This mechanism results in modulated phases across two length scales, yielding raft-like microdomains featuring nanoarrays at the nanoscale. By varying the density of lactose and molecular architecture of JGDs, the nanoarray type and size, shape, and spacing of the domains were controlled. Our findings offer insight into the potential primordial origins of rudimentary raft domains and highlight the crucial role of glycans within the glycocalyx.

Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials.

The diffusivity of water in aqueous cesium iodide solutions is larger than that in neat liquid water and vice versa for sodium chloride solutions. Such peculiar ion-specific behavior, called anomalous diffusion, is not reproduced in typical force field based molecular dynamics (MD) simulations due to inadequate treatment of ion-water interactions. Herein, this hurdle is tackled by using machine learned atomic potentials (MLPs) trained on data from density functional theory calculations. MLP based atomistic MD simulations of aqueous salt solutions reproduce experimentally determined thermodynamic, structural, dynamical, and transport properties, including their varied trends in water diffusivities across salt concentration. This enables an examination of their intermolecular structure to unravel the microscopic underpinnings of the differences in their transport properties. While both ions in CsI solutions contribute to the faster diffusion of water molecules, the competition between the heavy retardation by Na ions and the slight acceleration by Cl ions in NaCl solutions reduces their water diffusivity.

Why Dissolving Salt in Water Decreases Its Dielectric Permittivity.

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters.

Ultra-rare complement factor 8 coding variants in families with age-related macular degeneration.

Genome-wide association studies have uncovered 52 independent common and rare variants across 34 genetic loci, which influence susceptibility to age related macular degeneration (AMD). Of the 5 AMD-associated complement genes, complement factor H (CFH) and CFI exhibit a significant rare variant burden implicating a major contribution of the complement pathway to disease pathology. However, the efforts for developing AMD therapy have been challenging as of yet. Here, we report the identification of ultra-rare variants in complement factors 8A and 8B, two components of the terminal complement membrane attack complex (MAC), by whole exome sequencing of a cohort of AMD families. The identified C8 variants impact local interactions among proteins of C8 triplex in vitro, indicating their effect on MAC stability. Our results suggest that MAC, and not the early steps of the complement pathway, might be a more effective target for designing treatments for AMD.

Molecular Engineering of pH-Responsive Anchoring Systems onto Poly(ethylene glycol) Corona.

An adaptive surface that can sense and respond to environmental stimuli is integral to smart functional materials. Here, we report pH-responsive anchoring systems onto the poly(ethylene glycol) (PEG) corona of polymer vesicles. The hydrophobic anchor, pyrene, is reversibly inserted into the PEG corona through the reversible protonation of its covalently linked pH-sensing group. Depending on the pKa of the sensor, the pH-responsive region is engineered from acidic to neutral and basic conditions. The switchable electrostatic repulsion between the sensors contributes to the responsive anchoring behavior. Our findings provide a new responsive binding chemistry for the creation of smart nanomedicine and a nanoreactor.

Probing the Gold/Water Interface with Surface-Specific Spectroscopy.

Water is an integral component in electrochemistry, in the generation of the electric double layer, and in the propagation of the interfacial electric fields into the solution; however, probing the molecular-level structure of interfacial water near functioning electrode surfaces remains challenging. Due to the surface-specificity, sum-frequency-generation (SFG) spectroscopy offers an opportunity to investigate the structure of water near working electrochemical interfaces but probing the hydrogen-bonded structure of water at this buried electrode-electrolyte interface was thought to be impossible. Propagating the laser beams through the solvent leads to a large attenuation of the infrared light due to the absorption of water, and interrogating the interface by sending the laser beams through the electrode normally obscures the SFG spectra due to the large nonlinear response of conduction band electrons. Here, we show that the latter limitation is removed when the gold layer is thin. To demonstrate this, we prepared Au gradient films on CaF2 with a thickness between 0 and 8 nm. SFG spectra of the Au gradient films in contact with H2O and D2O demonstrate that resonant water SFG spectra can be obtained using Au films with a thickness of ∼2 nm or less. The measured spectra are distinctively different from the frequency-dependent Fresnel factors of the interface, suggesting that the features we observe in the OH stretching region indeed do not arise from the nonresonant response of the Au films. With the newfound ability to probe interfacial solvent structure at electrode/aqueous interfaces, we hope to provide insights into more efficient electrolyte composition and electrode design.

Adaptive insertion of a hydrophobic anchor into a poly(ethylene glycol) host for programmable surface functionalization.

Covalent and non-covalent molecular binding are two strategies to tailor surface properties and functions. However, the lack of responsiveness and requirement for specific binding groups makes spatiotemporal control challenging. Here, we report the adaptive insertion of a hydrophobic anchor into a poly(ethylene glycol) (PEG) host as a non-covalent binding strategy for surface functionalization. By using polycyclic aromatic hydrocarbons as the hydrophobic anchor, hydrophilic charged and non-charged functional modules were spontaneously loaded onto PEG corona in 2 min without the assistance of any catalysts and binding groups. The thermodynamically favourable insertion of the hydrophobic anchor can be reversed by pulling the functional module, enabling programmable surface functionalization. We anticipate that the adaptive molecular recognition between the hydrophobic anchor and the PEG host will challenge the hydrophilic understanding of PEG and enhance the progress in nanomedicine, advanced materials and nanotechnology.

Understanding the Dynamics of the Structural States of Cannabinoid Receptors and the Role of Different Modulators.

The cannabinoid receptors CB1R and CB2R are members of the G protein-coupled receptor (GPCR) family. These receptors have recently come to light as possible therapeutic targets for conditions affecting the central nervous system. However, because CB1R is known to have psychoactive side effects, its potential as a drug target is constrained. Therefore, targeting CB2R has become the primary focus of recent research. Using various molecular modeling studies, we analyzed the active, inactive, and intermediate states of both CBRs in this study. We conducted in-depth research on the binding properties of various groups of cannabinoid modulators, including agonists, antagonists, and inverse agonists, with all of the different conformational states of the CBRs. The binding effects of these modulators were studied on various CB structural features, including the movement of the transmembrane helices, the volume of the binding cavity, the internal fluids, and the important GPCR properties. Then, using in vitro experiments and computational modeling, we investigated how vitamin E functions as a lipid modulator to influence THC binding. This comparative examination of modulator binding to CBRs provides significant insight into the mechanisms of structural alterations and ligand affinity, which can directly help in the rational design of selective modulators that target either CB1R or CB2R.

Superhydrophilicity of α-alumina surfaces results from tight binding of interfacial waters to specific aluminols.

HYPOTHESIS: Understanding the microscopic driving force of water wetting is challenging and important for design of materials. The relations between structure, dynamics and hydrogen bonds of interfacial water can be investigated using molecular dynamics simulations. EXPERIMENTS AND SIMULATIONS: Contact angles at the alumina (0001) and (112‾0) surfaces are studied using both classical molecular dynamics simulations and experiments. To test the superhydrophilicity, the free energy cost of removing waters near the interfaces are calculated using the density fluctuations method. The strength of hydrogen bonds is determined by their lifetime and geometry. FINDINGS: Both surfaces are superhydrophilic and the (0001) surface is more hydrophilic. Interactions between surfaces and interfacial waters promote a templating effect whereby the latter are aligned in a pattern that follows the underlying lattice of the surfaces. Translational and rotational dynamics of interfacial water molecules are slower than in bulk water. Hydrogen bonds between water and both surfaces are asymmetric, water-to-aluminol ones are stronger than aluminol-to-water ones. Molecular dynamics simulations eliminate the impacts of surface contamination when measuring contact angles and the results reveal the microscopic origin of the macroscopic superhydrophilicity of alumina surfaces: strong water-to-aluminol hydrogen bonds.

Allosteric inhibitors of the main protease of SARS-CoV-2.

SARS-CoV-2 has raised the alarm to search for effective therapy for this virus. To date several vaccines have been approved but few available drugs reported recently still need approval from FDA. Remdesivir was approved for emergency use only. In this report, the SARS-CoV-2 3CLpro was expressed and purified. By using a FRET-based enzymatic assay, we have screened a library consisting of more than 300 different niclosamide derivatives and identified three molecules JMX0286, JMX0301, and JMX0941 as potent allosteric inhibitors against SARS-CoV-2 3CLpro, with IC50 values similar to that of known covalent inhibitor boceprevir. In a cell-based antiviral assay, these inhibitors can inhibit the virus growth with EC50 in the range of 2-3 µM. The mechanism of action of JMX0286, JMX0301, and JMX0941 were characterized by enzyme kinetics, affinity binding and protein-based substrate digestion. Molecular docking, molecular dynamics (MD) simulations and hydration studies suggested that JMX0286, JMX0301, JMX0941 bind specifically to an allosteric pocket of the SARS-CoV-2 3CL protease. This study provides three potent compounds for further studies.

Potential Pro-Inflammatory Effect of Vitamin E Analogs through Mitigation of Tetrahydrocannabinol (THC) Binding to the Cannabinoid 2 Receptor.

Vitamin E acetate, which is used as a diluent of tetrahydrocannabinol (THC), has been reported as the primary causative agent of e-cigarette, or vaping, product use-associated lung injury (EVALI). Here, we employ in vitro assays, docking, and molecular dynamics (MD) computer simulations to investigate the interaction of vitamin E with the membrane-bound cannabinoid 2 receptor (CB2R), and its role in modulating the binding affinity of THC to CB2R. From the MD simulations, we determined that vitamin E interacts with both CB2R and membrane phospholipids. Notably, the synchronized effect of these interactions likely facilitates vitamin E acting as a lipid modulator for the cannabinoid system. Furthermore, MD simulation and trajectory analysis show that when THC binds to CB2R in the presence of vitamin E, the binding cavity widens, facilitating the entry of water molecules into it, leading to a reduced interaction of THC with CB2R. Additionally, the interaction between THC and vitamin E in solution is stabilized by several H bonds, which can directly limit the interaction of free THCs with CB2R. Overall, both the MD simulations and the in vitro dissociation assay results indicate that THC binding to CB2R is reduced in the presence of vitamin E. Our study discusses the role of vitamin E in limiting the effect of THCs and its implications on the reported pathology of EVALI.

SPICA Force Field for Proteins and Peptides.

A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.

Dissolving salt is not equivalent to applying a pressure on water.

Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water's structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

Extracellular matrix dysfunction in Sorsby patient-derived retinal pigment epithelium.

Sorsby Fundus Dystrophy (SFD) is a rare form of macular degeneration that is clinically similar to age-related macular degeneration (AMD), and a histologic hallmark of SFD is a thick layer of extracellular deposits beneath the retinal pigment epithelium (RPE). Previous studies of SFD patient-induced pluripotent stem cell (iPSC) derived RPE differ as to whether these cultures recapitulate this key clinical feature by forming increased drusenoid deposits. The primary purpose of this study is to examine whether SFD patient-derived iPSC-RPE form basal deposits similar to what is found in affected family member SFD globes and to determine whether SFD iPSC RPE may be more oxidatively stressed. We performed a careful comparison of iPSC RPE from three control individuals, multiple iPSC clones from two SFD patients' iPSC RPE, and post-mortem eyes of affected SFD family members. We also examined the effect of CRISPR-Cas9 gene correction of the S204C TIMP3 mutation on RPE phenotype. Finally, targeted metabolomics with liquid chromatography and mass spectrometry analysis and stable isotope-labeled metabolite analysis were performed to determine whether SFD RPE are more oxidatively stressed. We found that SFD iPSC-RPE formed significantly more sub-RPE deposits (∼6-90 µm in height) compared to control RPE at 8 weeks. These deposits were similar in composition to the thick layer of sub-RPE deposits found in SFD family member globes by immunofluorescence staining and TEM imaging. S204C TIMP3 correction by CRISPR-Cas9 gene editing in SFD iPSC RPE cells resulted in significantly reduced basal laminar and sub-RPE calcium deposits. We detected a ∼18-fold increase in TIMP3 accumulation in the extracellular matrix (ECM) of SFD RPE, and targeted metabolomics showed that intracellular 4-hydroxyproline, a major breakdown product of collagen, is significantly elevated in SFD RPE, suggesting increased ECM turnover. Finally, SFD RPE cells have decreased intracellular reduced glutathione and were found to be more vulnerable to oxidative stress. Our findings suggest that elements of SFD pathology can be demonstrated in culture which may lead to insights into disease mechanisms.

Modeling Liquid Water by Climbing up Jacob's Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials.

Within the framework of Kohn-Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob's ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric-isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

Discovery of Novel Small-Molecule Inhibitors of SARS-CoV-2 Main Protease as Potential Leads for COVID-19 Treatment.

The main protease of SARS-CoV-2 virus, Mpro, is an essential element for viral replication, and inhibitors targeting Mpro are currently being investigated in many drug development programs as a possible treatment for COVID-19. An in vitro pilot screen of a highly focused collection of compounds was initiated to identify new lead scaffolds for Mpro. These efforts identified a number of hits. The most effective of these was compound SIMR-2418 having an inhibitory IC50 value of 20.7 µM. Molecular modeling studies were performed to understand the binding characteristics of the identified compounds. The presence of a cyclohexenone warhead group facilitated covalent binding with the Cys145 residue of Mpro. Our results highlight the challenges of targeting Mpro protease and pave the way toward the discovery of potent lead molecules.

Modeling the Ultrafast Electron Attachment Dynamics of Solvated Uracil.

Electron attachment to DNA by low energy electrons can lead to DNA damage, so a fundamental understanding of how electrons interact with the components of nucleic acids in solution is an open challenge. In solution, low energy electrons can generate presolvated electrons, epre-, which are efficiently scavanged by pyrimidine nucleobases to form transient negative ions, able to relax to either stable valence bound anions or undergo dissociative electron detachment or transfer to other parts of DNA/RNA leading to strand breakages. In order to understand the initial electron attachment dynamics, this paper presents a joint molecular dynamics and high-level electronic structure study into the behavior of the electronic states of the solvated uracil anion. Both the valence π* and nonvalence epre- states of the solvated uracil system are studied, and the effect of the solvent environment and the geometric structure of the uracil core are uncoupled to gain insight into the physical origin of the stabilization of the solvated uracil anion. Solvent reorganization is found to play a dominant role followed by relaxation of the uracil core.

Regulation and drug modulation of a voltage-gated sodium channel: Pivotal role of the S4-S5 linker in activation and slow inactivation.

Voltage-gated sodium (NaV) channels control excitable cell functions. While structural investigations have revealed conformation details of different functional states, the mechanisms of both activation and slow inactivation remain unclear. Here, we identify residue T140 in the S4-S5 linker of the bacterial voltage-gated sodium channel NaChBac as critical for channel activation and drug effects on inactivation. Mutations at T140 either attenuate activation or render the channel nonfunctional. Propofol, a clinical anesthetic known to inhibit NaChBac by promoting slow inactivation, binds to a pocket between the S4-S5 linker and S6 helix in a conformation-dependent manner. Using 19F-NMR to quantify site-specific binding by saturation transfer differences (STDs), we found strong STDs in inactivated, but not activated, NaChBac. Molecular dynamics simulations show a highly dynamic pocket in the activated conformation, limiting STD buildup. In contrast, drug binding to this pocket promotes and stabilizes the inactivated states. Our results provide direct experimental evidence showing distinctly different associations between the S4-S5 linker and S6 helix in activated and inactivated states. Specifically, an exchange occurs between interaction partners T140 and N234 of the same subunit in activation, and T140 and N225 of the domain-swapped subunit in slow inactivation. The drug action on slow inactivation of prokaryotic NaV channels seems to have a mechanism similar to the recently proposed "door-wedge" action of the isoleucine-phenylalanine-methionine (IFM) motif on the fast inactivation of eukaryotic NaV channels. Elucidating this gating mechanism points to a possible direction for conformation-dependent drug development.

Targeting SARS-CoV-2 M3CLpro by HCV NS3/4a Inhibitors: In Silico Modeling and In Vitro Screening.

Currently the entire human population is in the midst of a global pandemic caused by SARS-CoV-2 (Severe Acute Respiratory Syndrome CoronaVirus 2). This highly pathogenic virus has to date caused >71 million infections and >1.6 million deaths in >180 countries. Several vaccines and drugs are being studied as possible treatments or prophylactics of this viral infection. M3CLpro (coronavirus main cysteine protease) is a promising drug target as it has a significant role in viral replication. Here we use the X-ray crystal structure of M3CLpro in complex with boceprevir to study the dynamic changes of the protease upon ligand binding. The binding free energy was calculated for water molecules at different locations of the binding site, and molecular dynamics (MD) simulations were carried out for the M3CLpro/boceprevir complex, to thoroughly understand the chemical environment of the binding site. Several HCV NS3/4a protease inhibitors were tested in vitro against M3CLpro. Specifically, asunaprevir, narlaprevir, paritaprevir, simeprevir, and telaprevir all showed inhibitory effects on M3CLpro. Molecular docking and MD simulations were then performed to investigate the effects of these ligands on M3CLpro and to provide insights into the chemical environment of the ligand binding site. Our findings and observations are offered to help guide the design of possible potent protease inhibitors and aid in coping with the COVID-19 pandemic.

Enhanced Concanavalin†A Binding to Preorganized Mannose Nanoarrays in Glycodendrimersomes Revealed Multivalent Interactions.

The effect of the two-dimensional glycan display on glycan-lectin recognition remains poorly understood despite the importance of these interactions in a plethora of cellular processes, in (patho)physiology, as well as its potential for advanced therapeutics. Faced with this challenge we utilized glycodendrimersomes, a type of synthetic vesicles whose membrane mimics the surface of a cell and offers a means to probe the carbohydrate biological activity. These single-component vesicles were formed by the self-assembly of sequence-defined mannose-Janus dendrimers, which serve as surrogates for glycolipids. Using atomic force microscopy and molecular modeling we demonstrated that even mannose, a monosaccharide, was capable of organizing the sugar moieties into periodic nanoarrays without the need of the formation of liquid-ordered phases as assumed necessary for rafts. Kinetics studies of Concanavalin A binding revealed that those nanoarrays resulted in a new effective ligand yielding a ten-fold increase in the kinetic and thermodynamic constant of association.