Displaying publications 1 - 20 of 60 in total

Abstract:
Sort:
  1. Aznan AM, Abdullah Z, Lee VS, Tiekink ER
    Acta Crystallogr Sect E Struct Rep Online, 2014 Aug 1;70(Pt 8):58-61.
    PMID: 25249854 DOI: 10.1107/S1600536814012227
    The title compound, C12H11N3O2, is a second monoclinic polymorph (P21, with Z' = 4) of the previously reported monoclinic (P21/c, with Z' = 2) form [Akhmad Aznan et al. (2010 ▶). Acta Cryst. E66, o2400]. Four independent mol-ecules comprise the asymmetric unit, which have the common features of a syn disposition of the pyridine N atom and the toluene ring, and an intra-molecular amine-nitro N-H⋯O hydrogen bond. The differences between mol-ecules relate to the dihedral angles between the rings which range from 2.92 (19) to 26.24 (19)°. The geometry-optimized structure [B3LYP level of theory and 6-311 g+(d,p) basis set] has the same features except that the entire mol-ecule is planar. In the crystal, the three-dimensional architecture is consolidated by a combination of C-H⋯O, C-H⋯π, nitro-N-O⋯π and π-π inter-actions [inter-centroid distances = 3.649 (2)-3.916 (2) Å].
  2. Tan CW, Sam IC, Chong WL, Lee VS, Chan YF
    Antiviral Res, 2017 07;143:186-194.
    PMID: 28457855 DOI: 10.1016/j.antiviral.2017.04.017
    Zika virus (ZIKV) is an arthropod-borne flavivirus that causes newborn microcephaly and Guillian-Barré syndrome in adults. No therapeutics are available to treat ZIKV infection or other flaviviruses. In this study, we explored the inhibitory effect of glycosaminoglycans and analogues against ZIKV infection. Highly sulfated heparin, dextran sulfate and suramin significantly inhibited ZIKV infection in Vero cells. De-sulfated heparin analogues lose inhibitory effect, implying that sulfonate groups are critical for viral inhibition. Suramin, an FDA-approved anti-parasitic drug, inhibits ZIKV infection with 3-5 log10 PFU viral reduction with IC50value of ∼2.5-5 μg/ml (1.93 μM-3.85 μM). A time-of-drug-addition study revealed that suramin remains potent even when administrated at 1-24 hpi. Suramin inhibits ZIKV infection by preventing viral adsorption, entry and replication. Molecular dynamics simulation revealed stronger interaction of suramin with ZIKV NS3 helicase than with the envelope protein. Suramin warrants further investigation as a potential antiviral candidate for ZIKV infection. Heparan sulfate (HS) is a cellular attachment receptor for multiple flaviviruses. However, no direct ZIKV-heparin interaction was observed in heparin-binding analysis, and downregulate or removal of cellular HS with sodium chlorate or heparinase I/III did not inhibit ZIKV infection. This indicates that cell surface HS is not utilized by ZIKV as an attachment receptor.
  3. Kara J, Suwanhom P, Wattanapiromsakul C, Nualnoi T, Puripattanavong J, Khongkow P, et al.
    Arch Pharm (Weinheim), 2019 Jul;352(7):e1800310.
    PMID: 31125474 DOI: 10.1002/ardp.201800310
    Sixteen novel coumarin-based compounds are reported as potent acetylcholinesterase (AChE) inhibitors. The most active compound in this series, 5a (IC50 0.04 ± 0.01 µM), noncompetitively inhibited AChE with a higher potency than tacrine and galantamine. Compounds 5d, 5j, and 5 m showed a moderate antilipid peroxidation activity. The compounds showed cytotoxicity in the same range as the standard drugs in HEK-293 cells. Molecular docking demonstrated that 5a acted as a dual binding site inhibitor. The coumarin moiety occupied the peripheral anionic site and showed π-π interaction with Trp278. The tertiary amino group displayed significant cation-π interaction with Phe329. The aromatic group showed π-π interaction with Trp83 at the catalytic anionic site. The long chain of methylene lay along the gorge interacting with Phe330 via hydrophobic interaction. Molecular docking was applied to postulate the selectivity toward AChE of 5a in comparison with donepezil and tacrine. Structural insights into the selectivity of the coumarin derivatives toward huAChE were explored by molecular docking and 3D QSAR and molecular dynamics simulation for 20 ns. ADMET analysis suggested that the 2-(2-oxo-2H-chromen-4-yl)acetamides showed a good pharmacokinetic profile and no hepatotoxicity. These coumarin derivatives showed high potential for further development as anti-Alzheimer agents.
  4. Saoin S, Wisitponchai T, Intachai K, Chupradit K, Moonmuang S, Nangola S, et al.
    Asian Pac J Allergy Immunol, 2018 06;36(2):126-135.
    PMID: 28802032 DOI: 10.12932/AP-280217-0037
    BACKGROUND: AnkGAG1D4 is an artificial ankyrin repeat protein which recognizes the capsid protein (CA) of the human immunodeficiency virus type 1 (HIV-1) and exhibits the intracellular antiviral activity on the viral assembly process. Improving the binding affinity of AnkGAG1D4 would potentially enhance the AnkGAG1D4-mediated antiviral activity.

    OBJECTIVE: To augment the affinity of AnkGAG1D4 scaffold towards its CA target, through computational predictions and experimental designs.

    METHOD: Three dimensional structure of the binary complex formed by AnkGAG1D4 docked to the CA was used as a model for van der Waals (vdW) binding energy calculation. The results generated a simple guideline to select the amino acids for modifications. Following the predictions, modified AnkGAG1D4 proteins were produced and further evaluated for their CA-binding activity, using ELISA-modified method and bio-layer interferometry (BLI).

    RESULTS: Tyrosine at position 56 (Y56) in AnkGAG1D4 was experimentally identified as the most critical residue for CA binding. Rational substitutions of this residue diminished the binding affinity. However, vdW calculation preconized to substitute serine for tyrosine at position 45. Remarkably, the affinity for the viral CA was significantly enhanced in AnkGAG1D4-S45Y mutant, with no alteration of the target specificity.

    CONCLUSIONS: The S-to-Y mutation at position 45, based on the prediction of interacting amino acids and on vdW binding energy calculation, resulted in a significant enhancement of the affinity of AnkGAG1D4 ankyrin for its CA target. AnkGAG1D4-S45Y mutant represented the starting point for further construction of variants with even higher affinity towards the viral CA, and higher therapeutic potential in the future.

  5. Wisitponchai T, Shoombuatong W, Lee VS, Kitidee K, Tayapiwatana C
    BMC Bioinformatics, 2017 Apr 19;18(1):220.
    PMID: 28424069 DOI: 10.1186/s12859-017-1628-6
    BACKGROUND: Computational analysis of protein-protein interaction provided the crucial information to increase the binding affinity without a change in basic conformation. Several docking programs were used to predict the near-native poses of the protein-protein complex in 10 top-rankings. The universal criteria for discriminating the near-native pose are not available since there are several classes of recognition protein. Currently, the explicit criteria for identifying the near-native pose of ankyrin-protein complexes (APKs) have not been reported yet.

    RESULTS: In this study, we established an ensemble computational model for discriminating the near-native docking pose of APKs named "AnkPlex". A dataset of APKs was generated from seven X-ray APKs, which consisted of 3 internal domains, using the reliable docking tool ZDOCK. The dataset was composed of 669 and 44,334 near-native and non-near-native poses, respectively, and it was used to generate eleven informative features. Subsequently, a re-scoring rank was generated by AnkPlex using a combination of a decision tree algorithm and logistic regression. AnkPlex achieved superior efficiency with ≥1 near-native complexes in the 10 top-rankings for nine X-ray complexes compared to ZDOCK, which only obtained six X-ray complexes. In addition, feature analysis demonstrated that the van der Waals feature was the dominant near-native pose out of the potential ankyrin-protein docking poses.

    CONCLUSION: The AnkPlex model achieved a success at predicting near-native docking poses and led to the discovery of informative characteristics that could further improve our understanding of the ankyrin-protein complex. Our computational study could be useful for predicting the near-native poses of binding proteins and desired targets, especially for ankyrin-protein complexes. The AnkPlex web server is freely accessible at http://ankplex.ams.cmu.ac.th .

  6. Eid EEM, Alshehade SA, Almaiman AA, Kamran S, Lee VS, Alshawsh MA
    Biomedicines, 2023 Jul 04;11(7).
    PMID: 37509531 DOI: 10.3390/biomedicines11071891
    Leukemia, a condition characterized by the abnormal proliferation of blood cells, poses significant challenges in cancer treatment. Thymoquinone (TQ), a bioactive compound derived from black seed, has demonstrated anticancer properties, including telomerase inhibition and the induction of apoptosis. However, TQ's poor solubility and limited bioavailability hinder its clinical application. This study explored the use of Sulfobutylether-β-cyclodextrin (SBE-β-CD), a cyclodextrin derivative, to enhance the solubility and stability of TQ for leukemia treatment. SBE-β-CD offers low hemolytic activity and has been successfully employed in controlled drug release systems. The study investigated the formation of inclusion complexes between TQ and SBE-β-CD and evaluated their effects on leukemia cell growth and telomerase activity. The results indicated that the TQ/SBE-β-CD complex exhibited improved solubility and enhanced cytotoxic effects against K-562 leukemia cells compared to TQ alone, suggesting the potential of SBE-β-CD as a drug delivery system for TQ. The annexin V-FITC assay demonstrated increased apoptosis, while the qPCR quantification assay revealed reduced telomerase activity in leukemia cells treated with TQ/SBE-β-CD, supporting its anti-leukemic potential. The molecular docking analysis indicated a strong binding affinity between TQ and telomerase. However, further research is needed to optimize the apoptotic effects and minimize necrosis induction. In conclusion, TQ/SBE-β-CD shows promise as a novel strategy for leukemia treatment by inhibiting telomerase and enhancing the cytotoxic effects of TQ, offering a potential solution to overcome the limitations of TQ's poor solubility and bioavailability.
  7. Nutho B, Pengthaisong S, Tankrathok A, Lee VS, Ketudat Cairns JR, Rungrotmongkol T, et al.
    Biomolecules, 2020 Jun 15;10(6).
    PMID: 32549280 DOI: 10.3390/biom10060907
    β-Glucosidases and β-mannosidases hydrolyze substrates that differ only in the epimer of the nonreducing terminal sugar moiety, but most such enzymes show a strong preference for one activity or the other. Rice Os3BGlu7 and Os7BGlu26 β-glycosidases show a less strong preference, but Os3BGlu7 and Os7BGlu26 prefer glucosides and mannosides, respectively. Previous studies of crystal structures with glucoimidazole (GIm) and mannoimidazole (MIm) complexes and metadynamic simulations suggested that Os7BGlu26 hydrolyzes mannosides via the B2,5 transition state (TS) conformation preferred for mannosides and glucosides via their preferred 4H3/4E TS conformation. However, MIm is weakly bound by both enzymes. In the present study, we found that MIm was not bound in the active site of crystallized Os3BGlu7, but GIm was tightly bound in the -1 subsite in a 4H3/4E conformation via hydrogen bonds with the surrounding residues. One-microsecond molecular dynamics simulations showed that GIm was stably bound in the Os3BGlu7 active site and the glycone-binding site with little distortion. In contrast, MIm initialized in the B2,5 conformation rapidly relaxed to a E3/4H3 conformation and moved out into a position in the entrance of the active site, where it bound more stably despite making fewer interactions. The lack of MIm binding in the glycone site in protein crystals and simulations implies that the energy required to distort MIm to the B2,5 conformation for optimal active site residue interactions is sufficient to offset the energy of those interactions in Os3BGlu7. This balance between distortion and binding energy may also provide a rationale for glucosidase versus mannosidase specificity in plant β-glycosidases.
  8. Leong ST, Liew SY, Khaw KY, Ahmad Hassali H, Richomme P, Derbré S, et al.
    Bioorg Chem, 2023 Dec;141:106859.
    PMID: 37742494 DOI: 10.1016/j.bioorg.2023.106859
    A bio-assay guided fractionation strategy based on cholinesterase assay combined with 13C NMR-based dereplication was used to identify active metabolites from the bark of Mesua lepidota. Eight compounds were identified with the aid of the 13C NMR-based dereplication software, MixONat, i.e., sitosterol (1), stigmasterol (2), α-amyrin (3), friedelin (6), 3β-friedelinol (7), betulinic acid (9), lepidotol A (10) and lepidotol B (11). Further bio-assay guided isolation of active compounds afforded one xanthone, pyranojacareubin (12) and six coumarins; lepidotol A (10), lepidotol B (11), lepidotol E (13), lepidotin A (14), and lepidotin B (15), including a new Mammea coumarin, lepidotin C (16). All the metabolites showed strong to moderate butyrylcholinesterase (BChE) inhibition. Lepidotin B (15) exhibited the most potent inhibition towards BChE with a mix-mode inhibition profile and a Ki value of 1.03 µM. Molecular docking and molecular dynamics simulations have revealed that lepidotin B (15) forms stable interactions with key residues within five critical regions of BChE. These regions encompass residues Asp70 and Tyr332, the acyl hydrophobic pocket marked by Leu286, the catalytic triad represented by Ser198 and His438, the oxyanion hole (OH) constituted by Gly116 and Gly117, and the choline binding site featuring Trp82. To gauge the binding strength of lepidotin B (15) and to pinpoint pivotal residues at the binding interface, free energy calculations were conducted using the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach. This analysis not only predicted a favourable binding affinity for lepidotin B (15) but also facilitated the identification of significant residues crucial for the binding interaction.
  9. Lythell E, Suardíaz R, Hinchliffe P, Hanpaibool C, Visitsatthawong S, Oliveira ASF, et al.
    Chem Commun (Camb), 2020 Jun 23;56(50):6874-6877.
    PMID: 32432618 DOI: 10.1039/d0cc02520h
    MCR (mobile colistin resistance) enzymes catalyse phosphoethanolamine (PEA) addition to bacterial lipid A, threatening the "last-resort" antibiotic colistin. Molecular dynamics and density functional theory simulations indicate that monozinc MCR supports PEA transfer to the Thr285 acceptor, positioning MCR as a mono- rather than multinuclear member of the alkaline phosphatase superfamily.
  10. Saeed IM, Lee VS, Mazari SA, Si Ali B, Basirun WJ, Asghar A, et al.
    Chem Cent J, 2017;11:10.
    PMID: 28184241 DOI: 10.1186/s13065-016-0231-7
    Amine degradation is the main significant problems in amine-based post-combustion CO2 capture, causes foaming, increase in viscosity, corrosion, fouling as well as environmental issues. Therefore it is very important to develop the most efficient solvent with high thermal and chemical stability. This study investigated thermal degradation of aqueous 30% 2-aminoethylethanolamine (AEEA) using 316 stainless steel cylinders in the presence and absence of CO2 for 4 weeks. The degradation products were identified by gas chromatography mass spectrometry (GC/MS) and liquid chromatography-time-of-flight-mass spectrometry (LC-QTOF/MS). The results showed AEEA is stable in the absence of CO2, while in the presence of CO2 AEEA showed to be very unstable and numbers of degradation products were identified. 1-(2-Hydroxyethyl)-2-imidazolidinone (HEIA) was the most abundance degradation product. A possible mechanism for the thermal degradation of AEEA has been developed to explain the formation of degradation products. In addition, the reaction energy of formation of the most abundance degradation product HEIA was calculated using quantum mechanical calculation.
  11. Woon KL, Mustapa SAS, Mohd Jamel NS, Lee VS, Zakaria MZ, Ariffin A
    Chemphyschem, 2020 Sep 17.
    PMID: 32940952 DOI: 10.1002/cphc.202000612
    Material designs that use donor and acceptor units are often found in organic optoelectronic devices. Molecular level insight into the interactions between donors and acceptors are crucial for understanding how such interactions can modify the optical properties of the organic optoelectronic materials. In this paper, tris(4-(tert-butyl)phenyl)amine (pTPA) was synthesized as a donor in order to compare with unmodified triphenylamine (TPA) in a donor-acceptor system by having 2,4,6-triphenyl-1,3,5-triazine (TRZ) as an acceptor. Dimerization of donors and acceptors occurred in solvent when the concentration of solute is high. At 0 K, using a polarizable continuum model, the nitrogen atom of TPA is found to stack on top of the center of triazine of TRZ, whereas such alignment is offset in pTPA and TRZ. We attributed such alignment in TPA-TRZ as the result of attractive interactions between partial localization of 2pz electrons at the nitrogen atom of TPA and the π deficiency of triazine in TPA-TRZ. By taking into account random motions of the solvent effect at 300 K in quantum molecular dynamics and classical molecular dynamics simulations to interpret the marked difference in emission spectra between TPA-TRZ and pTPA-TRZ, it was revealed that the attractive interaction between pTPA and TRZ in toluene is weaker than TPA and TRZ. Because of the weaker attractive interaction between pTPA and TRZ in toluene, the dimers adopted numerous ground state conformations resulting in broad emission bands superimposed with multiple small Gaussian peaks. This is in contrast to TPA-TRZ which has only one dominant dimer conformation. This study demonstrates that the strength of intermolecular interactions between donors and acceptors should be taken into consideration in designing supramolecular structures.
  12. Lee VS, Sukumaran SD, Tan PK, Kuppusamy UR, Arumugam B
    Comput Biol Chem, 2021 Jun;92:107501.
    PMID: 33989998 DOI: 10.1016/j.compbiolchem.2021.107501
    Naturally occurring proteins are emerging as novel therapeutics in the protein-based biopharmaceutical industry for the treatment of diabetes and obesity. However, proteins are not suitable for oral delivery due to short half-life, reduced physical and chemical stability and low permeability across the membrane. Chemical modification has been identified as a formulation strategy to enhance the stability and bioavailability of protein drugs. The present study aims to study the effect of charge-specific modification of basic amino acids (Lys, Arg) and guanidination on the interaction of insulin with its receptor using molecular modelling. Our investigation revealed that the guanidination of insulin (Lys-NHC = NHNH2) enhanced and exerted stronger binding of the protein to its receptor through electrostatic interaction than native insulin (Lys-NH3+). Point mutations of Lys and Arg (R22, K29; R22K, K29; R22, K29R; R22K, K29R) were attempted and the effects on the interaction and stability between insulin/modified insulins and insulin receptor were also analyzed in this study. The findings from the study are expected to provide a better understanding of the possible mechanism of action of the modified protein at a molecular level before advancing to real experiments.
  13. Al-Nema M, Gaurav A, Lee VS
    Comput Biol Med, 2023 Apr 03;159:106869.
    PMID: 37071939 DOI: 10.1016/j.compbiomed.2023.106869
    In recent years, the PDE1B enzyme has become a desirable drug target for the treatment of psychological and neurological disorders, particularly schizophrenia disorder, due to the expression of PDE1B in brain regions involved in volitional behaviour, learning and memory. Although several inhibitors of PDE1 have been identified using different methods, none of these inhibitors has reached the market yet. Thus, searching for novel PDE1B inhibitors is considered a major scientific challenge. In this study, pharmacophore-based screening, ensemble docking and molecular dynamics simulations have been performed to identify a lead inhibitor of PDE1B with a new chemical scaffold. Five PDE1B crystal structures have been utilised in the docking study to improve the possibility of identifying an active compound compared to the use of a single crystal structure. Finally, the structure-activity- relationship was studied, and the structure of the lead molecule was modified to design novel inhibitors with a high affinity for PDE1B. As a result, two novel compounds have been designed that exhibited a higher affinity to PDE1B compared to the lead compound and the other designed compounds.
  14. Abdullah AA, Lee YK, Chin SP, Lim SK, Lee VS, Othman R, et al.
    Curr Med Chem, 2020;27(30):4945-5036.
    PMID: 30514185 DOI: 10.2174/0929867326666181204155336
    To date, there is still no approved anti-dengue agent to treat dengue infection in the market. Although the only licensed dengue vaccine, Dengvaxia is available, its protective efficacy against serotypes 1 and 2 of dengue virus was reported to be lower than serotypes 3 and 4. Moreover, according to WHO, the risk of being hospitalized and having severe dengue increased in seronegative individuals after they received Dengvaxia vaccination. Nevertheless, various studies had been carried out in search of dengue virus inhibitors. These studies focused on the structural (C, prM, E) and non-structural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5) of dengue virus as well as host factors as drug targets. Hence, this article provides an overall up-to-date review of the discovery of dengue virus inhibitors that are only targeting the structural and non-structural viral proteins as drug targets.
  15. Wongrattanakamon P, Lee VS, Nimmanpipug P, Jiranusornkul S
    Data Brief, 2016 Dec;9:35-42.
    PMID: 27626051 DOI: 10.1016/j.dib.2016.08.004
    The data is obtained from exploring the modulatory activities of bioflavonoids on P-glycoprotein function by ligand-based approaches. Multivariate Linear-QSAR models for predicting the induced/inhibitory activities of the flavonoids were created. Molecular descriptors were initially used as independent variables and a dependent variable was expressed as pFAR. The variables were then used in MLR analysis by stepwise regression calculation to build the linear QSAR data. The entire dataset consisted of 23 bioflavonoids was used as a training set. Regarding the obtained MLR QSAR model, R of 0.963, R (2)=0.927, [Formula: see text], SEE=0.197, F=33.849 and q (2)=0.927 were achieved. The true predictabilities of QSAR model were justified by evaluation with the external dataset (Table 4). The pFARs of representative flavonoids were predicted by MLR QSAR modelling. The data showed that internal and external validations may generate the same conclusion.
  16. Waris KH, Lee VS, Mohamad S
    Environ Sci Pollut Res Int, 2021 Sep;28(35):47785-47799.
    PMID: 34296410 DOI: 10.1007/s11356-021-15434-9
    The aim of this review is to highlight and provide an update on the current development of pesticide remediation methods, focusing on the utilization of different cyclodextrin (CD) molecules. Because of less environmental impact and non-toxic nature, CDs are beneficial for pesticide remediation, reducing environmental risk and health hazards. They are advantageous for the removal of pesticides from contaminated areas, as well as for better pesticide formulation and, posing significant effects on the hydrolysis or degradation of pesticides. The review focuses on the current trend and innovations regarding the methods and strategies employed for using CDs in designing pesticide remediation. Nowadays, in addition to the conventional experimental techniques, molecular simulation approaches are significantly contributing to the study of such phenomena and hence are recognized as a widely used tool.
  17. Nokinsee D, Shank L, Lee VS, Nimmanpipug P
    Enzyme Res, 2015;2015:262364.
    PMID: 26788364 DOI: 10.1155/2015/262364
    Tyrosinase is a key enzyme in melanogenesis. Generally, mushroom tyrosinase from A. bisporus had been used as a model in skin-whitening agent tests employed in the cosmetic industry. The recently obtained crystal structure of bacterial tyrosinase from B. megaterium has high similarity (33.5%) to the human enzyme and thus it was used as a template for constructing of the human model. Binding of tyrosinase to a series of its inhibitors was simulated by automated docking calculations. Docking and MD simulation results suggested that N81, N260, H263, and M280 are involved in the binding of inhibitors to mushroom tyrosinase. E195 and H208 are important residues in bacterial tyrosinase, while E230, S245, N249, H252, V262, and S265 bind to inhibitors and are important in forming pi interaction in human tyrosinase.
  18. Al-Nema M, Gaurav A, Lee VS
    Heliyon, 2020 Sep;6(9):e04856.
    PMID: 32984588 DOI: 10.1016/j.heliyon.2020.e04856
    Inhibition of phosphodiesterase 4 (PDE4) is a promising therapeutic approach for the treatment of inflammatory pulmonary disorders, i.e. asthma and chronic obstructive pulmonary disease. However, the treatment with non-selective PDE4 inhibitors is associated with side effects such as nausea and vomiting. Among the subtypes of PDE4 inhibited by these inhibitors, PDE4B is expressed in immune, inflammatory and airway smooth muscle cells, whereas, PDE4D is expressed in the area postrema and nucleus of the solitary tract. Thus, PDE4D inhibition is responsible for the emetic response. In this regard, a selective PDE4B inhibitor is expected to be a potential drug candidate for the treatment of inflammatory pulmonary disorders. Therefore, a shared feature pharmacophore model was developed and used as a query for the virtual screening of Maybridge and SPECS databases. A number of filters were applied to ensure only compounds with drug-like properties were selected. Accordingly, nine compounds have been identified as final hits, where HTS04529 showed the highest affinity and selectivity for PDE4B over PDE4D in molecular docking. The docked complexes of HTS04529 with PDE4B and PDE4D were subjected to molecular dynamics simulations for 100ns to assess their binding stability. The results showed that HTS04529 was bound tightly to PDE4B and formed a more stable complex with it than with PDE4D.
  19. Chong WL, Saparpakorn P, Sangma C, Lee VS, Hannongbua S
    Heliyon, 2023 Jan;9(1):e12667.
    PMID: 36618128 DOI: 10.1016/j.heliyon.2022.e12667
    SARS-CoV-2 virus continues to evolve and mutate causing most of the mutated variants resist to many of the therapeutic monoclonal antibodies (mAbs). Despite several mAbs retained neutralizing capability for Omicron BA.1 and BA.2, reduction in neutralization potency was reported. Hence, effort of searching for mAb that is broader in neutralization breadth without losing the neutralizing ability is continued. MW06 was reported with capability in neutralizing most of the variants of concern (VOC) and it binds to the conserved region (left flank) near epitope mAb sotrovimab (S309). In this study, binding affinity of mAb MW06 and its cocktail formulation with MW05 for receptor binding domain (RBD) SARS-CoV-2 virus was investigated under molecular dynamics simulations (MDs). Binding free energies computed by Molecular Mechanics Generalised Born Surface Area (MM-GBSA) algorithm predicted the binding affinity of MW06 for RBD BA.1 (-53 kcal/mol) as strong as RBD wildtype (-58 kcal/mol) while deterioration was observed for RBD BA.2 (-43 kcal/mol). Alike S309 and MW06, simulated cocktail mAb (MW05 and MW06)-RBD interactions suggested the neutralizing capability of the cocktail formulation for RBD BA.1 and BA.2 reduced. Meanwhile, residue pairs that favour the communication between the mAb and RBD have been identified by decomposing the free energy per pairwise residue basis. Apart from understanding the effects of mutation occurred in the RBD region on human angiotensin-converting enzyme 2 (hACE2) binding, impact of heavily mutated RBD on mAb-RBD interactions was investigated in this study as well. In addition to energetic profile obtained from MDs, plotting the dynamics cross-correlation map of the mAb-RBD complex under elastic network model (ENM) was aimed to understand the cross-correlations between residue fluctuations. It allows simple and rapid analysis on the motions or dynamics of the protein residues of mAbs and RBD in complex. Protein residues having correlated motions are normally part of the structural domains of the protein and their respective motions and protein function are related. Motion of mutated RBD residues and mAb residues was less correlated while their respective interactions energy computed to be higher. The combined techniques of MDs and ENM offered simplicity in understanding dynamics and energy contribution that explain binding affinity of mAb-RBD complexes.
Related Terms
Filters
Contact Us

Please provide feedback to Administrator (afdal@afpm.org.my)

External Links