<?xml version="1.0" encoding="utf-8" ?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.2 20190208//EN"
                  "JATS-archivearticle1.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="1.2" article-type="other">
<front>
<journal-meta>
<journal-id></journal-id>
<journal-title-group>
</journal-title-group>
<issn></issn>
<publisher>
<publisher-name></publisher-name>
</publisher>
</journal-meta>
<article-meta>
<permissions>
</permissions>
</article-meta>
</front>
<body>
<p><bold>Unraveling the Molecular Mechanism of Curcumin Inhibition
against White Spot Syndrome Virus VP26: An Integrated Computational
Study</bold></p>
<table-wrap>
  <table>
    <colgroup>
      <col width="18%" />
      <col width="62%" />
      <col width="20%" />
    </colgroup>
    <thead>
      <tr>
        <th><inline-graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image1.jpeg" />ajbms.knu.edu.af</th>
        <th><p><bold>Afghanistan Journal of Basic Medical
        Sciences</bold></p>
        <p>2026 Jan; 3(1): 38-52.</p></th>
        <th><graphic mimetype="image" mime-subtype="png" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image2.png" />
        <p>ISSN: 3005-6632</p></th>
      </tr>
    </thead>
    <tbody>
    </tbody>
  </table>
</table-wrap>
<p>Ahmad Wali Ataye <sup>1,2</sup>, Nasir Nazari <sup>3,4</sup>, Mahdi
Nowroozi <sup>3,4</sup>, Ghulam Aishan Sabawoon<sup>5</sup>, Nooh Amin
<sup>3,4</sup>, *Dawood Hossaini <sup>4</sup>, Aziz-ur-Rahman Niazi
<sup>1,6</sup></p>
<list list-type="order">
  <list-item>
    <p><italic>Department of Public Health, Faculty of Medicine, Afghan
    International Islamic University of Kabul, Afghanistan</italic></p>
  </list-item>
  <list-item>
    <p><italic>Department of Microbiology, Faculty of Allied, Kabul
    University of Medical Science, Kabul, Afghanistan.</italic></p>
  </list-item>
  <list-item>
    <p><italic>Medical Research and Technology Center, Khatam Al-Nabieen
    University, Kabul, Afghanistan</italic></p>
  </list-item>
  <list-item>
    <p><italic>Department of Paraclinical Sciences and Laboratory
    Management, Faculty of Medical Laboratory Technology, Khatam
    Al-Nabieen University, Kabul, Afghanistan</italic></p>
  </list-item>
  <list-item>
    <p><italic>Pharmacology department, Faculty of pharmacy, Shifa
    University, Kabul, Afghanistan</italic></p>
  </list-item>
  <list-item>
    <p><italic>Department of Public Health and Infectious Diseases,
    Faculty of Medicine, Herat University, Herat,
    Afghanistan</italic></p>
  </list-item>
</list>
<table-wrap>
  <table>
    <colgroup>
      <col width="29%" />
      <col width="71%" />
    </colgroup>
    <thead>
      <tr>
        <th><bold>A R ART I C L E I N F O</bold></th>
        <th><bold>A B S T R A C T</bold></th>
      </tr>
    </thead>
    <tbody>
      <tr>
        <td><p><bold>Type: Original Article</bold></p>
        <p>Received: 4 Dec, 2025</p>
        <p>Accepted: 29 Dec, 2025</p>
        <p><sup>*</sup>Corresponding Author:</p>
        <p>E-mails: <email>dawood.hossaini@knu.edu.af</email></p>
        <p><email>Dawoodhossaini75@gmail.com</email></p>
        <p><bold>To cite this article:</bold></p>
        <p>Ataye AW, Nazari N, Nowroozi M, Sabawoon GA, Amin N, Hossaini
        D, Niazi AR. Unraveling the molecular mechanism of curcumin
        inhibition against white spot syndrome virus VP26: an integrated
        computational study. Afghanistan Journal of Basic Medical
        Sciences. 2026 Jan;3(1):38–52.</p>
        <p>DOI:</p>
        <p><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.62134/khatamuni.143">https://doi.org/10.62134/khatamuni.143</ext-link></p></td>
        <td><p><bold>Background:</bold> The White Spot Syndrome Virus
        (WSSV) poses a significant threat to global shrimp aquaculture,
        necessitating the development of effective antiviral agents.
        Among the structural proteins of WSSV, VP26 plays a critical
        role in viral assembly and host-cell interactions, making it a
        promising target for therapeutic intervention. Curcumin, a
        bioactive compound derived from Curcuma longa, exhibits
        broad-spectrum antiviral potential, but its molecular mechanism
        of inhibition against WSSV remains unclear.</p>
        <p><bold>Methods:</bold> An integrated computational approach
        combining molecular docking and molecular dynamics (MD)
        simulations was employed to elucidate the inhibitory
        interactions between curcumin and VP26.</p>
        <p><bold>Results:</bold> Docking results revealed a favorable
        binding affinity of –6.53 kcal/mol, indicating a spontaneous and
        stable interaction predominantly stabilized by van der Waals
        forces and hydrogen bonding, particularly involving Arg136.
        Subsequent 100 ns MD simulations demonstrated that the
        VP26-curcumin complex maintained high structural stability, with
        consistent hydrogen bonding, short interatomic distances, and
        minimal deviation in the radius of gyration. Residue-specific
        flexibility analysis indicated localized increases in dynamics
        near the binding site, suggesting subtle conformational
        adaptation upon ligand binding. The MM/PBSA binding free energy
        (-93.46 kJ/mol) confirmed strong and stable complex
        formation.</p>
        <p><bold>Conclusion:</bold> Collectively, these findings provide
        atomistic insights into the binding mechanism of curcumin with
        VP26, supporting its potential as a natural antiviral inhibitor
        against WSSV and offering a foundation for the rational design
        of novel antiviral agents in aquaculture.</p>
        <p><bold>Keywords:</bold> White Spot Syndrome Virus, Curcumin,
        Molecular Docking, Molecular Dynamics Simulations</p></td>
      </tr>
    </tbody>
  </table>
</table-wrap>
<sec id="introduction">
  <title><bold>Introduction</bold></title>
  <p>The White Spot Syndrome Virus (WSSV) is a highly virulent pathogen
  and the etiological agent of White Spot Disease, a condition
  responsible for devastating economic losses in global shrimp
  aquaculture (1, 2). As a member of the Nimaviridae family, WSSV
  possesses one of the largest known viral genomes, a circular
  double-stranded DNA exceeding 300 kilobases that encodes for
  approximately 180 proteins (3-5). Among these, the structural proteins
  VP19, VP24, VP26, and VP28 are critically involved in viral assembly
  and host-cell interactions, making them promising targets for
  antiviral strategies (6).</p>
  <p>The search for effective therapeutics has increasingly turned
  towards natural products, which are a rich source of bioactive
  compounds with diverse pharmacological properties (7, 8). Curcumin,
  the primary curcuminoid derived from the rhizome of Curcuma longa L.
  (turmeric), is a prime example. This poly phenolic compound has
  garnered significant attention for its broad-spectrum biological
  activities, including potent anti-oxidant, anti-inflammatory, and
  anti-cancer effects (9, 10). Notably, its antiviral potential has been
  demonstrated against a range of viruses, suggesting a mechanism that
  may involve the disruption of viral entry or replication (11-13).
  Despite this promise, the precise molecular interactions between
  curcumin and specific viral targets remain largely unexplored.
  Experimental elucidation of these interactions can be resource
  intensive, highlighting the value of robust computational methods as a
  preliminary investigative tool (14). Integrated computational
  approaches, particularly molecular docking paired with molecular
  dynamics (MD) simulations, offer a powerful strategy to predict
  ligand-binding modes, characterize binding affinities, and evaluate
  the stability of resulting complexes under physiological conditions
  (15, 16). Molecular docking provides a static snapshot of potential
  binding poses, while MD simulations extend this insight by modeling
  the dynamic behavior and conformational stability of the
  protein-ligand complex over time. Given the critical role of the VP26
  protein in WSSV pathogens and the documented antiviral properties of
  curcumin, we hypothesize that curcumin may act as a potent inhibitor
  of VP26.</p>
  <p>To test this hypothesis, we employed an integrated in silico
  approach combining molecular docking and extensive MD simulations. Our
  objectives were to: 1) predict the binding affinity and precise
  binding pose of curcumin within the VP26 structure, 2) identify the
  key amino acid residues mediating this interaction, and 3) assess the
  stability and dynamic behavior of the curcumin-VP26 complex. The
  findings from this study provide crucial atomistic insights into the
  inhibitory mechanism of curcumin against WSSV, establishing a
  foundational framework for the future development of targeted
  antiviral agents in aquaculture.</p>
</sec>
<sec id="methods">
  <title><bold>Methods</bold></title>
  <sec id="protein-preparation">
    <title><italic><bold>Protein Preparation</bold></italic></title>
    <p>The three-dimensional structure of the (WSSV) VP26 protein was
    retrieved from the Protein Data Bank (PDB ID: 2EDM). Protein
    preparation was conducted using AutoDockTools 1.5.6 program (9). The
    preparation process involved removing all hetero atoms and
    crystallographic water molecules, adding hydrogen atoms, and
    assigning Gasteiger partial charges. The resulting structure was
    subsequently saved in the PDBQT format for docking analysis.</p>
  </sec>
  <sec id="section">
    <title></title>
  </sec>
  <sec id="ligand-preparation-and-binding-site-prediction">
    <title><italic><bold>Ligand Preparation and Binding Site
    Prediction</bold></italic></title>
    <p>The three-dimensional structure of curcumin (PubChem CID: 969516)
    was obtained from the PubChem database (17). The ligand geometry was
    optimized using the Gaussian software package (18) with the
    B3LYP/6-31G basis set to achieve a stable conformation.</p>
  </sec>
  <sec id="section-1">
    <title></title>
  </sec>
  <sec id="molecular-docking">
    <title><italic><bold>Molecular Docking</bold></italic></title>
    <p>Molecular docking studies were carried out to explore the binding
    affinity and interaction pattern between curcumin and the VP26
    protein using AutoDock 4.2 (8). A grid box with dimensions of 60 ×
    60 × 80 and a grid spacing of 0.375 Å was defined to encompass the
    predicted binding pocket. The docking procedure was performed with
    200 independent runs using the Lamarckian Genetic Algorithm (LGA).
    During docking, the protein was kept rigid, while all rotatable
    bonds in the ligand were allowed to move freely. The most favorable
    docking conformation, defined as the structure with the lowest
    binding energy in the most populated cluster, was selected for
    further analysis. Intermolecular interactions, including hydrogen
    bonds and van der Waals contacts, were visualized and analyzed using
    Ligplot+ program (19).</p>
  </sec>
  <sec id="section-2">
    <title></title>
  </sec>
  <sec id="molecular-dynamics-md-simulation">
    <title><italic><bold>Molecular Dynamics (MD)
    Simulation</bold></italic></title>
    <p>To evaluate the stability and dynamic behavior of the
    VP26-curcumin complex, MD simulations were performed using GROMACS
    2022.6 (20) with the AMBER99SB force field. Topology parameters for
    curcumin were generated using the ACPYPE Python interface (21). Two
    systems were simulated: the VP26–curcumin complex and the free VP26
    protein. Systems were solvated in a cubic periodic box filled with
    TIP3P water molecules, maintaining a minimum distance of 1.0 nm
    between the solute and the box boundaries. System neutrality was
    achieved by adding appropriate numbers of Na⁺ and Cl⁻ counterions,
    followed by adjustment to a physiological salt concentration of 0.15
    M. Energy minimization was performed using the steepest descent
    algorithm until the maximum force was reduced below 100 kJ/mol/nm.
    The minimized system was equilibrated in two stages. NVT
    equilibration (canonical ensemble) for 1 ns at 310 K using the
    Nose–Hoover thermostat. NPT equilibration (isothermal–isobaric
    ensemble) for 1 ns at 1 bar using the Parrinello–Rahman barostat. A
    production MD run was then performed for 100 ns with a 2-fs time
    step.</p>
  </sec>
</sec>
<sec id="results-and-discussion">
  <title><bold>Results and Discussion</bold></title>
  <p><italic><bold>Docking Analysis of the Curcumin-VP26
  Complex</bold></italic></p>
  <p>Molecular docking analysis using AutoDock predicted a stable
  binding mode for curcumin within a defined pocket of the VP26 protein.
  The calculated binding affinity for the most favorable pose was -6.53
  kcal/mol, corresponding to a predicted inhibition constant (Ki) of
  16.34 µM (Table 1). This energetically favorable score suggests a
  spontaneous and stable binding interaction under physiological
  conditions.</p>
  <p><bold>Table 1:</bold> Energetic components of the curcumin-VP26
  interactions obtained from molecular docking</p>
  <table-wrap>
    <table>
      <colgroup>
        <col width="19%" />
        <col width="18%" />
        <col width="25%" />
        <col width="18%" />
        <col width="21%" />
      </colgroup>
      <thead>
        <tr>
          <th><p><bold>Lowest binding energy</bold></p>
          <p><bold>(kcal/mol)</bold></p></th>
          <th><p><bold>Inhibition constant</bold></p>
          <p><bold>(Micro Molar)</bold></p></th>
          <th><p><bold>vdW + H- bond + dissolved</bold></p>
          <p><bold>energy (kcal/mol)</bold></p></th>
          <th><p><bold>electrostatic energy</bold></p>
          <p><bold>(kcal/mol)</bold></p></th>
          <th><p><bold>Torsional Free Energy</bold></p>
          <p><bold>(kcal/mol)</bold></p></th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <td>- 6.53</td>
          <td>16.34</td>
          <td>- 9.16</td>
          <td>- 0.36</td>
          <td>2.98</td>
        </tr>
      </tbody>
    </table>
  </table-wrap>
  <p>A detailed decomposition of the binding energy revealed that the
  interaction was primarily stabilized by van der Waals forces, hydrogen
  bonding, and desolvation effects, which collectively contributed -9.16
  kcal/mol. In contrast, the electrostatic contribution was relatively
  minor at -0.36 kcal/mol. The torsional free energy penalty, associated
  with the loss of conformational freedom upon binding, was +2.98
  kcal/mol.</p>
  <p>The binding orientation of curcumin, illustrated in Fig. 1, shows
  the ligand nestled within a hydrophobic cleft of the VP26 protein.
  Analysis of the specific molecular interactions identified several key
  amino acid residues forming van der Waals contacts with curcumin,
  including Tyr46, Asp47, Gln48, Met49, Pro72, Asp132, Thr134, Arg136,
  Asn186, and Gly188. A critical stabilizing feature of the complex is a
  conventional hydrogen bond formed between the guanidino group of
  Arg136 and the carbonyl oxygen atom of curcumin. This specific
  interaction, combined with the extensive network of non-bonded
  contacts, underscores the binding of curcumin to the VP26 binding site
  and provides a structural rationale for its potential inhibitory
  activity.</p>
  <graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image4.jpeg" />
  <p><bold>Fig. 1:</bold> Predicted binding pose of curcumin in the
  binding site of VP26 protein. The three-dimensional binding pose of
  curcumin within the surface representation of the VP26 binding pocket
  and two-dimensional ligand-protein interaction diagram generated by
  Ligplot<sup>+</sup>, detailing the key intermolecular contacts</p>
  <p>The binding affinity observed for curcumin against VP26 is
  comparable to other reports of curcumin’s interactions with viral and
  protein targets. For example, in silico studies targeting viral
  proteases and coat proteins of plant viruses reported curcumin docking
  energies in the range of −6.3 to −7.2 kcal/mol, indicating similar
  moderate binding strengths to viral structural proteins in silico
  (−7.2 kcal/mol for PVY coat protein) (22). Similar docking studies
  against human disease-related proteins also report curcumin binding
  energies in a range that supports its potential as a modest affinity
  inhibitor (e.g., docking to SARS-CoV-2 proteins and CDK2 in cancer
  studies) (23). These comparative studies suggest that curcumin tends
  to show consistent moderate binding affinities across diverse protein
  targets, reinforcing the hypothesis that curcumin can interact
  reliably with protein pockets through hydrophobic and hydrogen-bonding
  interactions. The dominant contribution of van der Waals and hydrogen
  bonding observed in our analysis is in agreement with these previous
  reports where hydrophobic contacts play central roles in stabilizing
  complexes of curcumin with viral domains and enzymes (24).</p>
  <sec id="molecular-dynamics-simulation-analysis">
    <title><italic><bold>Molecular Dynamics Simulation
    Analysis</bold></italic></title>
    <p>To validate the stability of the docked curcumin-VP26 complex and
    gain deeper insights into the molecular dynamics and structural
    consequences of binding, we performed 100 ns MD simulations for both
    the ligand-bound (VP26-Curcumin) and unbound (Free VP26) states.</p>
    <sec id="section-3">
      <title></title>
    </sec>
    <sec id="structural-and-dynamical-analysis">
      <title><bold>Structural and Dynamical Analysis</bold></title>
      <p>The structural integrity and global dynamics of VP26 in its
      free and curcumin-bound states were rigorously assessed by
      calculating the Root Mean Square Deviation (RMSD), Root Mean
      Square Fluctuation (RMSF), Radius of gyration (Rg), secondary
      structure composition, and through Principal Component Analysis
      (PCA). These metrics collectively provide a comprehensive view of
      the protein's stability, flexibility, and conformational
      landscape. The backbone RMSD analysis, depicted in Fig. 2,
      revealed the dynamic stability of both systems throughout the
      simulation. After an initial equilibration period, the Free VP26
      structure stabilized, maintaining an average RMSD of 0.446 ± 0.039
      nm.</p>
      <p><inline-graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image5.jpeg" /><bold>Fig.
      2:</bold> Backbone RMSD of trajectories of Free VP26 and the
      VP26-Curcumin complex over the 100 ns MD simulation</p>
      <p>The preservation of global compactness, as evidenced by nearly
      unchanged radius of gyration (Rg) values between Free VP26 and the
      VP26-Curcumin complex, further supports that curcumin binding does
      not disrupt the native fold of VP26. Similar result has been
      reported in MD studies of curcumin-related ligand complexes, where
      stable Rg profiles were observed for both unbound and ligand-bound
      forms throughout 100 ns trajectories, indicating that ligand
      engagement does not lead to major conformational expansion or
      unfolding but rather maintains overall structural integrity
      (25).</p>
      <p>A key observation emerges from the VP26-Curcumin complex, which
      settled into a distinct, stable state with a consistently higher
      average RMSD of 0.479 ± 0.033 nm (Table 2). This elevation in RMSD
      does not signify instability or unfolding, but rather indicates a
      subtle ligand-induced structural perturbation. Curcumin binding
      has shifted the protein's energy landscape, stabilizing a new
      conformational state that is slightly different from its native,
      unbound form. The well-defined plateau and low standard deviation
      in the complex's trajectory are critical, as they confirm that
      this altered conformation is itself stable and persistent. This
      represents the formation of a specific and well-defined
      ligand-protein complex, where VP26 has undergone a minor but
      stable structural adaptation to accommodate curcumin. This
      behavior is consistent with previous MD analyses demonstrating
      that ligand engagement can shift a protein’s energy landscape into
      a new conformational basin without compromising overall stability,
      leading to an alternative but stable structural sub-state sampled
      during the trajectory (25).</p>
      <p><bold>Table 2:</bold> Comparative summary of MD simulation
      parameters for the Free VP26 and VP26-Curcumin complex. Values
      represent the mean ± standard deviation calculated over the stable
      trajectory (last 20 ns)</p>
      <table-wrap>
        <table>
          <colgroup>
            <col width="24%" />
            <col width="26%" />
            <col width="25%" />
            <col width="25%" />
          </colgroup>
          <thead>
            <tr>
              <th><bold>System</bold></th>
              <th><bold>Mean RMSD (nm)</bold></th>
              <th><bold>Mean RMSF (nm)</bold></th>
              <th><bold>Mean Rg (nm)</bold></th>
            </tr>
          </thead>
          <tbody>
            <tr>
              <td>Free Vp26</td>
              <td>0.446±0.039</td>
              <td>0.139±0.088</td>
              <td>1.697±0.016</td>
            </tr>
            <tr>
              <td>Vp26 – Curcumin</td>
              <td>0.479±0.033</td>
              <td>0.143±0.089</td>
              <td>1.671±0.016</td>
            </tr>
          </tbody>
        </table>
      </table-wrap>
      <p>To dissect the local flexibility changes upon binding, we
      analyzed the per-residue RMSF (Fig. 3). While the overall mean
      RMSF values were similar (0.139 ± 0.088 nm for Free VP26 vs. 0.143
      ± 0.089 nm for the complex), a residue-specific analysis reveals
      critical, localized effects that provide deep insight into the
      binding mechanism.</p>
      <graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image6.jpeg" />
      <p><bold>Fig. 3:</bold> Per-residue RMSF of Free VP26 and the
      VP26-Curcumin complex</p>
      <p>Contrary to the typical expectation of rigidification at a
      binding site, we observed a significant increase in flexibility
      for several key residues directly involved in the interaction with
      curcumin. Most notably, Arg136, which forms a critical hydrogen
      bond with the ligand, shows a marked increase in fluctuation
      (0.170 nm in Free VP26 vs. 0.189 nm in the complex). The hydrogen
      bond with Arg136, while stable, appears to be dynamic in nature,
      pulling the residue into a new conformational sub-state that
      exhibits greater mobility than in the unbound protein.</p>
      <p>This localized increase in flexibility was not isolated. A
      pronounced enhancement in dynamics was also observed in the
      adjacent loop region spanning residues 60-75<bold>,</bold> which
      contained the binding site residue Pro72. Furthermore, the
      N-terminal region (residues 25<bold>-</bold>30) also becomes more
      flexible. These data paint a picture of a binding event that
      induces strain or dynamic allostery in its immediate vicinity.</p>
      <p>The potential impact of this binding induced dynamic adaptation
      on the global compactness of VP26 was assessed by monitoring the
      Rg, as shown in Fig. 4.</p>
      <p><inline-graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image7.jpeg" /><bold>Fig.
      4:</bold> Rg plots for Free VP26 and the VP26-Curcumin complex</p>
      <p>Such localized increases in fluctuation at or near ligand
      interaction sites are consistent with previous reports showing
      that ligand binding does not universally rigidify a protein but
      can instead induce region-specific changes in flexibility. For
      example, MD studies comparing ligand-bound and unbound protein
      structures have demonstrated that binding can increase mobility in
      certain loops or side chains adjacent to the binding site,
      reflecting a dynamic adaptation of the protein to accommodate the
      ligand rather than a simple rigidification of the binding
      interface (26). Furthermore, models of ligand-protein dynamics
      have shown that ligand engagement often leads to compensatory
      flexibility changes whereby regions directly interacting with the
      ligand become more dynamic while other regions may become less so,
      illustrating the complex and sometimes counterintuitive nature of
      protein flexibility upon binding (27). These observations support
      our interpretation that curcumin binding induces localized dynamic
      adaptation in VP26, particularly in loop regions and binding site
      residues, which may be part of the mechanism by which the protein
      accommodates ligand binding without global destabilization.</p>
      <p>The average Rg for the Free VP26 system was calculated to be
      1.697 ± 0.016 nm, while the VP26-Curcumin complex exhibited a very
      similar average Rg of 1.671 ± 0.016 nm. This minor difference
      suggests that the structural perturbation and localized
      flexibility changes do not translate into a major conformational
      expansion or compaction of the overall protein fold. The stability
      of the complex is further underscored by the identical and minimal
      standard deviation of the Rg, indicating a consistent and
      well-packed global structure throughout the simulation, despite
      the increased fluctuations at specific sites. The nearly identical
      Rg values observed for the Free VP26 and VP26–Curcumin systems
      indicate that ligand binding does not induce significant global
      expansion or compaction of the protein. In molecular dynamics
      simulations, Rg is a direct measure of protein compactness, and
      stable Rg values over time are interpreted as preservation of the
      overall tertiary structure despite local conformational
      fluctuations (28).</p>
      <p>To assess whether curcumin binding induced significant
      alterations in the protein's secondary structure, we analyzed the
      average percentage of each secondary structure element throughout
      the simulation (Fig. 5).</p>
      <graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image8.jpeg" />
      <p><bold>Fig. 5:</bold> Average percentage of the secondary
      structure composition of Free VP26 and the VP26-Curcumin complex
      throughout the last 20 ns MD simulation</p>
      <p>The overall composition was highly conserved between the Free
      VP26 and the VP26-Curcumin complex. B-Sheet (54.3% vs. 55.6%) and
      Coil (22.1% vs. 22.9%) remained the dominant structural elements
      in the free and bound states, respectively, indicating that no
      large-scale unfolding or refolding occurred. Minor fluctuations
      were observed in other elements: A-Helix content increased from
      2.4% to 3.0%, B-Bridge from 0.7% to 0.8%, and turn from 6.4% to
      6.5%. Conversely, a decrease was noted in Bend content, from 13.7%
      to 11.0%. Given the subtle nature of these changes, they do not
      point to a drastic structural reorganization. Instead, they likely
      reflect localized adjustments and a slight stabilization of
      certain elements, particularly B-sheets, upon ligand binding,
      which may contribute to the overall stability of the complex
      without fundamentally altering the protein's core
      architecture.</p>
      <p>To understand the effect of curcumin binding on the large-scale
      motions of the VP26, we performed Principal Component Analysis.
      The analysis of the eigenvalues showed that the first two
      principal eigenvectors (PC1 and PC2) represent the most dominant
      collective motions in both systems (Figure 6a). The projection of
      the simulation trajectories onto these two eigenvectors revealed
      the conformational space sampled by the protein (Fig. 6b and
      c).</p>
      <graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image9.jpeg" />
      <p><bold>Fig. 6:</bold> Principal Component Analysis (PCA) of the
      MD trajectories. (a) Eigenvalue of the first 10 eigenvectors.
      Projection of the trajectory frames onto the phase space defined
      by the first two principal eigenvectors for (b) Free VP26 and (c)
      VP26-Curcumin complex</p>
      <p>The Free VP26 system was characterized by a relatively broad
      distribution, indicating its inherent flexibility. The
      conformational ensemble of the VP26-Curcumin complex also samples
      a broad area, but with a subtle yet consistent shift and a slight
      narrowing of its density core compared to the free protein. This
      minor change in the conformational sampling suggests that the
      binding of curcumin induces a modest redistribution of the
      protein's essential motions, stabilizing a specific conformational
      sub-ensemble. While the effect is not dramatic, this subtle
      modulation of global dynamics, when combined with the local
      changes observed in RMSF, may contribute to the mechanism of
      curcumin's interaction with VP26. Similar PCA-based MD analyses in
      the literature have been used to differentiate functional states
      of proteins, identify metastable substates, and compare
      conformational ensembles between unbound and ligand-bound forms,
      highlighting how ligand binding can alter the accessible
      conformational landscape while preserving overall structural
      integrity (29). These observations support our interpretation that
      curcumin binding subtly modulates collective motions of VP26,
      redistributing the conformational sampling without introducing
      large-scale destabilization.</p>
    </sec>
  </sec>
  <sec id="binding-stability-and-interaction-analysis">
    <title><italic><bold>Binding Stability and Interaction
    Analysis</bold></italic></title>
    <p>To gain deeper insights into the stability of the curcumin-VP26
    complex and the specific molecular interactions that govern binding,
    we analyzed the contact frequency, minimum distance, hydrogen
    bonding, and binding free energy throughout the MD simulation
    trajectory.</p>
    <p>The stability of the curcumin-VP26 complex, initially predicted
    by molecular docking, was rigorously assessed during the 100 ns MD
    simulation. The number of atomic contacts between curcumin and VP26
    residues provides a measure of the extent of the interaction
    interface. As shown in Fig. 7a, the complex maintained a high and
    stable number of contacts throughout the simulation.</p>
    <graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image10.jpeg" />
    <p><bold>Fig. 7:</bold> Analysis of the atomic interactions and
    proximity in the curcumin-VP26 complex. <bold>(a)</bold> The
    time-evolution of the number of atomic contacts (cutoff distance 0.6
    nm) between curcumin and VP26 residues throughout the 100 ns
    molecular dynamics simulation <bold>(b)</bold> The minimum distance
    between the curcumin molecule and the VP26 binding site over the
    simulation time</p>
    <p>Quantitatively, the system exhibited an average of 1548 contacts,
    with a maximum of 2183 and a minimum of 215. The consistently high
    average and maximum values indicate a robust and dense network of
    interactions, suggesting the ligand is extensively engaged with the
    binding site residues. While the transient dips to a minimum of 215
    contacts reflect the dynamic nature of the binding interface, where
    side chains sample different conformations, the rapid recovery to a
    high number of contacts underscores the resilience and enduring
    nature of the association. Such analyses of contact persistence as a
    marker of binding stability are commonly used to validate docking
    poses in MD studies, where a stable interaction network correlates
    with maintained complex integrity over time (30).</p>
    <p>Complementing this, the minimum distance between the ligand and
    the protein binding pocket was analyzed as a critical metric for
    binding stability. As shown in Fig. 7b, the minimum distance between
    curcumin and the VP26 binding site remained consistently low
    throughout the simulation. The average minimum distance was
    calculated to be 0.21 nm, fluctuating within a very narrow range
    between 0.14 nm (min) and 0.28 nm (max). This consistently short
    proximity, far below a typical van der Waals cutoff, indicates that
    curcumin remained intimately associated and firmly positioned within
    its binding pocket without any significant dissociation events. The
    convergence of a high number of atomic contacts with a short and
    stable minimum distance provides strong, mutually reinforcing
    evidence for a tight and stable binding mode.</p>
    <p>Hydrogen bonding plays a pivotal role in conferring specificity
    and strength to protein-ligand interactions. The analysis of
    hydrogen bond formation (Fig. 8) throughout the simulation
    trajectory revealed that the curcumin-VP26 complex maintained
    persistent hydrogen bonding.</p>
    <graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image11.jpeg" />
    <p><bold>Fig. 8:</bold> Time-dependent evolution of the number of
    hydrogen bonds between curcumin and the VP26 protein during the 100
    ns molecular dynamics simulation</p>
    <p>The stability of the binding pose predicted by molecular docking
    was visually confirmed through a time evolution analysis of the
    simulation snapshots. As depicted in Fig. 9, representative
    structures extracted at 25 ns intervals (0, 25, 50, 75, and 100 ns)
    demonstrate that curcumin remains securely bound within the same
    hydrophobic pocket of VP26. The ligand's orientation and position
    are remarkably consistent from the initial frame to the final frame,
    with no significant translation or rotation observed. This visual
    evidence directly corroborates the quantitative data from the
    contact, and distance analyses, confirming that the initially docked
    pose is not an artifact but represents a genuine, stable binding
    mode.</p>
    <graphic mimetype="image" mime-subtype="jpeg" xlink:href="vertopal_8cad03f31d694698b241ff2b2222ba84/media/image12.jpeg" />
    <p><bold>Fig. 9:</bold> Structural snapshots of the curcumin-VP26
    complex at key time points (0, 25, 50, 75, and 100 ns) from the 100
    ns molecular dynamics simulation. VP26 is shown in cartoon
    representation (cyan), curcumin is depicted as red ball and
    sticks</p>
    <p>To quantitatively evaluate the binding affinity of curcumin for
    VP26, the Molecular Mechanics/Poisson-Boltzmann Surface Area
    (MM/PBSA) method was employed. The calculated binding free energies
    from the stable simulation trajectory are summarized in Table
    3<bold>.</bold></p>
    <p><bold>Table 3:</bold> MM/PBSA binding free energy components for
    the curcumin-VP26 complex (kJ/mol)</p>
    <table-wrap>
      <table>
        <colgroup>
          <col width="49%" />
          <col width="51%" />
        </colgroup>
        <thead>
          <tr>
            <th>Energy Component</th>
            <th>Mean ± Standard Deviation</th>
          </tr>
        </thead>
        <tbody>
          <tr>
            <td>van der Waals Energy</td>
            <td>-128.253 ± 27.688</td>
          </tr>
          <tr>
            <td>Electrostatic Energy</td>
            <td>-7.684 ± 10.672</td>
          </tr>
          <tr>
            <td>Polar Solvation Energy</td>
            <td>57.371 ± 17.188</td>
          </tr>
          <tr>
            <td>SASA Energy</td>
            <td>-14.891 ± 2.503</td>
          </tr>
          <tr>
            <td>Total Binding Energy</td>
            <td>-93.457 ± 23.540</td>
          </tr>
        </tbody>
      </table>
    </table-wrap>
    <p>The MM/PBSA results reveal a strongly favorable binding free
    energy of -93.46 kJ/mol for the curcumin-VP26 complex. Decomposition
    of the total energy into its components shows that the binding is
    primarily driven by a large, favorable van der Waals contribution
    (-128.25 kJ/mol), supplemented by a favorable non-polar solvation
    energy (-14.89 kJ/mol). This indicates that hydrophobic interactions
    and the close shape complementarity between curcumin and the VP26
    binding pocket are the dominant forces stabilizing the complex.
    While the electrostatic energy is favorable, its effect is largely
    counteracted by the unfavorable polar solvation penalty upon
    binding. This is a common phenomenon in biomolecular recognition,
    where the desolvation of charged/polar groups is energetically
    costly.</p>
  </sec>
  <sec id="conclusion">
    <title><bold>Conclusion</bold></title>
    <p>This study provides a molecular-level perspective on the
    potential interaction between curcumin and the VP26 protein of White
    Spot Syndrome Virus, highlighting curcumin’s capacity to associate
    stably with a functionally relevant region of the protein. The
    integrated computational framework suggests that such binding may
    induce localized conformational and dynamic constraints that could
    plausibly influence VP26-mediated processes during viral infection.
    Importantly, these observations should be interpreted as mechanistic
    hypotheses rather than direct evidence of antiviral activity, as
    they are derived exclusively from in silico analyses. Nevertheless,
    the findings contribute to a conceptual understanding of how natural
    polyphenolic compounds might target structural viral proteins and
    offer a rational starting point for subsequent experimental
    validation, structure-based optimization, and biological assessment
    in the context of shrimp aquaculture disease management.</p>
    <p><bold>Acknowledgments</bold></p>
    <p>The authors extend their heartfelt appreciation to the Khatam
    Al-Nabieen University Medical Research and Technology Center for its
    supportive scientific environment, institutional backing, and access
    to essential research resources that made this narrative review
    possible.</p>
    <p><bold>Founding</bold></p>
    <p>None.</p>
    <p><bold>Conflict of interest</bold></p>
    <p>The authors declare that there is no conflict of interests.</p>
    <p><bold>Reference</bold></p>
    <list list-type="order">
      <list-item>
        <p>Narasimman V, Ramachandran S. In silico analysis of low
        molecular weight sulfated chitosan from Sepia brevimana as
        potential inhibitors of white spot syndrome envelope proteins.
        Biomass Conversion and Biorefinery. 2024;14(16):18539-50.</p>
      </list-item>
      <list-item>
        <p>Wu C, Yang F. Localization studies of two white spot syndrome
        virus structural proteins VP51 and VP76. Virol J.
        2006;3(1):76.</p>
      </list-item>
      <list-item>
        <p>Rozenberg A, Brand P, Rivera N, Leese F, Schubart CD.
        Characterization of fossilized relatives of the White Spot
        Syndrome Virus in genomes of decapod crustaceans. BMC
        Evolutionary Biol. 2015;15(1):142.</p>
      </list-item>
      <list-item>
        <p>Cao X-T, Wu L-J, Xu F-L, Li X-C, Lan J-F. Pc Trim prevents
        early infection with white spot syndrome virus by inhibiting
        AP1-induced endocytosis. Cell Communication and Signaling.
        2023;21(1):104.</p>
      </list-item>
      <list-item>
        <p>Seok SH, Park JH, Cho SA, Baek MW, Lee HY, Kim DJ, et al.
        Cloning and sequencing of envelope proteins (VP19, VP28) and
        nucleocapsid proteins (VP15, VP35) of a white spot syndrome
        virus isolate from Korean shrimp. Dis Aquat Organ.
        2004;60(1):85-8.</p>
      </list-item>
      <list-item>
        <p>Osaka R, Kobayashi N, Shimada K, Ishii A, Oka N, Kondo K.
        VP26, a herpes simplex virus type 1 capsid protein, increases
        DNA methylation in COASY promoter region. Brain, Behavior, &amp;
        Immunity-Health. 2022;26:100545.</p>
      </list-item>
      <list-item>
        <p>Xie X, Yang F. White spot syndrome virus VP24 interacts with
        VP28 and is involved in virus infection. J Gen Virol.
        2006;87(7):1903-8.</p>
      </list-item>
      <list-item>
        <p>Dinesh S, Sudharsana S, Mohanapriya A, Itami T, Sudhakaran R.
        Molecular docking and simulation studies of Phyllanthus amarus
        phytocompounds against structural and nucleocapsid proteins of
        white spot syndrome virus. Biotech. 2017;7(5):353.</p>
      </list-item>
      <list-item>
        <p>Molla MHR, Aljahdali MO. Identification of phytochemical
        compounds to inhibit the matrix-like linker protein VP26 to
        block the assembles of white spot syndrome virus (WSSV) envelope
        and nucleocapsid protein of marine shrimp: In silico approach.
        Journal of King Saud University-Science. 2022;34(8):102346.</p>
      </list-item>
      <list-item>
        <p>Thomford NE, Senthebane DA, Rowe A, Munro D, Seele P, Maroyi
        A, et al. Natural products for drug discovery in the 21st
        century: innovations for novel drug discovery. Int J Mol Sci.
        2018;19(6):1578.</p>
      </list-item>
      <list-item>
        <p>Harvey AL. Natural products in drug discovery. Drug Discovery
        Today. 2008;13(19-20):894-901.</p>
      </list-item>
      <list-item>
        <p>Rezaei M, Soltani M, Alipoor E, Rezayat SM,
        Vasheghani-Farahani A, Yaseri M, et al. Effect of nano-curcumin
        supplementation on angina status, and traditional and novel
        cardiovascular risk factors in overweight or obese patients with
        coronary slow flow phenomenon: a randomized double-blind
        placebo-controlled clinical trial. BMC Nutr. 2024;10(1):73.</p>
      </list-item>
      <list-item>
        <p>Koohpeyma F, Khodaparast Z, Salehi S, Danesh S, Gheshlagh FM,
        Naseri A, et al. The ameliorative effects of curcumin
        nanomicelle on testicular damage in the mouse model of multiple
        sclerosis. BMC Complement Med Ther. 2024;24(1):200.</p>
      </list-item>
      <list-item>
        <p>Praditya D, Kirchhoff L, Brüning J, Rachmawati H, Steinmann
        J, Steinmann E. Anti-infective properties of the golden spice
        curcumin. Front Microbiol. 2019;10:912.</p>
      </list-item>
      <list-item>
        <p>Zhao P, Feng L, Jiang W, Wu P, Liu Y, Ren H, et al. Unveiling
        the emerging role of curcumin to alleviate ochratoxin A-induced
        muscle toxicity in grass carp (Ctenopharyngodon idella): in
        vitro and in vivo studies. J Anim Sci Biotechnol.
        2024;15(1):72.</p>
      </list-item>
      <list-item>
        <p>Nicoliche T, Bartolomeo CS, Lemes RMR, Pereira GC, Nunes TA,
        Oliveira RB, et al. Antiviral, anti-inflammatory and antioxidant
        effects of curcumin and curcuminoids in SH-SY5Y cells infected
        by SARS-CoV-2. Sci Rep. 2024;14(1):10696.</p>
      </list-item>
      <list-item>
        <p>Alkhathami AG, Alshahrani MY, Alshehri SA, Nasir N, Wahab S.
        Curcumin as a potential inhibitor of TGFβ3 computational
        insights for breast cancer therapy. Sci Rep.
        2025;15(1):2871.</p>
      </list-item>
      <list-item>
        <p>Soriano-Díaz I, Ortí E, Giussani A. On the Importance of
        Ligand-Centered Excited States in the Emission of Cyclometalated
        Ir(III) Complexes. Inorg Chem. 2021;60(17):13222-32.</p>
      </list-item>
      <list-item>
        <p>Patil R, Das S, Stanley A, Yadav L, Sudhakar A, Varma AK.
        Optimized hydrophobic interactions and hydrogen bonding at the
        target-ligand interface leads the pathways of drug-designing.
        PLoS One. 2010;5(8):e12029.</p>
      </list-item>
      <list-item>
        <p>Dondapati¹ RS, Agarwal R. Overview of GROMACS. Molecular
        Dynamics Simulation of Nanocomposites using BIOVIA Materials
        Studio, Lammps and Gromacs. 2025:76.</p>
      </list-item>
      <list-item>
        <p>Sousa da Silva AW, Vranken WF. ACPYPE - AnteChamber PYthon
        Parser interfacE. BMC Res Notes. 2012;5(1):367.</p>
      </list-item>
      <list-item>
        <p>Megahed AA, Barakat A, Embaby AA, Mohamed SH, Sadik AS,
        Elshaer NA. Antiviral activity of turmeric (Curcuma longa)
        against potato virus Y: in silico molecular docking analysis.
        Beni-Suef University J Basic Appl Sci. 2025;14(1):97.</p>
      </list-item>
      <list-item>
        <p>Suravajhala R, Parashar A, Choudhir G, Kumar A, Malik B,
        Nagaraj VA, et al. Molecular docking and dynamics studies of
        curcumin with COVID-19 proteins. Netw Model Anal Health Inform
        Bioinform. 2021;10(1):44.</p>
      </list-item>
      <list-item>
        <p>Srivastava AK, Kumar V, Roy BK. Insights from the molecular
        docking of curcumin to the virulent factors of Helicobacter
        pylori. Bioinformation. 2015;11(10):447-53.</p>
      </list-item>
      <list-item>
        <p>Rampogu S, Lee G, Park JS, Lee KW, Kim MO. Molecular Docking
        and Molecular Dynamics Simulations Discover Curcumin Analogue as
        a Plausible Dual Inhibitor for SARS-CoV-2. Int J Mol Sci.
        2022;23(3).</p>
      </list-item>
      <list-item>
        <p>Erman B. Effects of ligand binding upon flexibility of
        proteins. Proteins. 2015;83(5):805-8.</p>
      </list-item>
      <list-item>
        <p>Wankowicz SA, de Oliveira SH, Hogan DW, van den Bedem H,
        Fraser JS. Ligand binding remodels protein side-chain
        conformational heterogeneity. Elife. 2022;11.</p>
      </list-item>
      <list-item>
        <p>Hollenbeck JJ, McClain DL, Oakley MG. The role of helix
        stabilizing residues in GCN4 basic region folding and DNA
        binding. Protein Sci. 2002;11(11):2740-7.</p>
      </list-item>
      <list-item>
        <p>Yasuda T, Shigeta Y, Harada R. Efficient Conformational
        Sampling of Collective Motions of Proteins with Principal
        Component Analysis-Based Parallel Cascade Selection Molecular
        Dynamics. J Chem Inf Model. 2020;60(8):4021-9.</p>
      </list-item>
      <list-item>
        <p>Liu K, Watanabe E, Kokubo H. Exploring the stability of
        ligand binding modes to proteins by molecular dynamics
        simulations. J Comput Aided Mol Des. 2017;31(2):201-11.</p>
      </list-item>
    </list>
  </sec>
</sec>
</body>
<back>
</back>
</article>
