Search for: All records

Mutations in the new variants of SARS-CoV-2 spike protein modulates the dynamics of the neutralizing antibodies. Capturing such modulations from MD simulations and graph network model identifies the role of mutations in facilitating immune evasion. 

Infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) involves the attachment of the receptor-binding domain (RBD) of its spike proteins to the ACE2 receptors on the peripheral membrane of host cells. Binding is initiated by a down-to-up conformational change in the spike protein, the change that presents the RBD to the receptor. To date, computational and experimental studies that search for therapeutics have concentrated, for good reason, on the RBD. However, the RBD region is highly prone to mutations, and is therefore a hotspot for drug resistance. In contrast, we here focus on the correlations between the RBD and residues distant to it in the spike protein. This allows for a deeper understanding of the underlying molecular recognition events and prediction of the highest-effect key mutations in distant, allosteric sites, with implications for therapeutics. Also, these sites can appear in emerging mutants with possibly higher transmissibility and virulence, and preidentifying them can give clues for designing pan-coronavirus vaccines against future outbreaks. Our model, based on time-lagged independent component analysis (tICA) and protein graph connectivity network, is able to identify multiple residues that exhibit long-distance coupling with the RBD opening. Residues involved in the most ubiquitous D614G mutation and the A570D mutation of the highly contagious UK SARS-CoV-2 variant are predicted ab initio from our model. Conversely, broad-spectrum therapeutics like drugs and monoclonal antibodies can target these key distant-but-conserved regions of the spike protein. 

The molecular features that dictate interactions between functionalized nanoparticles and biomolecules are not well understood. This is in part because for highly charged nanoparticles in solution, establishing a clear connection between the molecular features of surface ligands and common experimental observables such as ζ potential requires going beyond the classical models based on continuum and mean field models. Motivated by these considerations, molecular dynamics simulations are used to probe the electrostatic properties of functionalized gold nanoparticles and their interaction with a charged peptide in salt solutions. Counterions are observed to screen the bare ligand charge to a significant degree even at the moderate salt concentration of 50 mM. As a result, the apparent charge density and ζ potential are largely insensitive to the bare ligand charge densities, which fall in the range of ligand densities typically measured experimentally for gold nanoparticles. While this screening effect was predicted by classical models such as the Manning condensation theory, the magnitudes of the apparent surface charge from microscopic simulations and mean-field models are significantly different. Moreover, our simulations found that the chemical features of the surface ligand ( e.g. , primary vs. quaternary amines, heterogeneous ligand lengths) modulate the interfacial ion and water distributions and therefore the interfacial potential. The importance of interfacial water is further highlighted by the observation that introducing a fraction of hydrophobic ligands enhances the strength of electrostatic binding of the charged peptide. Finally, the simulations highlight that the electric double layer is perturbed upon binding interactions. As a result, it is the bare charge density rather than the apparent charge density or ζ potential that better correlates with binding affinity of the nanoparticle to a charged peptide. Overall, our study highlights the importance of molecular features of the nanoparticle/water interface and underscores a set of design rules for the modulation of electrostatic driven interactions at nano/bio interfaces.