Results Of Molecular Dynamics Simulations Research Articles

Biomimetic management of bone healing stages: MOFs induce tunable degradability and enhanced angiogenesis-osteogenesis coupling

Currently, biodegradable biomaterials are extensively used in bone regeneration. Nevertheless, a significant challenge remains in accurately adjusting their degradation rates to match bone healing stages to achieve bone regeneration. Herein, we present a novel biodegradable bone grafts integrated with zeolite imidazole framework-8 (ZIF-8) nanocrystals, fabricated via microinjection molding. Apart from the reinforcing effect derived from the in situ fibrillation of PCL in PLGA, nanoscale ZIF-8 is inspiringly demonstrated to control the PLGA degradation. Both molecular dynamic simulation and experimental results indicate that Zn2+ derived from ZIF-8 nanocrystals acts as a Lewis acid to participate the degradation of PLGA matrix, thus successfully reconciles the contradiction between degradation of matrix and maintenance of reinforcing phase. Moreover, nanometric ZIF-8 modified bone grafts improve the angiogenesis-osteogenesis coupling. In vitro, it directly boosts HIF-1α and capillary-like network formation of endothelial cells (ECs), and RUNX2 and OCN expression in bone marrow mesenchymal stem cells (BMSCs). In vivo, bone grafts with a moderate concentration (0.5 %) of ZIF-8 nanoparticles aid vascularized bone regeneration by tunable biodegradation matching bone healing stages. Further, RNA sequencing shows PLGA/PCL@nZ bone grafts promote angiogenesis-osteogenesis by cooperatively regulating calcium signaling pathway in both ECs and BMSCs. Collectively, this study delineates a groundbreaking bone repair strategy with profound clinical application potential.

A systematic mutation analysis of 13 major SARS-CoV-2 variants

SARS-CoV-2 evolves constantly with various novel mutations. Due to their enhanced infectivity, transmissibility and immune evasion, a comprehensive understanding of the association between these mutations and the respective functional changes is crucial. However, previous mutation studies of major SARS-CoV-2 variants remain limited. Here, we performed systematic analyses of full-length amino acids mutation, phylogenetic features, protein physicochemical properties, molecular dynamics and immune escape as well as pseudotype virus infection assays among thirteen major SARS-CoV-2 variants. We found that Omicron exhibited the most abundant and complex mutation sites, higher indices of hydrophobicity and flexibility than other variants. The results of molecular dynamics simulation suggest that Omicron has the highest number of hydrogen bonds and strongest binding free energy between the S protein and ACE2 receptor. Furthermore, we revealed 10 immune escape sites in 13 major variants, some of them were reported previously, but four of which (i.e. 339/373/477/496) are first reported to be specific to Omicron, whereas 462 is specific to Epslion. The infectivity of these variants was confirmed by the pseudotype virus infection assays. Our findings may help us understand the functional consequences of the mutations within various variants and the underlying mechanisms of the immune escapes conferred by the S proteins.

Thermodynamic and transport properties of triangular-well fluids from molecular dynamics

Thermodynamic and transport properties of particles interacting via the triangular-well potential in the range of 1.4 ≤ λ ≤ 2.6 have been studied using perturbation theory and Molecular Dynamics simulations. We present thermodynamic properties such as vapour-liquid coexistence, vapour pressures, and density for different values of the attractive range potential. Good agreement is observed between the theoretical approach and computer simulations in a wide range of densities, temperatures, and pressures. Additionally, we show Molecular Dynamic simulation results of transport properties such as the self-diffusion coefficient as a function of both density and temperature. Finally, thermodynamic and transport properties of real fluids like oxygen, methane, hydrogen sulfide, and fluoromethane have been predicted using both approaches. Results of vapour-liquid coexistence, vapour pressures, and critical points were in excellent agreement with experimental data. The coupling of perturbation theory, Molecular Dynamics simulations, and the triangular-well pair potential offers a versatile and efficient method to explore various real fluid systems, which would include associating fluids, complex molecular liquids, and colloidal systems. This approach can advance our understanding of these systems using a simple approach and to contribute the development of more accurate and predictive models.

Correlation of rare earth coordination and spectral properties in Er3+ doped Na2O–Al2O3–SiO2 glasses with different Al2O3 concentrations by molecular dynamics simulations

The molecular structure of Er2O3 doped Na2O–Al2O3–SiO2 glasses with varying Na2O/Al2O3 ratios is explored via molecular dynamics (MD) simulations using the so-called inherent structure sampling method, which allows the calculation of a large number of local structures of low concentration, as needed to determine the surrounding of low concentration dopants. General structural parameters, including radial distribution functions, coordination numbers and interatomic distances of all network forming and network modifying ions are reported. However, in this work, special attention is devoted to the effect of Al2O3 concentration on the local surrounding of the doped Er3+ ions. It is shown that the Er atoms coordinate 5–6 oxygen ions in their first coordination sphere in the investigated glasses. The Er–O coordination number increases monotonically with increasing Al2O3 concentration and decreasing Na2O/Al2O3 ratio. It is found that the Er atoms are preferably connected to non-bridging oxygen atoms (NBO) in all glasses, even in the peraluminous composition. Additionally, the MD simulation results are compared to the glasses spectral properties that were already investigated in detail by Tanabe and Hanada. The increasing Er-NBO coordination number derived by MD simulations could be correlated with the increased peak splitting of the Er3+ absorption peaks reported by Tanabe and Hanada. Furthermore, a correlation between the Judd-Ofelt parameters published by Tanabe and Hanada and the Er3+ coordination in the glass structure is discussed. It is shown that the Er–O coordination increases with increasing Ω2 parameter as the Al2O3 concentration increases in the glass composition. A correlation of the average overall Er–O coordination number with the symmetry of the local Er site is proposed.

Topological isomerism design of Demulsifiers: A fusion of molecular dynamics simulation and experiment

Demulsifiers are crucial for enhancing crude oil extraction and transportation efficiency by facilitating the separation of oil–water emulsions. However, research conducted on the design and synthesis of demulsifiers based on the results of molecular dynamics simulation remains scarce. Herein, we designed eight block polyethers, and explored their interfacial behaviors at the oil–water interface via molecular dynamics simulation. Correspondingly, eight block polyether demulsifiers were synthesized, and their interface properties were tested. The simulation results were consistent with the experimental results, indicating that the energy of simulation systems could be efficiently reduced by designed molecules, and that the prepared demulsifiers with the same structure could achieve effective demulsification correspondingly. Based on the simulation and experimental results, a “hand-shaped” demulsification mechanism model was proposed. The demulsifiers with hydrophobic “palm” and hydrophilic “fingers” demonstrated a more efficient performance in demulsification. The hydrophobic “palm” was stable in the continuous oil phase, and the hydrophilic “fingers” could attract dispersed water droplets to aggregate, facilitating demulsification. The study also finds that longer hydrophilic “fingers” exert stronger molecular force. Combining molecular dynamics simulation and experiment significantly accelerates the development of new demulsifiers, highlighting the utility of this approach in the research on demulsifiers.

New analytical thermodynamic models developed for pure H2, CH4, CO2 and H2 containing mixtures based on molecular simulations

Pressure-volume-temperature-composition (PVTX) properties of gas mixtures are important to model and predict flow and compression performance in a range of industrial processes. It is time consuming and sometimes impractical to obtain these properties from experiments. Here, molecular dynamics (MD) simulations were conducted to compute the thermodynamic volumetric properties of pure H2, CH4, CO2, and binary mixtures of H2/CH4 and H2/CO2 at temperatures of 310.9–470 K, pressures up to 172 MPa. MD simulation results were used to develop five analytical equations explicit in pressure as a function of density or volume, temperature, and composition, for the pure gases and mixtures. Our simulated results matched very well with the predictions from the established GERG-2008 equation and data published by NIST. The simplified equations of state (EOS) developed in this work will be more computationally efficient than other EOS for mixed systems. This can be significant in reservoir simulations.

Mechanism of Cordyceps militaris in gouty nephropathy explored using network pharmacology and molecular docking technology

Cordyceps militaris Link. (C. militaris) with strong anti-hyperuricemia effects could be detected in mice urine and serum after oral administration as we previously proved. In this study, the mechanism and molecular effects of anti-hyperuricemia by C. militaris in gouty nephropathy (GN) mice were firstly investigated, and its key active structure was further explored. The results showed C. militaris could significantly reduce the uric acid (UA) level and related inflammatory factors. Meanwhile, C. militaris maintain the normal dynamic balance of UA by inhibiting hepatic xanthine oxidase (XO) and regulating the expression of renal transport proteins ABCG2, GLUT9 and URAT1. Then we employed network pharmacology to theoretically determine two active components in C. militaris, cordycepin and ergosterol, that act on GN. Moreover, the molecular docking and molecular dynamic simulation results showed cordycepin and ergosterol could stably bind to XO and cyclooxygenase-2 (COX-2) due to hydrogen bonding and hydrophobic interactions. This study could provide theoretical evidence for utilization of C. militaris in hyperuricemia-management functional foods.

Effect of voltage gradients on EK-PRB remediation: Experimental and molecular dynamics simulations

Electrokinetic-Permeable Reaction Barrier (EK-PRB) coupled remediation technology can effectively treat heavy metal-contaminated soil near coal mines. This study was conducted on cadmium (Cd), a widely present element in the soil of the mining area. To investigate the impact of the voltage gradient on the remediation effect of EK-PRB, the changes in current, power consumption, pH, and Cd concentration content during the macroscopic experiment were analyzed. A three-dimensional visualized kaolinite-heavy metal-water simulation system was constructed and combined with the Molecular Dynamics (MD) simulations to elucidate the migration mechanism and binding active sites of Cd on the kaolinite (001) crystalline surface at the microscopic scale. The results showed that the voltage gradient positively correlates with the current, power consumption, and Cd concentration during EK-PRB remediation, and the average removal efficiency increases non-linearly with increasing voltage gradient. Considering power consumption, average removal efficiency, and cost-effectiveness, the voltage range is between 1.5 and 3.0 V/cm, with 2.5 V/cm being the optimal value. The results of MD simulations and experiments correspond to each other. Cd2+ formed a highly stable adsorption structure in contrast to the Al–O sheet on the kaolinite (001) crystalline surface. The mean square displacement (MSD) curve of Cd2+ under the electric field exhibits anisotropy, the total diffusion coefficient DTotal increases and the Cd2+ migration rate accelerates. The electric field influences the microstructure of Cd2+ complexes. With the enhancement of the voltage gradient, the complexation between Cd2+ and water molecules is enhanced, and the interaction between Cd2+ and Cl− in solution is weakened.

Antigen specific VNAR screening in whitespotted bamboo shark (Chiloscyllium plagiosum) with next generation sequencing

IgNAR exhibits significant promise in the fields of cancer and anti-virus biotherapies. Notably, the variable regions of IgNAR (VNAR) possess comparable antigen binding affinity with much smaller molecular weight (∼12 kDa) compared to IgNAR. Antigen specific VNAR screening is a changeling work, which limits its application in medicine and therapy fields. Though phage display is a powerful tool for VNAR screening, it has a lot of drawbacks, such as small library coverage, low expression levels, unstable target protein, complicating and time-consuming procedures. Here we report VANR screening with next generation sequencing (NGS) could effectively overcome the limitations of phage display, and we successfully identified approximately 3000 BAFF-specific VNARs in Chiloscyllium plagiosum vaccinated with the BAFF antigen. The results of modelling and molecular dynamics simulation and ELISA assay demonstrated that one out of the top five abundant specific VNARs exhibited higher binding affinity to the BAFF antigen than those obtained through phage display screening. Our data indicates NGS would be an alternative way for VNAR screening with plenty of advantages.

Structural Pharmacology of TRPV4 Antagonists.

The nonselective calcium-permeable Transient Receptor Potential Cation Channel Subfamily V Member4 (TRPV4)channel regulates various physiological activities. Dysfunction of TRPV4 is linked to many severe diseases, including edema, pain, gastrointestinal disorders, lung diseases, and inherited neurodegeneration. Emerging TRPV4 antagonists show potential clinical benefits. However, the molecular mechanisms of TRPV4 antagonism remain poorly understood. Here, cryo-electron microscopy (cryo-EM) structures of human TRPV4 are presented in-complex with two potent antagonists, revealing the detailed binding pockets and regulatory mechanisms of TRPV4 gating. Both antagonists bind to the voltage-sensing-like domain (VSLD) and stabilize the channel in closed states. These two antagonists induce TRPV4 to undergo an apparent fourfold to twofold symmetry transition. Moreover, it is demonstrated that one of the antagonists binds to the VSLD extended pocket, which differs from the canonical VSLD pocket. Complemented with functional and molecular dynamics simulation results, this study provides crucial mechanistic insights into TRPV4 regulation by small-molecule antagonists, which may facilitate future drug discovery targeting TRPV4.

High‐Strength Wear‐Resistant Cu/Tungsten Carbide/Diamond Composites Fabricated by Powder Metallurgy

Copper‐based diamond composites are widely used for thermal management and wear‐resistant materials. In this work, Cu/diamond, Cu/tungsten carbide (WC)/diamond, and Cu–0.92Cr/WC/diamond composites are fabricated by high‐energy ball milling and rapid hot‐pressing sintering. Physical characteristics, compressive strength, and tribological properties of the studied composites are investigated. The compressive strength of Cu/7wt%WC/diamond composites is 102% higher than the Cu/diamond composites, reaching 415 MPa. Compression fractures extend from the interface due to compression cracks, consistent with the compression curves from molecular dynamics (MD) simulations. The crack propagates from the interface during the compression experiment, which is consistent with the result of MD simulation. The Cu/7wt%WC/diamond composites exhibit a wear rate of 0.6 × 10−6 mm3 (N m)−1 and a friction coefficient of 0.37 at 50 N. The addition of WC and diamond improves the compressive strength and wear resistance. These findings are helpful in the development of copper composites with high compression strength and wear resistance.

The crosslinking sites and molecular conformation of gelatin hydrogel modified by transglutaminase

This work aimed to investigate the covalent crosslinking sites induced by transglutaminase (TGase, EC 2.3.2.13) in gelatin hydrogel and their impact on the gelatin molecular conformation. LC-MS/MS results demonstrated that Lys644 and Gln1351 of α1 chain and Gln19 of α2 chain were the main sites. The results of molecular dynamic (MD) simulation indicated that at the moderate degree of covalent crosslinking (13.55%), Radius of gyration (Rg) values decreased from 11.5 to 10.0 nm at 4 °C and led to a more compact structure, therefore enhancing structural tightness and gel strength. As the degree of crosslinking rose, the number of crosslinking sites gradually increased, with a predominant distribution in the C-terminal pro-peptide of gelatin chain. At the high crosslinking degree (37.48%), Rg value decreased at 100 °C and gelatin chains were connected through a “head to tail” linkage configuration, resulting in a significant reduction in the triple helix-like structure. The structure maintained the gel-forming ability at high temperature while lowering the aggregation energy of water molecules from −2.359 × 105 to −2.455 × 105 kJ/mol, leading to thermal irreversibility of gelatin hydrogel.

Benchmarking docking, density functional theory and molecular dynamics studies to assess the aldose reductase inhibitory potential of Trigonella foenum-graecum compounds for managing diabetes-associated complications

Inhibition of aldose reductase (AR) could be a beneficial strategy for managing diabetes-associated complications. Trigonella foenum-graecum (TFG) is used around the globe as a traditional medicine for the management of diabetes. Our study aimed to assess the potential of TFG phytocompounds as inhibitors of AR in the context of diabetes-related complications. Our research work employed molecular docking, density functional theory (DFT) and molecular dynamics (MD) to evaluate the efficacy of TFG compounds. The study compared the predictive power of AutoDock and AutoDock Vina docking software and found that AutoDock Vina performs better in ranking and discriminating actives and decoys. The research identified five compounds as potential AR inhibitors from fifty-eight reported TFG phytoconstituents. Tigogenin and Gitogenin stood out as the most promising AR inhibitors. The electronic properties of the compounds were analysed through DFT studies and provided insights into their binding potential. Finally, the results of MD simulations indicated that Tigogenin and Gitogenin bound robustly with AR throughout the simulation period. This study predicted the AR inhibitory potential of TFG compounds for managing diabetes-associated complications and supports further drug development from TFG. The benchmarking approach used in this study improves the accuracy and dependability of bioactivity prediction.

Molecular Dynamics Simulations of Different Nanoparticles at Substrates.

We report the results of large-scale molecular dynamics simulations of adsorption nanoparticles on solid surfaces. The particles were modeled as stiff aggregates of spherical segments. Three types of particles were studied: rods, rectangles, and triangles built of the same number of segments. We show how the particle shape affects the adsorption, the structure of the surface layer, and the degree of the removal of particles from the solvent. The systems with different segment-segment and segment-surface interactions and different concentrations of particles were investigated. The ordered structures formed in adsorption monolayers were also analyzed. The results are consistent with experimental observations.

Investigation of cathinone analogs targeting human dopamine transporter using molecular modeling

In a step towards understanding the structure–property relationship among Synthetic Cathinones (SCs), a combined methodology based on Density Functional Theory (DFT), Administration, Distribution, Metabolism, Excretion, and Toxicity (ADMET) predictions, docking and molecular dynamics simulations have been applied to correlate physicochemical descriptors of various SCs to their biological activity. The results from DFT and molecular docking studies correlate well with each other explaining the biological activity trends of the studied SCs. Quantum mechanical descriptors viz. polarizability, electron affinity, ionization potential, chemical hardness, electronegativity, molecular electrostatic potential, and ion interaction studies unravel the distinguishingly reactive nature of Group D (pyrrolidine substituted) and Group E (methylenedioxy and pyrrolidine substituted) compounds. According to ADMET analysis, Group D and Group E molecules have a higher probability of permeating through the blood–brain barrier. Molecular docking results indicate that Phe76, Ala77, Asp79, Val152, Tyr156, Phe320, and Phe326 constitute the binding pocket residues of hDAT in which the most active ligands MDPV, MDPBP, and MDPPP are bound. Finally, to validate the derived quantum chemical descriptors and docking results, Molecular Dynamics (MD) simulations are performed with homology-modelled hDAT (human dopamine transporter). The MD simulation results revealed that the majority of SCs remain stable within the hDAT protein’s active sites via non-bonded interactions after 100 ns long simulations. The findings from DFT, ADMET analysis, molecular docking, and molecular dynamics simulation studies complement each other suggesting that pyrrolidine-substituted SCs (Group D and E), specifically, MPBP and PVN are proven potent SCs along with MDPV, validating various experimental observations. Communicated by Ramaswamy H. Sarma

The synergistic effect of body modification and surface coating for the moisture barrier of hydrophilic PU film

Polyurethane (PU) films have established a distinctive position among various protective materials due to their outstanding chemical and physical properties. However, their limited moisture barrier performance hinders their broader applications. To solve this difficult problem, PU film is first modified by adding reactive fillers (MgO) and physical barrier fillers (Modified graphene oxide, M-GO), and then Si-DLC coating is deposited on the surface of the film by hollow cathode plasma immersion ion implantation (HCPIII) method at near room temperature. MgO can react with the water molecules diffused into the film, while the addition of M-GO extends the diffusion path of water molecules. It is revealed that the addition of 5% MgO and 0.5% M-GO, along with the deposition of a 2.56 μm Si-doped DLC coating, leads to a remarkable 14.67-fold reduction in the water vapor transmittance rate (WVTR). The result of molecular dynamics simulations shows that the reactive sheet filler captures the water molecules and reduces the WVTR. In addition, the combination method of PU body modification and surface coating is suitable for polyurethane foam (PUF) to reduce the WVTR as a double coating. It's amazing that the WVTR of PUF significantly reduces by 15.29 times after double-coating. The combination of body modification and surface coating provides an effective route for significantly improving the moisture barrier property of hydrophilic polymer film and polymer foam.

Transport Dynamics of Water Molecules Confined between Lipid Membranes.

Water molecules confined between biological membranes exhibit a distinctive non-Gaussian displacement distribution, far different from that of bulk water. Here, we introduce a new transport equation for water molecules in the intermembrane space, quantitatively explaining molecular dynamics simulation results. We find that the unique transport dynamics of water molecules stems from the lateral diffusion coefficient fluctuation caused by their longitudinal motion in the direction perpendicular to the membranes. We also identify an interfacial region where water possesses distinct physical properties, which is unaffected by changes in the intermembrane separation.

Effect of temperature on the mechanical properties of Ni-based superalloys via molecular dynamics and crystal plasticity

Study on the influence of temperature on the mechanical properties across multiple scales has been a focus on the research of Hastelloy-X (HX) alloys for the application in high-temperature structure components. In this work, Molecular Dynamics (MD) and Crystal Plasticity (CP) are put together to solve it from atomic scale to mesoscopic scale. MD research indicates that the deformation of HX alloy occurs in two stages at temperature below 300 K: initially, as stacking fault deforms, stacking fault can transform into twinning with increasing strain. When the temperature exceeds 300 K, deformation primarily forms a stacking fault. The twinning deformation path transforms from intrinsic stacking fault to extrinsic stacking fault and then to twinning. A mesoscopic-scale CP model was developed using atomic-scale deformation mechanisms to bridge the gap between deformation mechanisms and experimental results. The CP results indicate a functional relationship between the strength of HX alloy and temperature. This relationship appears insensitive to crystal texture and grain shape. Incorporating grain morphology and texture into the model significantly impacts the strength response of calculating HX alloy. After the tensile deformation of HX alloy at 300 and 1173 K, the atomic scale deformation results characterized by transmission electron microscopy are aligned with the MD simulation results. The relationship between strength and temperature predicted by CP results has also been validated. A thorough investigation into the deformation behavior of HX alloys across different scales, employing MD and CP models, introduces a novel approach for predicting the mechanical properties of superalloys.

Discovery of a novel SHP2 allosteric inhibitor using virtual screening, FMO calculation, and molecular dynamic simulation.

SHP2 is a non-receptor protein tyrosine phosphatase to remove tyrosine phosphorylation. Functionally, SHP2 is an essential bridge to connect numerous oncogenic cell-signaling cascades including RAS-ERK, PI3K-AKT, JAK-STAT, and PD-1/PD-L1 pathways. This study aims to discover novel and potent SHP2 inhibitors using a hierarchical structure-based virtual screening strategy that combines molecular docking and the fragment molecular orbital method (FMO) for calculating binding affinity (referred to as the Dock-FMO protocol). For the SHP2 target, the FMO method prediction has a high correlation between the binding affinity of the protein-ligand interaction and experimental values (R2 = 0.55), demonstrating a significant advantage over the MM/PBSA (R2 = 0.02) and MM/GBSA (R2 = 0.15) methods. Therefore, we employed Dock-FMO virtual screening of ChemDiv database of ∼2,990,000 compounds to identify a novel SHP2 allosteric inhibitor bearing hydroxyimino acetamide scaffold. Experimental validation demonstrated that the new compound (E)-2-(hydroxyimino)-2-phenyl-N-(piperidin-4-ylmethyl)acetamide (7188-0011) effectively inhibited SHP2 in a dose-dependent manner. Molecular dynamics (MD) simulation analysis revealed the binding stability of compound 7188-0011 and the SHP2 protein, along with the key interacting residues in the allosteric binding site. Overall, our work has identified a novel and promising allosteric inhibitor that targets SHP2, providing a new starting point for further optimization to develop more potent inhibitors. All the molecular docking studies were employed to identify potential leads with Maestro v10.1. The protein-ligand binding affinities of potential leads were further predicted by FMO calculations at MP2/6-31G* level using GAMESS v2020 system. MD simulations were carried out with AmberTools18 by applying the FF14SB force field. MD trajectories were analyzed using VMD v1.9.3. MM/GB(PB)SA binding free energy analysis was carried out with the mmpbsa.py tool of AmberTools18. The docking and MD simulation results were visualized through PyMOL v2.5.0.

A new insight on a novel auxiliary collector 4-MBA synergize with BHA to enhance flotation of scheelite

The hydroxamic acid reagents are extremely pervasive in the collector of scheelite flotation. The high flotation cost is caused by its large dosage and high price. A novel auxiliary collector 4-methodybenzoicacid (4-MBA) was used in combination with benzohydroxamic acid (BHA) to optimize the flotation reagents system for scheelite in the study. A single 4-MBA has weak collection capacity for scheelite. Nevertheless, the combination of 4-MBA and BHA optimized the flotation system of scheelite by ensuring its recovery and reducing the reagent consumption. The results of solution surface tension and contact angle measurements showed that the surface tension of the mixed reagent solution is a relatively lower and the hydrophobicity of the scheelite surface is a relatively higher in the BHA and 4-MBA mixed reagent system. Additionally, a synergistic effect between BHA and 4-MBA on Pb2+ activated scheelite is revealed by the Zeta potential measurements, infrared spectroscopy (FT-IR) and X-ray photoelectron spectroscopy (XPS). The addition of 4-MBA promoted the physical and chemical adsorption of BHA. The results of molecular dynamics simulation (MDS) further explained that the collector adsorption layer of the surface of scheelite (112) is more compact under BHA and 4-MBA mixed reagent system, which prevents the vertical entry of water phase and improves the hydrophobicity of scheelite. Finally, the bench-scale flotation experiment verified that the addition of 4-MBA enhanced scheelite flotation with BHA, and decrease the flotation cost of scheelite.