Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Structural and dynamic changes in P-Rex1 upon activation by PIP 3 and inhibition by IP 4

PIP 3 -dependent Rac exchanger 1 (P-Rex1) is abundantly expressed in neutrophils and plays central roles in chemotaxis and cancer metastasis by serving as a guanine-nucleotide exchange factor (GEF) for Rac. The enzyme is synergistically activated by PIP 3 and heterotrimeric Gβγ subunits, but mechanistic details remain poorly understood. While investigating the regulation of P-Rex1 by PIP 3 , we discovered that Ins(1,3,4,5)P 4 (IP 4 ) inhibits P-Rex1 activity and induces large decreases in backbone dynamics in diverse regions of the protein. Cryo-electron microscopy analysis of the P-Rex1·IP 4 complex revealed a conformation wherein the pleckstrin homology (PH) domain occludes the active site of the Dbl homology (DH) domain. This configuration is stabilized by interactions between the first DEP domain (DEP1) and the DH domain and between the PH domain and a 4-helix bundle (4HB) subdomain that extends from the C-terminal domain of P-Rex1. Disruption of the DH–DEP1 interface in a DH/PH-DEP1 fragment enhanced activity and led to a more extended conformation in solution, whereas mutations that constrain the occluded conformation led to decreased GEF activity. Variants of full-length P-Rex1 in which the DH–DEP1 and PH–4HB interfaces were disturbed exhibited enhanced activity during chemokineinduced cell migration, confirming that the observed structure represents the autoinhibited state in living cells. Interactions with PIP 3 -containing liposomes led to disruption of these interfaces and increased dynamics protein-wide. Our results further suggest that inositol phosphates such as IP 4 help to inhibit basal P-Rex1 activity in neutrophils, similar to their inhibitory effects on phosphatidylinositol-3-kinase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Design and Optimization of the IMAGINE-X Diffractometer at NB1

An upgraded protein crystal diffractometer with Dynamic Nuclear Polarization (DNP) is proposed for installation during the next HFIR Beryllium Reflector Replacement. The diffractometer will replace the IMAGINE protein crystal diffraction instrument currently located at CG4D. The simulation design meets the needs as requested in the science requirements document.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Predicting Small Molecule Transfer Free Energies by Combining Molecular Dynamics Simulations and Deep Learning

Accurately predicting small molecule partitioning and hydrophobicity is critical in the drug discovery process. There are many heterogeneous chemical environments within a cell and entire human body. For example, drugs must be able to cross the hydrophobic cellular membrane to reach their intracellular targets, and hydrophobicity is an important driving force for drug–protein binding. Atomistic molecular dynamics (MD) simulations are routinely used to calculate free energies of small molecules binding to proteins, crossing lipid membranes, and solvation but are computationally expensive. Machine learning (ML) and empirical methods are also used throughout drug discovery but rely on experimental data, limiting the domain of applicability. We present atomistic MD simulations calculating 15,000 small molecule free energies of transfer from water to cyclohexane. This large data set is used to train ML models that predict the free energies of transfer. We show that a spatial graph neural network model achieves the highest accuracy, followed closely by a 3D-convolutional neural network, and shallow learning based on the chemical fingerprint is significantly less accurate. A mean absolute error of ~4 kJ/mol compared to the MD calculations was achieved for our best ML model. We also show that including data from the MD simulation improves the predictions, tests the transferability of each model to a diverse set of molecules, and show multitask learning improves the predictions. This work provides insight into the hydrophobicity of small molecules and ML cheminformatics modeling, and our data set will be useful for designing and testing future ML cheminformatics methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 2 (Lowest Temperature Extracted Protein Coordinate Trajectories)

Protein-only trajectories of the lowest-temperature window (310K) extracted from temperature replica-exchange molecular dynamics simulations of 23 different SARS CoV-2 systems, including S Protein ACE2-receptor binding-domain, MPro, PLPro, NSP3 ADRP (X-domain/macrodomain/phosphatase), NSP15 (endoribonuclease), NSP9, NSP10, NSP16, and N-protein N-terminus. Systems were prepared with charmm-gui and simulated using the GROMACS simulation software suite. Trajectories are provided in compressed dcd format with accompanying coordinate/topology files in pdb and psf formats. This data supplements the data release, DOI: 10.13139/OLCF/1650650 ('SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 1 (Simulation Initial Conditions and Parameters)')

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Towards predictive control of reversible nanoparticle assembly with solid-binding proteins

Although a broad range of ligand-functionalized nanoparticles and physico-chemical triggers have been exploited to create stimuli-responsive colloidal systems, little attention has been paid to the reversible assembly of unmodified nanoparticles with non-covalently bound proteins. Previously, we reported that a derivative of green fluorescent protein engineered with oppositely located silica-binding peptides mediates the repeated assembly and disassembly of 10-nm silica nanoparticles when pH is toggled between 7.5 and 8.5. We captured the subtle interplay between interparticle electrostatic repulsion and their protein-mediated short-range attraction with a multiscale model energetically benchmarked to collective system behavior captured by scattering experiments. Here, in this work, we show that both solution conditions (pH and ionic strength) and protein engineering (sequence and position of engineered silica-binding peptides) provide pathways for reversible control over growth and fragmentation, leading to clusters ranging in size from 25 nm protein-coated particles to micrometer-size aggregate. We further find that the higher electrolyte environment associated with successive cycles of base addition eventually eliminates reversibility. Our model accurately predicts these multiple length scales phenomena. The underpinning concepts provide design principles for the dynamic control of other protein- and particle-based nanocomposites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small molecule–binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase–1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule–binding proteins with tuned interaction energies is feasible entirely from computation.

Science & Technology - Other Topics↗
DOE OSTI2020

Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance

For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λ max value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λ max values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λ max of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λ max using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λ max values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λ max of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λ max values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Characterization of Gramicidin A in Triblock and Diblock Polymersomes and Hybrid Vesicles via Continuous Wave Electron Paramagnetic Resonance Spectroscopy

Studying membrane proteins in a native environment is crucial to understanding their structural and/or functional studies. Often, widely accepted mimetic systems have limitations that prevent the study of some membrane proteins. Micelles, bicelles, and liposomes are common biomimetic systems but have problems with membrane compatibility, limited lipid composition, and heterogeneity. To overcome these limitations, polymersomes and hybrid vesicles have become popular alternatives. Polymersomes form from amphiphilic triblock or diblock copolymers and are considered more robust than liposomes. Hybrid vesicles are a combination of lipids and block copolymers that form vesicles composed of a mixture of the two. These hybrid vesicles are appealing because they have the native lipid environment of bilayers but also the stability and customizability of polymersomes. Gramicidin A was incorporated into these polymersomes and characterized using continuous wave electron paramagnetic resonance (CW-EPR) and transmission electron microscopy (TEM). EPR spectroscopy is a powerful biophysical technique used to study the structure and dynamic properties of membrane proteins in their native environment. Spectroscopic studies of gramicidin A have been limited to liposomes; in this study, the membrane peptide is studied in both polymersomes and hybrid vesicles using CW-EPR spectroscopy. Lineshape analysis of spin-labeled gramicidin A revealed linewidth broadening, suggesting that the thicker polymersome membranes restrict the motion of the spin label more when compared to liposome membranes. Statement of Significance: Understanding membrane proteins’ structures and functions is critical in the study of many diseases. In order to study them in a native environment, membrane mimetics must be developed that can be suitable for obtaining superior biophysical data quality to characterize structural dynamics while maintaining their native functions and structures. Many currently widely accepted methods have limitations, such as a loss of native structure and function, heterogeneous vesicle formation, restricted lipid types for the vesicle formation for many proteins, and experimental artifacts, which leaves rooms for the development of new biomembrane mimetics. The triblock and diblock polymersomes and hybrid versicles utilized in this study may overcome these limitations and provide the stability and customizability of polymersomes, keeping the biocompatibility and functionality of liposomes for EPR studies of membrane proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Calcitonin Receptor N-Glycosylation Enhances Peptide Hormone Affinity by Controlling Receptor Dynamics

The class B G protein-coupled receptor (GPCR) calcitonin receptor (CTR) is a drug target for osteoporosis and diabetes. N-glycosylation of asparagine 130 in its extracellular domain (ECD) enhances calcitonin hormone affinity with the proximal GlcNAc residue mediating this effect through an unknown mechanism. Here, we present two crystal structures of salmon calcitonin-bound, GlcNAc-bearing CTR ECD at 1.78 and 2.85 Å resolutions and analyze the mechanism of the glycan effect. The N130 GlcNAc does not contact the hormone. Surprisingly, the structures are nearly identical to a structure of hormone-bound, N-glycan-free ECD, which suggested that the GlcNAc might affect CTR dynamics not observed in the static crystallographic snapshots. Hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations revealed that glycosylation stabilized a β-sheet adjacent to the N130 GlcNAc and the N-terminal α-helix near the peptide-binding site while increasing flexibility of the peptide-binding site turret loop. These changes due to N-glycosylation increased the ligand on-rate and decreased its off-rate. The glycan effect extended to RAMP-CTR amylin receptor complexes and was also conserved in the related CGRP receptor. These results reveal that N-glycosylation can modulate GPCR function by altering receptor dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Interplays of electron and nuclear motions along CO dissociation trajectory in myoglobin revealed by ultrafast X-rays and quantum dynamics calculations

Significance In this multidisciplinary study, we investigated local and global structural responses upon carbon monoxide (CO) dissociation from the heme iron(II) in myoglobin in solution using ultrafast X-ray absorption and scattering. CO dissociation was optically triggered and pulsed X-rays probed structural and electronic changes as a function of time, revealing an iron(II) spin state transition and nuclear motions propagating from the active site throughout the entire protein. Quantum mechanical computational results were used to assign the iron K-edge absorption features, and to map out electronic and nuclear trajectories during Fe–CO bond elongation and heme macrocycle doming. The study has important implications for heme–protein function and the interaction of small molecule substrates with the metal active centers in enzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Unraveling the Hsp70-ROS-autophagy axis in pentachlorophenol-challenged lung and liver epithelial cells

Pentachlorophenol (PCP) was extensively utilized as an organochlorine pesticide and wood preservative in the United States from the 1930s until the Environmental Protection Agency (EPA) imposed restrictions due to concerns about its toxicity and potential carcinogenic properties. Although it is no longer widely used, PCP remains a concern due to its environmental persistence and potential for long-term health effects. Significant occupational and environmental exposures have likely occurred, with the health and economic costs of PCP exposure potentially being substantial given its known toxicity. Notably, PCP exhibits rapid absorption through both the skin and respiratory system and has been shown to cause hepatotoxicity, developmental toxicity, immunotoxicity, irritation, and carcinogenicity in laboratory animal studies. PCP exposure induces oxidative stress, a key mechanism underlying its inflammatory and toxic effects, which can activate cellular stress responses including upregulation of heat shock protein 70 (Hsp70). Previous studies in lung and liver epithelial cells have shown that Hsp70 and oxidative stress play pivotal roles in triggering autophagy. This study establishes the critical role of the Hsp70-reactive oxygen species (ROS)-autophagy axis in regulating cellular responses to PCP exposure in human alveolar (A549) and liver carcinoma (HepG2) epithelial cells. Our research elucidated the molecular mechanisms underlying PCP's cellular effects, demonstrating that its exposure resulted in increased expression of autophagy-related proteins (Beclin-1, LC3B, ATG12, and ATG16), subunits of NADPH oxidase (NCF-1, NCF-2, NOX2, and Rac), and antioxidant proteins (SOD and GPx) in both lung and liver cell types. Notably, PCP augmented the interaction between Hsp70 and the autophagy regulator Beclin-1. Pretreatment with the ROS inhibitor N-acetylcysteine or Hsp70 knockdown markedly reversed PCP-induced responses. Our in-silico protein–protein docking analysis and molecular dynamics simulation studies revealed enhanced interactions and/or stable confirmations maintained throughout the simulations for TLR4-Hsp70 and Hsp70-Beclin-1 complexes in the presence of PCP. These findings provide a strong foundation for future studies, employing in vivo experimental models and human populations to identify promising targets for PCP-induced toxicity and cellular injury. As a result, these findings may have far-reaching implications for public health and environmental policy, ultimately leading to the identification of biomarkers and the development of more effective interventions for environmentally induced toxicity and diseases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A structured RNA motif locks Argonaute2:miR-122 onto the 5’ end of the HCV genome

microRNAs (miRNAs) form regulatory networks in metazoans. Viruses engage miRNA networks in numerous ways, with Flaviviridae members exploiting direct interactions of their RNA genomes with host miRNAs. For hepatitis C virus (HCV), binding of liver-abundant miR-122 stabilizes the viral RNA and regulates viral translation. Here, we investigate the structural basis for these activities, taking into consideration that miRNAs function in complex with Argonaute (Ago) proteins. The crystal structure of the Ago2:miR-122:HCV complex reveals a structured RNA motif that traps Ago2 on the viral RNA, masking its 5’ end from enzymatic attack. The trapped Ago2 can recruit host factor PCBP2, implicated in viral translation, while binding of a second Ago2:miR-122 competes with PCBP2, creating a potential molecular switch for translational control. Combined results reveal a viral RNA structure that modulates Ago2:miR-122 dynamics and repurposes host proteins to generate a functional analog of the mRNA cap-binding complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

TeraChem: A graphical processing unit-accelerated electronic structure package for large-scale ab initio molecular dynamics

TeraChem was born in 2008 with the goal of providing fast on-the-fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low-cost parallel computer architecture that became more accessible for general-purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Here, electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large-scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large-scale nonadiabatic dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Familial Alzheimer mutations stabilize synaptotoxic γ-secretase-substrate complexes

Mutations that cause familial Alzheimer’s disease (FAD) are found in amyloid precursor protein (APP) and presenilin, the catalytic component of γ-secretase, that together produce amyloid β-peptide (Aβ). Nevertheless, whether Aβ is the primary disease driver remains controversial. We report here that FAD mutations disrupt initial proteolytic events in the multistep processing of APP substrate C99 by γ-secretase. Cryoelectron microscopy reveals that a substrate mimetic traps γ-secretase during the transition state, and this structure aligns with activated enzyme-substrate complex captured by molecular dynamics simulations. In silico simulations and in cellulo fluorescence microscopy support stabilization of enzyme-substrate complexes by FAD mutations. Neuronal expression of C99 and/or presenilin-1 in Caenorhabditis elegans leads to synaptic loss only with FAD-mutant transgenes. Designed mutations that stabilize the enzyme-substrate complex and block Aβ production likewise led to synaptic loss. Collectively, these findings implicate the stalled process—not the products—of γ-secretase cleavage of substrates in FAD pathogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Elucidating interprotein energy transfer dynamics within the antenna network from purple bacteria

In photosynthesis, absorbed light energy transfers through a network of antenna proteins with near-unity quantum efficiency to reach the reaction center, which initiates the downstream biochemical reactions. While the energy transfer dynamics within individual antenna proteins have been extensively studied over the past decades, the dynamics between the proteins are poorly understood due to the heterogeneous organization of the network. Previously reported timescales averaged over such heterogeneity, obscuring individual interprotein energy transfer steps. Here, we isolated and interrogated interprotein energy transfer by embedding two variants of the primary antenna protein from purple bacteria, light-harvesting complex 2 (LH2), together into a near-native membrane disc, known as a nanodisc. We integrated ultrafast transient absorption spectroscopy, quantum dynamics simulations, and cryogenic electron microscopy to determine interprotein energy transfer timescales. By varying the diameter of the nanodiscs, we replicated a range of distances between the proteins. The closest distance possible between neighboring LH2, which is the most common in native membranes, is 25 Å and resulted in a timescale of 5.7 ps. Larger distances of 28 to 31 Å resulted in timescales of 10 to 14 ps. Corresponding simulations showed that the fast energy transfer steps between closely spaced LH2 increase transport distances by ∼15%. Overall, our results introduce a framework for well-controlled studies of interprotein energy transfer dynamics and suggest that protein pairs serve as the primary pathway for the efficient transport of solar energy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Surface Enhanced Raman Scattering Selectivity in Proteins Arises from Electron Capture and Resonant Enhancement of Radical Species

Plasmon-enhanced Raman scattering is a powerful approach to detecting and characterizing proteins in live and dynamic biological systems. However, the selective detection/enhancement of specific residues as well as spectral diffusion and fluctuations have complicated the interpretation of enhanced Raman spectra and images of biological matter. In this paper, we therefore demonstrate that the amino acid tryptophan (Trp) can capture an electron from an excited plasmon, which generates a radical anion that is resonantly enhanced: a visible excited electronic state slides into resonance upon charging. This surface enhanced resonance Raman scattering (SERRS) mechanism explains the persistence of Trp signatures in the SERS and TERS spectra of proteins. Evidence for this picture includes the observation of visible resonances in the UV-Vis extinction spectrum, changes in the ground state vibrational spectrum, and plasmon-resonance dependent behavior. DFT calculations support the experimental observations. The behavior observed from the free Trp molecule is shown to explain the SERS spectrum of the Trp-cage protein. In effect, resonant Raman scattering from radicals formed through plasmonic excitation represents an under-investigated mechanism that may be exploited for chemical sensing applications. This work was supported by the National Institutes of Health Award R01-GM109988. The authors thank Hannah Shafaat for helpful discussions and providing an azurin sample for comparison. The authors also thank Joel Harris for providing constructive comments on our manuscript.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Coarse-Grained Simulations for the Characterization and Optimization of Hybrid Protein–DNA Nanostructures

We present here the combination of experimental and computational modeling tools for the design and characterization of protein–DNA hybrid nanostructures. Our work incorporates several features in the design of these nanostructures: (1) modeling of the protein–DNA linker identity and length; (2) optimizing the design of protein–DNA cages to account for mechanical stresses; (3) probing the incorporation efficiency of protein–DNA conjugates into DNA nanostructures. The modeling tools were experimentally validated using structural characterization methods like cryo-TEM and AFM. Our method can be used for fitting low-resolution electron density maps when structural insights cannot be deciphered from experiments, as well as enable in-silico validation of nanostructured systems before their experimental realization. Finally, these tools will facilitate the design of complex hybrid protein–DNA nanostructures that seamlessly integrate the two different biomolecules.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The emergence and transmission dynamics of HIV-1 CRF07_BC in Mainland China

A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset (‘full1155’), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences (‘backbone_295’) and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989–August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.

60 APPLIED LIFE SCIENCES↗