Showing papers in "Journal of Physical Chemistry B in 2021"

A Data-Driven Hydrophobicity Scale for Predicting Liquid-Liquid Phase Separation of Proteins.

Abstract: An accurate model for macroscale disordered assemblies of biological macromolecules such as those formed in so-called membraneless organelles would greatly assist in studying their structure, function, and dynamics. Recent evidence has suggested that liquid-liquid phase separation (LLPS) underlies the formation of membraneless organelles. While the general mechanism of exchange of macromolecule/water for macromolecule/macromolecule interactions is known to be the driving force for LLPS, the specific interactions involved are not well understood. One way that protein-water and protein-protein interactions have been understood historically is via hydrophobicity scales. However, these scales are typically optimized for describing these relative interactions in certain cases, such as protein folding or insertion of proteins into membranes. To better describe the relative interactions of proteins that undergo LLPS, we have developed a new, data-driven hydrophobicity scale. To determine the new scale, we used coarse-grained molecular dynamics simulations using the hydrophobicity scale coarse-grained model, which relates the interactions between amino acids to their hydrophobicity. Hydrophobicity values were determined via the force-balance method on a library of proteins that includes unfolded, intrinsically disordered, and phase-separating proteins (PSP). The resulting hydrophobicity scale can better predict whether a given protein will undergo LLPS at physiological conditions by using coarse-grained molecular dynamics simulations than existing hydrophobicity scales. This new scale confirms the importance of π-π interactions between amino acids as important drivers of LLPS. This new hydrophobicity scale provides a convenient and compact description of protein-protein interactions for proteins that undergo LLPS and could be used to develop new models to describe interactions between PSP and other components, such as nucleic acids.

Biomolecular QM/MM Simulations: What Are Some of the "Burning Issues"?

Abstract: QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Che...

Dynamic Network Modeling of Allosteric Interactions and Communication Pathways in the SARS-CoV-2 Spike Trimer Mutants: Differential Modulation of Conformational Landscapes and Signal Transmission via Cascades of Regulatory Switches.

Abstract: The rapidly growing body of structural and biochemical studies of the SARS-CoV-2 spike glycoprotein has revealed a variety of distinct functional states with radically different arrangements of the receptor-binding domain, highlighting a remarkable function-driven conformational plasticity and adaptability of the spike proteins. In this study, we examined molecular mechanisms underlying conformational and dynamic changes in the SARS-CoV-2 spike mutant trimers through the lens of dynamic analysis of allosteric interaction networks and atomistic modeling of signal transmission. Using an integrated approach that combined coarse-grained molecular simulations, protein stability analysis, and perturbation-based modeling of residue interaction networks, we examined how mutations in the regulatory regions of the SARS-CoV-2 spike protein can differentially affect dynamics and allosteric signaling in distinct functional states. The results of this study revealed key functional regions and regulatory centers that govern collective dynamics, allosteric interactions, and control signal transmission in the SARS-CoV-2 spike proteins. We found that the experimentally confirmed regulatory hotspots that dictate dynamic switching between conformational states of the SARS-CoV-2 spike protein correspond to the key hinge sites and global mediating centers of the allosteric interaction networks. The results of this study provide a novel insight into allosteric regulatory mechanisms of SARS-CoV-2 spike proteins showing that mutations at the key regulatory positions can differentially modulate distribution of states and determine topography of signal communication pathways operating through state-specific cascades of control switch points. This analysis provides a plausible strategy for allosteric probing of the conformational equilibrium and therapeutic intervention by targeting specific hotspots of allosteric interactions and communications in the SARS-CoV-2 spike proteins.

Accelerating COVID-19 Research Using Molecular Dynamics Simulation.

Abstract: The COVID-19 pandemic has emerged as a global medico-socio-economic disaster. Given the lack of effective therapeutics against SARS-CoV-2, scientists are racing to disseminate suggestions for rapidly deployable therapeutic options, including drug repurposing and repositioning strategies. Molecular dynamics (MD) simulations have provided the opportunity to make rational scientific breakthroughs in a time of crisis. Advancements in these technologies in recent years have become an indispensable tool for scientists studying protein structure, function, dynamics, interactions, and drug discovery. Integrating the structural data obtained from high-resolution methods with MD simulations has helped in comprehending the process of infection and pathogenesis, as well as the SARS-CoV-2 maturation in host cells, in a short duration of time. It has also guided us to identify and prioritize drug targets and new chemical entities, and to repurpose drugs. Here, we discuss how MD simulation has been explored by the scientific community to accelerate and guide translational research on SARS-CoV-2 in the past year. We have also considered future research directions for researchers, where MD simulations can help fill the existing gaps in COVID-19 research.

Integrated Biophysical Modeling of the SARS-CoV-2 Spike Protein Binding and Allosteric Interactions with Antibodies.

Abstract: Structural and biochemical studies of the severe acute respiratory syndrome (SARS)-CoV-2 spike glycoproteins and complexes with highly potent antibodies have revealed multiple conformation-dependent epitopes highlighting conformational plasticity of spike proteins and capacity for eliciting specific binding and broad neutralization responses. In this study, we used coevolutionary analysis, molecular simulations, and perturbation-based hierarchical network modeling of the SARS-CoV-2 spike protein complexes with a panel of antibodies targeting distinct epitopes to explore molecular mechanisms underlying binding-induced modulation of dynamics and allosteric signaling in the spike proteins. Through coevolutionary analysis of the SARS-CoV-2 spike proteins, we identified highly coevolving hotspots and functional clusters that enable a functional cross-talk between distant allosteric regions in the SARS-CoV-2 spike complexes with antibodies. Coarse-grained and all-atom molecular dynamics simulations combined with mutational sensitivity mapping and perturbation-based profiling of the SARS-CoV-2 receptor-binding domain (RBD) complexes with CR3022 and CB6 antibodies enabled a detailed validation of the proposed approach and an extensive quantitative comparison with the experimental structural and deep mutagenesis scanning data. By combining in silico mutational scanning, perturbation-based modeling, and network analysis of the SARS-CoV-2 spike trimer complexes with H014, S309, S2M11, and S2E12 antibodies, we demonstrated that antibodies can incur specific and functionally relevant changes by modulating allosteric propensities and collective dynamics of the SARS-CoV-2 spike proteins. The results provide a novel insight into regulatory mechanisms of SARS-CoV-2 S proteins showing that antibody-escaping mutations can preferentially target structurally adaptable energy hotspots and allosteric effector centers that control functional movements and allosteric communication in the complexes.

Polymers with Dynamic Bonds: Adaptive Functional Materials for a Sustainable Future.

Abstract: Polymeric materials play critical role in many current technologies. Among them, adaptive polymeric materials with dynamic (reversible) bonds exhibit unique properties and provide exciting opportunities for various future technologies. Dynamic bonds enable structural rearrangements in polymer networks in specific conditions. Replacement of a few covalent bonds by dynamic bonds can enhance polymeric properties, e.g., strongly improve the toughness and the adhesive properties of polymers. Moreover, they provide recyclability and enable new properties, such as self-healing and shape memory effects. We briefly overview new developments in the field of polymers with dynamic bonds and current understanding of their dynamic properties. We further highlight several examples of unique properties of polymers with dynamic bonds and provide our perspectives for them to be used in many current and future applications.

Frustration in fuzzy protein complexes leads to interaction versatility

Abstract: Disordered proteins frequently serve as interaction hubs involving a constrained variety of partners. Complexes with different partners frequently exhibit distinct binding modes, involving regions that remain disordered in the bound state. While the conformational properties of disordered proteins are well-characterized in their free states, less is known about the molecular mechanisms by which specificity can be achieved not with one but with multiple partners. Using the energy landscape theory concept of protein frustration, we demonstrate that complexes of disordered proteins exhibit a high degree of local frustration, especically at the binding interface. These suboptimal interactions lead to the possibility of multiple bound substates, each displaying distinct frustration patterns, which are differently populated in complexes with different partners. These results explain how specificity of disordered proteins can be achieved without a single common bound conformation and how the confliict between different interactions can be used to control the binding to multiple partners.

Critical Interactions Between the SARS-CoV-2 Spike Glycoprotein and the Human ACE2 Receptor.

Abstract: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects human cells by binding its spike (S) glycoproteins to angiotensin-converting enzyme 2 (ACE2) receptors and causes the coronavirus disease 2019 (COVID-19). Therapeutic approaches to prevent SARS-CoV-2 infection are mostly focused on blocking S-ACE2 binding, but critical residues that stabilize this interaction are not well understood. By performing all-atom molecular dynamics (MD) simulations, we identified an extended network of salt bridges, hydrophobic and electrostatic interactions, and hydrogen bonds between the receptor-binding domain (RBD) of the S protein and ACE2. Mutagenesis of these residues on the RBD was not sufficient to destabilize binding but reduced the average work to unbind the S protein from ACE2. In particular, the hydrophobic end of RBD serves as the main anchor site and is the last to unbind from ACE2 under force. We propose that blocking the hydrophobic surface of RBD via neutralizing antibodies could prove to be an effective strategy to inhibit S-ACE2 interactions.

Coarse-Graining Organic Semiconductors: The Path to Multiscale Design.

Abstract: Four decades of molecular theory and computation have helped form the modern understanding of the physical chemistry of organic semiconductors. Whereas these efforts have historically centered arou...

Biomolecular Condensates: Sequence Determinants of Phase Separation, Microstructural Organization, Enzymatic Activity, and Material Properties

Abstract: This perspective article highlights recent progress and emerging challenges in understanding the formation and function of membraneless organelles (MLOs). A long-term goal in the MLO field is to identify the sequence-encoded rules that dictate the formation of compositionally controlled biomolecular condensates, which cells utilize to perform a wide variety of functions. The molecular organization of the different components within a condensate can vary significantly, ranging from a homogeneous mixture to core-shell droplet structures. We provide many examples to highlight the richness of the observed behavior and potential research directions for improving our mechanistic understanding. The tunable environment within condensates can, in principle, alter enzymatic activity significantly. We examine recent examples where this was demonstrated, including applications in synthetic biology. An important question about MLOs is the role of liquid-like material properties in biological function. We discuss the need for improved quantitative characterization tools and the development of sequence-structure-dynamics relationships.

Biophysics of Phase Separation of Disordered Proteins Is Governed by Balance between Short- And Long-Range Interactions.

Abstract: Intrinsically disordered proteins play a crucial role in cellular phase separation, yet the diverse molecular forces driving phase separation are not fully understood. It is of utmost importance to understand how peptide sequence, and particularly the balance between the peptides' short- and long-range interactions with other peptides, may affect the stability, structure, and dynamics of liquid-liquid phase separation in protein condensates. Here, using coarse-grained molecular dynamics simulations, we studied the liquid properties of the condensate in a series of polymers in which the ratio of short-range dispersion interactions to long-range electrostatic interactions varied. As the fraction of mutations that participate in short-range interactions increases at the expense of long-range electrostatic interactions, a significant decrease in the critical temperature of phase separation is observed. Nevertheless, sequences with a high fraction of short-range interactions exhibit stabilization, which suggests compensation for the loss of long-range electrostatic interactions. Decreased condensate stability is coupled with decreased translational diffusion of the polymers in the condensate, which may result in the loss of liquid characteristics in the presence of a high fraction of uncharged residues. The effect of exchanging long-range electrostatic interactions for short-range interactions can be explained by the kinetics of breaking intermolecular contacts with neighboring polymers and the kinetics of intramolecular fluctuations. While both time scales are coupled and increase as electrostatic interactions are lost, for sequences that are dominated by short-range interactions, the kinetics of intermolecular contact breakage significantly slows down. Our study supports the contention that different types of interactions can maintain protein condensates, however, long-range electrostatic interactions enhance its liquid-like behavior.

Water-in-Salt LiTFSI Aqueous Electrolytes. 1. Liquid Structure from Combined Molecular Dynamics Simulation and Experimental Studies.

Abstract: The concept of water-in-salt electrolytes was introduced recently, and these systems have been successfully applied to yield extended operation voltage and hence significantly improved energy density in aqueous Li-ion batteries. In the present work, results of X-ray scattering and Fourier-transform infrared spectra measurements over a wide range of temperatures and salt concentrations are reported for the LiTFSI (lithium bis(trifluoromethane sulfonyl)imide)-based water-in-salt electrolyte. Classical molecular dynamics simulations are validated against the experiments and used to gain additional information about the electrolyte structure. Based on our analyses, a new model for the liquid structure is proposed. Specifically, we demonstrate that at the highest LiTFSI concentration of 20 m the water network is disrupted, and the majority of water molecules exist in the form of isolated monomers, clusters, or small aggregates with chain-like configurations. On the other hand, TFSI- anions are connected to each other and form a network. This description is fundamentally different from those proposed in earlier studies of this system.

Introducing Memory in Coarse-Grained Molecular Simulations.

Abstract: Preserving the correct dynamics at the coarse-grained (CG) level is a pressing problem in the development of systematic CG models in soft matter simulation. Starting from the seminal idea of simple time-scale mapping, there have been many efforts over the years toward establishing a meticulous connection between the CG and fine-grained (FG) dynamics based on fundamental statistical mechanics approaches. One of the most successful attempts in this context has been the development of CG models based on the Mori-Zwanzig (MZ) theory, where the resulting equation of motion has the form of a generalized Langevin equation (GLE) and closely preserves the underlying FG dynamics. In this Review, we describe some of the recent studies in this regard. We focus on the construction and simulation of dynamically consistent systematic CG models based on the GLE, both in the simple Markovian limit and the non-Markovian case. Some recent studies of physical effects of memory are also discussed. The Review is aimed at summarizing recent developments in the field while highlighting the major challenges and possible future directions.

Pressure and Temperature Phase Diagram for Liquid-Liquid Phase Separation of the RNA-Binding Protein Fused in Sarcoma.

Abstract: Liquid-liquid phase separation (LLPS) of proteins and nucleic acids to form membraneless cellular compartments is considered to be involved in various biological functions. The RNA-binding protein fused in sarcoma (FUS) undergoes LLPS in vivo and in vitro. Here, we investigated the effects of pressure and temperature on the LLPS of FUS by high-pressure microscopy and high-pressure UV/vis spectroscopy. The phase-separated condensate of FUS was obliterated with increasing pressure but was observed again at a higher pressure. We generated a pressure-temperature phase diagram that describes the phase separation of FUS and provides a general understanding of the thermodynamic properties of self-assembly and phase separation of proteins. FUS has two types of condensed phases, observed at low pressure (LP-LLPS) and high pressure (HP-LLPS). The HP-LLPS state was more condensed and exhibited lower susceptibility to dissolution by 1,6-hexanediol and karyopherin-β2 than the LP-LLPS state. Moreover, molecular dynamic simulations revealed that electrostatic interactions were destabilized, whereas cation-π, π-π, and hydrophobic interactions were stabilized in HP-LLPS. When cation-π, π-π, and hydrophobic interactions were transiently stabilized in the cellular environment, the phase transition to HP-LLPS occurred; this might be correlated to the aberrant enrichment of cytoplasmic ribonucleoprotein granules, leading to amyotrophic lateral sclerosis.

Temperature-Responsive Peptide-Nucleotide Coacervates.

Abstract: Coacervates are a type of liquid-liquid phase separated (LLPS) droplets that can serve as models of membraneless organelles (MLOs) in living cells. Peptide-nucleotide coacervates have been widely used to mimic properties of ribonucleoprotein (RNP) granules, but the thermal stability and the role of base stacking is still poorly understood. Here, we report a systematic investigation of coacervates formed by five different nucleoside triphosphates (NTPs) with poly-l-lysine and poly-l-arginine as a function of temperature. All studied combinations exhibit an upper critical solution temperature (UCST), and a temperature-dependent critical salt concentration, originating from a significant nonelectrostatic contribution to the mixing free energy. Both the enthalpic and entropic parts of this nonelectrostatic interaction decrease in the order G/A/U/C/T, in accordance with nucleobase stacking free energies. Partitioning of two dyes proves that the local hydrophobicity inside the peptide-nucleotide coacervates is different for every nucleoside triphosphate. We derive a simple relation between the temperature and salt concentration at the critical point based on a mean-field model of phase separation. Finally, when different NTPs are mixed with one common oppositely charged peptide, hybrid coacervates were formed, characterized by a single intermediate UCST and critical salt concentration. NTPs with lower critical salt concentrations can remain condensed in mixed coacervates far beyond their original critical salt concentration. Our results show that NTP-based coacervates have a strong temperature sensitivity due to base stacking interactions and that mixing NTPs can significantly influence the stability of condensates and, by extension, their bioavailability.

Fast Prediction of Binding Affinities of the SARS-CoV-2 Spike Protein Mutant N501Y (UK Variant) with ACE2 and Miniprotein Drug Candidates.

Abstract: A recently identified variant of SARS-CoV-2 virus, known as the United Kingdom (UK) variant (lineage B.1.1.7), has an N501Y mutation on its spike protein. SARS-CoV-2 spike protein binds with angiotensin-converting enzyme 2 (ACE2), a key protein for the viral entry into the host cells. Here, we report an efficient computational approach, including the simple energy minimizations and binding free energy calculations, starting from an experimental structure of the binding complex along with experimental calibration of the calculated binding free energies, to rapidly and reliably predict the binding affinities of the N501Y mutant with human ACE2 (hACE2) and recently reported miniprotein and hACE2 decoy (CTC-445.2) drug candidates. It has been demonstrated that the N501Y mutation markedly increases the ACE2-spike protein binding affinity (Kd) from 22 to 0.44 nM, which could partially explain why the UK variant is more infectious. The miniproteins are predicted to have ∼10,000- to 100,000-fold diminished binding affinities with the N501Y mutant, creating a need for design of novel therapeutic candidates to overcome the N501Y mutation-induced drug resistance. The N501Y mutation is also predicted to decrease the binding affinity of a hACE2 decoy (CTC-445.2) binding with the spike protein by ∼200-fold. This convenient computational approach along with experimental calibration may be similarly used in the future to predict the binding affinities of potential new variants of the spike protein.

UMAP as a Dimensionality Reduction Tool for Molecular Dynamics Simulations of Biomacromolecules: A Comparison Study.

Abstract: Proteins are the molecular machines of life. The multitude of possible conformations that proteins can adopt determines their free-energy landscapes. However, the inherently high dimensionality of a protein free-energy landscape poses a challenge to deciphering how proteins perform their functions. For this reason, dimensionality reduction is an active field of research for molecular biologists. The uniform manifold approximation and projection (UMAP) is a dimensionality reduction method based on a fuzzy topological analysis of data. In the present study, the performance of UMAP is compared with that of other popular dimensionality reduction methods such as t-distributed stochastic neighbor embedding (t-SNE), principal component analysis (PCA), and time-structure independent components analysis (tICA) in the context of analyzing molecular dynamics simulations of the circadian clock protein VIVID. A good dimensionality reduction method should accurately represent the data structure on the projected components. The comparison of the raw high-dimensional data with the projections obtained using different dimensionality reduction methods based on various metrics showed that UMAP has superior performance when compared with linear reduction methods (PCA and tICA) and has competitive performance and scalable computational cost.

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

Abstract: 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.

Ionic Surfactants at Air/Water and Oil/Water Interfaces: A Comparison Based on Molecular Dynamics Simulations.

Abstract: Ionic surfactants are known to build up higher interfacial pressures at oil/water interfaces than at air/water interfaces for the same surfactant bulk concentration. Here, we systematically investi...

Shape-Changing Bottlebrush Polymers

Abstract: Bottlebrush polymers (BBPs), composed of relatively short polymeric side chains densely grafted on a polymer backbone, exhibit many unique characteristics and hold promise for a variety of applications This Perspective focuses on environmentally induced shape-changing behavior of BBPs at interface and in solution, particularly worm/star-globule shape transitions While BBPs with a single type of homopolymer or random copolymer side chains have been shown to undergo pronounced worm-to-globule shape changes in response to external stimuli, the collapsed brushes are unstable and prone to aggregation By introducing a second, solvophilic polymer into the side chains, either as a distinct type of side chain or as the outer block of block copolymer side chains, the collapsed brushes not only are stabilized but also create unimolecular micellar nanostructures, which can be used for, eg, encapsulation and delivery of substances The current challenges in the design, synthesis, and characterization of stimuli-responsive shape-changing BBPs are discussed

A Simple Explicit-Solvent Model of Polyampholyte Phase Behaviors and Its Ramifications for Dielectric Effects in Biomolecular Condensates.

Abstract: Biomolecular condensates such as membraneless organelles, underpinned by liquid-liquid phase separation (LLPS), are important for physiological function, with electrostatics, among other interaction types, being a prominent force in their assembly. Charge interactions of intrinsically disordered proteins (IDPs) and other biomolecules are sensitive to the aqueous dielectric environment. Because the relative permittivity of protein is significantly lower than that of water, the interior of an IDP condensate is expected to be a relatively low-dielectric regime, which aside from its possible functional effects on client molecules should facilitate stronger electrostatic interactions among the scaffold IDPs. To gain insight into this LLPS-induced dielectric heterogeneity, addressing in particular whether a low-dielectric condensed phase entails more favorable LLPS than that posited by assuming IDP electrostatic interactions are uniformly modulated by the higher dielectric constant of the pure solvent, we consider a simplified multiple-chain model of polyampholytes immersed in explicit solvents that are either polarizable or possess a permanent dipole. Notably, simulated phase behaviors of these systems exhibit only minor to moderate differences from those obtained using implicit-solvent models with a uniform relative permittivity equals to that of pure solvent. Buttressed by theoretical treatments developed here using random phase approximation and polymer field-theoretic simulations, these observations indicate a partial compensation of effects between favorable solvent-mediated interactions among the polyampholytes in the condensed phase and favorable polyampholyte-solvent interactions in the dilute phase, often netting only a minor enhancement of overall LLPS propensity from the very dielectric heterogeneity that arises from the LLPS itself. Further ramifications of this principle are discussed.

Presence of a SARS-CoV-2 Protein Enhances Amyloid Formation of Serum Amyloid A.

Abstract: A marker for the severeness and disease progress of COVID-19 is overexpression of serum amyloid A (SAA) to levels that in other diseases are associated with a risk for SAA amyloidosis. To understand whether SAA amyloidosis could also be a long-term risk of SARS-CoV-2 infections, we have used long all-atom molecular dynamic simulations to study the effect of a SARS-CoV-2 protein segment on SAA amyloid formation. Sampling over 40 μs, we find that the presence of the nine-residue segment SK9, located at the C-terminus of the envelope protein, increases the propensity for SAA fibril formation by three mechanisms: it reduces the stability of the lipid-transporting hexamer shifting the equilibrium toward monomers, it increases the frequency of aggregation-prone configurations in the resulting chains, and it raises the stability of SAA fibrils. Our results therefore suggest that SAA amyloidosis and related pathologies may be a long-term risk of SARS-CoV-2 infections.

Electrostatic Interactions Explain the Higher Binding Affinity of the CR3022 Antibody for SARS-CoV-2 than the 4A8 Antibody.

Abstract: A structural understanding of the mechanism by which antibodies bind SARS-CoV-2 at the atomic level is highly desirable as it can tell the development of more effective antibodies to treat Covid-19. Here, we use steered molecular dynamics (SMD) and coarse-grained simulations to estimate the binding affinity of the monoclonal antibodies CR3022 and 4A8 to the SARS-CoV-2 receptor-binding domain (RBD) and SARS-CoV-2 N-terminal domain (NTD). Consistent with experiments, our SMD and coarse-grained simulations both indicate that CR3022 has a higher affinity for SARS-CoV-2 RBD than 4A8 for the NTD, and the coarse-grained simulations indicate the former binds three times stronger to its respective epitope. This finding shows that CR3022 is a candidate for Covid-19 therapy and is likely a better choice than 4A8. Energetic decomposition of the interaction energies between these two complexes reveals that electrostatic interactions explain the difference in the observed binding affinity between the two complexes. This result could lead to a new approach for developing anti-Covid-19 antibodies in which good candidates must contain charged amino acids in the area of contact with the virus.

Molecular Understanding of Homogeneous Nucleation of CO2 Hydrates Using Transition Path Sampling.

Abstract: Carbon dioxide hydrate is a solid built from hydrogen-bond stabilized water cages that encapsulate individual CO2 molecules. As potential candidates for reducing greenhouse gases, hydrates have attracted attention from both the industry and scientific community. Under high pressure and low temperature, hydrates are formed spontaneously from a mixture of CO2 and water via nucleation and growth. Yet, for moderate undercooling, i.e., moderate supersaturation, studying hydrate formation with molecular simulations is very challenging due to the high nucleation barriers involved. We investigate the homogeneous nucleation mechanism of CO2 hydrate as a function of temperature using transition path sampling (TPS), which generates ensembles of unbiased dynamical trajectories across the high barrier between the liquid and solid states. The resulting path ensembles reveal that at high driving force (low temperature), amorphous structures are predominantly formed, with 4151062 cages being the most abundant. With increasing temperature, the nucleation mechanism changes, and 51262 becomes the most abundant cage type, giving rise to the crystalline sI structure. Reaction coordinate analysis can reveal the most important collective variable involved in the mechanism. With increasing temperature, we observe a shift from a single feature (size of the nucleus) to a 2-dimensional (size and cage type) variable as the salient ingredient of the reaction coordinate, and then back to only the nucleus size. This finding is in line with the underlying shift from an amorphous to a crystalline nucleation channel. Modeling such complex phase transformations using transition path sampling gives unbiased insight into the molecular mechanisms toward different polymorphs, and how these are determined by thermodynamics and kinetics. This study will be beneficial for researchers aiming to produce such hydrates with different polymorphic forms.

Challenges Encountered Applying Equilibrium and Nonequilibrium Binding Free Energy Calculations.

Abstract: Binding free energy calculations have become increasingly valuable to drive decision making in drug discovery projects. However, among other issues, inadequate sampling can reduce accuracy, limiting the value of the technique. In this paper, we apply absolute binding free energy calculations to ligands binding to T4 lysozyme L99A and HSP90 using equilibrium and nonequilibrium approaches. We highlight sampling problems encountered in these systems, such as slow side chain rearrangements and slow changes of water placement upon ligand binding. These same types of challenges are also likely to show up in other protein-ligand systems, and we propose some strategies to diagnose and test for such problems in alchemical free energy calculations. We also explore similarities and differences in how the equilibrium and the nonequilibrium approaches handle these problems. Our results show the large amount of work still to be done to make free energy calculations robust and reliable and provide insight for future research in this area.

Direct Observation of Cholesterol Dimers and Tetramers in Lipid Bilayers.

Abstract: Cholesterol is a ubiquitous component of mammalian cell membranes and affects membrane protein function. Although cholesterol-mediated formation of ordered membrane domains has been extensively studied, molecular-level structural information about cholesterol self-association has been absent. Here, we combine solid-state nuclear magnetic resonance (NMR) spectroscopy with all-atom molecular dynamics simulations to determine the oligomeric structure of cholesterol in phospholipid bilayers. Two-dimensional 13C-13C correlation spectra of differentially labeled cholesterol indicate that cholesterol self-associates in a face-to-face fashion at membrane concentrations from 17 to 44 mol %. 2D 13C and 19F spin-counting experiments allowed us to measure the average oligomeric number of these cholesterol clusters. At low cholesterol concentrations of ∼20%, the average cluster size is centered on dimers. At a high cholesterol concentration of 44%, which is representative of virus lipid envelopes and liquid-ordered domains of cell membranes, both dimers and tetramers are observed. The cholesterol dimers are found in both phase-separated membranes that contain sphingomyelin and in disordered and miscible membranes that are free of sphingomyelin. Molecular dynamics simulations support these experimental observations and moreover provide the lifetimes, stabilities, distributions, and structures of these nanoscopic cholesterol clusters. Taken together, these NMR and MD data strongly suggest that dimers are the basic structural unit of cholesterol in phospholipid bilayers. The direct observation of cholesterol dimers and tetramers provides a revised framework for studying cholesterol interactions with membrane proteins to regulate protein functions and for understanding the pathogenic role of cholesterol in diseases.

Transferable Ion Force Fields in Water from a Simultaneous Optimization of Ion Solvation and Ion-Ion Interaction.

Abstract: The poor performance of many existing nonpolarizable ion force fields is typically blamed on either the lack of explicit polarizability, the absence of charge transfer, or the use of unreduced Coulomb interactions. However, this analysis disregards the large and mostly unexplored parameter range offered by the Lennard-Jones potential. We use a global optimization procedure to develop water-model-transferable force fields for the ions K+, Na+, Cl-, and Br- in the complete parameter space of all Lennard-Jones interactions using standard mixing rules. No extra-thermodynamic assumption is necessary for the simultaneous optimization of the four ion pairs. After an optimization with respect to the experimental solvation free energy and activity, the force fields reproduce the concentration-dependent density, ionic conductivity, and dielectric constant with high accuracy. The force field is fully transferable between simple point charge/extended and transferable intermolecular potential water models. Our results show that a thermodynamically consistent force field for these ions needs only Lennard-Jones and standard Coulomb interactions.

Exploring the Minimum-Energy Pathways and Free-Energy Profiles of Enzymatic Reactions with QM/MM Calculations.

Abstract: Understanding molecular mechanisms of enzymatic reactions is of vital importance in biochemistry and biophysics. Here, we introduce new functions of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in the GENESIS program to compute the minimum-energy pathways (MEPs) and free-energy profiles of enzymatic reactions. For this purpose, an interface in GENESIS is developed to utilize a highly parallel electronic structure program, QSimulate-QM (https://qsimulate.com), calling it as a shared library from GENESIS. Second, algorithms to search the MEP are implemented, combining the string method (E et al. J. Chem. Phys. 2007, 126, 164103) with the energy minimization of the buffer MM region. The method implemented in GENESIS is applied to an enzyme, triosephosphate isomerase, which converts dihyroxyacetone phosphate to glyceraldehyde 3-phosphate in four proton-transfer processes. QM/MM-molecular dynamics simulations show performances of greater than 1 ns/day with the density functional tight binding (DFTB), and 10-30 ps/day with the hybrid density functional theory, B3LYP-D3. These performances allow us to compute not only MEP but also the potential of mean force (PMF) of the enzymatic reactions using the QM/MM calculations. The barrier height obtained as 13 kcal mol-1 with B3LYP-D3 in the QM/MM calculation is in agreement with the experimental results. The impact of conformational sampling in PMF calculations and the level of electronic structure calculations (DFTB vs B3LYP-D3) suggests reliable computational protocols for enzymatic reactions without high computational costs.

Binding of SARS-CoV-2 Fusion Peptide to Host Endosome and Plasma Membrane.

Abstract: During infection the SARS-CoV-2 virus fuses its viral envelope with cellular membranes of its human host. The viral spike (S) protein mediates both the initial contact with the host cell and the subsequent membrane fusion. Proteolytic cleavage of S at the S2' site exposes its fusion peptide (FP) as the new N-terminus. By binding to the host membrane, the FP anchors the virus to the host cell. The reorganization of S2 between virus and host then pulls the two membranes together. Here we use molecular dynamics (MD) simulations to study the two core functions of the SARS-CoV-2 FP: to attach quickly to cellular membranes and to form an anchor strong enough to withstand the mechanical force during membrane fusion. In eight 10 μs long MD simulations of FP in proximity to endosomal and plasma membranes, we find that FP binds spontaneously to the membranes and that binding proceeds predominantly by insertion of two short amphipathic helices into the membrane interface. Connected via a flexible linker, the two helices can bind the membrane independently, yet binding of one promotes the binding of the other by tethering it close to the target membrane. By simulating mechanical pulling forces acting on the C-terminus of the FP, we then show that the bound FP can bear forces up to 250 pN before detaching from the membrane. This detachment force is more than 10-fold higher than an estimate of the force required to pull host and viral membranes together for fusion. We identify a fully conserved disulfide bridge in the FP as a major factor for the high mechanical stability of the FP membrane anchor. We conclude, first, that the sequential binding of two short amphipathic helices allows the SARS-CoV-2 FP to insert quickly into the target membrane, before the virion is swept away after shedding the S1 domain connecting it to the host cell receptor. Second, we conclude that the double attachment and the conserved disulfide bridge establish the strong anchoring required for subsequent membrane fusion. Multiple distinct membrane-anchoring elements ensure high avidity and high mechanical strength of FP-membrane binding.

A Brief Guide to the Structure of High-Temperature Molten Salts and Key Aspects Making Them Different from Their Low-Temperature Relatives, the Ionic Liquids.

Abstract: High-temperature molten salt research is undergoing somewhat of a renaissance these days due to the apparent advantage of these systems in areas related to clean and sustainable energy harvesting and transfer In many ways, this is a mature field with decades if not already a century of outstanding work devoted to it Yet, much of this work was done with pioneering experimental and computational setups that lack the current day capabilities of synchrotrons and high-performance-computing systems resulting in deeply entrenched results in the literature that when carefully inspected may require revision Yet, in other cases, access to isotopically substituted ions make those pioneering studies very unique and prohibitively expensive to carry out nowadays There are many review articles on molten salts, some of them cited in this perspective, that are simply outstanding and we dare not try to outdo those Instead, having worked for almost a couple of decades already on their low-temperature relatives, the ionic liquids, this is the perspective article that some of the authors would have wanted to read when embarking on their research journey on high-temperature molten salts We hope that this will serve as a simple guide to those expanding from research on ionic liquids to molten salts and vice versa, particularly, when looking into their bulk structural features The article does not aim at being comprehensive but instead focuses on selected topics such as short- and intermediate-range order, the constraints on force field requirements, and other details that make the high- and low-temperature ionic melts in some ways similar but in others diametrically opposite