[Study on Migrants Found in Migration Solutions from Commercially Available Polyethylene Products].

Migrants found in migration solutions obtained from commercially available polyethylene products that may contain food were studied and analysed via liquid chromatography coupled with quadrupole time-of-flight mass spectrometry (LC-QTOF) for non-target screening and LC-MS/MS for quantifying 14 substances in migration solutions. Furthermore, an analytical approach based on the retention gap was developed for accurate separation techniques using LC-MS/MS. Irganox 1076 was detected at a maximum of 1.5 mg/kg, which was 1/4 of the Specific Migration Limit in the EU, in nine commercially available plastic bags tested. This is in accordance with European Regulation No 10/2011/EU. Furthermore, migration of Erucamide and Irgafos 168-oxide was confirmed.

Acceleration of generalized replica exchange with solute tempering simulations of large biological systems on massively parallel supercomputer.

Generalized replica exchange with solute tempering (gREST) is one of the enhanced sampling algorithms for proteins or other systems with rugged energy landscapes. Unlike the replica-exchange molecular dynamics (REMD) method, solvent temperatures are the same in all replicas, while solute temperatures are different and are exchanged frequently between replicas for exploring various solute structures. Here, we apply the gREST scheme to large biological systems containing over one million atoms using a large number of processors in a supercomputer. First, communication time on a multi-dimensional torus network is reduced by matching each replica to MPI processors optimally. This is applicable not only to gREST but also to other multi-copy algorithms. Second, energy evaluations, which are necessary for the multistate bennet acceptance ratio (MBAR) method for free energy estimations, are performed on-the-fly during the gREST simulations. Using these two advanced schemes, we observed 57.72 ns/day performance in 128-replica gREST calculations with 1.5 million atoms system using 16,384 nodes in Fugaku. These schemes implemented in the latest version of GENESIS software could open new possibilities to answer unresolved questions on large biomolecular complex systems with slow conformational dynamics.

[Surveillance of the Naturally Derived Benzoic Acid Levels in Fruits and Fruit Products, and Comparison of Analytical Methods].

Benzoic acid (BA) is typically found in natural food; therefore, naturally occurring BA must be distinguished from added BA preservatives. In this study, we investigated BA levels in 100 samples of fruit products and their fresh fruits as raw materials using dialysis and steam distillation approaches. BA was detected in the range (minimum-maximum) of 2.1-1380 µg/g and 2.2-1950 µg/g in dialysis and steam distillation, respectively. Steam distillation indicated higher BA levels than dialysis.

[Improvement of Nitrite Analysis Method in Food Products].

The conventional analysis method has problems with extraction efficiency, operability, and reproducibility. In this study, we attempted to solve these problems and improve the analytical method to obtain sufficient extraction efficiency and good operability and accuracy. The conventional method was able to get sufficient extraction in dried meat products, where the extraction efficiency of the conventional method was low, by increasing the concentration of sodium hydroxide solution at the time of homogenization. Suction filtration after adding the defoaming agent was added allowed for accurate volume adjustment. The turbidity of the extract caused by insufficient addition of zinc acetate solution was removed by increasing the amount of zinc acetate solution that was added. Turbidity caused by starch was removed by adding pancreatin. The RSD of the quantitative values was improved by adding sodium hydroxide solution and 80-90℃ water and immediately homogenizing. Furthermore, by changing the dilution factor of the extract solution in the colorimetric method, the inhibition of coloration by reducing substances was suppressed, and more accurate quantitative values could be obtained than with the conventional method. The recovery rate was 78.5-105% (RSD 0.7-5.8%), which was a good result. This method was considered to be a useful analytical method that can contribute to improving the inspection accuracy of nitrite ion analysis.

[Storage Stability of Nitrite in Nitrite Analysis of Meat Products].

For the analysis of nitrite ions in food, the stabilities of nitrite ions in meat products and their standard solutions were evaluated. Nitrite is easily oxidized or reduced; hence, products with standard solutions or colour retention agent must be carefully handled. To assess the stability and decreasing trend of nitrite, we examined the storage stability of standard solutions using calibration curves, the time course of nitrite in chopped meat products stored under different conditions, and the time course of nitrite in the sample solutions. Regarding calibration curves, the storage stability was determined for standard solutions that were prepared with ultrapure water at concentrations of 0.025 and 0.4 µg/mL and were stored at 5℃ for one year. The results revealed no changes in concentration of any solution over time, suggesting that no readjustments to the standard solution concentration were necessary before testing until one year after their preparation. Time course of nitrite in chopped meat products stored under different conditions showed a significant decrease in nitrite in refrigerated storage (5℃), whereas stability of nitrite was maintained for up to 1 day in frozen storage (-20℃) and for 14 days in frozen storage (-40℃). The time course of nitrite in the sample solutions showed that the quantitative values of nitrite in the extract remained unchanged within one week of extraction for the meat products tested in the study.

Multiple sub-state structures of SERCA2b reveal conformational overlap at transition steps during the catalytic cycle.

Sarco/endoplasmic reticulum Ca2+ ATPase (SERCA) pumps Ca2+ into the endoplasmic reticulum (ER). Herein, we present cryo-electron microscopy (EM) structures of three intermediates of SERCA2b: Ca2+-bound phosphorylated (E1P·2Ca2+) and Ca2+-unbound dephosphorylated (E2·Pi) intermediates and another between the E2P and E2·Pi states. Our cryo-EM analysis demonstrates that the E1P·2Ca2+ state exists in low abundance and preferentially transitions to an E2P-like structure by releasing Ca2+ and that the Ca2+ release gate subsequently undergoes stepwise closure during the dephosphorylation processes. Importantly, each intermediate adopts multiple sub-state structures including those like the next one in the catalytic series, indicating conformational overlap at transition steps, as further substantiated by atomistic molecular dynamic simulations of SERCA2b in a lipid bilayer. The present findings provide insight into how enzymes accelerate catalytic cycles.

Implementation of residue-level coarse-grained models in GENESIS for large-scale molecular dynamics simulations.

Residue-level coarse-grained (CG) models have become one of the most popular tools in biomolecular simulations in the trade-off between modeling accuracy and computational efficiency. To investigate large-scale biological phenomena in molecular dynamics (MD) simulations with CG models, unified treatments of proteins and nucleic acids, as well as efficient parallel computations, are indispensable. In the GENESIS MD software, we implement several residue-level CG models, covering structure-based and context-based potentials for both well-folded biomolecules and intrinsically disordered regions. An amino acid residue in protein is represented as a single CG particle centered at the Cα atom position, while a nucleotide in RNA or DNA is modeled with three beads. Then, a single CG particle represents around ten heavy atoms in both proteins and nucleic acids. The input data in CG MD simulations are treated as GROMACS-style input files generated from a newly developed toolbox, GENESIS-CG-tool. To optimize the performance in CG MD simulations, we utilize multiple neighbor lists, each of which is attached to a different nonbonded interaction potential in the cell-linked list method. We found that random number generations for Gaussian distributions in the Langevin thermostat are one of the bottlenecks in CG MD simulations. Therefore, we parallelize the computations with message-passing-interface (MPI) to improve the performance on PC clusters or supercomputers. We simulate Herpes simplex virus (HSV) type 2 B-capsid and chromatin models containing more than 1,000 nucleosomes in GENESIS as examples of large-scale biomolecular simulations with residue-level CG models. This framework extends accessible spatial and temporal scales by multi-scale simulations to study biologically relevant phenomena, such as genome-scale chromatin folding or phase-separated membrane-less condensations.

The inherent flexibility of receptor binding domains in SARS-CoV-2 spike protein.

Spike (S) protein is the primary antigenic target for neutralization and vaccine development for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It decorates the virus surface and undergoes large motions of its receptor binding domains (RBDs) to enter the host cell. Here, we observe Down, one-Up, one-Open, and two-Up-like structures in enhanced molecular dynamics simulations, and characterize the transition pathways via inter-domain interactions. Transient salt-bridges between RBDA and RBDC and the interaction with glycan at N343B support RBDA motions from Down to one-Up. Reduced interactions between RBDA and RBDB in one-Up induce RBDB motions toward two-Up. The simulations overall agree with cryo-electron microscopy structure distributions and FRET experiments and provide hidden functional structures, namely, intermediates along Down-to-one-Up transition with druggable cryptic pockets as well as one-Open with a maximum exposed RBD. The inherent flexibility of S-protein thus provides essential information for antiviral drug rational design or vaccine development.

[Study of Rapid HPLC Method for the Determination of Antioxidants in Foods, and Standard Solution Stability during Storage].

In this study, we simultaneously determined three antioxidants, butylated hydroxytoluene (BHT), butylated hydroxyanisole (BHA), and tert-butylhydroquinone (TBHQ), using HPLC equipped both with a photodiode array detector and a fluorescence detector in 25 minutes per sample. Due to the combined use of the two detectors, we could achieve improved target selectivity. Further, quantification at the specific wavelengths for each target substance particularly increased BHT detection sensitivity. This approach enabled us to avoid repeated measurements during daily inspections. Furthermore, detections were performed using LC-MS/MS instead of GC-MS to overcome the problem of helium gas shortage.In addition, we investigated antioxidant stability in standard solutions during storage. Although TBHQ was stable in methanol with ascorbic acid at -20℃, ascorbic acid storage has possibility to lead to decrease in BHT and BHA concentrations. We recognized that the mixture of BHT and BHA dissolved in methanol at 4℃ and that of BHT, BHA and TBHQ dissolved in methanol with ascorbic acid at -20℃ were suitable for about one year.

Structural and energetic analysis of metastable intermediate states in the E1P-E2P transition of Ca2+-ATPase.

Sarcoplasmic reticulum (SR) Ca2+-ATPase transports two Ca2+ ions from the cytoplasm to the SR lumen against a large concentration gradient. X-ray crystallography has revealed the atomic structures of the protein before and after the dissociation of Ca2+, while biochemical studies have suggested the existence of intermediate states in the transition between E1P⋅ADP⋅2Ca2+ and E2P. Here, we explore the pathway and free energy profile of the transition using atomistic molecular dynamics simulations with the mean-force string method and umbrella sampling. The simulations suggest that a series of structural changes accompany the ordered dissociation of ADP, the A-domain rotation, and the rearrangement of the transmembrane (TM) helices. The luminal gate then opens to release Ca2+ ions toward the SR lumen. Intermediate structures on the pathway are stabilized by transient sidechain interactions between the A- and P-domains. Lipid molecules between TM helices play a key role in the stabilization. Free energy profiles of the transition assuming different protonation states suggest rapid exchanges between Ca2+ ions and protons when the Ca2+ ions are released toward the SR lumen.

Optimized Hydrogen Mass Repartitioning Scheme Combined with Accurate Temperature/Pressure Evaluations for Thermodynamic and Kinetic Properties of Biological Systems.

In recent years, molecular dynamics (MD) simulations with larger time steps have been performed with the hydrogen-mass-repartitioning (HMR) scheme, where the mass of each hydrogen atom is increased to reduce high-frequency motion while the mass of a non-hydrogen atom bonded to a hydrogen atom is decreased to keep the total molecular mass unchanged. Here, we optimize the scaling factors in HMR and combine them with previously developed accurate temperature/pressure evaluations. The heterogeneous HMR scaling factors are useful to avoid the structural instability of amino acid residues having a five- or six-membered ring in MD simulations with larger time steps. It also reproduces kinetic properties, namely translational and rotational diffusions, if the HMR scaling factors are applied to only solute molecules. The integration scheme is tested for biological systems that include soluble/membrane proteins and lipid bilayers for about 200 µs MD simulations in total and give consistent results in MD simulations with both a small time step of 2.0 fs and a large, multiple time step integration with time steps of 3.5 fs (for fast motions) and 7.0 fs (for slower motions). We also confirm that the multiple time step integration scheme used in this study provides more accurate energy conservations than the RESPA/C1 and is comparable to the RESPA/C2 in NAMD. In summary, the current integration scheme combining the optimized HMR with accurate temperature/pressure evaluations can provide stable and reliable MD trajectories with a larger time step, which are computationally more than 2-fold efficient compared to the conventional methods.

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins.

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

[Attempt to Calculate the Efficiency of Methylation Using 1H-NMR for Analyzing Brominated Vegetable Oil].

We studied the efficiency of methylation for analyzing brominated vegetable oil (BVO). In this report, we investigated whether 1H-NMR is an applicable method for assessing the efficiency of methylation to analyze BVO. 1H-NMR sufficiently calculated the efficiency of methylation using each integral and the numbers of protons derived from the methyl group, which is characteristic in products, and the methine group, which is characteristic in unreacted substances. Additionally, the efficiency of methylation calculated via 1H-NMR was in good agreement with changes in the peak area of BVO fatty acid methyl esters (BVOFAMEs) after various heating times obtained from GC-FID analysis. Therefore, 1H-NMR is applicable for calculating the efficiency of methylation to analyze BVO.

Elucidation of interactions regulating conformational stability and dynamics of SARS-CoV-2 S-protein.

The ongoing COVID-19 pandemic caused by the new coronavirus, SARS-CoV-2, calls for urgent developments of vaccines and antiviral drugs. The spike protein of SARS-CoV-2 (S-protein), which consists of trimeric polypeptide chains with glycosylated residues on the surface, triggers the virus entry into a host cell. Extensive structural and functional studies on this protein have rapidly advanced our understanding of the S-protein structure at atomic resolutions, although most of these structural studies overlook the effect of glycans attached to the S-protein on the conformational stability and functional motions between the inactive down and active up forms. Here, we performed all-atom molecular dynamics simulations of both down and up forms of a fully glycosylated S-protein in solution as well as targeted molecular dynamics simulations between them to elucidate key interdomain interactions for stabilizing each form and inducing the large-scale conformational transitions. The residue-level interaction analysis of the simulation trajectories detects distinct amino acid residues and N-glycans as determinants on conformational stability of each form. During the conformational transitions between them, interdomain interactions mediated by glycosylated residues are switched to play key roles on the stabilization of another form. Electrostatic interactions, as well as hydrogen bonds between the three receptor binding domains, work as driving forces to initiate the conformational transitions toward the active form. This study sheds light on the mechanisms underlying conformational stability and functional motions of the S-protein, which are relevant for vaccine and antiviral drug developments.

Analyses of the brominated vegetable oil in soft drinks using gas chromatography-flame ionization detector and atmospheric pressure gas chromatography-quadrupole/time-of-flight mass spectrometry.

Herein, a new method for quantifying the brominated vegetable oil content in commercial soft drinks was developed, which accelerated the sample preparation process and improved analytical efficiency. First, simple and accurate chromatographic separation techniques were performed using a VF-5ht column for both GC-FID (quantitative) and APGC-QTOF (qualitative) analyses. The samples were subjected to chromatography on a reversed-phase solid-phase extraction cartridge. (PoraPak™RxnRP). Transesterification using a boron trifluoride methanol complex in methanol solution was performed. When validating this method, the analyte recovery percentages were between 82.2% and 99.9%, and the recovery and standard deviation of repeatability values were between 1.2% and 3.5%. Using an isotope library, the bromostearic acid methyl esters (9,10-dibromostearic acid methyl ester, 9,10,12,13-tetrabromostearic acid methyl ester, and hexabromostearic acid methyl ester) in the sample mixtures were qualitatively confirmed via APGC-QTOF. A novel aqueous ammonium adduct, which has not been previously reported, was also confirmed. These results indicated that this new method was simple, accurate, and also allowed for precise qualitative and quantitative confirmation as well as high reproducibility.

New parallel computing algorithm of molecular dynamics for extremely huge scale biological systems.

In this paper, we address high performance extreme-scale molecular dynamics (MD) algorithm in the GENESIS software to perform cellular-scale molecular dynamics (MD) simulations with more than 100,000 CPU cores. It includes (1) the new algorithm of real-space nonbonded interactions maximizing the performance on ARM CPU architecture, (2) reciprocal-space nonbonded interactions minimizing communicational cost, (3) accurate temperature/pressure evaluations that allows a large time step, and (4) effective parallel file inputs/outputs (I/O) for MD simulations of extremely huge systems. The largest system that contains 1.6 billion atoms was simulated using MD with a performance of 8.30 ns/day on Fugaku supercomputer. It extends the available size and time of MD simulations to answer unresolved questions of biomacromolecules in a living cell.

A singularity-free torsion angle potential for coarse-grained molecular dynamics simulations.

Conventional torsion angle potentials used in molecular dynamics (MD) have a singularity problem when three bonded particles are collinearly aligned. This problem is often encountered in coarse-grained (CG) simulations. Here, we propose a new form of the torsion angle potential, which introduces an angle-dependent modulating function. By carefully tuning the parameters for this modulating function, our method can eliminate the problematic angle-dependent singularity while being combined with existing models. As an example, we optimized the modulating function of the torsion angle potential for popular CG models of biomolecules based on the statistics over experimental structures deposited in the Protein Data Bank. By applying our method to designed and natural biomolecules, we show that the new torsion angle potential is able to eliminate the singularity problem while maintaining the structural features in the original models. Furthermore, by comparing our design with previous methods, we found that our new potential has advantages in computational efficiency and numerical stability. We strongly recommend the usage of our new potential in the CG simulations of flexible molecules.

[Determination of Nitrate Ion in Cheese and Sake].

The Japanese official analysis method for determination of nitrate ions in food products used as food additives is associated with various challenges. In some kinds of cheese, the extract becomes suspended. The volume of extracted solution is often not accurate owing to the presence of residues in the solution. Moreover, the determination with liquid chromatography-ultraviolet detection (HPLC-UV) is difficult owing to the influence of impurities. Sake usually does not contain lipids or proteins ; therefore, its analysis can be simplified by omitting the co-precipitation steps to remove them. In the present study, for cheese, the amount of sodium hydroxide solution that causes suspension was reduced, and the influence of residues was removed by adjusting the volume after suction filtration. Whereas, sake was diluted with water and centrifuged. Furthermore, solid-phase extraction (SPE) method using cartridge containing carbon molecular sieve to remove the influence of impurities on the chromatogram was successfully established. The recoveries of the nitrate ions were good outcomes of 91.3-99.6% (CV 0.9-4.5%) (n=5). The analysis range was 0.010-0.20 g/kg for cheese, 0.010-0.20 g/L for milk, and 0.010-0.10 g/kg for sake. The developed analysis methods are considered useful, because various challenges of the official analysis method can be solved and the operation are efficient.

Determination of nine preservatives in processed foods using a modified QuEChERS extraction and quantified by HPLC-PDA.

In this study, a new method was developed for simultaneously determining nine preservatives, that is, benzoic acid (BA), sorbic acid (SOA), dehydroacetic acid (DHA) and PHBAs (methyl p-hydroxybenzoate [PHBA-me], ethyl p-hydroxybenzoate [PHBA-et], isopropyl p-hydroxybenzoate [PHBA-ipro], propyl p-hydroxybenzoate [PHBA-npro]), isobutyl p-hydroxybenzoate [PHBA-ibut] and butyl p-hydroxybenzoate [PHBA-nbut]), in processed foods, employing liquid chromatography (LC). This procedure accelerated sample preparation and improved efficiency by employing modified quick, easy, cheap, effective, rugged and safe (QuEChERS) extraction without clean-up. Samples were prepared with 20 mL of acetonitrile/water (1:1) with the assistance of a ceramic stone. The extract solutions were diluted 10 times or according to the detection amount and then injected into an LC-PDA. This method showed good linearity, and the LOQs were 10 mg/kg for BA, SOA and DHA and 5 mg/kg for the PHBAs. When validating this method, the recoveries of the nine preservatives were in the range 77.0-99.6%, RSDr values were in the range 0.7-5.3% and those of RSDwr were in the range 2.3-8.4%. These results suggest that this new method is highly reproducible.