Assessing the role of residue Phe108 of cytochrome P450 3A4 in allosteric effects of midazolam metabolism.

Cytochrome P450 3A4 (CYP3A4) is involved in the metabolism of more drugs in clinical use than any other xenobiotic-metabolizing enzyme. CYP3A4-mediated drug metabolism is usually allosterically modulated by substrate concentration (homotropic allostery) and other drugs (heterotropic allostery), exhibiting unusual kinetic profiles and regiospecific metabolism. Recent studies suggest that residue Phe108 (F108) of CYP3A4 may have an important role in drug metabolism. In this work, residue mutations were coupled with well-tempered metadynamics simulations to assess the importance of F108 in the allosteric effects of midazolam metabolism. Comparing the simulation results of the wild-type and mutation systems, we identify that the π-π interaction and steric effect between the F108 side chain and midazolam is favorable for the stable binding of substrate in the active site. F108 also plays an important role in the transition of substrate binding mode, which mainly induces the transition of substrate binding mode by forming π-π interactions with multiple aromatic rings of the substrate. Moreover, the side chain of F108 is closely related to the radius and depth of the 2a and 2f channels, and F108 may further regulate drug metabolism by affecting the pathway, orientation, or time of substrate entry into the CYP3A4 active site or product egress from the active site. Altogether, we suggest that F108 affects drug metabolism and regulatory mechanisms by affecting substrate binding stability, binding mode transition, and channel characteristics of CYP3A4. Our findings could promote the understanding of complicated allosteric mechanisms in CYP3A4-mediated drug metabolism, and the knowledge could be used for drug development and disease treatment.

VARIDT 2.0: structural variability of drug transporter.

The structural variability data of drug transporter (DT) are key for research on precision medicine and rational drug use. However, these valuable data are not sufficiently covered by the available databases. In this study, a major update of VARIDT (a database previously constructed to provide DTs' variability data) was thus described. First, the experimentally resolved structures of all DTs reported in the original VARIDT were discovered from PubMed and Protein Data Bank. Second, the structural variability data of each DT were collected by literature review, which included: (a) mutation-induced spatial variations in folded state, (b) difference among DT structures of human and model organisms, (c) outward/inward-facing DT conformations and (d) xenobiotics-driven alterations in the 3D complexes. Third, for those DTs without experimentally resolved structural variabilities, homology modeling was further applied as well-established protocol to enrich such valuable data. As a result, 145 mutation-induced spatial variations of 42 DTs, 1622 inter-species structures originating from 292 DTs, 118 outward/inward-facing conformations belonging to 59 DTs, and 822 xenobiotics-regulated structures in complex with 57 DTs were updated to VARIDT (https://idrblab.org/varidt/ and http://varidt.idrblab.net/). All in all, the newly collected structural variabilities will be indispensable for explaining drug sensitivity/selectivity, bridging preclinical research with clinical trial, revealing the mechanism underlying drug-drug interaction, and so on.

Novel Anthraquinone-Based Benzenesulfonamide Derivatives and Their Analogues as Potent Human Carbonic Anhydrase Inhibitors with Antitumor Activity: Synthesis, Biological Evaluation, and In Silico Analysis.

In this study, we designed two series of novel anthraquinone-based benzenesulfonamide derivatives and their analogues as potential carbonic anhydrase inhibitors (CAIs) and evaluated their inhibitory activities against off-target human carbonic anhydrase II (hCA II) isoform and tumor-associated human carbonic anhydrase IX (hCA IX) isoform. Most of these compounds exhibited good inhibitory activities against hCA II and IX. The compounds that exhibited the best hCA inhibition were further studied against the MDA-MB-231, MCF-7, and HepG2 cell lines under hypoxic and normoxic conditions. Additionally, the compounds exhibiting the best antitumor activity were subjected to apoptosis and mitochondrial membrane potential assays, which revealed a significant increase in the percentage of apoptotic cells and a notable decrease in cell viability. Molecular docking studies were performed to demonstrate the presence of numerous hydrogen bonds and hydrophobic interactions between the compounds and the active site of hCA. Absorption, distribution, metabolism, excretion (ADME) predictions showed that all of the compounds had good pharmacokinetic and physicochemical properties.

The Potential of Cyclodextrins as Inhibitors for the BM2 Protein: An In Silico Investigation.

The influenza BM2 transmembrane domain (BM2TM), an acid-activated proton channel, is an attractive antiviral target due to its essential roles during influenza virus replication, whereas no effective inhibitors have been reported for BM2. In this study, we draw inspiration from the properties of cyclodextrins (CDs) and hypothesize that CDs of appropriate sizes may possess the potential to act as inhibitors of the BM2TM proton channel. To explore this possibility, molecular dynamics simulations were employed to assess their inhibitory capabilities. Our findings reveal that CD4, CD5, and CD6 are capable of binding to the BM2TM proton channel, resulting in disrupted water networks and reduced hydrogen bond occupancy between H19 and the solvent within the BM2TM channel necessary for proton conduction. Notably, CD4 completely obstructs the BM2TM water channel. Based on these observations, we propose that CD4, CD5, and CD6 individually contribute to diminishing the proton transfer efficiency of the BM2 protein, and CD4 demonstrates promising potential as an inhibitor for the BM2 proton channel.

Insights into the binding mechanism between α-TOH and CYP4F2: A homology modeling, molecular docking, and molecular dynamics simulation study.

α-Tocopherol (α-TOH) is a potent antioxidant. The concentrations of α-TOH in plasma are closely related to human health. α-TOH can be regulated by the metabolism of cytochrome P450 4F2 (CYP4F2). However, the atomic-level basis for this regulation process remains elusive. Here, we successfully constructed the structure of CYP4F2 by homology modeling and obtained the α-TOH-CYP4F2 complex models using molecular docking. Three parallel 500 ns molecular dynamics simulations were performed on each complex model to investigate the details of the interaction between α-TOH and CYP4F2. MM-GBSA method combined with principal component analysis shows that 8 key residues establish a hydrophobic cavity stabilizing α-TOH in the pocket of CYP4F2 and S423 forms an important hydrogen bond with α-TOH anchoring α-TOH in the favorable position for ω-hydroxylation. Based on our simulation results and the experimental facts, we designed mutation simulation experiments to clarify the important role of two key residues (S423 and V433) in the binding of α-TOH with CYP4F2. The results show that the mutations directly or indirectly change the binding mode of α-TOH and decrease its binding affinity with CYP4F2, which is unfavorable for ω-hydroxylation. Our results could enrich the information on structure-function relationships of CYP4F2 and provide valuable insights into the regulatory mechanism of CYP4F2 on the metabolism of α-TOH.

The Important Role of Transporter Structures in Drug Disposition, Efficacy, and Toxicity.

The ATP-binding cassette (ABC) and solute carrier (SLC) transporters are critical determinants of drug disposition, clinical efficacy, and toxicity as they specifically mediate the influx and efflux of various substrates and drugs. ABC transporters can modulate the pharmacokinetics of many drugs via mediating the translocation of drugs across biologic membranes. SLC transporters are important drug targets involved in the uptake of a broad range of compounds across the membrane. However, high-resolution experimental structures have been reported for a very limited number of transporters, which limits the study of their physiologic functions. In this review, we collected structural information on ABC and SLC transporters and described the application of computational methods in structure prediction. Taking P-glycoprotein (ABCB1) and serotonin transporter (SLC6A4) as examples, we assessed the pivotal role of structure in transport mechanisms, details of ligand-receptor interactions, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms. The data collected contributes toward safer and more effective pharmacological treatments. SIGNIFICANCE STATEMENT: The experimental structure of ATP-binding cassette and solute carrier transporters was collected, and the application of computational methods in structure prediction was described. P-glycoprotein and serotonin transporter were used as examples to reveal the pivotal role of structure in transport mechanisms, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms.

Active Binding Modes of Caffeine with Cytochrome P450 1A2 Determine Its Metabolite Profiles.

Caffeine is a very common kind of nervous stimulant, and it is primarily metabolized by Cytochrome P450 1A2 (CYP1A2) in the human body. Over the years, determining the interactions between caffeine and CYP1A2 has been a tough issue. The active binding modes and the catalytic regioselectivity of the metabolism between CYP1A2 and caffeine remain unclear. Here, to investigate the interactions between CYP1A2 and caffeine, we constructed the all-sequence CYP1A2-caffeine-membrane system using a multiple template approach. According to our simulation results, four active binding modes between CYP1A2 and caffeine that correspond to the four metabolic sites of caffeine are determined. What is more, a pre-reaction state for the CYP1A2-catalyzed reaction at caffeine's N3 site is identified. A more preponderant active binding mode might be the reason why the N3 site of caffeine becomes the primary metabolic site. Our findings could enhance our knowledge of the interactions between CYP1A2 and caffeine and help us better understand the regioselectivity of the metabolism between CYP1A2 and caffeine.

Dissecting the Effect of Temperature on Hyperthermophilic Pf2001 Esterase Dimerization by Molecular Dynamics.

Pf2001 esterase (Pf2001) from Pyrococcus furiosus has hyperthermophilic properties and exerts a biocatalytic function in a dimeric state. Crystal structures revealed that the structural rearrangement of the cap domain is responsible for the Pf2001 dimer formation. However, the details of the cap domain remodeling and the effects of temperature on the dimerization process remain elusive at the molecular level, taking into account that experimental methods are difficult to capture the dynamic process of dimerization to some extent. Herein, four dimer models based on the monomeric crystal structure (PDB ID: 5G59) were constructed to investigate the conformational transition details and temperature effects in the dimerization by conventional molecular dynamics and accelerated molecular dynamics simulations. Our simulation results indicate that the monomer undergoes a conformational change into a "preparatory state" at high temperatures, which is more favorable for its transformation into a stable dimer. The subsequent free energy landscape analysis further identifies four intermediate states (from separated state to dimeric state) and discloses that a more accessible α-helix driven by stronger hydrophobic interactions induces a rearrangement of the cap domain, displaying a "tic-tac-toe" activation feature that is important for stabilizing the dimer interface and facilitating the formation of hydrophobic pockets. In addition, the electrostatic potential surface analysis illustrates that the weaker electrostatic repulsion (Lys and Arg) in the dimer interface at high temperatures is also a key factor for dimer stabilization. Altogether, our results can provide molecular-level insight into the dimer formation process of hyperthermophilic esterase and would be useful to understand the enzymatic specificity of α/ß-hydrolase.

A direct approach toward investigating DNA-ligand interactions via surface-enhanced Raman spectroscopy combined with molecular dynamics simulations.

Small molecules that interfere with DNA replication can trigger genomic instability, which makes these molecules valuable in the search for anticancer drugs. Thus, interactions between DNA and its ligands at the molecular level are of great significance. In the present study, a new method based on surface-enhanced Raman spectroscopy (SERS) combined with molecular dynamics simulations has been proposed for analyzing the interactions between DNA and its ligands. The SERS signals of DNA hairpins (ST: d(CGACCAACGTGTCGCCTGGTCG), AP1: d(CGCACAACGTGTCGCCTGTGCG)), pure argininamide, and their complexes, were obtained, and the characteristic peak sites of the DNA secondary structure and argininamide ligand-binding region were analyzed. Molecular dynamics calculations predicted that argininamide binds to the 8C and 9G bases of AP1 via hydrogen bonding. Our method successfully detected the changes of SERS fingerprint peaks of hydrogen bonds and bases between argininamide and DNA hairpin bases, and their binding sites and action modes were consistent with the predicted results of the molecular dynamics simulations. This SERS technology combined with the molecular dynamics simulation detection platform provides a general analysis tool, with the advantage of effective, rapid, and sensitive detection. This platform can obtain sufficient molecular level conformational information to provide avenues for rapid drug screening and promote progress in several fields, including targeted drug design.

How does multiple substrate binding lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan? A theoretical study.

CYP2D6 is one of the most important metalloenzymes involved in the biodegradation of many drug molecules in the human body. It has been found that multiple substrate binding can lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan (DM), but the corresponding theoretical mechanism is rarely reported. Therefore, we chose DM as the probe and performed molecular dynamics simulations and quantum mechanical calculations on CYP2D6-DM systems to investigate the mechanism of how the multiple substrate binding leads to the substrate inhibition of CYP2D6 metabolizing substrates. According to our results, three gate residues (Arg221, Val374, and Phe483) for the catalytic pocket are determined. We also found that the multiple substrate binding can lead to substrate inhibition by reducing the stability of CYP2D6 binding DM and increasing the reactive activation energy of the rate-determining step. Our findings would help to understand the substrate inhibition of CYP2D6 metabolizing the DM and enrich the knowledge of the drug-drug interactions for the cytochrome P450 superfamily.

Differences of Atomic-Level Interactions between Midazolam and Two CYP Isoforms 3A4 and 3A5.

CYP 3A4 and CYP 3A5 are two important members of the human cytochrome P450 family. Although their overall structures are similar, the local structures of the active site are different, which directly leads to obvious individual differences in drug metabolic efficacy and toxicity. In this work, midazolam (MDZ) was selected as the probe substrate, and its interaction with two proteins, CYP 3A4 and CYP 3A5, was studied by molecular dynamics simulation (MD) along with the calculation of the binding free energy. The results show that two protein-substrate complexes have some similarities in enzyme-substrate binding; that is, in both complexes, Ser119 forms a high occupancy hydrogen bond with MDZ, which plays a key role in the stability of the interaction between MDZ and the enzymes. However, the complex formed by CYP 3A4 and MDZ is more stable, which may be attributed to the sandwich structure formed by the fluorophenyl group of the substrate with Leu216 and Leu482. Our study interprets the binding differences between two isoform-substrate complexes and reveals a structure-function relationship from the atomic perspective, which is expected to provide a theoretical basis for accurately measuring the effectiveness and toxicity of drugs for individuals in the era of precision medicine.

Molecular Insights into the Heterotropic Allosteric Mechanism in Cytochrome P450 3A4-Mediated Midazolam Metabolism.

Cytochrome P450 3A4 (CYP3A4) is the main P450 enzyme for drug metabolism and drug-drug interactions (DDIs), as it is involved in the metabolic process of approximately 50% of drugs. A detailed mechanistic elucidation of DDIs mediated by CYP3A4 is commonly believed to be critical for drug optimization and rational use. Here, two typical probes, midazolam (MDZ, substrate) and testosterone (TST, allosteric effector), are used to investigate the molecular mechanism of CYP3A4-mediated heterotropic allosteric interactions, through conventional molecular dynamics (cMD) and well-tempered metadynamics (WT-MTD) simulations. Distance monitoring shows that TST can stably bind in two potential peripheral sites (Site 1 and Site 2) of CYP3A4. The binding of TST at these two sites can induce conformational changes in CYP3A4 flexible loops on the basis of conformational analysis, thereby promoting the transition of the MDZ binding mode and affecting the ratio of MDZ metabolites. According to the results of the residue interaction network, multiple allosteric communication pathways are identified that can provide vivid and applicable insights into the heterotropic allostery of TST on MDZ metabolism. Comparing the regulatory effects and the communication pathways, the allosteric effect caused by TST binding in Site 2 seems to be more pronounced than in Site 1. Our findings could provide a deeper understanding of CYP3A4-mediated heterotropic allostery at the atomic level and would be helpful for rational drug use as well as the design of new allosteric modulators.

The regioselectivity of the interaction between dextromethorphan and CYP2D6.

CYP2D6 is an important enzyme of the cytochrome P450 superfamily, and catalyzes nearly 25% of the drugs sold in the market. For decades, the interactions and metabolism between CYP2D6 and substrates have been a hot topic. However, the key factors of the catalytic regioselectivity for CYP2D6 still remain controversial. Here, we construct four systems to explore the interaction between dextromethorphan (DM) and CYP2D6. A new binding mode of CYP2D6 is defined, and two key residues (residue Asp301 and residue Glu216) are discovered working simultaneously to stabilize the DM at the reactive site by forming water bridge hydrogen bonds when CYP2D6 binds DM. Our results also indicate that the substrate concentration could mediate the binding mode between the substrate and CYP2D6 by decreasing the volume of the catalytic pocket, which is not conducive to the O-demethylation of DM but benefits the N-demethylation of DM. These results could shed light on the process of CYP2D6 binding to the substrate, and help to better understand the regioselectivity of CYP2D6 catalyzing the substrates.

Investigation of the molecular and mechanistic basis for the regioselective metabolism of midazolam by cytochrome P450 3A4.

Cytochrome P450 3A4 (CYP3A4) is the most important P450 enzyme for drug metabolism and drug-drug interaction, due to it being responsible for the biotransformation of approximately 50% of clinically used drugs. Advance knowledge of the molecular and mechanistic basis of CYP3A4 regioselective metabolism is beneficial for understanding the production of metabolites, and may allow personalized metabolic pathways or designing pathway-specific therapeutics. In this work, we focus on investigating the ligand-receptor interactions, substrate conformational transition, and key factors regulating the specificity of metabolic pathways using midazolam (MDZ) as a probe. Here, three types of substrate-binding conformations related to the diversity of MDZ metabolites are identified. The results also suggest that an allosteric site for MDZ is located near the F'-helix, A-anchor, and C-terminal loop of CYP3A4. The presence of an effector in the allosteric site can accelerate the conformational transition of the substrate via modulating a "sandwich" structure, and may affect the proportion of metabolites at high substrate concentration. We hope that the results can improve the understanding of the CYP3A4 structure and function, and provide a new perspective for drug development.

Multiple Molecular Dynamics Simulations and Free-Energy Predictions Uncover the Susceptibility of Variants of HIV-1 Protease against Inhibitors Darunavir and KNI-1657.

HIV-1 protease (PR) is considered to be the main targets of anti-AIDS drug design because of its role in the proteolytic processing of viral polyproteins. However, the emergence of drug-resistant HIV has become a major problem in the therapy of HIV-1-infected patients. Focused on the complexes of wild type (WT) PR and two mutant PRs (V32I/L33F/I54M/V82I and V32I/L33F/I54M/I84 V) with inhibitors Darunavir (DRV) and KNI-1657 (KNI), respectively, we have conducted research on the conformational dynamics and the resistance mechanism caused by residue mutations through multiple molecular dynamics (MD) simulations combined with an energy (MM-PBSA and solvated interaction energy (SIE)) prediction. The results indicate that mutated residues of PR alter the distance between flap regions and catalytic sites, the volume of the inner catalytic site, and the curling degree of the flap tips, thereby affecting DRV and KNI inhibitor binding to PR. These mutated residues reduced the binding affinity of the two mutant PRs to DRV, resulting in drug resistance, whereas the two mutant PRs increase the binding affinity with KNI, indicating they enhance the sensitivity to KNI. Compared with the WT PR, the changes in van der Waals interaction and electrostatic interaction in the two variant PRs play a vital part in the binding of PR with DRV and KNI. These results may supply valuable guidance for the design of anti-AIDS drugs targeting PR.

How DNA affects the hyperthermophilic protein Ape10b2 for oligomerization: an investigation using multiple short molecular dynamics simulations.

Alba2 is a hyperthermophilic DNA-binding protein, and DNA plays a crucial role in the Alba2 oligomerization process. It is a pity that there is limited research in terms of how DNA affects the conformational change of Alba2 in oligomerization. Herein, we complement the crystal structure of the Ape10b2 (belongs to Alba2)-dsDNA complex (PDB ID: 3U6Y) and employ multiple short molecular dynamics (MSMD) simulations to illuminate the influence of DNA on Ape10b2 at four temperatures (300, 343, 363, and 373 K). Our results indicate that DNA could cause the conformational changes of two important regions (loop1 and loop5), which may be beneficial for protein oligomerization. The results of hydrogen bond analysis show that the increasing number of hydrogen bonds between two monomers of Ape10b2 may also be a favorable factor for oligomerization. In addition, Ape10b2 can stabilize DNA by electrostatic interactions with an increase in temperature, and five residues (Arg40, Arg42, Asn43, Asn45, and Arg46) play a stabilizing role during protein binding to DNA. Our findings could help in understanding the favorable factors leading to protein oligomerization, which contributes to enzyme engineering research from an industrial perspective.

A theoretical study on the signal transduction process of bacterial photoreceptor PpSB1 based on the Markov state model.

Light-oxygen-voltage (LOV) domains are blue light sensors and play an important role in signal transduction in many organisms. Generally, LOV domains use chromophores to absorb photons, and then photochemical reactions will occur to convert light energy into chemical energy and transduce it to the main chain of proteins. These reactions can cause conformational rearrangement of proteins, and thus leading to signal transduction. Therefore, it is important to study the signal transduction process of LOV domains for understanding the control mechanism of cellular functions. However, how small photochemical changes in the active sites of the LOV domains lead to large conformational rearrangements of proteins, which in turn lead to signal transduction, has been puzzling us for a long time. Currently, the LOV domains are mainly studied in plants. The signal transduction mechanism of LOV domains in bacteria is still unclear. In this work, the Markov state model (MSM) combined with molecular dynamics (MD) simulations was applied to investigate the signal transduction process of the LOV protein from pseudomonas putida (PpSB1-LOV). The present work will play an important role in understanding the signal transduction mechanism of PpSB1-LOV domains, which may provide theoretical basis for the design and improvement of LOV-based optogenetic tools.

Studies of Interaction Mechanism between Pyrido [3,4-d] Pyrimidine Inhibitors and Mps1.

Monopolar spindle 1 (Mps1), a dual-specific kinase, is related to the proper execution of chromosome biorientation and mitotic checkpoint signaling. The overexpression of Mps1 promotes the occurrence of cancer or the survival of aneuploid cancer cells, in other words, the reduction of Mps1 will severely reduce the viability of human cancer cells. Therefore, Mps1 is a potential target for cancer treatment. Recently, a series of novel pyrido [3,4-d] pyrimidine derivatives targeting Mps1 with high biological activity were synthesized. The crystal structure of Mps1 in complex with pyrido [3,4-d] pyrimidine derivatives was also reported, but there were no specific mechanism studies for this series of small molecule inhibitors. In this study, complexes binding modes were probed by molecular docking and further validated by molecular dynamics simulations and the molecular mechanics/generalized Born surface area (MM/GBSA) method. The results indicated that the van der Waals interactions and the nonpolar solvation energies were responsible to the basis for favorable binding free energies, all inhibitors interacted with residues I531, V539, M602, C604, N606, I607, L654, I663, and P673 of Mps1. By analyzing the hydrogen bonds, we found the residues G605 and K529 in Mps1 formed stable hydrogen bonds with compounds, it was more conducive to activities of Mps1 inhibitors. According to the above analysis, we further designed five new compounds. We found that compounds IV and V were better potential Mps1 inhibitors through docking and ADMET prediction. The obtained new insights not only were helpful in understanding the binding mode of inhibitors in Mps1, but also provided important references for further rational design of Mps1 inhibitors.

Exploring the Distinct Binding and Activation Mechanisms for Different CagA Oncoproteins and SHP2 by Molecular Dynamics Simulations.

CagA is a major virulence factor of Helicobacter pylori. H. pylori CagA is geographically subclassified into East Asian CagA and Western CagA, which are characterized by the presence of a EPIYA-D or EPIYA-C segment. The East Asian CagA is more closely associated with gastric cancer than the Western CagA. In this study, molecular dynamic (MD) simulations were performed to investigate the binding details of SHP2 and EPIYA segments, and to explore the allosteric regulation mechanism of SHP2. Our results show that the EPIYA-D has a stronger binding affinity to the N-SH2 domain of SHP2 than EPIYA-C. In addition, a single EPIYA-D binding to N-SH2 domain of SHP2 can cause a deflection of the key helix B, and the deflected helix B could squeeze the N-SH2 and PTP domains to break the autoinhibition pocket of SHP2. However, a single EPIYA-C binding to the N-SH2 domain of SHP2 cannot break the autoinhibition of SHP2 because the secondary structure of the key helix B is destroyed. However, the tandem EPIYA-C not only increases its binding affinity to SHP2, but also does not significantly break the secondary structure of the key helix B. Our study can help us better understand the mechanism of gastric cancer caused by Helicobacter pylori infection.

Multiple Molecular Dynamics Simulations of the Inhibitor GRL-02031 Complex with Wild Type and Mutant HIV-1 Protease Reveal the Binding and Drug-Resistance Mechanism.

Human immunodeficiency virus type 1 (HIV-1) protease is regarded as a fascinating target for drug development against HIV infection. However, mutations causing drug resistance severely limit the efficiency of the recently marketed drugs in the treatment of HIV replication. To elucidate the binding mechanism of HIV-1 protease with promising inhibitor GRL-02031 and further to probe the resistance mechanism associated with mutations (I47V, L76V, V82A, and N88D) to the inhibitor, we applied multiple molecular dynamics (MMD) simulations along with energy analysis by the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methodology on specific HIV-1 protease with GRL-0231 complexes. On the basis of detail analysis of the simulations, we revealed key characteristics that constitute the drug resistance of four mutation HIV-1 proteases toward GRL-02031: substitution of the side chain in these four mutation residues leads to a change in the distances between the flaps and catalytic sites, thereby reducing the affinity for GRL-02031 with these four mutation proteases, even though the L76V and N88D residues cannot directly contact GRL-02031. The results of energy analysis according to the MM-PBSA and SIE methods further indicated that hydrophobic interaction was considered to be the prime driving force for inhibitor GRL-02031 binding to protease and the decrease in van der Waals interactions between inhibitor GRL-02031 and mutant proteases as the primary cause of the drug resistance. Analyses of the hydrogen bonds and atomic interactions further provided detailed explanations for the resistance of these four mutation proteases toward inhibitor GRL-02031. The present study provides potential guidance on the structure-based inhibitors' design targeting HIV-1 protease.