Novel umami peptides from two Termitomyces mushrooms and molecular docking to the taste receptor T1R1/T1R3

Wild edible Termitomyces mushrooms are popular in Southwest China and umami is important flavor qualities of edible mushrooms.This study aimed to understand the umami taste of Termitomyces intermedius and Termitomyces aff.bulborhizus.Ten umami peptides from aqueous extracts were separated using a Sephadex G-15 gel filtration chromatography.The intense umami fraction was evaluated by both sensory evaluation and electronic tongue.They were identified as KLNDAQAPK,DSTDEKFLR,VGKGAHLSGEH,MLKKKKLA,SLGFGGPPGY,TVATFSSSTKPDD,AMDDDEADLLLLAM,VEDEDEKPKEK,SPEEKKEEET and PEGADKPNK.Seven peptides,except VEDEDEKPKEK,SPEEKKEEET and PEGADKPNK were selectively synthesized to verify their taste characteristics.All these 10 peptides had umami or salt taste.The 10 peptides were conducted by molecular docking to study their interaction with identified peptides and the umami taste receptor T1R1/T1R3.All these 10 peptides perfectly docked the active residues in the T1R3 subunit.Our results provide theoretical basis for the umami taste and address the umami mechanism of two wild edible Termitomyces mushrooms.

Distinct molecular targets of ProEGCG from EGCG and superior inhibition of angiogenesis signaling pathways for treatment of endometriosis

Endometriosis is a common chronic gynecological disease with endometrial cell implantation outside the uterus.Angiogenesis is a major pathophysiology in endometriosis.Our previous studies have demonstrated that the prodrug of epigallocatechin gallate(ProEGCG)exhibits superior anti-endometriotic and anti-angiogenic effects compared to epigallocatechin gallate(EGCG).However,their direct binding targets and underlying mechanisms for the differential effects remain unknown.In this study,we demonstrated that oral ProEGCG can be effective in preventing and treating endometriosis.Additionally,1D and 2D Proteome Integral Solubility Alteration assay-based chemical proteomics identified metadherin(MTDH)and PX domain containing serine/threonine kinase-like(PXK)as novel binding targets of EGCG and ProEGCG,respectively.Computational simulation and BioLayer interferometry were used to confirm their binding affinity.Our results showed that MTDH-EGCG inhibited protein kinase B(Akt)-mediated angiogenesis,while PXK-ProEGCG inhibited epidermal growth factor(EGF)-mediated angiogenesis via the EGF/hypoxia-inducible factor(HIF-1a)/vascular endothelial growth factor(VEGF)pathway.In vitro and in vivo knockdown assays and microvascular network imaging further confirmed the involvement of these signaling pathways.Moreover,our study demonstrated that ProEGCG has superior therapeutic effects than EGCG by targeting distinct signal transduction pathways and may act as a novel antiangiogenic therapy for endometriosis.

Molecular Mechanism and Molecular Design of Lubricating Oil Antioxidants

To overcome the limitations of traditional experimental“trial and error”methods in lubricant additive design,a new molecular design method based on molecular structure parameters is established here.The molecular mechanism of the antioxidant reaction of hindered phenol,diphenylamine,and alkyl sulfide are studied via molecular simulations.Calculation results show that the strong electron-donating ability and high hydrogen-donating activity of the antioxidant molecule and the low hydrogen-abstracting activity of free radicals formed after dehydrogenation are the internal molecular causes of the shielding of phenol and diphenylamine from scavenging peroxy free radicals,and the strong electron-donating ability is the internal molecular cause of the high activity of thioether in decomposing alkyl hydrogen peroxide.Based on this antioxidant molecular mechanism,a molecular design rule of antioxidant is proposed,namely“high EHOMO,large Q(S),low bond dissociation energy BDE(O—H)and BDE(N—H)”.Two new antioxidants,PAS-I and PAS-II,are designed and prepared by chemical bonding of hindered phenol,diphenylamine,and sulfur atoms.Experimental results show that these antioxidants both have excellent antioxidant effects in lubricating oil,and that PAS-II is the superior antioxidant,consistent with theoretical predictions.

Molecular phylogeny and taxonomy of Phlomoides(Lamiaceae subfamily Lamioideae)in China:Insights from molecular and morphological data

Phlomoides,with 150-170 species,is the second largest and perhaps most taxonomically challenging genus within the subfamily Lamioideae(Lamiaceae).With about 60 species,China is one of three major biodiversity centers of Phlomoides.Although some Phlomoides species from China have been included in previous molecular phylogenetic studies,a robust and broad phylogeny of this lineage has yet to be completed.Moreover,given the myriad new additions to the genus,the existing infrageneric classification needs to be evaluated and revised.Here,we combine molecular and morphological data to investigate relationships within Phlomoides,with a focus on Chinese species.We observed that plastid DNA sequences can resolve relationships within Phlomoides better than nuclear ribosomal internal and external transcribed spacer regions(nrITS and nrETS).Molecular phylogenetic analyses confirm the monophyly of Phlomoides,but most previously defined infrageneric groups are not monophyletic.In addition,morphological analysis demonstrates the significant taxonomic value of eight characters to the genus.Based on our molecular phylogenetic analyses and morphological data,we establish a novel section Notochaete within Phlomoides,and propose three new combinations as well as three new synonyms.This study presents the first molecular phylogenetic analyses of Phlomoides in which taxa representative of the entire genus are included,and highlights the phylogenetic and taxonomic value of several morphological characters from species of Phlomoides from China.Our study suggests that a taxonomic revision and reclassification for the entire genus is necessary in the future.

Absorption characteristics,model,and molecular mechanism of hydrogen sulfide in morpholine acetate aqueous solution

The solubility of H_(2)S was measured in solutions of N-butyl-N-methylmorpholine acetate([Bmmorp][Ac])containing 20%-40%(mass)water at experimental temperatures ranged from 298.15 to 328.15 K and pressures up to 320 k Pa.The total solubility of H_(2)S increased with higher temperatures,lower pressures,and reduced water content.The reaction equilibrium thermodynamic model was used to correlate the solubility data.The results indicate that the chemical reaction equilibrium constant decrease with increasing water content and temperature,whereas Henry constant increase with increasing water content and temperature.Compared with other ionic liquids,H_(2)S exhibits a higher physical absorption enthalpy and a lower chemical absorption enthalpy in[Bmmorp][Ac]aqueous solution.This suggests that[Bmmorp][Ac]has a strong physical affinity for H_(2)S and low energy requirement for desorption.Quantum chemical methods were used to investigate the molecular mechanism of H_(2)S absorption in ionic liquids.The interaction energy analysis revealed that the binding of H_(2)S with the ionic liquid in a1:2 ratio is more stable.Detailed analyses by the methods of the interaction region indicator and the atoms in molecules were conducted to the interactions between H_(2)S and the ionic liquid.

Comprehensive understanding of glioblastoma molecular phenotypes:classification,characteristics,and transition

Among central nervous system-associated malignancies,glioblastoma(GBM)is the most common and has the highest mortality rate.The high heterogeneity of GBM cell types and the complex tumor microenvironment frequently lead to tumor recurrence and sudden relapse in patients treated with temozolomide.In precision medicine,research on GBM treatment is increasingly focusing on molecular subtyping to precisely characterize the cellular and molecular heterogeneity,as well as the refractory nature of GBM toward therapy.Deep understanding of the different molecular expression patterns of GBM subtypes is critical.Researchers have recently proposed tetra fractional or tripartite methods for detecting GBM molecular subtypes.The various molecular subtypes of GBM show significant differences in gene expression patterns and biological behaviors.These subtypes also exhibit high plasticity in their regulatory pathways,oncogene expression,tumor microenvironment alterations,and differential responses to standard therapy.Herein,we summarize the current molecular typing scheme of GBM and the major molecular/genetic characteristics of each subtype.Furthermore,we review the mesenchymal transition mechanisms of GBM under various regulators.

Balanced Fracturing and Cold-welding of Magnesium during Ball Milling Assisted by Carbon Coating:Experimental and Molecular Dynamic Simulation

The lignite-derived carbon from self-protection pyrolysis was employed to balance the fracturing and cold-welding of magnesium during ball milling.Particle size analysis indicates that the introduction of lignite-derived carbon can effectively reduce the particle size of Mg while the introduction of graphite does no help.Besides,the effect of lignite-derived carbon on crystallite size reduction of Mg is also better than graphite.A moderate cold-welding phenomenon was observed after ball-milling Mg with the lignite-derived carbon,suggesting less Mg is wasted on the milling vials and balls.Molecular dynamic simulations reveal that the balanced fracturing and cold-welding of magnesium during ball milling is mainly attributed to the special structure of the lignite-derived carbon:graphitized short-range ordered stacking function as dry lubricant and irregular shape/sharp edge function as milling aid.The preliminary findings in current study are expected to offer implications for designing efficient Mg-based hydrogen storage materials.

B-COPNA resin formation from ethylene tar light fractions:Process development and mechanical exploration by molecular simulation

An efficient utilization strategy of ethylene tar(ET),the main by-product of the ethylene cracking unit,is urgently required to meet demands for modern petrochemical industry.On the other hand,condensed polynuclear aromatic resin of moderate condensation degree(B-COPNA)is a widely used carbon material due to its superb processability,the production of which is,however,seriously limited by the high cost of raw materials.Under such context,an interesting strategy was proposed in this study for producing B-COPNA resin using crosslinked light fractions of ethylene tar(ETLF,boiling point<260℃)facilitated by molecular simulation.1,4-Benzenedimethanol(PXG)was first selected as the crosslinking agent according to the findings of molecular simulation.The effects of operating conditions,including reactions temperature,crosslinking agent,and catalyst content on the softening point and yield of B-COPNA resin products were then investigated to optimize the process.The reaction mechanism of resin production was studied by analyzing the molecular structure and transition state of ETLF and crosslinking agents.It was shown that PXG exhibited a superior capacity of withdrawing electrons and a higher electrophilic reactivity than other crosslinking agents.In addition to the highest yield and greatest heat properties,PXG-prepared resin contained the most condensed aromatics.The corresponding optimized conditions of resin preparation were 180℃,1:1.9(PXG:ETLF),and 3%(mass)of catalyst content with a resin yield of 78.57%.It was the electrophilic substitution reaction that occurred between the ETLF and crosslinking agent molecules that were responsible for the resin formation,according to the experimental characterization and molecular simulation.Hence,it was confirmed that the proposed strategy and demonstrated process can achieve a clean and high value-added utilization of ETLF via B-COPNA resin preparation,bringing huge economic value to the current petrochemical industry.

Molecular dynamics simulations on the interactions between nucleic acids and a phospholipid bilayer

Recently,lipid nanoparticles(LNPs)have been extensively investigated as non-viral carriers of nucleic acid vaccines due to their high transport efficiency,safety,and straightforward production and scalability.However,the molecular mechanism underlying the interactions between nucleic acids and phospholipid bilayers within LNPs remains elusive.In this study,we employed the all-atom molecular dynamics simulation to investigate the interactions between single-stranded nucleic acids and a phospholipid bilayer.Our findings revealed that hydrophilic bases,specifically G in single-stranded RNA(ssRNA)and single-stranded DNA(ssDNA),displayed a higher propensity to form hydrogen bonds with phospholipid head groups.Notably,ssRNA exhibited stronger binding energy than ssDNA.Furthermore,divalent ions,particularly Ca2+,facilitated the binding of ssRNA to phospholipids due to their higher binding energy and lower dissociation rate from phospholipids.Overall,our study provides valuable insights into the molecular mechanisms underlying nucleic acidphospholipid interactions,with potential implications for the nucleic acids in biotherapies,particularly in the context of lipid carriers.

Collective Molecular Machines: Multidimensionality and Reconfigurability

Molecular machines are key to cellular activity where they are involved in converting chemical and light energy into efficient mechanical work.During the last 60 years,designing molecular structures capable of generating unidirectional mechanical motion at the nanoscale has been the topic of intense research.Effective progress has been made,attributed to advances in various fields such as supramolecular chemistry,biology and nanotechnology,and informatics.However,individual molecular machines are only capable of producing nanometer work and generally have only a single functionality.In order to address these problems,collective behaviors realized by integrating several or more of these individual mechanical units in space and time have become a new paradigm.In this review,we comprehensively discuss recent developments in the collective behaviors of molecular machines.In particular,collective behavior is divided into two paradigms.One is the appropriate integration of molecular machines to efficiently amplify molecular motions and deformations to construct novel functional materials.The other is the construction of swarming modes at the supramolecular level to perform nanoscale or microscale operations.We discuss design strategies for both modes and focus on the modulation of features and properties.Subsequently,in order to address existing challenges,the idea of transferring experience gained in the field of micro/nano robotics is presented,offering prospects for future developments in the collective behavior of molecular machines.

Molecular dynamics simulation of the flow mechanism of shear-thinning fluids in a microchannel

Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.

Effect of Molecular Weight on Thermoelectric Performance of P3HT Analogues with 2-Propoxyethyl Side Chains

By replacing hexyl chains in poly(3-hexylthiophene)(P3HT)with 2-propoxyethyls,four poly(3-(2-propoxyethyl)thiophene)(P3POET)homopolymers with comparable polydispersity indexes(PDI)and regioregularities were prepared herein in addition with step increment of about 7 kDa on numberaverage molecular weight(M_(n))from around 11 to 32 kDa(accordingly denoted as P11k,P18k,P25k,and P32k).When doped in film by FeCl_(3)at the optimized conditions,the maximum power factor(PF_(max))increases greatly from 4.3μW·m^(-1)·K^(-2)for P11k to 8.8μW·m^(-1)·K^(-2)for P18k,and further to 9.7μW·m^(-1)·K^(-2)for P25k,followed by a slight decrease to 9.2μW·m^(-1)·K^(-2)for P32k.The close Seebeck coefficients(S)at PF_(max)are observed in these doped polymer films due to their consistent frontier orbital energy levels and Fermi levels.The main contribution to this PF_(max)evolution thus comes from the corresponding conductivities(σ).Theσvariation of the doped films can be rationally correlated with their microstructure evolution.

Molecular features of gastroenteropancreatic neuroendocrine carcinoma: A comparative analysis with lung neuroendocrine carcinoma and digestive adenocarcinomas

Objective: There is an ongoing debate about whether the management of gastroenteropancreatic(GEP)neuroendocrine carcinoma(NEC) should follow the guidelines of small-cell lung cancer(SCLC). We aim to identify the genetic differences of GEPNEC and its counterpart.Methods: We recruited GEPNEC patients as the main cohort, with lung NEC and digestive adenocarcinomas as comparative cohorts. All patients undergone next-generation sequencing(NGS). Different gene alterations were compared and analyzed between GEPNEC and lung NEC(LNEC), GEPNEC and adenocarcinoma to yield the remarkable genes.Results: We recruited 257 patients, including 99 GEPNEC, 57 LNEC, and 101 digestive adenocarcinomas.Among the mutations, KRAS, RB1, TERT, IL7R, and CTNNB1 were found to have different gene alterations between GEPNEC and LNEC samples. Specific genes for each site were revealed: gastric NEC(TERT amplification),colorectal NEC(KRAS mutation), and bile tract NEC(ARID1A mutation). The gene disparities between small-cell NEC(SCNEC) and large-cell NEC(LCNEC) were KEAP1 and CDH1. Digestive adenocarcinoma was also compared with GEPNEC and suggested RB1, APC, and KRAS as significant genes. The TP53/RB1 mutation pattern was associated with first-line effectiveness. Putative targetable genes and biomarkers in GEPNEC were identified in22.2% of the patients, and they had longer progression-free survival(PFS) upon targetable treatment [12.5 months vs. 3.0 months, HR=0.40(0.21-0.75), P=0.006].Conclusions: This work demonstrated striking gene distinctions in GEPNEC compared with LNEC and adenocarcinoma and their clinical utility.

Seed Storability in Rice: Physiological Foundations, Molecular Mechanisms, and Applications in Breeding

Long-term storage of crop seeds is critical for the conservation of germplasm resources, ensuring food supply, and supporting sustainable production. Rice, as a major food staple, has a substantial stock for consumption and production worldwide. However, its food value and seed viability tend to decline during storage. Understanding the physiological responses and molecular mechanisms of aging tolerance forms the basis for enhancing seed storability in rice. This review outlines the latest progress in influential factors, evaluation methods, and identification indices of seed storability. It also discusses the physiological consequences, molecular mechanisms, and strategies for breeding aging-tolerant rice in detail. Finally, it highlights challenges in seed storability research that require future attention. This review offers a theoretical foundation and research direction for uncovering the mechanisms behind seed storability and breeding aging-tolerant rice.

Molecular simulation study on the evolution process of hydrate residual structures into hydrate

The clathrate hydrate memory effect is a fascinating phenomenon with potential applications in carbon capture,utilization and storage(CCUS),gas separation,and gas storage as it can accelerate the secondary formation of clathrate hydrate.However,the underlying mechanism of this effect remains unclear.To gain a better understanding of the mechanism,we conducted molecular dynamic simulations to simulate the initial formation and reformation processes of methane hydrate.In this work,we showed the evolution process of hydrate residual structures into hydrate cages.The simulation results indicate that the residual structures are closely related to the existence of hydrate memory effect,and the higher the contribution of hydrate dissociated water to the hydrate nucleation process,the faster the hydrate nucleation.After hydrate dissociation,the locally ordered structures still exist after hydrate dissociation and can promote the formation of cluster structures,thus accelerating hydrate nucleation.Additionally,the nucleation process of hydrate and the formation process of clusters are inseparable.The size of clusters composed of cup-cage structures is critical for hydrate nucleation.The residence time at high temperature after hydrate decomposition will affect the strength of the hydrate memory effect.Our simulation results provide microscopic insights into the occurrence of the hydrate memory effect and shed light on the hydrate reformation process at the molecular scale.

Temperature effect on nanotwinned Ni under nanoindentation using molecular dynamic simulation

Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.

Microwave irradiation-induced alterations in physicochemical properties and methane adsorption capability of coals:An experimental study using carbon molecular sieve

In order to comprehend the applicability of microwave irradiation for recovering coalbed methane,it is necessary to evaluate the microwave irradiation-induced alterations in coals with varying levels of metamorphism.In this work,the carbon molecular sieve combined with KMnO_(4)oxidation was selected to fabricate carbon molecular sieve with diverse oxidation degrees,which can serve as model substances toward coals.Afterwards,the microwave irradiation dependences of pores,functional groups,and highpressure methane adsorption characteristics of model substances were studied.The results indicated that microwave irradiation causes rearrangement of oxygen-containing functional groups,which could block the micropores with a size of 0.40-0.60 nm in carbon molecular sieve;meanwhile,naphthalene and phenanthrene generated by macro-molecular structure pyrolysis due to microwave irradiation could block the micropores with a size of 0.70-0.90 nm.These alterations in micropore structure weaken the saturated methane adsorption capacity of oxidized carbon molecular sieve by 2.91%-23.28%,suggesting that microwave irradiation could promote methane desorption.Moreover,the increased mesopores found for oxidized carbon molecular sieve after microwave irradiation could benefit CH4 diffusion.In summary,the oxidized carbon molecular sieve can act as model substances toward coals with different ranks.Additionally,microwave irradiation is a promising technology to enhance coalbed methane recovery.

Morphological comparison and molecular marker screening of three Skeletonema species found in Changjiang(Yangtze)River Basin

In our recent investigations of diatom diversity,we studied three species,namely,Skeletonema costatum,Skeletonema subsalsum,and Skeletonema potamos.Although they have been found frequently in Changjiang(Yangtze)River Basin,their morphological and molecular identification is difficult in taxonomy.Therefore,to integrate morphological and molecular biological approaches,we compared systematically their morphological characters and performed phylogenetic analysis.Twelve strains of Skeletonema were collected and isolated from Shanghai and Jiangsu,China,and their morphological characteristics were examined by light microscopy(LM)and the scanning electron microscopy(SEM).Based on morphological comparison,we determined that S.potamos is easy to distinguish from the other two species.The heavily silicified areolae,undulated or cleft distal ends of terminal fultoportula processes(TFPPs),absence of basal pores of fultoportula processes(FPPs),the rootlike protrusions of FPPs,and no interlocking connection are the stable characteristics that can be used to identify S.potamos.However,there are only two features that can distinguish S.costatum from S.subsalsum,namely the location of terminal rimoportulae(TRPs)and the distal shape of TFPPs.In addition,we amplified and sequenced nine common genetic markers from the strains,from which 101 sequences were obtained,constructed phylogenetic trees based on the nine genes and evaluated that seven genes can be used to identify S.potamos,and revealed that S.subsalsum is the closest known relative of S.costatum,and only ATP synthetase beta-subunit gene(atp B)is able to distinguish them from each other,which strongly support that it is an effective molecular marker for Skeletonema.This work provided a theoretical basis for the taxonomic study of Skeletonema.

Analysis of CH_(4) and H_(2) Adsorption on Heterogeneous Shale Surfaces Using aMolecular Dynamics Approach

Determining the adsorption of shale gas on complex surfaces remains a challenge in molecular simulation studies.Difficulties essentially stem from the need to create a realistic shale structure model in terms of mineral heterogeneityand multiplicity.Moreover,precise characterization of the competitive adsorption of hydrogen andmethane in shale generally requires the experimental determination of the related adsorptive capacity.In thisstudy,the adsorption of adsorbates,methane(CH_(4)),and hydrogen(H_(2))on heterogeneous shale surface modelsof Kaolinite,Orthoclase,Muscovite,Mica,C_(60),and Butane has been simulated in the frame of a moleculardynamic’s numerical technique.The results show that these behaviors are influenced by pressure and potentialenergy.On increasing the pressure from 500 to 2000 psi,the sorption effect for CH_(4)significantly increasesbut shows a decline at a certain stage(if compared to H_(2)).The research findings also indicate that raw shalehas a higher capacity to adsorb CH_(4)compared to hydrogen.However,in shale,this difference is negligible.

Optimization of Cooperative RelayingMolecular Communications for Nanomedical Applications

Recently,nano-systems based on molecular communications via diffusion(MCvD)have been implemented in a variety of nanomedical applications,most notably in targeted drug delivery system(TDDS)scenarios.Furthermore,because the MCvD is unreliable and there exists molecular noise and inter symbol interference(ISI),cooperative nano-relays can acquire the reliability for drug delivery to targeted diseased cells,especially if the separation distance between the nano transmitter and nano receiver is increased.In this work,we propose an approach for optimizing the performance of the nano system using cooperative molecular communications with a nano relay scheme,while accounting for blood flow effects in terms of drift velocity.The fractions of the molecular drug that should be allocated to the nano transmitter and nano relay positioning are computed using a collaborative optimization problem solved by theModified Central Force Optimization(MCFO)algorithm.Unlike the previous work,the probability of bit error is expressed in a closed-form expression.It is used as an objective function to determine the optimal velocity of the drug molecules and the detection threshold at the nano receiver.The simulation results show that the probability of bit error can be dramatically reduced by optimizing the drift velocity,detection threshold,location of the nano-relay in the proposed nano system,and molecular drug budget.