Neuropilin-1 Protein May Serve as a Receptor for SARS-CoV-2 Infection: Evidence from Molecular Dynamics Simulations.

The binding of the virus to host cells is the first step in viral infection. Human cell angiotensin converting enzyme 2 (ACE2) is the most popular receptor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), while other receptors have recently been observed in experiments. Neuropilin-1 protein (NRP1) is one of them, but the mechanism of its binding to the wild type (WT) and different variants of the virus remain unclear at the atomic level. In this work, all-atom umbrella sampling simulations were performed to clarify the binding mechanism of NRP1 to the spike protein fragments 679-685 of the WT, Delta, and Omicron BA.1 variants. We found that the Delta variant binds most strongly to NRP1, while the affinity for Omicron BA.1 slightly decreases for NRP1 compared to that of WT, and the van der Waals interaction plays a key role in stabilizing the studied complexes. The change in the protonation state of the His amino acid results in different binding free energies between variants. Consistent with the experiment, decreasing the pH was shown to increase the binding affinity of the virus to NRP1. Our results indicate that Delta and Omicron mutations not only affect fusogenicity but also affect NRP1 binding. In addition, we argue that viral evolution does not further improve NRP1 binding affinity which remains in the µM range but may increase immune evasion.

Binding of SARS-CoV-2 Nonstructural Protein 1 to 40S Ribosome Inhibits mRNA Translation.

Experimental evidence has established that SARS-CoV-2 NSP1 acts as a factor that restricts cellular gene expression and impedes mRNA translation within the ribosome's 40S subunit. However, the precise molecular mechanisms underlying this phenomenon have remained elusive. To elucidate this issue, we employed a combination of all-atom steered molecular dynamics and coarse-grained alchemical simulations to explore the binding affinity of mRNA to the 40S ribosome, both in the presence and absence of SARS-CoV-2 NSP1. Our investigations revealed that the binding of SARS-CoV-2 NSP1 to the 40S ribosome leads to a significant enhancement in the binding affinity of mRNA. This observation, which aligns with experimental findings, strongly suggests that SARS-CoV-2 NSP1 has the capability to inhibit mRNA translation. Furthermore, we identified electrostatic interactions between mRNA and the 40S ribosome as the primary driving force behind mRNA translation. Notably, water molecules were found to play a pivotal role in stabilizing the mRNA-40S ribosome complex, underscoring their significance in this process. We successfully pinpointed the specific SARS-CoV-2 NSP1 residues that play a critical role in triggering the translation arrest.

A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta.

Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-ß (Aß42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the ß-content of Aß42; however, it increases the tendency of the monomeric Aß42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aß42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aß42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.

Serotonin Promotes Vesicular Association and Fusion by Modifying Lipid Bilayers.

The primary event in chemical neurotransmission involves the fusion of a membrane-limited vesicle at the plasma membrane and the subsequent release of its chemical neurotransmitter cargo. The cargo itself is not known to have any effect on the fusion event. However, amphiphilic monoamine neurotransmitters (e.g., serotonin and dopamine) are known to strongly interact with lipid bilayers and to affect their mechanical properties, which can in principle impact membrane-mediated processes. Here, we probe whether serotonin can enhance the association and fusion of artificial lipid vesicles in vitro. We employ fluorescence correlation spectroscopy and total internal reflection fluorescence microscopy to measure the attachment and fusion of vesicles whose lipid compositions mimic the major lipid components of synaptic vesicles. We find that the association between vesicles and supported lipid bilayers is strongly enhanced in a serotonin dose-dependent manner, and this drives an increase in the rate of spontaneous fusion. Molecular dynamics simulations and fluorescence spectroscopy data show that serotonin insertion increases the water content of the hydrophobic part of the bilayer. This suggests that the enhanced membrane association is likely driven by an energetically favorable drying transition. Other monoamines, such as dopamine and norepinephrine, but not other related species, such as tryptophan, show similar effects on membrane association. Our results reveal a lipid bilayer-mediated mechanism by which monoamines can themselves modulate vesicle fusion, potentially adding to the control toolbox for the tightly regulated process of neurotransmission in vivo.

Exploring Structural Insights of Aß42 and α-Synuclein Monomers and Heterodimer: A Comparative Study Using Implicit and Explicit Solvent Simulations.

Protein misfolding, aggregation, and fibril formation play a central role in the development of severe neurological disorders, including Alzheimer's and Parkinson's diseases. The structural stability of mature fibrils in these diseases is of great importance, as organisms struggle to effectively eliminate amyloid plaques. To address this issue, it is crucial to investigate the early stages of fibril formation when monomers aggregate into small, toxic, and soluble oligomers. However, these structures are inherently disordered, making them challenging to study through experimental approaches. Recently, it has been shown experimentally that amyloid-ß 42 (Aß42) and α-synuclein (α-Syn) can coassemble. This has motivated us to investigate the interaction between their monomers as a first step toward exploring the possibility of forming heterodimeric complexes. In particular, our study involves the utilization of various Amber and CHARMM force-fields, employing both implicit and explicit solvent models in replica exchange and conventional simulation modes. This comprehensive approach allowed us to assess the strengths and weaknesses of these solvent models and force fields in comparison to experimental and theoretical findings, ensuring the highest level of robustness. Our investigations revealed that Aß42 and α-Syn monomers can indeed form stable heterodimers, and the resulting heterodimeric model exhibits stronger interactions compared to the Aß42 dimer. The binding of α-Syn to Aß42 reduces the propensity of Aß42 to adopt fibril-prone conformations and induces significant changes in its conformational properties. Notably, in AMBER-FB15 and CHARMM36m force fields with the use of explicit solvent, the presence of Aß42 significantly increases the ß-content of α-Syn, consistent with the experiments showing that Aß42 triggers α-Syn aggregation. Our analysis clearly shows that although the use of implicit solvent resulted in too large compactness of monomeric α-Syn, structural properties of monomeric Aß42 and the heterodimer were preserved in explicit-solvent simulations. We anticipate that our study sheds light on the interaction between α-Syn and Aß42 proteins, thus providing the atom-level model required to assess the initial stage of aggregation mechanisms related to Alzheimer's and Parkinson's diseases.

SARS-CoV-2 Omicron Subvariants Do Not Differ Much in Binding Affinity to Human ACE2: A Molecular Dynamics Study.

The emergence of the variant of concern Omicron (B.1.1.529) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exacerbates the COVID-19 pandemic due to its high contagious ability. Studies have shown that the Omicron binds human ACE2 more strongly than the wild type. The prevalence of Omicron in new cases of COVID-19 promotes novel lineages with improved receptor binding affinity and immune evasion. To shed light on this open problem, in this work, we investigated the binding free energy of the receptor binding domain of the Omicron lineages BA.2, BA.2.3.20, BA.3, BA4/BA5, BA.2.75, BA.2.75.2, BA.4.6, XBB.1, XBB.1.5, BJ.1, BN.1, BQ.1.1, and CH.1.1 to human ACE2 using all-atom molecular dynamics simulation and the molecular mechanics Poisson-Boltzmann surface area method. The results show that these lineages have increased binding affinity compared to the BA.1 lineage, and BA.2.75 and BA.2.75.2 subvariants bind ACE2 more strongly than others. However, in general, the binding affinities of the Omicron lineages do not differ significantly from each other. The electrostatic force dominates over the van der Waals force in the interaction between Omicron lineages and human cells. Based on our results, we argue that viral evolution does not further improve the affinity of SARS-CoV-2 for ACE2 but may increase immune evasion.

Synonymous Mutations Can Alter Protein Dimerization Through Localized Interface Misfolding Involving Self-entanglements.

Synonymous mutations in messenger RNAs (mRNAs) can reduce protein-protein binding substantially without changing the protein's amino acid sequence. Here, we use coarse-grain simulations of protein synthesis, post-translational dynamics, and dimerization to understand how synonymous mutations can influence the dimerization of two E. coli homodimers, oligoribonuclease and ribonuclease T. We synthesize each protein from its wildtype, fastest- and slowest-translating synonymous mRNAs in silico and calculate the ensemble-averaged interaction energy between the resulting dimers. We find synonymous mutations alter oligoribonuclease's dimer properties. Relative to wildtype, the dimer interaction energy becomes 4% and 10% stronger, respectively, when translated from its fastest- and slowest-translating mRNAs. Ribonuclease T dimerization, however, is insensitive to synonymous mutations. The structural and kinetic origin of these changes are misfolded states containing non-covalent lasso-entanglements, many of which structurally perturb the dimer interface, and whose probability of occurrence depends on translation speed. These entangled states are kinetic traps that persist for long time scales. Entanglements cause altered dimerization energies for oligoribonuclease, as there is a large association (odds ratio: 52) between the co-occurrence of non-native self-entanglements and weak-binding dimer conformations. Simulated at all-atom resolution, these entangled structures persist for long timescales, indicating the conclusions are independent of model resolution. Finally, we show that regions of the protein we predict to have changes in entanglement are also structurally perturbed during refolding, as detected by limited-proteolysis mass spectrometry. Thus, non-native changes in entanglement at dimer interfaces is a mechanism through which oligomer structure and stability can be altered.

Deciphering the free energy landscapes of SARS-CoV-2 wild type and Omicron variant interacting with human ACE2.

The binding of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein to the host cell receptor angiotensin-converting enzyme 2 (ACE2) is the first step in human viral infection. Therefore, understanding the mechanism of interaction between RBD and ACE2 at the molecular level is critical for the prevention of COVID-19, as more variants of concern, such as Omicron, appear. Recently, atomic force microscopy has been applied to characterize the free energy landscape of the RBD-ACE2 complex, including estimation of the distance between the transition state and the bound state, xu. Here, using a coarse-grained model and replica-exchange umbrella sampling, we studied the free energy landscape of both the wild type and Omicron subvariants BA.1 and XBB.1.5 interacting with ACE2. In agreement with experiment, we find that the wild type and Omicron subvariants have similar xu values, but Omicron binds ACE2 more strongly than the wild type, having a lower dissociation constant KD.

Identification of potential anti-hyperglycemic compounds in Cordyceps militaris ethyl acetate extract: in vitro and in silico studies.

Cordyceps militaris has been long known for valuable health benefits by folk experience and was recently reported with diabetes-tackling evidences, thus deserving extending efforts on screening for component-activity relationship. In this study, experiments were carried out to find the evidence, justification, and input for computations on the potential against diabetes-related protein structures: PDB-4W93, PDB-3W37, and PDB-4A3A. Liquid chromatography identified 14 bioactive compounds in the ethyl acetate extract (1-14) and quantified the contents of cordycepin (0.11%) and adenosine (0.01%). Bioassays revealed the overall potential of the extract against α-amylase (IC50 = 6.443 ± 0.364 mg.mL-1) and α-glucosidase (IC50 = 2.580 ± 0.194 mg.mL-1). A combination of different computational platforms was used to select the most promising candidates for applications as anti-diabetic bio-inhibitors, i.e. 1 (ground state: -888.49715 a.u.; dipole moment 3.779 Debye; DS¯ -12.3 kcal.mol-1; polarizability 34.7 Å3; logP - 1.30), 10 (ground state: -688.52406 a.u.; dipole moment 5.487 Debye; DS¯ -12.6 kcal.mol-1; polarizability 24.9 Å3; logP - 3.39), and 12 (ground state: -1460.07276 a.u.; dipole moment 3.976 Debye; DS¯ -12.5 kcal.mol-1; polarizability 52.4 Å3; logP - 4.39). The results encourage further experimental tests on cordycepin (1), mannitol (10), and adenosylribose (12) to validate their in-practice diabetes-related activities, thus conducive to hypoglycemic applications.Communicated by Ramaswamy H. Sarma.

Amyloid-ß Tetramers and Divalent Cations at the Membrane/Water Interface: Simple Models Support a Functional Role.

Charge polarization at the membrane interface is a fundamental process in biology. Despite the lower concentration compared to the abundant monovalent ions, the relative abundance of divalent cations (Ca2+, Mg2+, Zn2+, Fe2+, Cu2+) in particular spaces, such as the neuron synapse, raised many questions on the possible effects of free multivalent ions and of the required protection of membranes by the eventual defects caused by the free forms of the cations. In this work, we first applied a recent realistic model of divalent cations to a well-investigated model of a polar lipid bilayer, di-myristoyl phosphatidyl choline (DMPC). The full atomistic model allows a fairly good description of changes in the hydration of charged and polar groups upon the association of cations to lipid atoms. The lipid-bound configurations were analyzed in detail. In parallel, amyloid-ß 1-42 (Aß42) peptides assembled into tetramers were modeled at the surface of the same bilayer. Two of the protein tetramers' models were loaded with four Cu2+ ions, the latter bound as in DMPC-free Aß42 oligomers. The two Cu-bound models differ in the binding topology: one with each Cu ion binding each of the monomers in the tetramer; one with pairs of Cu ions linking two monomers into dimers, forming tetramers as dimers of dimers. The models here described provide hints on the possible role of Cu ions in synaptic plasticity and of Aß42 oligomers in storing the same ions away from lipids. The release of structurally disordered peptides in the synapse can be a mechanism to recover ion homeostasis and lipid membranes from changes in the divalent cation concentration.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

How soluble misfolded proteins bypass chaperones at the molecular level.

Subpopulations of soluble, misfolded proteins can bypass chaperones within cells. The extent of this phenomenon and how it happens at the molecular level are unknown. Through a meta-analysis of the experimental literature we find that in all quantitative protein refolding studies there is always a subpopulation of soluble but misfolded protein that does not fold in the presence of one or more chaperones, and can take days or longer to do so. Thus, some misfolded subpopulations commonly bypass chaperones. Using multi-scale simulation models we observe that the misfolded structures that bypass various chaperones can do so because their structures are highly native like, leading to a situation where chaperones do not distinguish between the folded and near-native-misfolded states. More broadly, these results provide a mechanism by which long-time scale changes in protein structure and function can persist in cells because some misfolded states can bypass components of the proteostasis machinery.

Is Posttranslational Folding More Efficient Than Refolding from a Denatured State: A Computational Study.

The folding of proteins into their native conformation is a complex process that has been extensively studied over the past half-century. The ribosome, the molecular machine responsible for protein synthesis, is known to interact with nascent proteins, adding further complexity to the protein folding landscape. Consequently, it is unclear whether the folding pathways of proteins are conserved on and off the ribosome. The main question remains: to what extent does the ribosome help proteins fold? To address this question, we used coarse-grained molecular dynamics simulations to compare the mechanisms by which the proteins dihydrofolate reductase, type III chloramphenicol acetyltransferase, and d-alanine-d-alanine ligase B fold during and after vectorial synthesis on the ribosome to folding from the full-length unfolded state in bulk solution. Our results reveal that the influence of the ribosome on protein folding mechanisms varies depending on the size and complexity of the protein. Specifically, for a small protein with a simple fold, the ribosome facilitates efficient folding by helping the nascent protein avoid misfolded conformations. However, for larger and more complex proteins, the ribosome does not promote folding and may contribute to the formation of intermediate misfolded states cotranslationally. These misfolded states persist posttranslationally and do not convert to the native state during the 6 µs runtime of our coarse-grain simulations. Overall, our study highlights the complex interplay between the ribosome and protein folding and provides insight into the mechanisms of protein folding on and off the ribosome.

Unusual Robustness of Neurotransmitter Vesicle Membranes against Serotonin-Induced Perturbations.

Nature confines hundreds of millimolar of amphiphilic neurotransmitters, such as serotonin, in synaptic vesicles. This appears to be a puzzle, as the mechanical properties of lipid bilayer membranes of individual major polar lipid constituents of synaptic vesicles [phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylserine (PS)] are significantly affected by serotonin, sometimes even at few millimolar concentrations. These properties are measured by atomic force microscopy, and their results are corroborated by molecular dynamics simulations. Complementary 2H solid-state NMR measurements also show that the lipid acyl chain order parameters are strongly affected by serotonin. The resolution of the puzzle lies in the remarkably different properties displayed by the mixture of these lipids, at molar ratios mimicking those of natural vesicles (PC:PE:PS:Cholesterol = 3:5:2:5). Bilayers constituting of these lipids are minimally perturbed by serotonin, and show only a graded response at physiological concentrations (>100 mM). Significantly, the cholesterol (up to 33% molar ratio) plays only a minor role in dictating these mechanical perturbations, with PC:PE:PS:Cholesterol = 3:5:2:5 and 3:5:2:0 showing similar perturbations. We infer that nature uses an emergent mechanical property of a specific mixture of lipids, all individually vulnerable to serotonin, to appropriately respond to physiological serotonin levels.

Identifying inhibitors of NSP16-NSP10 of SARS-CoV-2 from large databases.

The COVID-19 pandemic, which has already claimed millions of lives, continues to pose a serious threat to human health, requiring the development of new effective drugs. Non-structural proteins of SARS-CoV-2 play an important role in viral replication and infection. Among them, NSP16 (non-structured protein 16) and its cofactor NSP10 (non-structured protein 10) perform C2'-O methylation at the 5' end of the viral RNA, which promotes efficient virus replication. Therefore, the NSP16-NSP10 complex becomes an attractive target for drug development. Using a multi-step virtual screening protocol which includes Lipinski's rule, docking, steered molecular dynamics and umbrella sampling, we searched for potential inhibitors from the PubChem and anti-HIV databases. It has been shown that CID 135566620 compound from PubChem is the best candidate with an inhibition constant in the sub-µM range. The Van der Waals interaction was found to be more important than the electrostatic interaction in the binding affinity of this compound to NSP16-NSP10. Further in vitro and in vivo studies are needed to test the activity of the identified compound against COVID-19.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulation of cancer cell membrane perforated by shockwave induced bubble collapse.

It has been widely accepted that cancer cells are softer than their normal counterparts. This motivates us to propose, as a proof-of-concept, a method for the efficient delivery of therapeutic agents into cancer cells, while normal cells are less affected. The basic idea of this method is to use a water jet generated by the collapse of the bubble under shockwaves to perforate pores in the cell membrane. Given a combination of shockwave and bubble parameters, the cancer membrane is more susceptible to bending, stretching, and perforating than the normal membrane because the bending modulus of the cancer cell membrane is smaller than that of the normal cell membrane. Therefore, the therapeutic agent delivery into cancer cells is easier than in normal cells. Adopting two well-studied models of the normal and cancer membranes, we perform shockwave induced bubble collapse molecular dynamics simulations to investigate the difference in the response of two membranes over a range of shockwave impulse 15-30 mPa s and bubble diameter 4-10 nm. The simulation shows that the presence of bubbles is essential for generating a water jet, which is required for perforation; otherwise, pores are not formed. Given a set of shockwave impulse and bubble parameters, the pore area in the cancer membrane is always larger than that in the normal membrane. However, a too strong shockwave and/or too large bubble results in too fast disruption of membranes, and pore areas are similar between two membrane types. The pore closure time in the cancer membrane is slower than that in the normal membrane. The implications of our results for applications in real cells are discussed in some details. Our simulation may be useful for encouraging future experimental work on novel approaches for cancer treatment.

Protein aggregation rate depends on mechanical stability of fibrillar structure.

The formation of the fibrillar structure of amyloid proteins/peptides is believed to be associated with neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. Since the rate of aggregation can influence neurotoxicity, finding the key factors that control this rate is of paramount importance. It was recently found that the rate of protein aggregation is related to the mechanical stability of the fibrillar structure such that the higher the mechanical stability, the faster the fibril is formed. However, this conclusion was supported by a limited dataset. In this work, we expand the previous study to a larger dataset, including the wild type of Aß42 peptide and its 20 mutants, the aggregation rate of which was measured experimentally. By using all-atom steered molecular dynamics (SMD) simulations, we can assess the mechanical stability of the fibril structure, which is characterized by the rupture force, pulling work, and unbinding free energy barrier. Our result confirms that mechanical stability is indeed related to the aggregation rate. Since the estimation of the aggregation rate using all-atom simulations is almost forbidden by the current computational capabilities, our result is useful for predicting it based on information obtained from fast SMD simulations for fibrils.

Antibody-nanobody combination increases their neutralizing activity against SARS-CoV-2 and nanobody H11-H4 is effective against Alpha, Kappa and Delta variants.

The global spread of COVID-19 is devastating health systems and economies worldwide. While the use of vaccines has yielded encouraging results, the emergence of new variants of SARS-CoV-2 shows that combating COVID-19 remains a big challenge. One of the most promising treatments is the use of not only antibodies, but also nanobodies. Recent experimental studies revealed that the combination of antibody and nanobody can significantly improve their neutralizing ability through binding to the SARS-CoV-2 spike protein, but the molecular mechanisms underlying this observation remain largely unknown. In this work, we investigated the binding affinity of the CR3022 antibody and H11-H4 nanobody to the SARS-CoV-2 receptor binding domain (RBD) using molecular modeling. Both all-atom steered molecular dynamics simulations and coarse-grained umbrella sampling showed that, consistent with the experiment, CR3022 associates with RBD more strongly than H11-H4. We predict that the combination of CR3022 and H11-H4 considerably increases their binding affinity to the spike protein. The electrostatic interaction was found to control the association strength of CR3022, but the van der Waals interaction dominates in the case of H11-H4. However, our study for a larger set of nanobodies and antibodies showed that the relative role of these interactions depends on the specific complex. Importantly, we showed Beta, Gamma, Lambda, and Mu variants reduce the H11-H4 activity while Alpha, Kappa and Delta variants increase its neutralizing ability, which is in line with experiment reporting that the nanobody elicited from the llama is very promising for fighting against the Delta variant.

SARS-CoV-2 Omicron Variant Binds to Human Cells More Strongly than the Wild Type: Evidence from Molecular Dynamics Simulation.

The emergence of the variant of concern Omicron (B.1.1.529) of the severe acute respiratory syndrome coronavirus 2 has aggravated the Covid-19 pandemic due to its very contagious ability. The high infection rate may be due to the high binding affinity of Omicron to human cells, but both experimental and computational studies have yielded conflicting results on this issue. Some studies have shown that the Omicron variant binds to human angiotensin-converting enzyme 2 (hACE2) more strongly than the wild type (WT), but other studies have reported comparable binding affinities. To shed light on this open problem, in this work, we calculated the binding free energy of the receptor binding domain (RBD) of the WT and Omicron spike protein to hACE2 using all-atom molecular dynamics simulation and the molecular mechanics Poisson-Boltzmann surface area method. We showed that Omicron binds to human cells more strongly than the WT due to increased RBD charge, which enhances electrostatic interaction with negatively charged hACE2. N440K, T478K, E484A, Q493R, and Q498R mutations in the RBD have been found to play a critical role in the stability of the RBD-hACE2 complex. The effect of homogeneous and heterogeneous models of glycans coating the viral RBD and the peptidyl domain of hACE2 was examined. Although the total binding free energy is not sensitive to the glycan model, the distribution of per-residue interaction energies depends on it. In addition, glycans have a little effect on the binding affinity of the WT RBD to hACE2.

Determination of Multidirectional Pathways for Ligand Release from the Receptor: A New Approach Based on Differential Evolution.

Steered molecular dynamics (SMD) simulation is a powerful method in computer-aided drug design as it can be used to access the relative binding affinity with high precision but with low computational cost. The success of SMD depends on the choice of the direction along which the ligand is pulled from the receptor-binding site. In most simulations, the unidirectional pathway was used, but in some cases, this choice resulted in the ligand colliding with the complex surface of the exit tunnel. To overcome this difficulty, several variants of SMD with multidirectional pulling have been proposed, but they are not completely devoid of disadvantages. Here, we have proposed to determine the direction of pulling with a simple scoring function that minimizes the receptor-ligand interaction, and an optimization algorithm called differential evolution is used for energy minimization. The effectiveness of our protocol was demonstrated by finding expulsion pathways of Huperzine A and camphor from the binding site of Torpedo California acetylcholinesterase and P450cam proteins, respectively, and comparing them with the previous results obtained using memetic sampling and random acceleration molecular dynamics. In addition, by applying this protocol to a set of ligands bound with LSD1 (lysine specific demethylase 1), we obtained a much higher correlation between the work of pulling force and experimental data on the inhibition constant IC50 compared to that obtained using the unidirectional approach based on minimal steric hindrance.