An Effective MM/GBSA Protocol for Absolute Binding Free Energy Calculations: A Case Study on SARS-CoV-2 Spike Protein and the Human ACE2 Receptor

The binding free energy calculation of protein–ligand complexes is necessary for research into virus–host interactions and the relevant applications in drug discovery. However, many current computational methods of such calculations are either inefficient or inaccurate in practice. Utilizing implicit solvent models in the molecular mechanics generalized Born surface area (MM/GBSA) framework allows for efficient calculations without significant loss of accuracy. Here, GBNSR6, a new flavor of the generalized Born model, is employed in the MM/GBSA framework for measuring the binding affinity between SARS-CoV-2 spike protein and the human ACE2 receptor. A computational protocol is developed based on the widely studied Ras–Raf complex, which has similar binding free energy to SARS-CoV-2/ACE2. Two options for representing the dielectric boundary of the complexes are evaluated: one based on the standard Bondi radii and the other based on a newly developed set of atomic radii (OPT1), optimized specifically for protein–ligand binding. Predictions based on the two radii sets provide upper and lower bounds on the experimental references: −14.7(ΔGbindBondi)<−10.6(ΔGbindExp.)<−4.1(ΔGbindOPT1) kcal/mol. The consensus estimates of the two bounds show quantitative agreement with the experiment values. This work also presents a novel truncation method and computational strategies for efficient entropy calculations with normal mode analysis. Interestingly, it is observed that a significant decrease in the number of snapshots does not affect the accuracy of entropy calculation, while it does lower computation time appreciably. The proposed MM/GBSA protocol can be used to study the binding mechanism of new variants of SARS-CoV-2, as well as other relevant structures.     

Zhang–Zhang Polynomials of Multiple Zigzag Chains Revisited: A Connection with the John–Sachs Theorem

Multiple zigzag chains Zm,n of length n and width m constitute an important class of regular graphene flakes of rectangular shape. The physical and chemical properties of these basic pericondensed benzenoids can be related to their various topological invariants, conveniently encoded as the coefficients of a combinatorial polynomial, usually referred to as the ZZ polynomial of multiple zigzag chains Zm,n. The current study reports a novel method for determination of these ZZ polynomials based on a hypothesized extension to John–Sachs theorem, used previously to enumerate Kekulé structures of various benzenoid hydrocarbons. We show that the ZZ polynomial of the Zm,n multiple zigzag chain can be conveniently expressed as a determinant of a Toeplitz (or almost Toeplitz) matrix of size m2×m2 consisting of simple hypergeometric polynomials. The presented analysis can be extended to generalized multiple zigzag chains Zkm,n, i.e., derivatives of Zm,n with a single attached polyacene chain of length k. All presented formulas are accompanied by formal proofs. The developed theoretical machinery is applied for predicting aromaticity distribution patterns in large and infinite multiple zigzag chains Zm,n and for computing the distribution of spin densities in biradical states of finite multiple zigzag chains Zm,n.     

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.     

Cross-Talk between Overlap Interactions in Biomolecules: A Case Study of the β-Turn Motif

Noncovalent interactions play a pivotal role in regulating protein conformation, stability and dynamics. Among the quantum mechanical (QM) overlap-based noncovalent interactions, n→π* is the best understood with studies ranging from small molecules to β-turns of model proteins such as GB1. However, these investigations do not explore the interplay between multiple overlap interactions in contributing to local structure and stability. In this work, we identify and characterize all noncovalent overlap interactions in the β-turn, an important secondary structural element that facilitates the folding of a polypeptide chain. Invoking a QM framework of natural bond orbitals, we demonstrate the role of several additional interactions such as n→σ* and π→π* that are energetically comparable to or larger than n→π*. We find that these interactions are sensitive to changes in the side chain of the residues in the β-turn of GB1, suggesting that the n→π* may not be the only component in dictating β-turn conformation and stability. Furthermore, a database search of n→σ* and π→π* in the PDB reveals that they are prevalent in most proteins and have significant interaction energies (∼1 kcal/mol). This indicates that all overlap interactions must be taken into account to obtain a comprehensive picture of their contributions to protein structure and energetics. Lastly, based on the extent of QM overlaps and interaction energies, we propose geometric criteria using which these additional interactions can be efficiently tracked in broad database searches.     

Carbon Atoms Speaking Out: How the Geometric Sensitivity of 13C Chemical Shifts Leads to Understanding the Colour Tuning of Phycocyanobilin in Cph1 and AnPixJ

We present a combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics–statistical approach for the interpretation of nuclear magnetic resonance (NMR) chemical shift patterns in phycocyanobilin (PCB). These were originally associated with colour tuning upon photoproduct formation in red/green-absorbing cyanobacteriochrome AnPixJg2 and red/far-red-absorbing phytochrome Cph1Δ2. We pursue an indirect approach without computation of the absorption frequencies since the molecular geometry of cofactor and protein are not accurately known. Instead, we resort to a heuristic determination of the conjugation length in PCB through the experimental NMR chemical shift patterns, supported by quantum chemical calculations. We have found a characteristic correlation pattern of 13C chemical shifts to specific bond orders within the π-conjugated system, which rests on the relative position of carbon atoms with respect to electron-withdrawing groups and the polarisation of covalent bonds. We propose the inversion of this regioselective relationship using multivariate statistics and to apply it to the known experimental NMR chemical shifts in order to predict changes in the bond alternation pattern. Therefrom the extent of electronic conjugation, and eventually the change in absorption frequency, can be derived. In the process, the consultation of explicit mesomeric formulae plays an important role to qualitatively account for possible conjugation scenarios of the chromophore. While we are able to consistently associate the NMR chemical shifts with hypsochromic and bathochromic shifts in the Pg and Pfr, our approach represents an alternative method to increase the explanatory power of NMR spectroscopic data in proteins.     

Kinetic Study on Alpha-Form Crystallization of Mixed-Acid Triacylglycerols POP,PPO, and Their Mixture

The crystallization behavior of the metastable α form of triacylglycerols (TAGs) plays a critical role as a precursor for the crystallization of more stable β′ and β forms for various applications in food and pharmaceutical products. However, precise analysis of the crystallization kinetics of α has not been performed, likely due to its rapid and complex behavior. This paper presents the observation results of the initial stages of the isothermal crystallization kinetics of α forms of 1,3-dipalmitoyl-2-oleoyl-glycerol (POP), 1,2-dipalmitoyl-3-oleoyl-rac-glycerol (rac-PPO), and molecular compound (MC) crystals of a POP/rac-PPO (1/1) mixture (MCPOP/PPO) using synchrotron radiation time-resolved X-ray diffraction and polarized optical microscopy. In all the TAGs, α crystals with a worm-like morphology started to grow rapidly in the first stage. Then, the α crystals slowly transformed into more stable forms in different manners for different TAG samples. In POP, the conversion was simple, as the α-2 form transformed into γ-3, whereas in rac-PPO, the lamellar distance values of the α-2 form continuously decreased with time and changed into the α-3 form. In the MCPOP/PPO crystals, in contrast, separate crystallization of α-2 of a rac-PPO fraction initially occurred, followed by the crystallization of α-2 of POP, and the two α forms merged into α-2 of MCPOP/PPO. This separate crystallization was caused by large differences in the crystallization kinetics of the α forms of POP and rac-PPO.     

Understanding Topological Insulators in Real Space

A real space understanding of the Su–Schrieffer–Heeger model of polyacetylene is introduced thanks to delocalization indices defined within the quantum theory of atoms in molecules. This approach enables to go beyond the analysis of electron localization usually enabled by topological insulator indices—such as IPR—enabling to differentiate between trivial and topological insulator phases. The approach is based on analyzing the electron delocalization between second neighbors, thus highlighting the relevance of the sublattices induced by chiral symmetry. Moreover, the second neighbor delocalization index, δi,i+2, also enables to identify the presence of chirality and when it is broken by doping or by eliminating atom pairs (as in the case of odd number of atoms chains). Hints to identify bulk behavior thanks to δ1,3 are also provided. Overall, we present a very simple, orbital invariant visualization tool that should help the analysis of chirality (independently of the crystallinity of the system) as well as spreading the concepts of topological behavior thanks to its relationship with well-known chemical concepts.     

Electrochemical Resistive DNA Biosensor for the Detection of HPV Type 16

In this work, a low-cost and rapid electrochemical resistive DNA biosensor based on the current relaxation method is described. A DNA probe, complementary to the specific human papillomavirus type 16 (HPV-16) sequence, was immobilized onto a screen-printed gold electrode. DNA hybridization was detected by applying a potential step of 30 mV to the system, composed of an external capacitor and the modified electrode DNA/gold, for 750 µs and then relaxed back to the OCP, at which point the voltage and current discharging curves are registered for 25 ms. From the discharging curves, the potential and current relaxation were evaluated, and by using Ohm’s law, the charge transfer resistance through the DNA-modified electrode was calculated. The presence of a complementary sequence was detected by the change in resistance when the ssDNA is transformed in dsDNA due to the hybridization event. The target DNA concentration was detected in the range of 5 to 20 nM. The results showed a good fit to the regression equation ΔRtotal(Ω)=2.99 × [DNA]+81.55, and a detection limit of 2.39 nM was obtained. As the sensing approach uses a direct current, the electronic architecture of the biosensor is simple and allows for the separation of faradic and nonfaradaic contributions. The simple electrochemical resistive biosensor reported here is a good candidate for the point-of-care diagnosis of HPV at a low cost and in a short detection time.     

Ca+ Ions Solvated in Helium Clusters

We present a combined experimental and theoretical investigation on Ca+ ions in helium droplets, HeNCa+. The clusters have been formed in the laboratory by means of electron-impact ionization of Ca-doped helium nanodroplets. Energies and structures of such complexes have been computed using various approaches such as path integral Monte Carlo, diffusion Monte Carlo and basin-hopping methods. The potential energy functions employed in these calculations consist of analytical expressions following an improved Lennard-Jones formula whose parameters are fine-tuned by exploiting ab initio estimations. Ion yields of HeNCa+ -obtained via high-resolution mass spectrometry- generally decrease with N with a more pronounced drop between N=17 and N=25, the computed quantum HeNCa+ evaporation energies resembling this behavior. The analysis of the energies and structures reveals that covering Ca+ with 17 He atoms leads to a cluster with one of the smallest energies per atom. As new atoms are added, they continue to fill the first shell at the expense of reducing its stability, until N=25, which corresponds to the maximum number of atoms in that shell. Behavior of the evaporation energies and radial densities suggests liquid-like cluster structures.     

TiCl4 Dissolved in Ionic Liquid Mixtures from Аb Initio Molecular Dynamics Simulations

To gain a deeper understanding of the TiCl4 solvation effects in multi-component ionic liquids, we performed ab initio molecular dynamics simulations of 1-butyl-3-methylimidazolium [C4C1Im]+, tetrafluoroborate [BF4]−, chloride [Cl]− both with and without water and titanium tetrachloride TiCl4. Complex interactions between cations and anions are observed in all investigated systems. By further addition of water and TiCl4 this complex interaction network is extended. Observations of the radial distribution functions and number integrals show that water and TiCl4 not only compete with each other to interact mainly with [Cl]−, which strongly influences the cation-[BF4]− interaction, but also interact with each other, which leads to the fact that in certain systems the cation-anion interaction is enhanced. Further investigations of the Voronoi polyhedra analysis have demonstrated that water has a greater impact on the nanosegregated system than TiCl4 which is also due to the fact of the shear amount of water relative to all other components and its higher mobility compared to TiCl4. Overall, the polar network of the IL mixture collapses by including water and TiCl4. In the case of [Cl]− chloride enters the water continuum, while [BF4]− remains largely unaffected, which deeply affects the interaction of the ionic liquid (IL) network.     

Insight into the Nucleation Mechanism of p-Methoxybenzoic Acid in Ethanol-Water System from Metastable Zone Width

The metastable zone width (MSZW) of p-methoxybenzoic acid (PMBA) in an ethanol-water system was measured using the polythermal method. The nucleation order m obtained by the Nývlt’s model indicates the nucleation of PMBA following a progressive nucleation mechanism at low saturation temperature (m = 3.18–7.50) and an instantaneous nucleation mechanism at high saturation temperature (m = 1.46–2.55). Then, combined with the metastable zone experiment and the Sangwal model, we found that the MSZW and the interfacial energy reached the maximum when the mass fraction of ethanol was 0.8, which resulted in the smallest crystal product size. Meanwhile, the maximum rcrit and ΔGcrit obtained based on the modified Sangwal model indicating the PMBA needs to overcome a higher nucleation barrier in the ethanol mass fraction of 0.8. Finally, we proposed a preferential strategy for adjusting MSZW by correlating the interfacial energy with the change in ethanol mass fraction, saturation temperature, and cooling rate, respectively.     

Significance of Shape Factor in Heat Transfer Performance of Molybdenum-Disulfide Nanofluid in Multiple Flow Situations; A Comparative Fractional Study

In this modern era, nanofluids are considered one of the advanced kinds of heat transferring fluids due to their enhanced thermal features. The present study is conducted to investigate that how the suspension of molybdenum-disulfide (MoS2) nanoparticles boosts the thermal performance of a Casson-type fluid. Sodium alginate (NaAlg) based nanofluid is contained inside a vertical channel of width d and it exhibits a flow due to the movement of the left wall. The walls are nested in a permeable medium, and a uniform magnetic field and radiation flux are also involved in determining flow patterns and thermal behavior of the nanofluid. Depending on velocity boundary conditions, the flow phenomenon is examined for three different situations. To evaluate the influence of shape factor, MoS2 nanoparticles of blade, cylinder, platelet, and brick shapes are considered. The mathematical modeling is performed in the form of non-integer order operators, and a double fractional analysis is carried out by separately solving Caputo-Fabrizio and Atangana-Baleanu operators based fractional models. The system of coupled PDEs is converted to ODEs by operating the Laplace transformation, and Zakian’s algorithm is applied to approximate the Laplace inversion numerically. The solutions of flow and energy equations are presented in terms of graphical illustrations and tables to discuss important physical aspects of the observed problem. Moreover, a detailed inspection on shear stress and Nusselt number is carried out to get a deep insight into skin friction and heat transfer mechanisms. It is analyzed that the suspension of MoS2 nanoparticles leads to ameliorating the heat transfer rate up to 9.5%. To serve the purpose of achieving maximum heat transfer rate and reduced skin friction, the Atangana-Baleanu operator based fractional model is more effective. Furthermore, it is perceived that velocity and energy functions of the nanofluid exhibit significant variations because of the different shapes of nanoparticles.     

Studies of the Formation of Inclusion Complexes Derivatives of Cinnamon Acid with α-Cyclodextrin in a Wide Range of Temperatures Using Conductometric Methods

The electrical conductivities of aqueous solutions of sodium salts of trans-4-hydroxycinnamic acid (trans-p-coumaric acid), trans-3,4-dihydroxycinnamic acid (trans-caffeic acid), trans-4-hydroxy-3-methoxycinnamic acid, (trans-ferulic acid) and trans-3-phenylacrylic acid (trans-cinnamic acid) with α-cyclodextrin were measured in the temperature range of 288.15 K–318.15 K. For the first time in the literature, using the limiting molar conductivity (Λmo) obtained from conductivity measurements, the values of the complexation constants (K_f) of the salts of phenolic acid derivatives with α-cyclodextrin were determined using a modified low concentration chemical model (IcCM). An attempt was also made to analyze the individual thermodynamic functions ΔGo, ΔHo and ΔSo describing the complexation process as a function of temperature changes. The obtained results show that the process of formation of inclusion complexes is exothermic and is spontaneous.     

Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative

Human immunodeficiency virus (HIV) is a life life-threatening and serious infection caused by a virus that attacks CD4+ T-cells, which fight against infections and make a person susceptible to other diseases. It is a global public health problem with no cure; therefore, it is highly important to study and understand the intricate phenomena of HIV. In this article, we focus on the numerical study of the path-tracking damped oscillatory behavior of a model for the HIV infection of CD4+ T-cells. We formulate fractional dynamics of HIV with a source term for the supply of new CD4⁺ T-cells depending on the viral load via the Caputo–Fabrizio derivative. In the formulation of fractional HIV dynamics, we replaced the constant source term for the supply of new CD4+ T-cells from the thymus with a variable source term depending on the concentration of the viral load, and introduced a term that describes the incidence of the HIV infection of CD4+ T-cells. We present a novel numerical scheme for fractional view analysis of the proposed model to highlight the solution pathway of HIV. We inspect the periodic and chaotic behavior of HIV for the given values of input factors using numerical simulations.     

Development of a Fast and Robust UHPLC Method for Apixaban In-Process Control Analysis

In-process control (IPC) is an important task during chemical syntheses in pharmaceutical industry. Despite the fact that each chemical reaction is unique, the most common analytical technique used for IPC analysis is high performance liquid chromatography (HPLC). Today, the so-called “Quality by Design” (QbD) principle is often being applied rather than “Trial and Error” approach for HPLC method development. The QbD approach requires only for a very few experimental measurements to find the appropriate stationary phase and optimal chromatographic conditions such as the composition of mobile phase, gradient steepness or time (tG), temperature (T), and mobile phase pH. In this study, the applicability of a multifactorial liquid chromatographic optimization software was studied in an extended knowledge space. Using state-of-the-art ultra-high performance liquid chromatography (UHPLC), the analysis time can significantly be shortened. By using UHPLC, it is possible to analyse the composition of the reaction mixture within few minutes. In this work, a mixture of route of synthesis of apixaban was analysed on short narrow bore column (50 × 2.1 mm, packed with sub-2 µm particles) resulting in short analysis time. The aim of the study was to cover a relatively narrow range of method parameters (tG, T, pH) in order to find a robust working point (zone). The results of the virtual (modeled) robustness testing were systematically compared to experimental measurements and Design of Experiments (DoE) based predictions.     

Two Intercalation Mechanisms of Oxazole Yellow Dimer (YOYO-1) into DNA

The oxazole yellow dye, YOYO-1 (a symmetric homodimer), is a commonly used molecule for staining DNA. We applied the brightness analysis to study the intercalation of YOYO-1 into the DNA. We distinguished two binding modes of the dye to dsDNA: mono-intercalation and bis-intercalation. Bis-intercalation consists of two consecutive mono-intercalation steps, characterised by two distinct equilibrium constants (with the average number of base pair per binding site equals 3.5): K1=3.36±0.43×107M−1 and K2=1.90±0.61×105M−1, respectively. Mono-intercalation dominates at high concentrations of YOYO-1. Bis-intercalation occurs at low concentrations.     

Magnetic-Field-Orientation Dependent Thermal Entanglement of a Spin-1 Heisenberg Dimer: The Case Study of Dinuclear Nickel Complex with an Uniaxial Single-Ion Anisotropy

The bipartite entanglement in pure and mixed states of a quantum spin-1 Heisenberg dimer with exchange and uniaxial single-ion anisotropies is quantified through the negativity in a presence of the external magnetic field. At zero temperature the negativity shows a marked stepwise dependence on a magnetic field with two abrupt jumps and plateaus, which can be attributed to the quantum antiferromagnetic and quantum ferrimagnetic ground states. The magnetic-field-driven phase transition between the quantum antiferromagnetic and quantum ferrimagnetic ground states manifests itself at nonzero temperatures by a local minimum of the negativity, which is followed by a peculiar field-induced rise of the negativity observable in a range of moderately strong magnetic fields. The rising temperature generally smears out abrupt jumps and plateaus of the negativity, which cannot be distinguished in the relevant dependencies above a certain temperature. It is shown that the thermal entanglement is most persistent against rising temperature at the magnetic field, for which an energy gap between a ground state and a first excited state is highest. Besides, temperature variations of the negativity of the spin-1 Heisenberg dimer with an easy-axis single-ion anisotropy may exhibit a singular point-kink, at which the negativity has discontinuity in its first derivative. The homodinuclear nickel complex [Ni2(Medpt)2(μ-ox)(H2O)2](ClO4)2·2H2O provides a suitable experimental platform of the antiferromagnetic spin-1 Heisenberg dimer, which allowed us to estimate a strength of the bipartite entanglement between two exchange-coupled Ni2+ magnetic ions on the grounds of the interaction constants reported previously from the fitting procedure of the magnetization data. It is verified that the negativity of this dinuclear compound is highly magnetic-field-orientation dependent due to presence of a relatively strong uniaxial single-ion anisotropy.     

Catechin Photolysis Suppression by Aluminum Chloride under Alkaline Conditions and Assessment with Liquid Chromatography–Mass Spectrometry

Tea is rich in catechins and aluminum. In this study, the process of catechin photolysis was applied as a model for examining the effects of aluminum chloride (AlCl₃) on the structural changes of catechin and the alteration of aluminum complexes under blue light irradiation (BLI) at pH 8 using liquid chromatography and mass spectrometry techniques. Additionally, the effects of anions on catechin upon the addition of AlCl₃ and treatment with BLI were also studied. In this study, when 1 mM catechin was treated with BLI, a superoxide anion radical (O2•−) was generated in an air-saturated aqueous solution, in addition to forming a dimeric catechin (proanthocyanidin) via a photon-induced redox reaction. The relative percentage of catechin was found to be 59.0 and 95.7 for catechin treated with BLI and catechin upon the addition of 1 mM AlCl₃ treated with BLI, respectively. It suggested that catechin treated with BLI could be suppressed by AlCl₃, while AlCl₃ did not form a complex with catechin in the photolytic system. However, under the same conditions, it was also found that the addition of AlCl₃ inhibited the photolytic formation of O2•−, and reduced the generation of proanthocyanidin, suggesting that the disconnection of proanthocyanidin was achieved by AlCl₃ acting as a catalyst under treatment with BLI. The influence of 1 mM fluoride (F−) and 1 mM oxalate (C2O42−) ions on the photolysis of 1 mM catechin upon the addition of 1 mM AlCl₃ and treatment with BLI was found to be insignificant, implying that, during the photolysis of catechin, the Al species were either neutral or negatively charged and the aluminum species did not form a complex with anions in the photolytic system. Therefore, aluminum, which is an amphoteric species, has an inherent potential to stabilize the photolysis of catechin in an alkaline conditions, while suppressing the O2•− and proanthocyanidin generation via aluminum ion catalysis in the catechin/Al system under treatment with BLI.     

Determination of Polycyclic Aromatic Hydrocarbons and Their Methylated Derivatives in Sewage Sludge from Northeastern China: Occurrence,Profiles and Toxicity Evaluation

This paper assesses the occurrence, distribution, source, and toxicity of polycyclic aromatic hydrocarbons (PAHs), and their methylated form (Me-PAHs) in sewage sludge from 10 WWTPs in Northeastern China was noted. The concentrations of ∑PAHs, ∑Me-PAHs ranged from 567 to 5040 and 48.1 to 479ng.g−1dw, which is greater than the safety limit for sludge in agriculture in China. High and low molecular weight 4 and 2-ring PAHs and Me-PAHs in sludge were prevalent. The flux of sludge PAHs and Me-PAHs released from ten WWTPs, in Heilongjiang province, was calculated to be over 100 kg/year. Principal component analysis (PCA), diagnostic ratios and positive matrix factorization (PMF) determined a similar mixed pyrogenic and petrogenic source of sewage sludge. The average values of Benzo[a]pyrene was below the safe value of 600ng.g−1 dependent on an incremental lifetime cancer risk ILCR of 10⁻⁶. Sludge is an important source for the transfer of pollutants into the environment, such as PAHs and Me-PAHs. Consequently, greater consideration should be given to its widespread occurrence.     

The Flow of Blood-Based Hybrid Nanofluids with Couple Stresses by the Convergent and Divergent Channel for the Applications of Drug Delivery

This research work aims to scrutinize the mathematical model for the hybrid nanofluid flow in a converging and diverging channel. Titanium dioxide and silver TiO2 and Ag are considered as solid nanoparticles while blood is considered a base solvent. The couple-stress fluid model is essentially use to describe the blood flow. Therefore, the couple-stress term was used in the recent study with the existence of a magnetic field and a Darcy–Forchheiner porous medium. The heat absorption/omission and radiation terms were also included in the energy equation for the sustainability of drug delivery. An endeavor was made to link the recent study with the applications of drug delivery. It has already been revealed by the available literature that the combination of TiO2 with any other metal can destroy cancer cells more effectively than TiO2 separately. Both the walls are stretchable/shrinkable, whereas flow is caused by a source or sink with α as a converging/diverging parameter. Governing equations were altered into the system of non-linear coupled equations by using the similarity variables. The homotopy analysis method (HAM) was applied to obtain the preferred solution. The influences of the modeled parameters have been calculated and displayed. The confrontation of wall shear stress and hybrid nanofluid flow increased as the couple stress parameter rose, which indicates an improvement in the stability of the base fluid (blood). The percentage (%) increase in the heat transfer rate with the variation of nanoparticle volume fraction was also calculated numerically and discussed theoretically.