Evaluating binding and interaction of selected pesticides with serum albumin proteins by Raman, 1H NMR, mass spectrometry and molecular dynamics simulation.

Addressing the acute pesticide poisoning and toxicity to humans, is a global challenge of top priority. Serum albumin is the most abundant plasma protein, capable of binding with herbicide and pesticide residues. This study reports multifaceted approaches for in-depth and robust investigation of the molecular interactions of selected pesticides, including propanil (PPL), bromoxynil (BXL), metolachlor (MLR) and glyphosate (GPE) with bovine serum albumin (BSA) proteins using experimental (Raman and FTIR spectroscopy, native mass spectrometry and high field 1H NMR), molecular dynamics (MD) simulation and principal component analysis (PCA). The binding of pesticides with BSA resulted in BSA amide I and amide II Raman spectral shifts. PCA of Raman spectra of serum-pesticide complexes showed the grouping of pesticides on the score plot based on the similarities and differences in pesticides' chemical structures. Native mass spectrometry results revealed strong adduct formation of the pesticides with the protein. The observed changes in chemical shifts, peak broadening or peak disappearance of characteristic proton signals of the pesticides, indicated altered chemical environments due to binding BSA-pesticides interactions. The results of MD simulation conducted for over 500 ns revealed strong pesticides interaction with LEU197, LEU218, LEU237, TRP213, SER286 and ILE289 residues to the site I of BSA. Free energy landscapes provided insights into the conformational changes in BSA on the binding of pesticides. Overall, the experimental and computational results are in consonant and indicate the binding of pesticides into the site I and site II (sub-domain IIA) of the BSA via hydrogen bonding, non-covalent and hydrophobic interactions.Communicated by Ramaswamy H. Sarma.

Structural dynamics and functional analysis of Saprolegnia parasitica chitin synthases 5 in a phospholipid bilayer.

Saprolegnia parasitica is an oomycete responsible for a fish disease called saprolegniosis, which poses an economic and environmental burden on aquaculture production. In Saprolegnia, CHS5 of S. parasitica (SpCHS5) contains an N-terminal domain, a catalytic domain of the glycosyltransferase -2 family containing a GT-A fold, and a C-terminal transmembrane domain. No three-dimensional structure of SpCHS5 is reported yet disclosing the structural details of this protein. We have developed a structural model of full-length SpCHS5 and validated it by molecular dynamics simulation technique. From the 1 microsecond simulations, we retrieved the stable RoseTTAFold model SpCHS5 protein to explain characteristics and structural features. Furthermore, from the analysis of the movement of chitin in the protein cavity, we assumed that ARG 482, GLN 527, PHE 529, PHE 530, LEU 540, SER 541, TYR 544, ASN 634, THR 641, TYR 645, THR 641, ASN 772 residues as a main cavity lining site. In SMD analysis, we investigated the opening of the transmembrane cavity required for chitin translocation. The pulling of chitin from the internal cavity to the extracellular region was observed through steered molecular dynamics simulations. A comparison of the initial and final structures of chitin complex showed that there's a transmembrane cavity opening in the simulations. Overall, this present work will help us understand the structural and functional basis of CHS5 and design inhibitors against SpCHS5.Communicated by Ramaswamy H. Sarma.

Antiviral peptides inhibiting the main protease of SARS-CoV-2 investigated by computational screening and in vitro protease assay.

The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plays an important role in viral replication and transcription and received great attention as a vital target for drug/peptide development. Therapeutic agents such as small-molecule drugs or peptides that interact with the Cys-His present in the catalytic site of Mpro are an efficient way to inhibit the protease. Although several emergency-approved vaccines showed good efficacy and drastically dropped the infection rate, evolving variants are still infecting and killing millions of people globally. While a small-molecule drug (Paxlovid) received emergency approval, small-molecule drugs have low target specificity and higher toxicity. Besides small-molecule drugs, peptide therapeutics are thus gaining increasing popularity as they are easy to synthesize and highly selective and have limited side effects. In this study, we investigated the therapeutic value of 67 peptides targeting Mpro using molecular docking. Subsequently, molecular dynamics (MD) simulations were implemented on eight protein-peptide complexes to obtain molecular-level information on the interaction between these peptides and the Mpro active site, which revealed that temporin L, indolicidin, and lymphocytic choriomeningitis virus (LCMV) GP1 are the best candidates in terms of stability, interaction, and structural compactness. These peptides were synthesized using the solid-phase peptide synthesis protocol, purified by reversed-phase high-performance liquid chromatography (RP-HPLC), and authenticated by mass spectrometry (MS). The in vitro fluorometric Mpro activity assay was used to validate the computational results, where temporin L and indolicidin were observed to be very active against SARS-CoV-2 Mpro with IC50 values of 38.80 and 87.23 µM, respectively. A liquid chromatography-MS (LC-MS) assay was developed, and the IC50 value of temporin L was measured at 23.8 µM. The solution-state nuclear magnetic resonance (NMR) structure of temporin L was determined in the absence of sodium dodecyl sulfate (SDS) micelles and was compared to previous temporin structures. This combined investigation provides critical insights and assists us to further develop peptide inhibitors of SARS-CoV-2 Mpro through structural guided investigation.

Elucidating the Structure, Dynamics, and Interaction of a Choline Chloride and Citric Acid Based Eutectic System by Spectroscopic and Molecular Modeling Investigations.

Eutectic solvent systems are versatile solvents that have found widespread use in numerous applications. Traditional solvents are homogeneous, having only one component, and their chemistry is relatively simple, with some exceptions. On the other hand, deep eutectic solvents (DESs) comprise binary components, generally a donor and an acceptor in hydrogen bonding with varying ratios. The interaction chemistry among the donor and acceptor involved in hydrogen bonding in DESs is complicated. Although numerous research is focused on the synthesis and application of DESs, few studies are reported to elucidate the complex structure and dynamic and interaction behavior of DESs. In this study, we employed calorimetry, vibrational spectroscopy techniques including FTIR and Raman, and nuclear magnetic resonance to derive insight into the structural feature and noncovalent contact of choline chloride (ChCl) and citric acid (CA) while they formed DESs. The 1:1 ChCl/CA eutectic system showed phase transitions and melting peaks with the most pronounced peak at 156.22 °C, suggesting the DESs melting at a lower temperature than the melting temperatures of ChCl and CA. In addition to IR and Raman findings, 1H NMR investigations demonstrate hydrogen bonding intermolecular interactions between ChCl and CA, supporting the formation of 1:1 ChCl/CA DESs based on the deshielded chemical shifts of the proton for Ch. The interaction of the chloride anion with the methyl protons (H4) and methylene protons (H3) of ChCl as well as the strong hydrogen bonding interactions between the hydroxyl hydrogen (H1) of ChCl with one of CA's carbonyl oxygens both supported the formation of conformer E. In addition, molecular dynamics followed by the density functional theory (DFT) was employed to visualize the structure and interaction of DESs using the ωB97XD theory and 6-311++G (d,p) basis set. Both experimental and theoretical IR, Raman, and structural analyses provided evidence of the formation of DESs by possessing hydrogen bonds. These multifaceted experimental and computational investigations provide details of structural and intermolecular interactions of ChCl/CA DESs.

Therapeutic Potential of Antiviral Peptides against the NS2B/NS3 Protease of Zika Virus.

The NS2B/NS3 protease is highly conserved among various proteases of the Zika virus, making it an important therapeutic target for developing broad-spectrum antiviral drugs. The NS2B/NS3 protease is a crucial enzyme in the replication cycle of Zika virus and plays a significant role in viral maturation and assembly. Inhibiting the activity of this protease can potentially prevent viral replication, making it an attractive target for developing therapies against Zika virus infection. This work screens 429 antiviral peptides in comparison with substrate peptide against the NS2B/NS3 of Zika virus using molecular docking and molecular dynamics (MD) simulation. Based on the docking screening, MD simulation conducted for the best four peptides including AVP0239, AVP0642, AVP0660, and AVP2044, could be effective against NS2B/NS3. These results were compared with the control substrate peptide. Further analysis indicates that AVP0642 and AVP2044 are the most promising candidates. The interaction analysis showed that the catalytic site residues including His51, Asp75, Ser135 and other non-catalytic residues such as Asp129, Asp83, and Asp79 contribute substantial interactions. Hydrogen bonds (41%) and hydrophobic interactions (33%) are observed as the prominent non-covalent interaction prompting the peptide-protein complex formation. Furthermore, the structure-activity relationship (SAR) illustrates that positively charged (Lys, Arg) residues in the peptides dominate the interactions. This study provides the basis for developing novel peptide-based protease inhibitors for Zika virus.

Structure and dynamics of whole-sequence homology model of ORF3a protein of SARS-CoV-2: An insight from microsecond molecular dynamics simulations.

The ORF3a is a large accessory protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which plays an important role in virulence and viral replication; especially in inflammasome activation and apoptosis. However,, the existing cryo-EM structure of SARS-CoV-2 ORF3a is incomplete, . making it challenging to understand its structural and functional features. The aim of this study is to investigate the dynamic behaviors of the full-sequence homology model of ORF3a and compare it with the cryo-EM structure using microsecond molecular dynamics simulations. The previous studies indicated that the unresolved residues of the cryo-EM structure are not only involved in the pathogenesis of the SARS-CoV-2 but also exhibit a significant antigenicity. The dynamics scenario of homology model revealed higher RMSD, Rg, and SASA values with stable pattern when compared to the cryo-EM structure. Moreover, the RMSF analysis demonstrated higher fluctuations at specific positions (1-43, 97-110, 172-180, 219-243) in the model structure, whereas the cryo-EM structure displayed lower overall drift (except 1-43) in comparison to the model structure.Secondary structural features indicated that a significant unfolding in the transmembrane domains and ß-strand at positions 166 to 172, affecting the stability and compactness of the cryo-EM structure , whereas the model exhibited noticeable unfolding in transmembrane domains and small-coiled regions in the N-terminal. , The results from molecular docking and steered molecular dynamics investigations showed the model structure had a greater number of non-bonding interactions, leading to enhanced stability when compared to the cryo-EM structure. Consequently, higher forces were necessary for unbinding of the baricitinib and ruxolitinib inhibitors from the model structure.. Our findings can help better understanding of the significance of unresolved residues at the molecular level. Additionally, this information can guide researchers for experimental endeavors aimed at completing the full-sequence structure of the ORF3a.Communicated by Ramaswamy H. Sarma.

Unraveling the impact of ORF3a Q57H mutation on SARS-CoV-2: insights from molecular dynamics.

ORF3a is a conserved accessory protein of SARS-CoV-2, linked to viral infection and pathogenesis, with acquired mutations at various locations. Previous studies have shown that the occurrence of the Q57H mutation is higher in comparison to other positions in ORF3a. This mutation is known to induce conformational changes, yet the extent of structural alteration and its role in the viral adaptation process remain unknown. Here we performed molecular dynamics (MD) simulations of wt-ORF3a, Q57H, and Q57A mutants to analyze structural changes caused by mutations compared to the native protein. The MD analysis revealed that Q57H and Q57A mutants show significant structural changes in the dimer conformation than the wt-ORF3a. This dimer conformer narrows down the ion channel cavity, which reduces Na + or K + permeability leading to decrease the antigenic response that can help the virus to escape the host immune system. Non-bonding interaction analysis shows the Q57H mutant has more interacting residues, resulting in more stability within dimer conformation than the wt-ORF3a and Q57A. Moreover, both mutant dimers (Q57H and Q57A) form a novel salt-bridge interaction at the same position between A:Asp142 and B:Lys61, whereas such an interaction is absent in the wt-ORF3a dimer. We have also noticed that the TM3 domain's flexibility in Q57H is increased because of strong inter-domain interactions of TM1 and TM2 within the dimer conformation. These unusual interactions and flexibility of Q57H mutant can have significant impacts on the SARS-CoV-2 adaptations, virulence, transmission, and immune system evasion. Our findings are consistent with the previous experimental data and provided details information on the structural perturbation in ORF3a caused by mutations, which can help better understand the structural change at the molecular level as well as the reason for the high virulence properties of this variant.Communicated by Ramaswamy H. Sarma.

Structural and functional effects of the L84S mutant in the SARS-COV-2 ORF8 dimer based on microsecond molecular dynamics study.

The L84S mutation has been observed frequently in the ORF8 protein of SARS-CoV-2, which is an accessory protein involved in various important functions such as virus propagation, pathogenesis, and evading the immune response. However, the specific effects of this mutation on the dimeric structure of ORF8 and its impacts on interactions with host components and immune responses are not well understood. In this study, we performed one microsecond molecular dynamics (MD) simulation and analyzed the dimeric behavior of the L84S and L84A mutants in comparison to the native protein. The MD simulations revealed that both mutations caused changes in the conformation of the ORF8 dimer, influenced protein folding mechanisms, and affected the overall structural stability. In particular, the 73YIDI76 motif has found to be significantly affected by the L84S mutation, leading to structural flexibility in the region connecting the C-terminal ß4 and ß5 strands. This flexibility might be responsible for virus immune modulation.  The free energy landscape (FEL) and principle component analysis (PCA) have also supported our investigation. Overall, the L84S and L84A mutations affect the ORF8 dimeric interfaces by reducing the frequency of protein-protein interacting residues (Arg52, Lys53, Arg98, Ile104, Arg115, Val117, Asp119, Phe120, and Ile121) in the ORF8 dimer.  Our findings provide detail insights for further research in designing structure-based therapeutics against the SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

In silico peptide-based therapeutics against human colorectal cancer by the activation of TLR5 signaling pathways.

OBJECTIVE: Colorectal cancer (CRC) is the third leading cause of cancer-related deaths in both men and women. Toll-like receptor 5 (TLR5), an autoimmune signaling receptor that plays a role in cancer, can be exploited for the suppression of human colon cancer. Salmonella flagellin protein, a novel agonist of TLR5 activating downstream signaling, could be a basis for designing anticancer peptides. METHODS: The three-dimensional crystal structure of TLR5 (PDB ID: 3J0A, Resolution = 26.0 Å) was optimized using the AMBER force field in the YASARA suit. In silico enzymatic digestion tool, PeptideCutter, was used to identify peptides from Salmonella flagellin, an agonist against human TLR5. The 3D structure of the peptides was generated using PEP-FOLD3. These peptides were screened against human TLR5 using shape complementarity principles based on the binding affinity and interactions with the active residue of TLR5 monomer, and the selected peptides were further validated by molecular dynamic (MD) simulation. RESULTS: In this study, we generated 42 peptides from Salmonella flagellin protein by in silico protein digestion. Then, based on a new hidden Markov model sub-optimal conformation sampling approach as well as the size of the fragments, we select 38 effective peptides from these 42 cleavages. These peptides were screened against the monomeric Xray structure of human TLR5 using shape complementarity principles. Based on the binding affinity and interactions with the active residue of TLR5 monomer (residues 294 and 366 of TLR5), nine top-scored peptides were selected for the initial molecular dynamic (MD) simulation. Among these peptides, Clv10, Clv17, and Clv28 showed high stability and less flexibility during MD simulation. A 1 µs MD simulation was performed on TLR5-Clv10, TLR-Clv17, and TLR5-Clv28 complexes to further analyze the stability, conformational changes, and binding mode (Clv10, Clv17, and Clv28). During this MD study, the peptides showed high salt bridges and ionic interactions with residue ASP294 and residue ASP366 throughout the simulation and remained in the concave of the human TLR5 monomer. The RMSD and Rg values showed that the peptide-protein complexes become stable after 200 ns of contraction and extraction. CONCLUSION: These findings can facilitate the rational design of selected peptides as an agonist of TLR5, which have antitumor activity, suppress colorectal cancer tumors, and can be used as promising candidates and novel agonists of TLR5.

Large scale peptide screening against main protease of SARS CoV-2.

The COVID-19 pandemic has been a public health emergency, with deadly forms constantly emerging around the world, highlighting the dire need for highly effective antiviral therapeutics. Peptide therapeutics show significant potential for this viral disease due to their efficiency, safety, and specificity. Here, two thousand seven hundred eight antibacterial peptides were screened computationally targeting the Main protease (Mpro) of SARS CoV-2. Six top-ranked peptides according to their binding scores, binding pose were investigated by molecular dynamics to explore the interaction and binding behavior of peptide-Mpro complexes. The structural and energetic characteristics of Mpro-DRAMP01760 and Mpro-DRAMP01808 complexes fluctuated less during a 250 ns MD simulation. In addition, three peptides (DRAMP01760, DRAMP01808, and DRAMP01342) bind strongly to Mpro protein, according to the free energy landscape and principal component analysis. Peptide helicity and secondary structure analysis are in agreement with our findings. Interaction analysis of protein-peptide complexes demonstrated that Mpro's residue CYS145, HIS41, PRO168, GLU166, GLN189, ASN142, MET49, and THR26 play significant contributions in peptide-protein attachment. Binding free energy analysis (MM-PBSA) demonstrated the energy profile of interacting residues of Mpro in peptide-Mpro complexes. To summarize, the peptides DRAMP01808 and DRAMP01760 may be highly Mpro specific, resulting disruption in a viral replication and transcription. The results of this research are expected to assist future research toward the development of antiviral peptide-based therapeutics for Covid-19 treatment.

A review on structural, non-structural, and accessory proteins of SARS-CoV-2: Highlighting drug target sites.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is a highly transmittable and pathogenic human coronavirus that first emerged in China in December 2019. The unprecedented outbreak of SARS-CoV-2 devastated human health within a short time leading to a global public health emergency. A detailed understanding of the viral proteins including their structural characteristics and virulence mechanism on human health is very crucial for developing vaccines and therapeutics. To date, over 1800 structures of non-structural, structural, and accessory proteins of SARS-CoV-2 are determined by cryo-electron microscopy, X-ray crystallography, and NMR spectroscopy. Designing therapeutics to target the viral proteins has several benefits since they could be highly specific against the virus while maintaining minimal detrimental effects on humans. However, for ongoing and future research on SARS-CoV-2, summarizing all the viral proteins and their detailed structural information is crucial. In this review, we compile comprehensive information on viral structural, non-structural, and accessory proteins structures with their binding and catalytic sites, different domain and motifs, and potential drug target sites to assist chemists, biologists, and clinicians finding necessary details for fundamental and therapeutic research.

Potential antiviral peptides against the nucleoprotein of SARS-CoV-2.

Nucleoprotein is a conserved structural protein of SARS-CoV-2, which is involved in several functions, including replication, packaging, and transcription. In this research, 21 antiviral peptides that are known to have inhibitory function against nucleoprotein in several other viruses, were screened computationally against the nucleoprotein of SARS-CoV-2. The complexes of five best performing peptides (AVP1142, AVP1145, AVP1148, AVP1150, AVP1155) with nucleoprotein were selected for subsequent screening via 5 ns molecular dynamics (MD) simulation. Two peptides, namely AVP1145 and AVP1155, came out as promising candidates and hence were selected for 200 ns MD simulation for further validation, incorporating a DMPC-based membrane environment. In the long MD simulation, both AVP1155 and AVP1145 utilized multiple residues-mainly aromatic, acidic, and nonpolar residues-as interacting points to remain in contact with the nucleoprotein and formed predominantly hydrogen bonds along with hydrophobic and electrostatic interactions. However, AVP1155 proved to be superior to AVP1145 when its complex with nucleoprotein was analyzed in terms of root-mean-square deviation, root-mean-square fluctuation, radius of gyration, solvent accessible surface area and free energy landscape. In a nutshell, the findings of this research may guide future studies in the development of selective peptide inhibitors of SARS-CoV-2 nucleoprotein. Supplementary Information: The online version contains supplementary material available at 10.1007/s11696-022-02514-4.

Isolation and In Silico Prediction of Potential Drug-like Compounds with a New Dimeric Prenylated Quinolone Alkaloid from Zanthoxylum rhetsa (Roxb.) Root Extracts Targeted against SARS-CoV-2 (Mpro).

A new dimeric prenylated quinolone alkaloid, named 2,11-didemethoxy-vepridimerine A, was isolated from the root bark of Zanthoxylum rhetsa, together with twelve known compounds. The structure of the new compound was elucidated on the basis of spectroscopic investigations (NMR and Mass). The interaction of the isolated compounds with the main protease of SARS-CoV-2 (Mpro) was evaluated using molecular docking followed by MD simulations. The result suggests that 2,11-didemethoxy-vepridimerine A, the new compound, has the highest negative binding affinity against the Mpro with a free energy of binding of -8.5 Kcal/mol, indicating interaction with the Mpro. This interaction was further validated by 100 ns MD simulation. This implies that the isolated new compound, which can be employed as a lead compound for an Mpro-targeting drug discovery program, may be able to block the action of Mpro.

Anatomy of Protein Electrospray Ionization Mass Spectra by Superconducting Tunnel Junction Mass and Energy Spectrometry.

Cryogenic superconducting tunnel junction (STJ) detectors have the advantage of single-particle sensitivity, high quantum efficiency, low noise, and the ability to detect the time and relative impact energy of deposited ions. This makes them attractive for use in mass spectrometry (MS) and as a form of energy spectrometry. STJ cryodetectors have been coupled to time-of-flight (TOF) mass spectrometers equipped with a matrix-assisted laser desorption ionization (MALDI) source and to an electrospray ionization (ESI) TOF mass spectrometer. Here, a lab-made linear quadrupole ion trap (LIT) mass spectrometer system was coupled to an ESI source and a 16-channel Nb-STJ array with improved readout electronics. The goal was to investigate fundamentals of ESI-generated protein ions by further exploiting the advantage of resolving these ions in a third dimension of the relative energy deposited into the STJs. The proteins equine cytochrome c, bovine carbonic anhydrase, bovine serum albumin, and murine immunoglobulin G were studied using this ESI-LIT-STJ-MS instrument. Multiply charged monomers, multimers, and fragments from metastable ions were resolved from monomer peaks by differences in ion deposition energy even when these ions have the same mass-to-charge ratio as the corresponding monomer. The determination of a fragment mass from metastable decomposition is accomplished without knowing the charge state of the fragment. The average charge state of the multimers is reduced with each addition of a protein which is presumed to be a direct reflection of the surface area available for charging. Multiply charged in-source fragments have also been observed and distinguished in the mass spectrum of carbonic anhydrase by using the differences in the energy deposited in the STJs.

Insights into the leaves of Ceriscoides campanulata: Natural proanthocyanidins alleviate diabetes, inflammation, and esophageal squamous cell cancer via in vitro and in silico models.

Fourteen flavones (1-14) including twelve polymethoxylated flavones, two A-type proanthocyanidins (oligomeric flavonoids) (15, 16), one benzoyl glucoside (17), one triterpenoid (18), and one phenylpropanoid (19) were isolated from the leaves of the South Asian medicinal plant Ceriscoides campanulata (Roxb.) Tirveng (Rubiaceae). The structures of the compounds were identified based on their spectroscopic and spectrometric data and in comparison with literature data. Isolated compounds were tested in vitro against inflammatory enzymes (COX-2, iNOS), pro-inflammatory cytokines (IL-1ß, IL-6, TNF-α), esophageal squamous carcinoma cell line (TE13), and carbohydrate digestion enzymes (α-amylase, α-glucosidase). Proanthocyanidins 15 and 16 significantly attenuated the LPS-induced inflammatory response of COX-2, iNOS, IL-1ß, IL-6, TNF-α in RAW 264.7 cells. Proanthocyanidins also satisfactorily inhibited the regrowth (64%), migration (51%), and formation of tumor-spheres (48%) in ESCC cell line TE13 at 50% toxic concentration. Compounds 15 and 16 showed the most potent effect against mammalian α-amylase (IC50 8.4 ± 0.3 µM and 3.5 ± 0.0 µM, respectively) compared to reference standard acarbose (IC50 5.9 ± 0.1 µM). As yeast α-glucosidase inhibitors, compounds 15 and 16 also displayed significant activities (IC50 6.2 ± 0.3 and 4.7 ± 0.1 µM, respectively), while compounds 1-6 displayed weaker α-glucosidase inhibitory activities, ranging from 49 to 142 µM, compared to acarbose (IC50 665 ± 42 µM). In an anticholinesterase assay, compounds 1, 2, 6 (IC50 51 ± 2, 53 ± 7, 64 ± 5 µM, respectively), and 4 (IC50 44 ± 1 µM) showed moderate inhibitory activities against acetylcholinesterase and butyrylcholinesterase, respectively. Furthermore, molecular docking and molecular dynamic simulation analyses of compounds 15 and 16 were performed against human pancreatic α-amylase and human lysosomal acid α-glucosidase to elucidate the interactions of these compounds in the respective enzymes' active sites.

Structure and dynamics of membrane protein in SARS-CoV-2.

SARS-CoV-2 membrane (M) protein performs a variety of critical functions in virus infection cycle. However, the expression and purification of membrane protein structure is difficult despite tremendous progress. In this study, the 3 D structure is modeled followed by intensive validation and molecular dynamics simulation. The lack of suitable homologous templates (>30% sequence identities) leads us to construct the membrane protein models using template-free modeling (de novo or ab initio) approach with Robetta and trRosetta servers. Comparing with other model structures, it is evident that trRosetta (TM-score: 0.64; TM region RMSD: 2 Å) can provide the best model than Robetta (TM-score: 0.61; TM region RMSD: 3.3 Å) and I-TASSER (TM-score: 0.45; TM region RMSD: 6.5 Å). 100 ns molecular dynamics simulations are performed on the model structures by incorporating membrane environment. Moreover, secondary structure elements and principal component analysis (PCA) have also been performed on MD simulation data. Finally, trRosetta model is utilized for interpretation and visualization of interacting residues during protein-protein interactions. The common interacting residues including Phe103, Arg107, Met109, Trp110, Arg131, and Glu135 in the C-terminal domain of M protein are identified in membrane-spike and membrane-nucleocapsid protein complexes. The active site residues are also predicted for potential drug and peptide binding. Overall, this study might be helpful to design drugs and peptides against the modeled membrane protein of SARS-CoV-2 to accelerate further investigation. Communicated by Ramaswamy H. Sarma.

Cysteine focused covalent inhibitors against the main protease of SARS-CoV-2.

In viral replication and transcription, the main protease (Mpro) of SARS-CoV-2 plays an important role and appears to be a vital target for drug design. In Mpro, there is a Cys-His catalytic dyad, and ligands that interact with the Cys145 assumed to be an effective approach to inhibit the Mpro. In this study, approximately 1400 cysteine-focused ligands were screened to identify the best candidates that can act as potent inhibitors against Mpro. Our results show that the selected ligands strongly interact with the key Cys145 and His41 residues. Covalent docking was performed for the selected candidates containing the acrylonitrile group, which can form a covalent bond with Cys145. All atoms molecular dynamics (MD) simulation was performed on the selected four inhibitors including L1, L2, L3 and L4 to validate the docking interactions. Our results were also compared with a control ligand, α-ketoamide (11r). Principal component analysis on structural and energy data obtained from the MD trajectories shows that L1, L3, L4 and α-ketoamide (11r) have structural similarity with the apo-form of the Mpro. Quantitative structure-activity relationship method was employed for pattern recognition of the best ligands, which discloses that ligands containing acrylonitrile and amide warheads can show better performance. ADMET analysis displays that our selected candidates appear to be safer inhibitors. Our combined studies suggest that the best cysteine focused ligands can help to design an effective lead drug for COVID-19 treatment. Communicated by Ramaswamy H. Sarma.

A student led computational screening of peptide inhibitors against main protease of SARS-CoV-2.

The main protease of SARS-CoV-2 is a promising drug target due to its functional role as a catalytic dyad in mediating proteolysis during the viral life cycle. In this study, experimentally proven 14 HIV protease peptides were screened against the main protease of SARS-CoV-2. Fourteen middle and high school "student researchers" were trained on relevant computational tools, provided with necessary biological and chemical background and scientific article writing. They performed the primary screening via molecular docking and the best performing complexes were subjected to molecular dynamics simulations. Molecular docking revealed that HIP82 and HIP1079 can bind with the catalytic residues, however after molecular dynamics simulation only HIP1079 retained its interaction with the catalytic sites. The student researchers were also trained to write scientific article and were involved with drafting of the manuscript. This project provided the student researchers an insight into multi-disciplinary research in biology and chemistry, inspired them about practical approaches of computational chemistry in solving a real-world problem like a global pandemic. This project also serves as an example to introduce scientific inquiry, research methodology, critical thinking, scientific writing, and communication for high school students.

Identification of potent inhibitors against transmembrane serine protease 2 for developing therapeutics against SARS-CoV-2.

In viral binding and entry, the Spike(S) protein of SARS-CoV-2 uses transmembrane serine protease 2 (TMPRSS2) for priming to cleavage themselves. In this study, we have screened 'drug-like' 7476 ligands and found that over thirty ligands can effectively inhibit the TMPRSS-2 better than the control ligand. Finally, the three best drug agents L1, L2, and L6 were selected according to their average binding affinities and fitting score. These ligands interact with Asp435, Cys437, Ser436, Trp461, and Cys465 amino acid residues. The three best candidates and a reported drug Nafamostat mesylate (NAM) were selected to run 250 ns molecular dynamics (MD) simulations. Various properties of ligand-protein interactions obtained from MD simulation such as bonds, angle, dihedral, planarity, coulomb, and van der Waals (VdW) were used for principal component analysis (PCA) calculation. PCA discloses the evidence of the structural similarities to the corresponding complexes of L1, L2, and L6 with the complex of TMPRSS2(TM) and Nafamostat mesylate (TM-NAM). Moreover, Quantitative structure-activity relationship (QSAR) pattern recognition was generated using PCA for the investigation of structural similarities among the selected ligands. Multiple Linear Regression (MLR) model was built to predict the binding energy compared to the binding energy obtained from molecular docking. The MLR regression model reveals an accuracy of 80% for the prediction of the binding energy of ligands. ADMET analysis demonstrates that these drug agents are appeared to be safer inhibitors. These three ligands can be used as potential inhibitors against the TMPRSS2.Communicated by Ramaswamy H. Sarma.

Designing potent inhibitors against the multidrug resistance P-glycoprotein.

The multidrug transporter P-glycoprotein is an ATP binding cassette (ABC) exporter responsible for resistance to tumor cells during chemotherapy. This study was designed with computational approaches aimed at identifying the best potent inhibitors of P-glycoprotein. Although many compounds have been suggested to inhibit P-glycoprotein, however, their information on bioavailability, selectivity, ADMET properties, and molecular interactions has not been revealed. Molecular docking, ADMET analysis, molecular dynamics, Principal component analysis (PCA), and binding free energy calculations were performed. Two compounds D1 and D2 showed the best docking score against P-glycoprotein and both compounds have 4-thiazolidinone derivatives containing indolin-3 one moiety are novel anti-tumor compounds. ADMET calculation analysis predicted D1 and D2 to have acceptable pharmacokinetic properties. The MD simulation discloses that D1-P-glycoprotein and D2-P-glycoprotein complexes are in stable conformation as apo-form. Hydrophobic amino acid such as phenylalanine plays significant on the interactions of inhibitors. Principal component analysis shows that both complexes are relatively similar variables as apo-form except planarity and Columbo energy profile. In addition, Quantitative Structural Activity Relationship (QSAR) of the ligand candidates were subjected to the principal component analysis (PCA) for pattern recognition. Partial-least-square (PLS) regression analysis was further utilized to model drug candidates' QSAR for subsequent prediction of the binding energy of validated drug candidates. PCA revealed groupings of the drug candidates based on the similarity or differences in drug candidates QSAR. Moreover, the developed PLS regression accurately predicted the values of the binding energy of drug candidates, with low residual error of prediction.Communicated by Ramaswamy H. Sarma.