ABSTRACT
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.
Subject(s)
Molecular Dynamics Simulation , Protein Biosynthesis , RNA, Messenger , Ribosome Subunits, Small, Eukaryotic , SARS-CoV-2 , Viral Nonstructural Proteins , Viral Nonstructural Proteins/metabolism , Viral Nonstructural Proteins/chemistry , Viral Nonstructural Proteins/genetics , RNA, Messenger/metabolism , RNA, Messenger/genetics , RNA, Messenger/chemistry , Ribosome Subunits, Small, Eukaryotic/metabolism , Ribosome Subunits, Small, Eukaryotic/chemistry , SARS-CoV-2/metabolism , SARS-CoV-2/chemistry , SARS-CoV-2/genetics , Protein Binding , Humans , Static ElectricityABSTRACT
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.
Subject(s)
Molecular Dynamics Simulation , Neuropilin-1 , Protein Binding , SARS-CoV-2 , Spike Glycoprotein, Coronavirus , Neuropilin-1/metabolism , Neuropilin-1/chemistry , Humans , SARS-CoV-2/chemistry , SARS-CoV-2/metabolism , Spike Glycoprotein, Coronavirus/chemistry , Spike Glycoprotein, Coronavirus/metabolism , COVID-19/virology , COVID-19/metabolism , Angiotensin-Converting Enzyme 2/metabolism , Angiotensin-Converting Enzyme 2/chemistryABSTRACT
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.
Subject(s)
Amyloid beta-Peptides , Disulfides , Intrinsically Disordered Proteins , Molecular Dynamics Simulation , Amyloid beta-Peptides/chemistry , Amyloid beta-Peptides/metabolism , Disulfides/chemistry , Intrinsically Disordered Proteins/chemistry , Humans , Protein Structure, Secondary , Peptide Fragments/chemistry , Protein ConformationABSTRACT
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.
Subject(s)
Lipid Bilayers , Molecular Dynamics Simulation , Serotonin , Lipid Bilayers/chemistry , Lipid Bilayers/metabolism , Serotonin/chemistry , Serotonin/metabolism , Membrane Fusion , Synaptic Vesicles/metabolism , Synaptic Vesicles/chemistry , Spectrometry, Fluorescence , Hydrophobic and Hydrophilic InteractionsABSTRACT
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.
Subject(s)
Amyloid beta-Peptides , Molecular Dynamics Simulation , Solvents , alpha-Synuclein , Humans , alpha-Synuclein/chemistry , alpha-Synuclein/metabolism , Amyloid beta-Peptides/chemistry , Amyloid beta-Peptides/metabolism , Peptide Fragments/chemistry , Peptide Fragments/metabolism , Protein Multimerization , Solvents/chemistryABSTRACT
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.
Subject(s)
Angiotensin-Converting Enzyme 2 , Molecular Dynamics Simulation , SARS-CoV-2 , Humans , COVID-19ABSTRACT
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.
Subject(s)
COVID-19 , SARS-CoV-2 , Spike Glycoprotein, Coronavirus , Humans , Angiotensin-Converting Enzyme 2 , Mutation , Protein BindingABSTRACT
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.
Subject(s)
Cell Membrane , Escherichia coli , Exoribonucleases , Protein Multimerization , Silent Mutation , Escherichia coli/enzymology , Exoribonucleases/chemistry , Exoribonucleases/genetics , Kinetics , Protein Folding , Protein Multimerization/genetics , Cell Membrane/enzymologyABSTRACT
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.
ABSTRACT
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.
Subject(s)
Lecithins , Lipid Bilayers , Cations, Divalent , Membranes , WaterABSTRACT
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.
Subject(s)
COVID-19 , SARS-CoV-2 , Humans , Pandemics , Spike Glycoprotein, Coronavirus , AntibodiesABSTRACT
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.
Subject(s)
Molecular Chaperones , Protein Folding , Molecular Chaperones/metabolismABSTRACT
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.
Subject(s)
Protein Biosynthesis , Protein Folding , Ribosomes/metabolism , Proteins/metabolism , Molecular Dynamics SimulationABSTRACT
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.
Subject(s)
Phosphatidylethanolamines , Serotonin , Phosphatidylethanolamines/chemistry , Lipid Bilayers/chemistry , Phosphatidylcholines/chemistry , Phosphatidylserines/chemistry , Cholesterol/chemistry , Phospholipids/chemistryABSTRACT
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.
ABSTRACT
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.
Subject(s)
Molecular Dynamics Simulation , Neoplasms , Cell Membrane , Membranes , WaterABSTRACT
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.
Subject(s)
Alzheimer Disease , Protein Aggregates , Alzheimer Disease/metabolism , Amyloid/chemistry , Amyloid beta-Peptides/chemistry , Humans , Mechanical Phenomena , Molecular Dynamics Simulation , Peptide Fragments/chemistry , Protein StabilityABSTRACT
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.
Subject(s)
COVID-19 , Single-Domain Antibodies , Antibodies, Monoclonal , Antibodies, Neutralizing , Antibodies, Viral , COVID-19/prevention & control , Humans , SARS-CoV-2 , Spike Glycoprotein, CoronavirusABSTRACT
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.
Subject(s)
COVID-19 , SARS-CoV-2 , Humans , Molecular Dynamics Simulation , Pandemics , Peptidyl-Dipeptidase A/chemistry , Protein Binding , Spike Glycoprotein, Coronavirus/chemistryABSTRACT
The amyloid cascade hypothesis states that senile plaques, composed of amyloid ß (Aß) fibrils, play a key role in Alzheimer's disease (AD). However, recent experiments have shown that Aß oligomers are more toxic to neurons than highly ordered fibrils. The molecular mechanism underlying this observation remains largely unknown. One of the possible scenarios for neurotoxicity is that Aß peptides create pores in the lipid membrane that allow Ca2+ ions to enter cells, resulting in a signal of cell apoptosis. Hence, one might think that oligomers are more toxic due to their higher ability to create ion channels than fibrils. In this work, we study the effect of Aß42 dodecamer and fibrils on a neuronal membrane, which is similar to that observed in AD patients, using all-atom molecular dynamics simulations. Due to short simulation times, we cannot observe the formation of pores, but useful insight on the early events of this process has been obtained. Namely, we showed that dodecamer distorts the lipid membrane to a greater extent than fibrils, which may indicate that ion channels can be more easily formed in the presence of oligomers. Based on this result, we anticipate that oligomers are more toxic than mature fibrils, as observed experimentally. Moreover, the Aß-membrane interaction was found to be governed by the repulsive electrostatic interaction between Aß and the ganglioside GM1 lipid. We calculated the bending and compressibility modulus of the membrane in the absence of Aß and obtained good agreement with the experiment. We predict that the dodecamer will increase the compressibility modulus but has little effect on the bending modulus. Due to the weak interaction with the membrane, fibrils insignificantly change the membrane elastic properties.