Functionalized Fullerene for Inhibition of SARS-CoV-2 Variants.

As virus outbreaks continue to pose a challenge, a nonspecific viral inhibitor can provide significant benefits, especially against respiratory viruses. Polyglycerol sulfates recently emerge as promising agents that mediate interactions between cells and viruses through electrostatics, leading to virus inhibition. Similarly, hydrophobic C60 fullerene can prevent virus infection via interactions with hydrophobic cavities of surface proteins. Here, two strategies are combined to inhibit infection of SARS-CoV-2 variants in vitro. Effective inhibitory concentrations in the millimolar range highlight the significance of bare fullerene's hydrophobic moiety and electrostatic interactions of polysulfates with surface proteins of SARS-CoV-2. Furthermore, microscale thermophoresis measurements support that fullerene linear polyglycerol sulfates interact with the SARS-CoV-2 virus via its spike protein, and highlight importance of electrostatic interactions within it. All-atom molecular dynamics simulations reveal that the fullerene binding site is situated close to the receptor binding domain, within 4 nm of polyglycerol sulfate binding sites, feasibly allowing both portions of the material to interact simultaneously.

Charge Matters: Mutations in Omicron Variant Favor Binding to Cells.

Evidence is strengthening to suggest that the novel SARS-CoV-2 mutant Omicron, with its more than 60 mutations, will spread and dominate worldwide. Although the mutations in the spike protein are known, the molecular basis for why the additional mutations in the spike protein that have not previously occurred account for Omicron's higher infection potential, is not understood. We propose, based on chemical rational and molecular dynamics simulations, that the elevated occurrence of positively charged amino acids in certain domains of the spike protein (Delta: +4; Omicron: +5 vs. wild type) increases binding to cellular polyanionic receptors, such as heparan sulfate due to multivalent charge-charge interactions. This observation is a starting point for targeted drug development.

Role of entropy in determining the phase behavior of protein solutions induced by multivalent ions.

Recent experiments have reported lower critical solution temperature (LCST) phase behavior of aqueous solutions of proteins induced by multivalent ions, where the solution phase separates upon heating. This phenomenon is linked to complex hydration effects that result in a net entropy gain upon phase separation. To decipher the underlying molecular mechanism, we use all-atom molecular dynamics simulations along with the two-phase thermodynamic method for entropy calculation. Based on simulations of a single BSA protein in various salt solutions (NaCl, CaCl2, MgCl2, and YCl3) at temperatures (T) ranging 283-323 K, we find that the cation-protein binding affinity increases with T, reflecting its thermodynamic driving force to be entropic in origin. We show that in the cation binding process, many tightly bound water molecules from the solvation shells of a cation and the protein are released to the bulk, resulting in entropy gain. To rationalize the LCST behavior, we calculate the ζ-potential that shows charge inversion of the protein for solutions containing multivalent ions. The ζ-potential increases with T. Performing simulations of two BSA proteins, we demonstrate that the protein-protein binding is mediated by multiple cation bridges and involves similar dehydration effects that cause a large entropy gain which more than compensates for rotational and translational entropy losses of the proteins. Thus, the LCST behavior is entropy-driven, but the associated solvation effects are markedly different from hydrophobic hydration. Our findings have direct implications for tuning the phase behavior of biological and soft-matter systems, e.g., protein condensation and crystallization.

Ultra-high permeable phenine nanotube membranes for water desalination.

Nanopore desalination technology hinges on high water-permeable membranes which, at the same time, block ions efficiently. In this study, we consider a recently synthesized [Science363, 151-155 (2019)] phenine nanotube (PNT) for water desalination applications. Using both equilibrium and non-equilibrium molecular dynamics simulations, we show that the PNT membrane completely rejects salts, but permeates water at a rate which is an order-of-magnitude higher than that of all the membranes used for water filtration. We provide the microscopic mechanisms of salt rejection and fast water-transport by calculating the free-energy landscapes and electrostatic potential profiles. A collective diffusion model accurately predicts the water permeability obtained from the simulations over a wide range of pressure gradients. We propose a method to calculate the osmotic water permeability from the equilibrium simulation data and find that it is very high for the PNT membrane. These remarkable properties of PNT can be applied in various nanofluidic applications, such as ion-selective channels, ionic transistors, sensing, molecular sieving, and blue energy harvesting.

Concerted Interactions between Multiple gp41 Trimers and the Target Cell Lipidome May Be Required for HIV-1 Entry.

The HIV-1 envelope glycoprotein gp41 mediates the fusion between viral and host cell membranes leading to virus entry and target cell infection. Despite years of research, important aspects of this process such as the number of gp41 trimers involved and how they orchestrate the rearrangement of the lipids in the apposed membranes along the fusion pathway remain obscure. To elucidate these molecular underpinnings, we performed coarse-grained molecular dynamics simulations of HIV-1 virions pinned to the CD4 T cell membrane by different numbers of gp41 trimers. We built realistic cell and viral membranes by mimicking their respective lipid compositions. We found that a single gp41 was inadequate for mediating fusion. Lipid mixing between membranes, indicating the onset of fusion, was efficient when three or more gp41 trimers pinned the membranes. The gp41 trimers interacted strongly with many different lipids in the host cell membrane, triggering lipid configurational rearrangements, exchange, and mixing. Simpler membranes, comprising fewer lipid types, displayed strong resistance to fusion, revealing the crucial role of the lipidomes in HIV-1 entry. Performing simulations at different temperatures, we estimated the free energy barrier to lipid mixing, and hence membrane stalk formation, with three and four tethering gp41 trimers to be ∼6.2 kcal/mol, a >4-fold reduction over estimates without gp41. Together, these findings present molecular-level, quantitative insights into the early stages of gp41-mediated HIV-1 entry. Preventing the requisite gp41 molecules from tethering the membranes or altering membrane lipid compositions may be potential intervention strategies.

Polysulfates Block SARS-CoV-2 Uptake through Electrostatic Interactions*.

Here we report that negatively charged polysulfates can bind to the spike protein of SARS-CoV-2 via electrostatic interactions. Using a plaque reduction assay, we compare inhibition of SARS-CoV-2 by heparin, pentosan sulfate, linear polyglycerol sulfate (LPGS) and hyperbranched polyglycerol sulfate (HPGS). Highly sulfated LPGS is the optimal inhibitor, with an IC50 of 67 µg mL-1 (approx. 1.6 µm). This synthetic polysulfate exhibits more than 60-fold higher virus inhibitory activity than heparin (IC50 : 4084 µg mL-1 ), along with much lower anticoagulant activity. Furthermore, in molecular dynamics simulations, we verified that LPGS can bind more strongly to the spike protein than heparin, and that LPGS can interact even more with the spike protein of the new N501Y and E484K variants. Our study demonstrates that the entry of SARS-CoV-2 into host cells can be blocked via electrostatic interactions, therefore LPGS can serve as a blueprint for the design of novel viral inhibitors of SARS-CoV-2.

Understanding enhanced mechanical stability of DNA in the presence of intercalated anticancer drug: Implications for DNA associated processes.

Most of the anticancer drugs bind to double-stranded DNA (dsDNA) by intercalative-binding mode. Although experimental studies have become available recently, a molecular-level understanding of the interactions between the drug and dsDNA that lead to the stability of the intercalated drug is lacking. Of particular interest are the modifications of the mechanical properties of dsDNA observed in experiments. The latter could affect many biological functions, such as DNA transcription and replication. Here, we probe, via all-atom molecular dynamics (MD) simulations, the change in the mechanical properties of intercalated drug-DNA complexes for two intercalators, daunomycin and ethidium. We find that, upon drug intercalation, the stretch modulus of DNA increases significantly, whereas its persistence length and bending modulus decrease. Steered MD simulations reveal that it requires higher forces to stretch the intercalated dsDNA complexes than the normal dsDNA. Adopting various pulling protocols to study force-induced DNA melting, we find that the dissociation of dsDNA becomes difficult in the presence of intercalators. The results obtained here provide a plausible mechanism of function of the anticancer drugs, i.e., via altering the mechanical properties of DNA. We also discuss long-time consequences of using these drugs, which require further in vivo investigations.

Rational design of amphiphilic fluorinated peptides: evaluation of self-assembly properties and hydrogel formation.

Advanced peptide-based nanomaterials composed of self-assembling peptides (SAPs) are of emerging interest in pharmaceutical and biomedical applications. The introduction of fluorine into peptides, in fact, offers unique opportunities to tune their biophysical properties and intermolecular interactions. In particular, the degree of fluorination plays a crucial role in peptide engineering as it can be used to control the characteristics of fluorine-specific interactions and, thus, peptide conformation and self-assembly. Here, we designed and explored a series of amphipathic peptides by incorporating the fluorinated amino acids (2S)-4-monofluoroethylglycine (MfeGly), (2S)-4,4-difluoroethylglycine (DfeGly) and (2S)-4,4,4-trifluoroethylglycine (TfeGly) as hydrophobic components. This approach enabled studying the impact of fluorination on secondary structure formation and peptide self-assembly on a systematic basis. We show that the interplay between polarity and hydrophobicity, both induced differentially by varying degrees of side chain fluorination, does affect peptide folding significantly. A greater degree of fluorination promotes peptide fibrillation and subsequent formation of physical hydrogels in physiological conditions. Molecular simulations revealed the key role played by electrostatically driven intra-chain and inter-chain contact pairs that are modulated by side chain fluorination and give insights into the different self-organization behaviour of selected peptides. Our study provides a systematic report about the distinct features of fluorinated oligomeric peptides with potential applications as peptide-based biomaterials.

Fine-tuning the DNA conductance by intercalation of drug molecules.

In this work we study the structure-transport property relationships of small ligand intercalated DNA molecules using a multiscale modeling approach where extensive ab initio calculations are performed on numerous MD-simulated configurations of dsDNA and dsDNA intercalated with two different intercalators, ethidium and daunomycin. DNA conductance is found to increase by one order of magnitude upon drug intercalation due to the local unwinding of the DNA base pairs adjacent to the intercalated sites, which leads to modifications of the density of states in the near-Fermi-energy region of the ligand-DNA complex. Our study suggests that the intercalators can be used to enhance or tune the DNA conductance, which opens new possibilities for their potential applications in nanoelectronics.

What do we know about DNA mechanics so far?

The DNA molecule, apart from carrying the genetic information, plays a crucial role in a variety of biological processes and finds applications in drug design, nanotechnology and nanoelectronics. The molecule undergoes significant structural transitions under the influence of forces due to physiological and non-physiological environments. Here, we summarize the insights gained from simulations and single-molecule experiments on the structural transitions and mechanics of DNA under force, as well as its elastic properties, in various environmental conditions, and discuss appealing future directions.

Unfolding Dynamics of Ubiquitin from Constant Force MD Simulation: Entropy-Enthalpy Interplay Shapes the Free-Energy Landscape.

Force probe methods are routinely used to study conformational transitions of biomolecules at single-molecule level. In contrast to simple kinetics, some proteins show complex response to mechanical perturbations that is manifested in terms of unusual force-dependent kinetics. Here, we study, via fully atomistic molecular dynamics simulations, constant force-induced unfolding of ubiquitin protein. Our simulations reveal a crossover at an intermediate force (about 400 pN) in the unfolding rate versus force curve. We find by calculation of multidimensional free-energy landscape (FEL) of the protein that the complex unfolding kinetics is intimately related to the force-dependent modifications in the FEL. Pearson correlation coefficient analysis allowed us to identify two appropriate order parameters describing the unfolding transition. The crossover in the rate can be explained in terms of an interplay between entropy and enthalpy with relative importance changing from low force to high force. We rationalize the results by using multidimensional transition-state theory.

Translocation of Bioactive Molecules through Carbon Nanotubes Embedded in the Lipid Membrane.

One of the major challenges of nanomedicine and gene therapy is the effective translocation of drugs and genes across cell membranes. In this study, we describe a systematic procedure that could be useful for efficient drug and gene delivery into the cell. Using fully atomistic molecular dynamics (MD) simulations, we show that molecules of various shapes, sizes, and chemistries can be spontaneously encapsulated in a single-walled carbon nanotube (SWCNT) embedded in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer, as we have exemplified with dendrimers, asiRNA, ssDNA, and ubiquitin protein. We compute the free energy gain by the molecules upon their entry inside the SWCNT channel to quantify the stability of these molecules inside the channel as well as to understand the spontaneity of the process. The free energy profiles suggest that all molecules can enter the channel without facing any energy barrier but experience a strong energy barrier (≫kBT) to translocate across the channel. We propose a theoretical model for the estimation of encapsulation and translocation times of the molecules. Whereas the model predicts the encapsulation time to be of the order of few nanoseconds, which match reasonably well with those obtained from the simulations, it predicts the translocation time to be astronomically large for each molecule considered in this study. This eliminates the possibility of passive diffusion of the molecules through the CNT-nanopore spanning across the membrane. To counter this, we put forward a mechanical method of ejecting the encapsulated molecules by pushing them with other free-floating SWCNTs of diameter smaller than the pore diameter. The feasibility of the proposed method is also demonstrated by performing MD simulations. The generic strategy described here should work for other molecules as well and hence could be potentially useful for drug- and gene-delivery applications.

The SPL7013 dendrimer destabilizes the HIV-1 gp120-CD4 complex.

The poly (l-lysine)-based SPL7013 dendrimer with naphthalene disulphonate surface groups blocks the entry of HIV-1 into target cells and is in clinical trials for development as a topical microbicide. Its mechanism of action against R5 HIV-1, the HIV-1 variant implicated in transmission across individuals, remains poorly understood. Using docking and fully atomistic MD simulations, we find that SPL7013 binds tightly to R5 gp120 in the gp120-CD4 complex but weakly to gp120 alone. Further, the binding, although to multiple regions of gp120, does not occlude the CD4 binding site on gp120, suggesting that SPL7013 does not prevent the binding of R5 gp120 to CD4. Using MD simulations to compute binding energies of several docked structures, we find that SPL7013 binding to gp120 significantly weakens the gp120-CD4 complex. Finally, we use steered molecular dynamics (SMD) to study the kinetics of the dissociation of the gp120-CD4 complex in the absence of the dendrimer and with the dendrimer bound in each of the several stable configurations to gp120. We find that SPL7013 significantly lowers the force required to rupture the gp120-CD4 complex and accelerates its dissociation. Taken together, our findings suggest that SPL7013 compromises the stability of the R5 gp120-CD4 complex, potentially preventing the accrual of the requisite number of gp120-CD4 complexes across the virus-cell interface, thereby blocking virus entry.