Nickel-catalyzed asymmetric hydrogenation for the preparation of α-substituted propionic acids.

Transition metal-catalyzed asymmetric hydrogenation is one of the most efficient methods for the preparation of chiral α-substituted propionic acids. However, research on this method, employing cleaner earth-abundant metal catalysts, is still insufficient in both academic and industrial contexts. Herein, we report an efficient nickel-catalyzed asymmetric hydrogenation of α-substituted acrylic acids affording the corresponding chiral α-substituted propionic acids with up to 99.4% ee (enantiomeric excess) and 10,000 S/C (substrate/catalyst). In particular, this method can be used to obtain (R)-dihydroartemisinic acid with 99.8:0.2 dr (diastereomeric ratio) and 5000 S/C, which is an essential intermediate for the preparation of the antimalarial drug Artemisinin. The reaction mechanism has been investigated via experiments and DFT (Density Functional Theory) calculations, which indicate that the protonolysis of the C-Ni bond of the key intermediate via an intramolecular proton transfer from the carboxylic acid group of the substrate, is the rate-determining step.

A firm-push-to-open and light-push-to-lock strategy for a general chemical platform to develop activatable dual-modality NIR-II probes.

Activatable near-infrared (NIR) imaging in the NIR-II range is crucial for deep tissue bioanalyte tracking. However, designing such probes remains challenging due to the limited availability of general chemical strategies. Here, we introduced a foundational platform for activatable probes, using analyte-triggered smart modulation of the π-conjugation system of a NIR-II-emitting rhodamine hybrid. By tuning the nucleophilicity of the ortho-carboxy moiety, we achieved an electronic effect termed "firm-push-to-open and light-push-to-lock," which enables complete spirocyclization of the probe before sensing and allows for efficient zwitterion formation when the light-pushing aniline carbamate trigger is transformed into a firm-pushing aniline. This platform produces dual-modality NIR-II imaging probes with ~50-fold fluorogenic and activatable photoacoustic signals in live mice, surpassing reported probes with generally below 10-fold activatable signals. Demonstrating generality, we successfully designed probes for hydrogen peroxide (H2O2) and hydrogen sulfide (H2S). We envision a widespread adoption of the chemical platform for designing activatable NIR-II probes across diverse applications.

Unveiling Conformational States of CDK6 Caused by Binding of Vcyclin Protein and Inhibitor by Combining Gaussian Accelerated Molecular Dynamics and Deep Learning.

CDK6 plays a key role in the regulation of the cell cycle and is considered a crucial target for cancer therapy. In this work, conformational transitions of CDK6 were identified by using Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy landscapes (FELs). DL finds that the binding pocket as well as the T-loop binding to the Vcyclin protein are involved in obvious differences of conformation contacts. This result suggests that the binding pocket of inhibitors (LQQ and AP9) and the binding interface of CDK6 to the Vcyclin protein play a key role in the function of CDK6. The analyses of FELs reveal that the binding pocket and the T-loop of CDK6 have disordered states. The results from principal component analysis (PCA) indicate that the binding of the Vcyclin protein affects the fluctuation behavior of the T-loop in CDK6. Our QM/MM-GBSA calculations suggest that the binding ability of LQQ to CDK6 is stronger than AP9 with or without the binding of the Vcyclin protein. Interaction networks of inhibitors with CDK6 were analyzed and the results reveal that LQQ contributes more hydrogen binding interactions (HBIs) and hot interaction spots with CDK6. In addition, the binding pocket endures flexibility changes from opening to closing states and the Vcyclin protein plays an important role in the stabilizing conformation of the T-loop. We anticipate that this work could provide useful information for further understanding the function of CDK6 and developing new promising inhibitors targeting CDK6.

Molecular Mechanism of Phosphorylation-Mediated Impacts on the Conformation Dynamics of GTP-Bound KRAS Probed by GaMD Trajectory-Based Deep Learning.

The phosphorylation of different sites produces a significant effect on the conformational dynamics of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations were combined with deep learning (DL) to explore the molecular mechanism of the phosphorylation-mediated effect on conformational dynamics of the GTP-bound KRAS. The DL finds that the switch domains are involved in obvious differences in conformation contacts and suggests that the switch domains play a key role in the function of KRAS. The analyses of free energy landscapes (FELs) reveal that the phosphorylation of pY32, pY64, and pY137 leads to more disordered states of the switch domains than the wild-type (WT) KRAS and induces conformational transformations between the closed and open states. The results from principal component analysis (PCA) indicate that principal motions PC1 and PC2 are responsible for the closed and open states of the phosphorylated KRAS. Interaction networks were analyzed and the results verify that the phosphorylation alters interactions of GTP and magnesium ion Mg2+ with the switch domains. It is concluded that the phosphorylation pY32, pY64, and pY137 tune the activity of KRAS through changing conformational dynamics and interactions of the switch domains. We anticipated that this work could provide theoretical aids for deeply understanding the function of KRAS.

Binding Mechanism of Inhibitors to BRD4 and BRD9 Decoded by Multiple Independent Molecular Dynamics Simulations and Deep Learning.

Bromodomain 4 and 9 (BRD4 and BRD9) have been regarded as important targets of drug designs in regard to the treatment of multiple diseases. In our current study, molecular dynamics (MD) simulations, deep learning (DL) and binding free energy calculations are integrated to probe the binding modes of three inhibitors (H1B, JQ1 and TVU) to BRD4 and BRD9. The MD trajectory-based DL successfully identify significant functional function domains, such as BC-loop and ZA-loop. The information from the post-processing analysis of MD simulations indicates that inhibitor binding highly influences the structural flexibility and dynamic behavior of BRD4 and BRD9. The results of the MM-GBSA calculations not only suggest that the binding ability of H1B, JQ1 and TVU to BRD9 are stronger than to BRD4, but they also verify that van der Walls interactions are the primary forces responsible for inhibitor binding. The hot spots of BRD4 and BRD9 revealed by residue-based free energy estimation provide target sites of drug design in regard to BRD4 and BRD9. This work is anticipated to provide useful theoretical aids for the development of selective inhibitors over BRD family members.

3D-QSAR pharmacophore modeling, virtual screening, molecular docking, MD simulations, in vitro and in vivo studies to identify potential anti-hyperplasia drugs.

Psoriasis is a common immune-mediated skin condition characterized by aberrant keratinocytes and cell proliferation. The purpose of this study was to explore the FDA-approved drugs by 3D-QSAR pharmacophore model and evaluate their efficiency by in-silico, in vitro, and in vivo psoriasis animal model. A 3D-QSAR pharmacophore model was developed by utilizing HypoGen algorithm using the structural features of 48 diaryl derivatives with diverse molecular patterns. The model was validated by a test set of 27 compounds, by cost analysis method, and Fischer's randomization test. The correlation coefficient of the best model (Hypo2) was 0.9601 for the training set while it was 0.805 for the test set. The selected model was taken as a 3D query for the virtual screening of over 3000 FDA-approved drugs. Compounds mapped with the pharmacophore model were further screened through molecular docking. The hits that showed the best docking results were screened through in silico skin toxicity approach. Top five hits were selected for the MD simulation studies. Based on MD simulations results, the best two hit molecules, that is, ebastine (Ebs) and mebeverine (Mbv) were selected for in vitro and in vivo antioxidant studies performed in mice. TNF-α and COX pro-inflammatory mediators, biochemical assays, histopathological analyses, and immunohistochemistry observations confirmed the anti-inflammatory response of the selected drugs. Based on these findings, it appeared that Ebs can effectively treat psoriasis-like skin lesions and down-regulate inflammatory responses which was consistent with docking predictions and could potentially be employed for further research on inflammation-related skin illnesses such as psoriasis.

MCPNET: Development of an interpretable deep learning model based on multiple conformations of the compound for predicting developmental toxicity.

The development of deep learning models for predicting toxicological endpoints has shown great promise, but one of the challenges in the field is the accuracy and interpretability of these models. The bioactive conformation of a compound plays a critical role for it to bind in the target. It is a big issue to figure out the bioactive conformation in deep learning without the co-crystal structure or highly precise molecular simulations. In this study, we developed a deep learning framework of Multi-Conformation Point Network (MCPNET) to construct classification and regression models, respectively, based on electrostatic potential distributions on vdW surfaces around multiple conformations of the compound using a dataset of compounds with developmental toxicity in zebrafish embryo. MCPNET applied 3D multi-conformational surface point cloud to extract the molecular features for model training, which may be critical for capturing the structural diversity of compounds. The models achieved an accuracy of 85 % on the classification task and R2 of 0.66 on the regression task, outperforming traditional machine learning models and other deep learning models. The key feature of our model is its interpretability with the component visualization to identify the factors contributing to the prediction and to understand the compound action mechanism. MCPNET may predict the conformation quietly close to the bioactive conformation of a compound by attention-based multi-conformation pooling mechanism. Our results demonstrated the potential of deep learning based on 3D molecular representations in accurately predicting developmental toxicity. The source code is publicly available at https://github.com/Superlit-CC/MCPNET.

Skin-Conformable Flexible and Stretchable Ultrasound Transducer for Wearable Imaging.

Ultrasound imaging offers a noninvasive, radiation-free method for visualizing internal tissues and organs, with deep penetration capabilities. This has established it as a crucial tool for physicians in diagnosing internal tissue pathologies and monitoring human conditions. Nonetheless, conventional ultrasound probes are often characterized by their rigidity and bulkiness. Designing a transducer that can seamlessly adapt to the contours and dynamics of soft, curved human skin presents significant technical hurdles. We present a novel flexible and stretchable ultrasound transducer (FSUT) designed for adaptability to large-curvature surfaces while preserving superior imaging quality. Central to this breakthrough is the innovative use of screen-printed silver nanowires (AgNWs) coupled with a composite elastic substrate, together ensuring robust and stable electrical and mechanical connections. Standard tensile and fatigue tests verify its durability. The mechanical, electrical, and acoustic properties of FSUTs are characterized using standard methods, with large tensile strains (≥110%), high flexibility ( R ≥ 1.4 mm), and lightweight ( ≤ 1.58 g) to meet the needs of wearable devices. Center frequency and -6-dB bandwidth are approximately 5.3 MHz and 66.47%, respectively. Images of the commercial anechoic cyst phantom yielded an axial and lateral resolution (depths of 10-70 mm) of approximately 0.31 and 0.46, and 0.34 and 0.84 mm, respectively. The complex curved surface imaging capabilities of FSUT were tested on agar-gelatin-based breast cyst phantoms under different curvatures. Finally, ultrasound images of the thyroid, brachial, and carotid arteries were also obtained from volunteer wearing FSUT.

Shotgun lipidomics of human subretinal fluids under rod-dominant retina reveals cone-dominated lipids.

Subretinal fluid (SRF) accumulates between photoreceptor outer segments and retinal pigment epithelium during rhegmatogenous retinal detachment. Biomolecular components such as lipids originate from cells surrounding the SRF. Knowledge of the composition of these molecules in SRF potentially provides mechanistic insight into the physiologic transfer of lipids between retinal tissue compartments. Using mass spectrometry and tandem mass spectrometry analysis on an electrospray ionization quadrupole-time-of-flight mass spectrometer, we identified a total of 115 lipid molecular species of 11 subclasses and 9 classes in two samples from two patients with rhegmatogenous retinal detachment. These included 47 glycerophosphocholines, 6 glycerophosphoethanolamines, 1 glycerophosphoinositol, 18 sphingomyelins, 9 cholesteryl esters, free cholesterol, 3 ceramides, 22 triacylglycerols and 8 free fatty acids. Glycerophosphocholines were of the highest intensity. By minimizing the formation of different adduct forms or clustering ions of different adducts, we determined the relative intensity of lipid molecular species within the same subclasses. The profiles were compared with those of retinal cells available in the published literature. The glycerophosphocholine profile of SRF was similar to that of cone outer segments, suggesting that outer segment degradation products are constitutively released into the interphotoreceptor matrix, appearing in SRF during detachment. This hypothesis was supported by the retinal distributions of corresponding lipid synthases' mRNA expression obtained from an online resource based on publicly available single-cell sequencing data. In contrast, based on lipid profiles and relevant gene expression in this study, the sources of free cholesterol and cholesteryl esters in SRF appeared more ambiguous, possibly reflecting that outer retina takes up plasma lipoproteins. Further studies to identify and quantify lipids in SRF will help better understand etiology of diseases relevant to outer retina.

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters.

Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.

The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes.

Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.

Insights into the selectivity of a brain-penetrant CDK4/6 vs CDK1/2 inhibitor for glioblastoma used in multiple replica molecular dynamics simulations.

Cyclin dependent kinases (CDKs) play an important role in cell cycle regulation and their dysfunction is associated with many cancers. That is why CDKs have been attractive targets for the treatment of cancer. Glioblastoma is a cancer caused by the aberrant expression of CDK4/6, so exploring the mechanism of the selection of CDK4/6 toward inhibitors relative to the other family members CDK1/2 is essential. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), molecular mechanics Poisson-Boltzmann/Generalized Born surface area (MM-PB/GBSA) and other methods were integrated to decipher the selectively binding mechanism of the inhibitor N1J to CDK4/6 and CDK1/2. Molecular electrostatic potential (MESP) analysis provides an explanation for the N1J selectivity. Residue-based free energy decomposition reveals that most of the hot residues are located at the same location of CDKs proteins, but the different types of residues in different proteins cause changes in binding energy, which is considered as a potential developmental direction to improve the selectivity of inhibitors to CDK4/6. These results provide insights into the source of inhibitor and CDK4/6 selectivity for the future development of more selective inhibitors.Communicated by Ramaswamy H. Sarma.

Multilayer meta-matching: translating phenotypic prediction models from multiple datasets to small data.

Resting-state functional connectivity (RSFC) is widely used to predict phenotypic traits in individuals. Large sample sizes can significantly improve prediction accuracies. However, for studies of certain clinical populations or focused neuroscience inquiries, small-scale datasets often remain a necessity. We have previously proposed a "meta-matching" approach to translate prediction models from large datasets to predict new phenotypes in small datasets. We demonstrated large improvement of meta-matching over classical kernel ridge regression (KRR) when translating models from a single source dataset (UK Biobank) to the Human Connectome Project Young Adults (HCP-YA) dataset. In the current study, we propose two meta-matching variants ("meta-matching with dataset stacking" and "multilayer meta-matching") to translate models from multiple source datasets across disparate sample sizes to predict new phenotypes in small target datasets. We evaluate both approaches by translating models trained from five source datasets (with sample sizes ranging from 862 participants to 36,834 participants) to predict phenotypes in the HCP-YA and HCP-Aging datasets. We find that multilayer meta-matching modestly outperforms meta-matching with dataset stacking. Both meta-matching variants perform better than the original "meta-matching with stacking" approach trained only on the UK Biobank. All meta-matching variants outperform classical KRR and transfer learning by a large margin. In fact, KRR is better than classical transfer learning when less than 50 participants are available for finetuning, suggesting the difficulty of classical transfer learning in the very small sample regime. The multilayer meta-matching model is publicly available at GITHUB_LINK.

Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study.

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continues to spread globally, and rapid viral evolution and the emergence of new variants pose challenges to pandemic control. During infection, the spike protein of SARS-CoV-2 interacts with the human ACE2 protein via its receptor binding domain (RBD), and it is known that engineered forms of ACE2 can compete with wild-type (WT) ACE2 for binding to inhibit infection. Here, we conducted multiple replica molecular dynamics (MRMD) simulations to study the mechanisms of the engineered ACE2 variants 3N39 and 3N94 and provide directions for optimization. Our findings reveal that engineered ACE2 is notably more efficacious in systems that show weaker binding to WT ACE2 (i.e., WT and BA.1 RBD), but also faces immune escape as the virus evolves. Moreover, by modifying residue types near the binding interface, engineered ACE2 alters the electrostatic potential distribution and reconfigures the hydrogen bonding network, which results in modified binding to the RBD. However, this structural rearrangement does not occur in all RBD variants. In addition, we identified potentially engineerable beneficial residues and potentially engineerable detrimental residues in both ACE2 and RBD. Functional conservation can thus enable the optimization of these residues and improve the binding competitiveness of engineered ACE2, which therefore provides additional immune escape prevention. Finally, we conclude that these findings have implications for understanding the mechanisms responsible for engineered ACE2 and can help us to develop engineered ACE2 proteins that show superior performance.

Effect of Glass Fiber-Reinforced Plastic Waste on the Mechanical Properties of Concrete and Evaluation of Its Feasibility for Reuse in Concrete Applications.

The disposal of glass fiber-reinforced plastic (GFRP) waste has become an urgent issue in both the engineering and environmental fields. In this study, the feasibility of reusing mechanically recycled GFRP in concrete was evaluated. Secondary screening of the recycled material was conducted to obtain different types of products, and the recycled GFRP (rGFRP) was characterized. Subsequently, the effect of rGFRP on concrete performance was evaluated using different dosages (0%, 1%, 3%, 5%) of rGFRP powder and rGFRP cluster (with different sizes and fiber contents) to replace fine aggregate in concrete preparation. The experimental results indicated that the addition of rGFRP powder has no significant impact on the mechanical properties of concrete, while the addition of a small amount of rGFRP cluster slightly improves the compressive strength and splitting tensile strength of concrete. Additionally, the short fibers in rGFRP improve the failure mode of concrete, and increased fiber content and longer fiber length demonstrate a more pronounced reinforcing effect. The challenges and potential directions for future research in the realm of reusing rGFRP in concrete are discussed at the end. A systematic process for reusing GFRP waste in concrete is proposed to address the primary challenges and provide guidance for its practical engineering application.

Exploring conformation changes of Janus kinase 2 pseudokinase mediated by mutations through Gaussian accelerated molecular dynamics and principal component analysis.

The pseudokinase domain (JH2) of the protein tyrosine kinase (Janus kinase 2, JAK2) regulates the activity of a tyrosine kinase domain (JH1) in JAK2, which is further affected by mutations in the JH2. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by construction of free energy landscapes (FELs) and principal component analysis (PCA) were performed to study effect of two mutations V617F and V617F/E596A on the conformations of the ATP-bound JH2. The dynamic analyses reveal that mutations affect the structural flexibility and correlated motions of the JH2, meanwhile also change the dynamics behavior of the P-loop and αC-helix of the JH2. The information from FELs unveils that mutations induce less energy states than the free JH2 and the WT one. The analyses of interaction networks uncover that mutations affect the salt bridge interactions of ATP with K581, K677 and R715 and alter hydrogen bonding interactions (HBIs) of ATP with the JH2. The changes in conformations of the JH2 and ATP-JH2 interaction networks caused by mutations in turn generate effect on the activity regulations of the JH2 on the JH1. This work is expected to provide significant theoretical helps for deeply understanding the function of the JH2 and drug design toward JAK2.Communicated by Ramaswamy H. Sarma.

The impact of different service states of tunnel lighting on traffic safety.

The service states of tunnel lighting will directly affect the lighting conditions, which affect traffic safety. Therefore, it is imperative to evaluate and predict traffic safety accurately in different lighting states. In this research, three hundred experimental scenarios of the service states of tunnel lighting were designed and implemented to evaluate the impact of different service states of tunnel lighting on traffic safety. The evaluation was achieved through a visual identification experiment in a physical tunnel. The experimental results show higher simulated vehicle speeds pose a greater threat to traffic safety. The severity of lighting attenuation contributes to an increased risk to traffic safety. An increase in the number of luminaires failure also poses a greater threat to traffic safety. The newly proposed traffic safety factor was employed to evaluate traffic safety quantitatively in road tunnels. To improve the accuracy and comprehensiveness of the traffic safety factor prediction in different lighting service states, an advanced neural network prediction system was developed. The prediction system was constructed using the Sparrow Search Algorithm (SSA) to optimize Extreme Learning Machine (ELM) neural network, and the dataset from the experiment was used for the prediction model. The SSA-ELM neural network model is a reliable model that can predict the traffic safety factor comprehensively and accurately. The recommended threshold value for the traffic safety factor is 0.6. When the value falls below 0.6, it shows that the service states of tunnel lighting pose a threat to traffic safety in the tunnel. These findings can provide insights into the safe and energy-efficient maintenance of road tunnels.

Synthesis of 2,3-Fused Quinazolinones via the Radical Cascade Pathway and Reaction Mechanistic Studies by DFT Calculations.

An efficient radical cascade cyclization of unactivated alkenes toward the synthesis of a series of ring-fused quinazolinones has been developed in moderate to excellent yields using commercially available ethers, alkanes, and alcohols, respectively, under a base-free condition in a short time without a transition metal as catalyst. Notably, the transformations can be carried out with the advantages of a broad substrate scope and high atomic economy. Density functional theory calculations and wavefunction analyses were performed to elucidate the radical reaction mechanism.

Copper-Catalyzed Chemoselective Asymmetric Hydrogenation of C=O Bonds of Exocyclic α,ß-Unsaturated Pentanones.

A highly chemoselective earth-abundant transition metal copper catalyzed asymmetric hydrogenation of C=O bonds of exocyclic α,ß-unsaturated pentanones was realized using H2 . The desired products were obtained with up to 99 % yield and 96 % ee (enantiomeric excess) (99 % ee, after recrystallization). The corresponding chiral exocyclic allylic pentanol products can be converted into several bioactive molecules. The hydrogenation mechanism was investigated via deuterium-labelling experiments and control experiments, which indicate that the keto-enol isomerization rate of the substrate is faster than that of the hydrogenation and also show that the Cu-H complex can only catalyze chemoselectively the asymmetric reduction of the carbonyl group. Computational results indicate that the multiple attractive dispersion interactions (MADI effect) between the catalyst with bulky substituents and substrate play important roles which stabilize the transition states and reduce the generation of by-products.

Identification Mechanism of BACE1 on Inhibitors Probed by Using Multiple Separate Molecular Dynamics Simulations and Comparative Calculations of Binding Free Energies.

ß-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward the treatment of AD.