Surveillance of Adenovirus and Norovirus Contaminants in the Water and Shellfish of Major Oyster Breeding Farms and Fishing Ports in Taiwan

The enteric viruses, including adenovirus (AdVs) and norovirus (NoVs), in shellfish is a significant food safety risk. This study investigated the prevalence, seasonal occurrence, genetic diversity, and quantification of AdVs and NoVs in the water and cultured shellfish samples at the four major coastal oyster breeding farms (COBF), five major fishing ports (FP), and their markets in Taiwan. The AdVs/NoVs in the water and shellfish samples were isolated by the membrane filtration and direct elution methods. The RNA of NoVs was reverse-transcribed into complementary DNA through reverse transcription reaction. Further NoVs and AdVs were detected using nested PCR. A higher detection rate was recorded in the low-temperature period than high-temperature. Detection difference was noted between nested PCR and qPCR outcomes for AdVs. The total detection rate of AdVs was higher in the water samples (COBF-40.6%, FP 20%) than the shellfish samples (COBF-11.7% and FP 6.3%). The AdVs load in the water and shellfish samples ranged from 1.23 × 103 to 1.00 × 106 copies/L and 3.57 × 103 to 4.27 × 104 copies/100g, respectively. The total detection of NoVs was highest in the water samples of the FP and their market shellfish samples (11.1% and 3.2%, respectively). Genotyping and phylogenetic analysis were identified as the prevalent AdVs and NoVs genotypes in the water and shellfish samples: A species HAdVs serotype 12; F species HAdVs serotype 41; and C species PAdVs serotype 5 (NoVs GI.2, GI.3 and GII.2). No significant differences were observed between the presence of AdVs, and all of the water quality parameters evaluated (heterotrophic plate count, water temperature, turbidity, pH, salinity, and dissolved oxygen). The virus contamination occurs mainly due to the direct discharge of domestic sewage, livestock farm, and fishing market wastewater into the coastal environment. Thus, this study suggested framing better estuarine management to prevent AdVs/NoVs transmission in water and cultured/distributed shellfish

Combinatorial Virtual Library Screening Study of Transforming Growth Factor-β2–Chondroitin Sulfate System

Transforming growth factor-beta (TGF-β), a member of the TGF-β cytokine superfamily, is known to bind to sulfated glycosaminoglycans (GAGs), but the nature of this interaction remains unclear. In a recent study, we found that preterm human milk TGF-β2 is sequestered by chondroitin sulfate (CS) in its proteoglycan form. To understand the molecular basis of the TGF-β2–CS interaction, we utilized the computational combinatorial virtual library screening (CVLS) approach in tandem with molecular dynamics (MD) simulations. All possible CS oligosaccharides were generated in a combinatorial manner to give 24 di- (CS02), 192 tetra- (CS04), and 1536 hexa- (CS06) saccharides. This library of 1752 CS oligosaccharides was first screened against TGF-β2 using the dual filter CVLS algorithm in which the GOLDScore and root-mean-square-difference (RMSD) between the best bound poses were used as surrogate markers for in silico affinity and in silico specificity. CVLS predicted that both the chain length and level of sulfation are critical for the high affinity and high specificity recognition of TGF-β2. Interestingly, CVLS led to identification of two distinct sites of GAG binding on TGF-β2. CVLS also deduced the preferred composition of the high specificity hexasaccharides, which were further assessed in all-atom explicit solvent MD simulations. The MD results confirmed that both sites of binding form stable GAG–protein complexes. More specifically, the highly selective CS chains were found to engage the TGF-β2 monomer with high affinity. Overall, this work present key principles of recognition with regard to the TGF-β2–CS system. In the process, it led to the generation of the in silico library of all possible CS oligosaccharides, which can be used for advanced studies on other protein–CS systems. Finally, the study led to the identification of unique CS sequences that are predicted to selectively recognize TGF-β2 and may out-compete common natural CS biopolymers

In-Depth Molecular Dynamics Study of All Possible Chondroitin Sulfate Disaccharides Reveals Key Insight into Structural Heterogeneity and Dynamism

GAGs exhibit a high level of conformational and configurational diversity, which remains untapped in terms of the recognition and modulation of proteins. Although GAGs are suggested to bind to more than 800 biologically important proteins, very few therapeutics have been designed or discovered so far. A key challenge is the inability to identify, understand and predict distinct topologies accessed by GAGs, which may help design novel protein-binding GAG sequences. Recent studies on chondroitin sulfate (CS), a key member of the GAG family, pinpointing its role in multiple biological functions led us to study the conformational dynamism of CS building blocks using molecular dynamics (MD). In the present study, we used the all-atom GLYCAM06 force field for the first time to explore the conformational space of all possible CS building blocks. Each of the 16 disaccharides was solvated in a TIP3P water box with an appropriate number of counter ions followed by equilibration and a production run. We analyzed the MD trajectories for torsional space, inter- and intra-molecular H-bonding, bridging water, conformational spread and energy landscapes. An in-house phi and psi probability density analysis showed that 1→3-linked sequences were more flexible than 1→4-linked sequences. More specifically, phi and psi regions for 1→4-linked sequences were held within a narrower range because of intra-molecular H-bonding between the GalNAc O5 atom and GlcA O3 atom, irrespective of sulfation pattern. In contrast, no such intra-molecular interaction arose for 1→3-linked sequences. Further, the stability of 1→4-linked sequences also arose from inter-molecular interactions involving bridged water molecules. The energy landscape for both classes of CS disaccharides demonstrated increased ruggedness as the level of sulfation increased. The results show that CS building blocks present distinct conformational dynamism that offers the high possibility of unique electrostatic surfaces for protein recognition. The fundamental results presented here will support the development of algorithms that help to design longer CS chains for protein recognition

A molecular dynamics-based algorithm for evaluating the glycosaminoglycan mimicking potential of synthetic, homogenous, sulfated small molecules

<div><p>Glycosaminoglycans (GAGs) are key natural biopolymers that exhibit a range of biological functions including growth and differentiation. Despite this multiplicity of function, natural GAG sequences have not yielded drugs because of problems of heterogeneity and synthesis. Recently, several homogenous non-saccharide glycosaminoglycan mimetics (NSGMs) have been reported as agents displaying major therapeutic promise. Yet, it remains unclear whether sulfated NSGMs structurally mimic sulfated GAGs. To address this, we developed a three-step molecular dynamics (MD)-based algorithm to compare sulfated NSGMs with GAGs. In the first step of this algorithm, parameters related to the range of conformations sampled by the two highly sulfated molecules as free entities in water were compared. The second step compared identity of binding site geometries and the final step evaluated comparable dynamics and interactions in the protein-bound state. Using a test case of interactions with fibroblast growth factor-related proteins, we show that this three-step algorithm effectively predicts the GAG structure mimicking property of NSGMs. Specifically, we show that two unique dimeric NSGMs mimic hexameric GAG sequences in the protein-bound state. In contrast, closely related monomeric and trimeric NSGMs do not mimic GAG in either the free or bound states. These results correspond well with the functional properties of NSGMs. The results show for the first time that appropriately designed sulfated NSGMs can be good structural mimetics of GAGs and the incorporation of a MD-based strategy at the NSGM library screening stage can identify promising mimetics of targeted GAG sequences.</p></div

The consistency of intermolecular hydrogen bonds across the MD simulation (final 10 ns) for FGF2–FGFR1 residues.

<p>The occurrence of inter-molecular hydrogen bonds between FGF2–FGFR1 residues (shown on y-axis) and small molecule ligands (NSGMs and HS06) are shown for each frame of the final 10 ns of the MD run. A) FGF2–FGFR1–G2.1 complex; B) FGF2–FGFR1–G2.2 complex; C) FGF2–FGFR1–HSO6(²S_O) complex; D) FGF2–FGFR1–HS06(¹C₄) complex.</p

The consistency of intermolecular hydrogen bonds across the MD simulation (final 10 ns) for FGF2 residues.

<p>The occurrence of inter-molecular hydrogen bonds between FGF2 residues (shown on y-axis) and small molecule ligands (NSGMs and HS06) are shown for each frame of the final 10 ns of the MD run. A) FGF2–G2.1 complex; B) FGF2–G2.2 complex; C) FGF2–HSO6(²S_O) complex; D) FGF2–HS06(¹C₄) complex.</p

Comparison of NSGMs and HS06 in the protein bound state using overall average inter-molecular hydrogen bond occupancy (A) and total binding free energy (B).

<p>A) shows the stability of bound NSGMs and HS06 based on average residue-level, inter-molecular hydrogen bond occupancy; B) shows the stability of the bound NSGMs and HS06 based on total binding free energy ΔG (in simulated kcal/mol, error bars shows standard deviation). See text for details.</p

MD-based structural equivalence of NSGMs to HS06 in free solution using minimum volume enclosing ellipsoid (MVEE) as one of the comparable parameters.

<p>A) shows MVEE of each MD frame for all four NSGMs; B) shows MVEE for HS06 sequence containing IdoA2S in ²S_O and ¹C₄ conformations; and C) compiles the results across the MD simulation in a box plot.</p