The 2019 novel coronavirus (SARS-CoV-2) epidemic, which was first reported in December 2019 in Wuhan, China, was declared a pandemic by the World Health Organization in March 2020. Genetically, SARS-CoV-2 is closely related to SARS-CoV, which caused a global epidemic with 8096 confirmed cases in more than 25 countries from 2002 to 2003. Given the significant morbidity and mortality rate, the current pandemic poses a danger to all of humanity, prompting us to understand the activity of SARS-CoV-2 at the atomic level. Experimental studies have revealed that spike proteins of both SARS-CoV-2 and SARS-CoV bind to angiotensin-converting enzyme 2 (ACE2) before entering the cell for replication. However, the binding affinities reported by different groups seem to contradict each other. Wrapp et al. ( Science 2020 , 367 , 1260–1263) showed that the spike protein of SARS-CoV-2 binds to the ACE2 peptidase domain (ACE2-PD) more strongly than does SARS-CoV, and this fact may be associated with a greater severity of the new virus. However, Walls et al. ( Cell 2020 , 181 , 281–292) reported that SARS-CoV-2 exhibits a higher binding affinity, but the difference between the two variants is relatively small. To understand the binding mechnism and experimental results, we investigated how the receptor binding domain (RBD) of SARS-CoV (SARS-CoV-RBD) and SARS-CoV-2 (SARS-CoV-2-RBD) interacts with a human ACE2-PD using molecular modeling. We applied a coarse-grained model to calculate the dissociation constant and found that SARS-CoV-2 displays a 2-fold higher binding affinity. Using steered all-atom molecular dynamics simulations, we demonstrate that, like a coarse-grained simulation, SARS-CoV-2-RBD was associated with ACE2-PD more strongly than was SARS-CoV-RBD, as evidenced by a higher rupture force and larger pulling work. We show that the binding affinity of both viruses to ACE2 is driven by electrostatic interactions.

Despite years of intensive research, little is known about oligomeric structures present during Alzheimer’s disease (AD). Excess of amyloid beta (Aβ) peptides and their aggregation are the basis of the amyloid cascade hypothesis, which attempts to explain the causes of AD. Because of the intrinsically disordered nature of Aβ monomers and the high aggregation rate of oligomers, their structures are almost impossible to resolve using experimental methods. For this reason, we used a physics-based coarse-grained force field to extensively search for the conformational space of the Aβ42 tetramer, which is believed to be the smallest stable Aβ oligomer and the most toxic one. The resulting structures were subsequently optimized, tested for stability, and compared with the proposed experimental fibril models, using molecular dynamics simulations in two popular all-atom force fields. Our results show that the Aβ42 tetramer can form polymorphic stable structures, which may explain different pathways of Aβ aggregation. The models obtained comprise the outer and core chains and, therefore, are significantly different from the structure of mature fibrils. We found that interaction with water is the reason why the tetramer is more compact and less dry inside than fibrils. Physicochemical properties of the proposed all-atom structures are consistent with the available experimental observations and theoretical expectations. Therefore, we provide possible models for further study and design of higher order oligomers.

The outbreak of a new coronavirus SARS-CoV-2 (severe acute respiratory syndrome–coronavirus 2) has caused a global COVID-19 (coronavirus disease 2019) pandemic, resulting in millions of infections and thousands of deaths around the world. There is currently no drug or vaccine for COVID-19, but it has been revealed that some commercially available drugs are promising, at least for treating symptoms. Among them, remdesivir, which can block the activity of RNA-dependent RNA polymerase (RdRp) in old SARS-CoV and MERS-CoV viruses, has been prescribed to COVID-19 patients in many countries. A recent experiment showed that remdesivir binds to SARS-CoV-2 with an inhibition constant of μM, but the exact target has not been reported. In this work, combining molecular docking, steered molecular dynamics, and umbrella sampling, we examined its binding affinity to two targets including the main protease (Mpro), also known as 3C-like protease, and RdRp. We showed that remdesivir binds to Mpro slightly weaker than to RdRp, and the corresponding inhibition constants, consistent with the experiment, fall to the μM range. The binding mechanisms of remdesivir to two targets differ in that the electrostatic interaction is the main force in stabilizing the RdRp–remdesivir complex, while the van der Waals interaction dominates in the Mpro–remdesivir case. Our result indicates that remdesivir can target not only RdRp but also Mpro, which can be invoked to explain why this drug is effective in treating COVID-19. We have identified residues of the target protein that make the most important contribution to binding affinity, and this information is useful for drug development for this disease.

No abstract

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.

The outbreak of a new coronavirus SARS-CoV-2 (severe acute respiratory syndrome–<br>coronavirus 2) has caused a global CoVid-19 (coronavirus disease 2019) pandemic, resulting in millions of infections and thousands of deaths around the world. There is currently no drug or vaccine for CoVid-19, but it has been revealed that some commercially available drugs are promising, at least for treating symptoms. Among them, Remdesivir, which can block the activity of RNA-dependent RNA polymerase (RdRp) in old SARS-CoV and MERS-CoV viruses, has been prescribed to CoVid-19 patients in many countries. A recent experiment showed that Remdesivir binds to SARS-CoV-2 with an inhibition constant of μM, but the exact target has not been reported. In this work, combining molecular docking, steered molecular dynamics and umbrella sampling we examined its binding affinity to two targets including the main protease (Mpro), also known as 3C-like protease, and RdRp. We showed that Remdesivir binds to Mpro slightly weaker than to RdRp and the corresponding inhibition constants, consistent with the experiment, fall to the μM range. The binding mechanisms of<br>Remdesivir to two targets differ in that electrostatic interaction is the main force in stabilizing the RdRp-Remdesivir complex, while the van der Waals interaction dominates in the MproRemdesivir case. Our result indicates that Remdesivir can target not only RdRp but also Mpro, which can be invoked to explain why this drug is effective in treating Covid-19. We have identified residues of the target protein that make the most important contribution to binding affinity, and this information is useful for drug development for this disease. <br>

Amyloid-β (Aβ) peptides form assemblies that are pathological hallmarks of Alzheimer's disease. Aβ oligomers are soluble, mobile, and toxic forms of the peptide that act in the extracellular space before assembling into protofibrils and fibrils. Therefore, oligomers play an important role in the mechanism of Alzheimer's disease. Since it is difficult to determine by experiment the atomic structures of oligomers, which accumulate fast and are polymorphic, computer simulation is a useful tool to investigate elusive oligomers' structures. In this work, we report extended all-atom molecular dynamics simulations, both canonical and replica exchange, of Aβ(1−42) trimer starting from two different initial conformations: (i) the pose produced by the best docking of a monomer aside of a dimer (simulation 1), representing oligomers freshly formed by assembling monomers, and (ii) a configuration extracted from an experimental mature fibril structure (simulation 2), representing settled oligomers in equilibrium with extended fibrils. We showed that in simulation 1, regions with small β-barrels are populated, indicating the chance of spontaneous formation of domains resembling channel-like structures. These structural domains are alternative to those more representative of mature fibrils (simulation 2), the latter showing a stable bundle of C-termini that is not sampled in simulation 1. Moreover, trimer of Aβ(1−42) can form internal pores that are large enough to be accessed by water molecules and Ca 2+ ions.

Divalent cations have a strong impact on the properties of phospholipid membranes, where amyloid-β peptides exert effects related to possible functional or pathological roles. In this work, we use an atomistic computational model of dimyristoyl-phosphatidylcholine (DMPC) membrane bilayers. We perturb this model with a simple model of divalent cations (Mg2+) and with a single amyloid-β (Aβ) peptide of 42 residues, both with and without a single Cu2+ ion bound to the N-terminus. In agreement with the experimental results reported in the literature, the model confirms that divalent cations locally destabilize the DMPC membrane bilayer and, for the first time, that the monomeric form of Aβ helps in avoiding the interactions between divalent cations and DMPC, preventing significant effects on the DMPC bilayer properties. These results are discussed in the frame of a protective role of the diluted Aβ peptide floating around phospholipid membranes.