Atomistic De-novo Inhibitor Generation-Guided Drug Repurposing for SARS-CoV-2 Spike Protein with Free-Energy Validation by Well-Tempered Metadynamics

Unraveling the folding-assisted unbinding mechanism of TCF with its binding partner β-catenin

This study utilizes molecular dynamics simulations aided with multiple walker parallel bias metadynamics to investigate the TCF unbinding mechanism from the β-catenin interface. The results, consistent with experimental binding affinity calculations, unveil a folding-assisted unbinding mechanism.

Insights from in silico study of receptor energetics of SARS-CoV-2 variants

The emergence of new variants of the novel coronavirus SARS-CoV-2 with increased infectivity, superior virulence, high transmissibility, and unmatched immune escape has demonstrated the adaptability and evolutionary fitness of the virus. The subject of relative order of the binding affinity of SARS-CoV-2 variants with the human ACE2 (hACE2) receptor is hotly debated and its resolution has implications for drug design and development. In this work, we have investigated the energetics of the binding of receptor binding domain (RBD) of SARS-CoV-2 variants of concern (VOCs): Beta (B.1.351), Delta (B.1.617.2), Omicron (B.1.1.529), variant of interest (VOI): Kappa (B.1.617.1), and Delta Plus (B.1.617.2.1) variant with the human ACE2 receptor by using the umbrella sampling (US) method. Our work indicates that Delta and Delta Plus variants have greater values of the US binding free energy than Wild-type (WT), whereas Beta, Kappa, and Omicron variants have lower values. Further analysis of hydrogen bonding, salt bridges, non-bonded interaction energy, and contact surface area at the RBD-hACE2 interface establish Delta as the variant with the highest binding affinity among these variants.

Targeting RdRp of SARS-CoV-2 with De Novo Molecule Generation

Viruses are known for their extremely high mutation rates, allowing them to evade both the human immune system and many forms of standard medicine. Despite this, the RNA dependent RNA polymerase (RdRp) of the RNA viruses has been largely conserved, and any significant mutation of this protein is unlikely. The recent COVID-19 pandemic presents a need for therapeutics. We have designed a de novo drug design algorithm that generates strong binding ligands from scratch, based on only the structure of the target protein's receptor. In this paper, we applied our method to target SARS-CoV-2 RdRp and generated several de novo molecules. We then chose some drug molecules based on the structural similarity to some of our strongest binding de novo molecules. Subsequently, we showed, using rigorous all-atom explicit-water free energy calculations in near-microsecond time scales using state-of-the-art well-tempered metadynamics simulations, that some of our de novo generated ligands bind more strongly to RdRp than the recent FDA approved drug remdesivir in its active form, remdesivir triphosphate (RTP). We elucidated the binding mechanism for some of the top binders and compared it with RTP. We believe that this work will be useful both by presenting lead structures for RdRp inhibition and by delivering key insights into the residues of the protein potentially involved in the binding/unbinding of these small molecule drugs, leading to more targeted studies in the future.

Mechanism of the Early Catalytic Events in the Collagenolysis by Matrix Metalloproteinase-1

Metalloproteinase-1 (MMP-1) catalyzed collagen degradation is essential for a wide variety of normal physiological processes, while at the same time contributing to several diseases in humans. Therefore, a comprehensive understanding of this process is of great importance. Although crystallographic and spectroscopic studies provided fundamental information about the structure and function of MMP-1, the precise mechanism of collagen degradation especially considering the complex and flexible structure of the substrate, remains poorly understood. In addition, how the protein environment dynamically reorganizes at the atomic scale into a catalytically active state capable of collagen hydrolysis remains unknown. In this study, we applied experimentally-guided multiscale molecular modeling methods including classical molecular dynamics (MD), well-tempered (WT) classical metadynamics (MetD), combined quantum mechanics/molecular mechanics (QM/MM) MD and QM/MM MetD simulations to explore and characterize the early catalytic events of MMP-1 collagenolysis. Importantly the study provided a complete atomic and dynamic description of the transition from the open to the closed form of the MMP-1•THP complex. Notably, the formation of catalytically active Michaelis complex competent for collagen cleavage was characterized. The study identified the changes in the coordination state of the catalytic zinc(II) associated with the conformational transformation and the formation of catalytically productive ES complex. Our results confirm the essential role of the MMP-1 catalytic domain's α-helices (hA, hB and hC) and the linker region in the transition to the catalytically competent ES complex. Overall, the results provide unique mechanistic insight into the conformational transformations and associated changes in the coordination state of the catalytic zinc(II) that would be important for the design of effective MMP-1 inhibitors.

Structures of the SARS-CoV-2 spike glycoprotein and applications for novel drug development

COVID-19 caused by SARS-CoV-2 has raised a health crisis worldwide. The high morbidity and mortality associated with COVID-19 and the lack of effective drugs or vaccines for SARS-CoV-2 emphasize the urgent need for standard treatment and prophylaxis of COVID-19. The receptor-binding domain (RBD) of the glycosylated spike protein (S protein) is capable of binding to human angiotensin-converting enzyme 2 (hACE2) and initiating membrane fusion and virus entry. Hence, it is rational to inhibit the RBD activity of the S protein by blocking the RBD interaction with hACE2, which makes the glycosylated S protein a potential target for designing and developing antiviral agents. In this study, the molecular features of the S protein of SARS-CoV-2 are highlighted, such as the structures, functions, and interactions of the S protein and ACE2. Additionally, computational tools developed for the treatment of COVID-19 are provided, for example, algorithms, databases, and relevant programs. Finally, recent advances in the novel development of antivirals against the S protein are summarized, including screening of natural products, drug repurposing and rational design. This study is expected to provide novel insights for the efficient discovery of promising drug candidates against the S protein and contribute to the development of broad-spectrum anti-coronavirus drugs to fight against SARS-CoV-2.

Point mutations in SARS-CoV-2 variants induce long-range dynamical perturbations in neutralizing antibodies

Monoclonal antibodies are emerging as a viable treatment for the coronavirus disease 19 (COVID-19). However, newly evolved variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can reduce the efficacy of currently available antibodies and can diminish vaccine-induced immunity. Here, we demonstrate that the microscopic dynamics of neutralizing monoclonal antibodies can be profoundly modified by the mutations present in the spike proteins of the SARS-COV-2 variants currently circulating in the world population. The dynamical perturbations within the antibody structure, which alter the thermodynamics of antigen recognition, are diverse and can depend both on the nature of the antibody and on the spatial location of the spike mutation. The correlation between the motion of the antibody and that of the spike receptor binding domain (RBD) can also be changed, modulating binding affinity. Using protein-graph-connectivity networks, we delineated the mutant-induced modifications in the information-flow along allosteric pathway throughout the antibody. Changes in the collective dynamics were spatially distributed both locally and across long-range distances within the antibody. On the receptor side, we identified an anchor-like structural element that prevents the detachment of the antibodies; individual mutations there can significantly affect the antibody binding propensity. Our study provides insight into how virus neutralization by monoclonal antibodies can be impacted by local mutations in the epitope via a change in dynamics. This realization adds a new layer of sophistication to the efforts for rational design of monoclonal antibodies against new variants of SARS-CoV2, taking the allostery in the antibody into consideration.

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.

Copper(II) import and reduction are dependent on His-Met clusters in the extracellular amino terminus of human copper transporter-1

Copper(I) is an essential metal for all life forms. Though Cu(II) is the most abundant and stable state, its reduction to Cu(I) via an unclear mechanism is prerequisite for its bioutilization. In eukaryotes, the copper transporter-1 (CTR1) is the primary high-affinity copper importer, although its mechanism and role in Cu(II) reduction remain uncharacterized. Here we show that extracellular amino-terminus of human CTR1 contains two methionine-histidine clusters and neighboring aspartates that distinctly bind Cu(I) and Cu(II) preceding its import. We determined that hCTR1 localizes at the basolateral membrane of polarized MDCK-II cells and that its endocytosis to Common-Recycling-Endosomes is regulated by reduction of Cu(II) to Cu(I) and subsequent Cu(I) coordination by the methionine cluster. We demonstrate the transient binding of both Cu(II) and Cu(I) during the reduction process is facilitated by aspartates that also act as another crucial determinant of hCTR1 endocytosis. Mutating the first Methionine cluster (⁷Met-Gly-Met⁹) and Asp¹³ abrogated copper uptake and endocytosis upon copper treatment. This phenotype could be reverted by treating the cells with reduced and nonreoxidizable Cu(I). We show that histidine clusters, on other hand, bind Cu(II) and are crucial for hCTR1 functioning at limiting copper. Finally, we show that two N-terminal His-Met-Asp clusters exhibit functional complementarity, as the second cluster is sufficient to preserve copper-induced CTR1 endocytosis upon complete deletion of the first cluster. We propose a novel and detailed mechanism by which the two His-Met-Asp residues of hCTR1 amino-terminus not only bind copper, but also maintain its reduced state, crucial for intracellular uptake.

Distant residues modulate conformational opening in SARS-CoV-2 spike protein

The novel coronavirus (SARS-CoV-2) pandemic resulted in the largest public health crisis in recent times. Significant drug design effort against SARS-CoV-2 is focused on the receptor-binding domain (RBD) of the spike protein, although this region is highly prone to mutations causing therapeutic resistance. We applied deep data analysis methods on all-atom molecular dynamics simulations to identify key non-RBD residues that play a crucial role in spike−receptor binding and infection. Because the non-RBD residues are typically conserved across multiple coronaviruses, they can be targeted by broad-spectrum antibodies and drugs to treat infections from new strains that might appear during future epidemics.

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.

Atomistic De-novo Inhibitor Generation-Guided Drug Repurposing for SARS-CoV-2 Spike Protein with Free-Energy Validation by Well-Tempered Metadynamics

Computational drug design is increasingly becoming important with new and unforeseen diseases like COVID-19. In this study, we present a new computational de novo drug design and repurposing method and applied it to find plausible drug candidates for the receptor binding domain (RBD) of SARS-CoV-2 (COVID-19). Our study comprises three steps: atom-by-atom generation of new molecules around a receptor, structural similarity mapping to existing approved and investigational drugs, and validation of their binding strengths to the viral spike proteins based on rigorous all-atom, explicit-water well-tempered metadynamics free energy calculations. By choosing the receptor binding domain of the viral spike protein, we showed that some of our new molecules and some of the repurposable drugs have stronger binding to RBD than hACE2. To validate our approach, we also calculated the free energy of hACE2 and RBD, and found it to be in an excellent agreement with experiments. These pool of drugs will allow strategic repurposing against COVID-19 for a particular prevailing conditions.