Multitarget Potential Drug Candidates for High-Grade Gliomas Identified by Multiple Reaction Monitoring Coupled with In Silico Drug Repurposing.

High-grade gliomas (HGGs) are extremely aggressive primary brain tumors with high mortality rates. Despite notable progress achieved by clinical research and biomarkers emerging from proteomics studies, efficacious drugs and therapeutic targets are limited. This study used targeted proteomics, in silico molecular docking, and simulation-based drug repurposing to identify potential drug candidates for HGGs. Importantly, we performed multiple reaction monitoring (MRM) on differentially expressed proteins with putative roles in the development and progression of HGGs based on our previous work and the published literature. Furthermore, in silico molecular docking-based drug repurposing was performed with a customized library of FDA-approved drugs to identify multitarget-directed ligands. The top drug candidates such as Pazopanib, Icotinib, Entrectinib, Regorafenib, and Cabozantinib were explored for their drug-likeness properties using the SwissADME. Pazopanib exhibited binding affinities with a maximum number of proteins and was considered for molecular dynamic simulations and cell toxicity assays. HGG cell lines showed enhanced cytotoxicity and cell proliferation inhibition with Pazopanib and Temozolomide combinatorial treatment compared to Temozolomide alone. To the best of our knowledge, this is the first study combining MRM with molecular docking and simulation-based drug repurposing to identify potential drug candidates for HGG. While the present study identified five multitarget-directed potential drug candidates, future clinical studies in larger cohorts are crucial to evaluate the efficacy of these molecular candidates. The research strategy and methodology used in the present study offer new avenues for innovation in drug discovery and development which may prove useful, particularly for cancers with low cure rates.

Mutations in V84I & A184V of blaTEM cluster plays a pivotal role in the dynamics of Ω-loop leading to genesis of IR-TEM.

The Enterobacteriaceae family exhibits resistance to antibiotics by producing ß-Lactamase. Mutations in ß-Lactamase, have led to the generation of inhibitor resistant variants known as IR-TEM. In the present study, IR-TEM ß-Lactamase of Enterobacter hormaechei and Enterobacter asburiae was compared. To investigate the mechanism behind the conferred mutation, we studied its interaction with Clavulanic acid, (ß-Lactamase inhibitor) with different lineages of TEM and IR-TEM. We found that Clavulanic acid quickly left the binding pockets of both variants using molecular dynamics (MD) simulations. Interestingly, mutations at the V84I and A184V positions were found to drastically influence the protein dynamics. Mutating the residues at V84I and A184V positions by computational mutagenesis in Enterobacter hormaechei, it was observed that the residues on the Ω-loop as well as a few downstream residues were primarily involved in generating resistance towards inhibitors by conferring increased flexibility to the loop. This further strongly supports the notion that inhibitor resistance in ß-Lactamase is conferred through allosteric regulation, wherein mutations in positions 84 and 184 may play an important role in regulating the movement of the Ω-loop. These two positions determine the lineage pattern between two clusters in TEM-1 and TEM-116. To date, no reports have been made regarding the importance of these mutations and their dynamics in Ω-loop.Communicated by Ramaswamy H. Sarma.

Delineating the Structure-Dynamics-Binding Differences among BA.1, BA.4/5, and BF.7 SARS-CoV-2 Variants through Atomistic Simulations: Correlation with Structural and Epidemiological Features.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an RNA virus possessing a spike (S) protein that facilitates the entry of the virus into human cells. The emergence of highly transmissible and fit SARS-CoV-2 variants has been driven by the positive selection of mutations within the S-protein. Notable among these variants are alpha, beta, gamma, delta, and omicron (BA.1), with the latter contributing to significant global health challenges and impacting populations worldwide. Recently, a novel subvariant of BA.1, named BF.7, has surfaced, purportedly exhibiting elevated transmissibility and infectivity rates. In order to comprehend and compare the transmissibility and disease progression characteristics of distinct SARS-CoV-2 variants, we performed an extensive comparative analysis utilizing all-atom molecular dynamics (MD) simulations (in triplicate) to investigate the structural, dynamic, and binding features of BA.1, BA.4/5, and BF.7. Our simulation findings, energetic analysis, and assessment of physicochemical properties collectively illuminate the dominance of the BA.1 variant over the others, a trend that is further substantiated by the sustained global prevalence of BA.1 relative to BA.4/5 and BF.7. Additionally, our simulation results align well with the reported cryoelectron microscopy (cryo-EM) structural data and epidemiological characteristics obtained from the Global Initiative on Sharing All Influenza Data (GISAID). This study presents a comprehensive comparative elucidation of the critical structural, dynamic, and binding attributes of these variants, providing insights into the predominance of BA.1 and its propensity to continuously generate numerous novel subvariants.

Identification of inhibitors against SARS-CoV-2 variants of concern using virtual screening and metadynamics-based enhanced sampling.

Among the variants of SARS-CoV-2, some are more infectious than the Wild-type. Interestingly, these mutations enable the virus to evade the therapeutic efforts. Hence, there is a need for candidate drug molecules that can potently bind with all the variants. We have adopted a strategy combining virtual screening, molecular docking followed by rigorous sampling by metadynamics simulations to find candidate molecules. From our results we found four highly potent drug candidates that can bind to the Spike-RBD of all the variants of the virus. Additionally, we also found that certain signature residues on the RBM region commonly bind to each of these inhibitors. Thus, our study not only gives information on the chemical compounds, but also residues on the proteins which could be targeted for future drug and vaccine development studies.

MoS2 nanosheets effectively bind to the receptor binding domain of the SARS-CoV-2 spike protein and destabilize the spike-human ACE2 receptor interactions.

The use of nanotechnology is becoming increasingly significant as a tool that can provide a range of options for the identification, inactivation, and therapy of coronavirus disease 2019 (COVID-19). The potential of nanoparticles as an alternative therapeutic agent to inactivate SARS-CoV-2 is continually being investigated. Herein, we have explored the interaction of 2D molybdenum disulfide (MoS2) nanosheets with the SARS-CoV-2 spike protein, human ACE2 receptor and the complex formed between them through molecular docking and atomistic simulations. The results indicated that MoS2 nanosheets can effectively bind to the receptor binding domain (RBD) of the spike protein with good docking energies. It is interesting to note that this also applied to the extensively glycosylated spike protein and its variations, Kappa and Delta. A significant loss of secondary structures was observed when MoS2 nanosheets interacted with the RBD of the spike protein. The nanosheets interacted strongly with the proteins through a number of hydrogen bonds and van der Waals interactions. Moreover, the binding of the MoS2 nanosheets at different locations of the RBD or ACE2 in the spike-RBD/ACE2 complex resulted in significant conformational changes. Detailed energetics and solvent accessibility calculations revealed that, when present at the interface, MoS2 nanosheets can be a potential inhibitory agent. The findings were supported by de-wetting calculations, indicating strong adherence of the RBD of spike protein on the MoS2 nanosheet and de-stability of the spike-ACE2 interaction. Thus, the findings clearly demonstrate the antiviral potential of 2D MoS2 nanosheets, prompting its further exploration for combating COVID-19.

How Does Temperature Affect the Dynamics of SARS-CoV-2 M Proteins? Insights from Molecular Dynamics Simulations.

Enveloped viruses, in general, have several transmembrane proteins and glycoproteins, which assist the virus in entry and attachment onto the host cells. These proteins also play a significant role in determining the shape and size of the newly formed virus particles. The lipid membrane and the embedded proteins affect each other in non-trivial ways during the course of the viral life cycle. Unraveling the nature of the protein-protein and protein-lipid interactions, under various environmental and physiological conditions, could therefore prove to be crucial in development of therapeutics. Here, we study the M protein of SARS-CoV-2 to understand the effect of temperature on the properties of the protein-membrane system. The membrane-embedded dimeric M proteins were studied using atomistic and coarse-grained molecular dynamics simulations at temperatures ranging between 10 and 50 °C. While temperature-induced fluctuations are expected to be monotonic, we observe a steady rise in the protein dynamics up to 40 °C, beyond which it surprisingly reverts back to the low-temperature behavior. Detailed investigation reveals disordering of the membrane lipids in the presence of the protein, which induces additional curvature around the transmembrane region. Coarse-grained simulations indicate temperature-dependent aggregation of M protein dimers. Our study clearly indicates that the dynamics of membrane lipids and integral M protein of SARS-CoV-2 enables it to better associate and aggregate only at a certain temperature range (i.e., ~ 30-40 °C). This can have important implications in the protein aggregation and subsequent viral budding/fission processes.

Molecular insights into the differential dynamics of SARS-CoV-2 variants of concern.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has affected the lives and livelihood of millions of individuals around the world. It has mutated several times after its first inception, with an estimated two mutations occurring every month. Although we have been successful in developing vaccines against the virus, the emergence of variants has enabled it to escape therapy. Few of the generated variants are also reported to be more infectious than the wild-type (WT). In this study, we analyze the attributes of all RBD/ACE2 complexes for the reported VOCs, namely, Alpha, Beta, Gamma, and Delta through computer simulations. Results indicate differences in orientation and binding energies of the VOCs from the WT. Overall, it was observed that electrostatic interactions play a major role in the binding of the complexes. Detailed residue level energetics revealed that the most prominent changes in interaction energies were seen particularly at the mutated residues which were present at RBD/ACE2 interface. We found that the Delta variant is one of the most tightly bound variants of SARS-CoV-2 with dynamics similar to WT. The high binding affinity of RBD towards ACE2 is indicative of an increase in viral transmission and infectivity. The details presented in our study provide additional information for the design and development of effective therapeutic strategies for the emerging variants of the virus in the future.

Scanning the RBD-ACE2 molecular interactions in Omicron variant.

The emergence of new SARS-CoV-2 variants poses a threat to the human population where it is difficult to assess the severity of a particular variant of the virus. Spike protein and specifically its receptor binding domain (RBD) which makes direct interaction with the ACE2 receptor of the human has shown prominent amino acid substitutions in most of the Variants of Concern. Here, by using all-atom molecular dynamics simulations we compare the interaction of Wild-type RBD/ACE2 receptor complex with that of the latest Omicron variant of the virus. We observed a very interesting diversification of the charge, dynamics and energetics of the protein complex formed upon mutations. These results would help us in understanding the molecular basis of binding of the Omicron variant with that of SARS-CoV-2 Wild-type.

Computational studies reveal Fluorine based quinolines to be potent inhibitors for proteins involved in SARS-CoV-2 assembly.

World is witnessing one of the worst pandemics of this century caused by SARS-CoV-2 virus which has affected millions of individuals. Despite rapid efforts to develop vaccines and drugs for COVID-19, the disease is still not under control. Chloroquine (CQ) and Hydroxychloroquine (HCQ) are two very promising inhibitors which have shown positive effect in combating the disease in preliminary experimental studies, but their use was reduced due to severe side-effects. Here, we performed a theoretical investigation of the same by studying the binding of the molecules with SARS-COV-2 Spike protein, the complex formed by Spike and ACE2 human receptor and a human serine protease TMPRSS2 which aids in cleavage of the Spike protein to initiate the viral activation in the body. Both the molecules had shown very good docking energies in the range of -6kcal/mol. Subsequently, we did a high throughput screening for other potential quinoline candidates which could be used as inhibitors. From the large pool of ligand candidates, we shortlisted the top three ligands (binding energy -8kcal/mol). We tested the stability of the docked complexes by running Molecular Dynamics (MD) simulations where we observed the stability of the quinoline analogues with the Spike-ACE2 and TMPRSS2 nevertheless the quinolines were not stable with the Spike protein alone. Thus, although the inhibitors bond very well with the protein molecules their intrinsic binding affinity depends on the protein dynamics. Moreover, the quinolines were stable when bound to electronegative pockets of Spike-ACE2 or TMPRSS2 but not with Viral Spike protein. We also observed that a Fluoride based compound: 3-[3-(Trifluoromethyl)phenyl]quinoline helps the inhibitor to bind with both Spike-ACE2 and TMPRSS2 with equal probability. The molecular details presented in this study would be very useful for developing quinoline based drugs for COVID-19 treatment.

Accelerating COVID-19 Research Using Molecular Dynamics Simulation.

The COVID-19 pandemic has emerged as a global medico-socio-economic disaster. Given the lack of effective therapeutics against SARS-CoV-2, scientists are racing to disseminate suggestions for rapidly deployable therapeutic options, including drug repurposing and repositioning strategies. Molecular dynamics (MD) simulations have provided the opportunity to make rational scientific breakthroughs in a time of crisis. Advancements in these technologies in recent years have become an indispensable tool for scientists studying protein structure, function, dynamics, interactions, and drug discovery. Integrating the structural data obtained from high-resolution methods with MD simulations has helped in comprehending the process of infection and pathogenesis, as well as the SARS-CoV-2 maturation in host cells, in a short duration of time. It has also guided us to identify and prioritize drug targets and new chemical entities, and to repurpose drugs. Here, we discuss how MD simulation has been explored by the scientific community to accelerate and guide translational research on SARS-CoV-2 in the past year. We have also considered future research directions for researchers, where MD simulations can help fill the existing gaps in COVID-19 research.

Investigation of the Effect of Temperature on the Structure of SARS-CoV-2 Spike Protein by Molecular Dynamics Simulations.

Statistical and epidemiological data imply temperature sensitivity of the SARS-CoV-2 coronavirus. However, the molecular level understanding of the virus structure at different temperature is still not clear. Spike protein is the outermost structural protein of the SARS-CoV-2 virus which interacts with the Angiotensin Converting Enzyme 2 (ACE2), a human receptor, and enters the respiratory system. In this study, we performed an all atom molecular dynamics simulation to study the effect of temperature on the structure of the Spike protein. After 200 ns of simulation at different temperatures, we came across some interesting phenomena exhibited by the protein. We found that the solvent exposed domain of Spike protein, namely S1, is more mobile than the transmembrane domain, S2. Structural studies implied the presence of several charged residues on the surface of N-terminal Domain of S1 which are optimally oriented at 10-30°C. Bioinformatics analyses indicated that it is capable of binding to other human receptors and should not be disregarded. Additionally, we found that receptor binding motif (RBM), present on the receptor binding domain (RBD) of S1, begins to close around temperature of 40°C and attains a completely closed conformation at 50°C. We also found that the presence of glycan moieties did not influence the observed protein dynamics. Nevertheless, the closed conformation disables its ability to bind to ACE2, due to the burying of its receptor binding residues. Our results clearly show that there are active and inactive states of the protein at different temperatures. This would not only prove beneficial for understanding the fundamental nature of the virus, but would be also useful in the development of vaccines and therapeutics.

Identification of Factors Promoting HBV Capsid Self-Assembly by Assembly-Promoting Antivirals.

Around 270 million individuals currently live with hepatitis B virus (HBV) infection. Heteroaryldihydropyrimidines (HAPs) are a family of antivirals that target the HBV capsid protein and induce aberrant self-assembly. The capsids formed resemble the native capsid structure but are unable to propagate the virus progeny because of a lack of RNA/DNA. Under normal conditions, self-assembly is initiated by the viral genome. The mode of action of HAPs, however, remains largely unknown. In this work, using molecular dynamics simulations, we attempted to understand the action of HAP by comparing the dynamics of capsid proteins with and without HAPs. We found that the inhibitor is more stable in higher oligomers. It retains its stability in the hexamer throughout 1 µs of simulation. Our results also show that the inhibitor might help in stabilizing the C-terminus, the HBc 149-183 arginine-rich domain of the capsid protein. The C-termini of dimers interact with each other, assisted by the HAP inhibitor. During capsid assembly, the termini are supposed to directly interact with the viral genome, thereby suggesting that the viral genome might work in a similar way to stabilize the capsid protein. Our results may help in understanding the underlying molecular mechanism of HBV capsid self-assembly, which should be crucial for exploring new drug targets and structure-based drug design.

Mechanism of p27 Unfolding for CDK2 Reactivation.

Cell-cycle regulatory protein, CDK2 is active when bound to its complementary partner protein, CyclinA or E. Recent discovery of the Kip/Cip family of proteins has indicated that the activity of CDK2 is also regulated by a member protein, p27. Although, the mechanism of CDK2 inhibition by p27 binding is known from crystal structure, little is known about the mechanism of CDK2 reactivation. Here, we execute classical and accelerated molecular dynamics simulations of unphosphorylated- and phosphorylated-p27 bound CDK2/CyclinA to unravel the CDK2 reactivation mechanism at molecular-to-atomic detail. Results suggest that the phosphorylation of p27 Y88 residue (pY88-p27) first disrupts the p27/CDK2 hybrid ß-sheet and subsequently ejects the p27 310 helix from CDK2 catalytic cleft. The unbinding of p27 from CDK2/CyclinA complex, thus, follows a two-step unfolding mechanism, where the 310 helix ejection constitutes the rate-limiting step. Interestingly, the unfolding of p27 leaves CDK2/CyclinA in an active state, where the prerequisite CDK2-CyclinA interfacial contacts were regained and ATP achieved its native position for smooth transfer of phosphate. Our findings match very well with NMR chemical shift data that indicated the flip-out of p27 310 helix from CDK2 pocket and kinetic experiments that exhibited significant kinase activity of CDK2 when saturated with pY88-p27.

Why are the truncated cyclin Es more effective CDK2 activators than the full-length isoforms?

Cell cycle regulating enzymes, CDKs, become activated upon association with their regulatory proteins, cyclins. The G1 cyclin, cyclin E, is overexpressed and present in low molecular weight (LMW) isoforms in breast cancer cells and tumor tissues. In vivo and in vitro studies have shown that these LMW isoforms of cyclin E hyperactivate CDK2 and accelerate the G1-S phase of cell division. The molecular basis of CDK2 hyperactivation due to LMW cyclin E isoforms in cancer cells is, however, unknown. Here, we employ a computational approach, combining homology modeling, bioinformatics analyses, molecular dynamics (MD) simulations, and principal component analyses to unravel the key structural features of CDK2-bound full-length and LMW isoforms of cyclin E1 and correlate those features to their differential activity. Results suggest that the missing N- and C-terminal regions of the cyclin E LMW isoforms constitute the Nuclear Localization Sequence (NLS) and PEST domains and are intrinsically disordered. These regions, when present in the full-length cyclin E/CDK2 complex, weaken the cyclin-CDK interface packing due to the loss of a large number of key interface interactions. Such weakening is manifested in the decreased contact area and increased solvent accessibility at the interface and also by the absence of concerted motions between the two partner proteins in the full-length complex. More effective packing and interactions between CDK2 and LMW cyclin E isoforms, however, produce more efficient protein-protein complexes that accelerate the cell division processes in cancer cells, where these cyclin E isoforms are overexpressed.

Molecular basis of differential selectivity of cyclobutyl-substituted imidazole inhibitors against CDKs: insights for rational drug design.

Cyclin-dependent kinases (CDKs) belong to the CMGC subfamily of protein kinases and play crucial roles in eukaryotic cell division cycle. At least seven different CDKs have been reported to be implicated in the cell cycle regulation in vertebrates. These CDKs are highly homologous and contain a conserved catalytic core. This makes the design of inhibitors specific for a particular CDK difficult. There is, however, growing need for CDK5 specific inhibitors to treat various neurodegenerative diseases. Recently, cis-substituted cyclobutyl-4-aminoimidazole inhibitors have been identified as potent CDK5 inhibitors that gave up to 30-fold selectivity over CDK2. Available IC50 values also indicate a higher potency of this class of inhibitors over commercially available drugs, such as roscovitine. To understand the molecular basis of higher potency and selectivity of these inhibitors, here, we present molecular dynamics simulation results of CDK5/p25 and CDK2/CyclinE complexed with a series of cyclobutyl-substituted imidazole inhibitors and roscovitine. The atomic details of the stereospecificity and selectivity of these inhibitors are obtained from energetics and binding characteristics to the CDK binding pocket. The study not only complements the experimental findings, but also provides a wealth of detailed information that could help the structure-based drug designing processes.