Binding modes of three inhibitors 8CA, F8A and I4A to A-FABP studied based on molecular dynamics simulation

Revealing the Role of the Arg and Lys in Shifting Paradigm from BTK Selective Inhibition to the BTK/HCK Dual Inhibition - Delving into the Inhibitory Activity of KIN-8194 against BTK, and HCK in the Treatment of Mutated BTKCys481 Waldenström Macroglobulinemia: A Computational Approach

BACKGROUND

Despite the early success of Bruton's tyrosine kinase (BTK) inhibitors in the treatment of Waldenström macroglobulinemia (WM), these single-target drug therapies have limitations in their clinical applications, such as drug resistance. Several alternative strategies have been developed, including the use of dual inhibitors, to maximize the therapeutic potential of these drugs.

OBJECTIVE

Recently, the pharmacological activity of KIN-8194 was repurposed to serve as a 'dual-target' inhibitor of BTK and Hematopoietic Cell Kinase (HCK). However, the structural dual inhibitory mechanism remains unexplored, hence the aim of this study.

METHODS

Conducting predictive pharmacokinetic profiling of KIN-8194, as well as demonstrating a comparative structural mechanism of inhibition against the above-mentioned enzymes.

RESULTS

Our results revealed favourable binding affinities of -20.17 kcal/mol, and -35.82 kcal/mol for KIN-8194 towards HCK and BTK, respectively. Catalytic residues Arg137/174 and Lys42/170 in BTK and Arg303 and Lys75/173/244/247 in HCK were identified as crucial mediators of the dual binding mechanism of KIN-8194, corroborated by high per-residue energy contributions and consistent high-affinity interactions of these residues. Prediction of the pharmacokinetics and physicochemical properties of KIN-8194 further established its inhibitory potential, evidenced by the favourable absorption, metabolism, excretion, and minimal toxicity properties. Structurally, KIN-8194 impacted the stability, flexibility, solvent-accessible surface area, and rigidity of BTK and HCK, characterized by various alterations observed in the bound and unbound structures, which proved enough to disrupt their biological function.

CONCLUSION

These structural insights provided a baseline for the understanding of the dual inhibitory activity of KIN- 8194. Establishing the cruciality of the interactions between the KIN-8194 and Arg and Lys residues could guide the structure-based design of novel dual BTK/HCK inhibitors with improved therapeutic activities.

Computationally Designed Anti-LuxP DNA Aptamer Suppressed Flagellar Assembly- and Quorum Sensing-Related Gene Expression in Vibrio parahaemolyticus

(1) Background: Quorum sensing (QS) is the chemical communication between bacteria that sense chemical signals in the bacterial population to control phenotypic changes through the regulation of gene expression. The inhibition of QS has various potential applications, particularly in the prevention of bacterial infection. QS can be inhibited by targeting the LuxP, a periplasmic receptor protein that is involved in the sensing of the QS signaling molecule known as the autoinducer 2 (AI-2). The sensing of AI-2 by LuxP transduces the chemical information through the inner membrane sensor kinase LuxQ protein and activates the QS cascade. (2) Methods: An in silico approach was applied to design DNA aptamers against LuxP in this study. A method combining molecular docking and molecular dynamics simulations was used to select the oligonucleotides that bind to LuxP, which were then further characterized using isothermal titration calorimetry. Subsequently, the bioactivity of the selected aptamer was examined through comparative transcriptome analysis. (3) Results: Two aptamer candidates were identified from the ITC, which have the lowest dissociation constants (Kd) of 0.2 and 0.5 micromolar. The aptamer with the lowest Kd demonstrated QS suppression and down-regulated the flagellar-assembly-related gene expression. (4) Conclusions: This study developed an in silico approach to design an aptamer that possesses anti-QS properties.

Co-Binding of JQ1 and Venetoclax Exhibited Synergetic Inhibitory Effect for Cancer Therapy; Potential Line of Treatment for the Waldenström Macroglobulinemia Lymphoma

In recent times, the development of combination therapy has been a focal point in drug discovery. This article explores the potential synergistic effect of co-administration of Bcl2 inhibitor Venetoclax and BET inhibitor JQ1. We envisioned that the 'dual-site'-binding of Bcl2 has significant prospects and paves the way for the next round of rational design of potent Waldenström macroglobulinemia (WM) therapy. The preferential binding mechanisms of the multi-catalytic sites of the Bcl2 enzyme have been a subject of debate in the literature. This study conducted a systematic procedure to explore the preferred binding modes and the structural effects of co-binding at each catalytic active site. Interestingly, a mutual enhanced binding effect was observed - Venetoclax increased the binding affinity of JQ1 by 11.5 %, while JQ1 boosted the binding affinity of Venetoclax by 16.3 % when compared with individual inhibition of each drug. This synergistic binding effect has significantly increased protein stability, with substantial correlated movements and multiple van der Waals interactions. The structural and thermodynamic insights unveiled in this report would assist the future design of improved combined therapy against WM.

Multi-catalytic Sites Inhibition of Bcl2 Induces Expanding of Hydrophobic Groove: A New Avenue Towards Waldenström Macroglobulinemia Therapy

B-cell lymphoma 2 (Bcl2) is a key protein regulator of apoptosis. The hydrophobic groove in Bcl2 is a unique structural feature to this class of enzymes and found to have a profound impact on protein overall structure, function, and dynamics. Dynamics of the hydrophobic groove is an essential determinant of the catalytic activity of Bcl2, an implicated protein in Waldenström macroglobulinemia (WM). The mobility of α3-α4 helices around the catalytic site of the protein remains crucial to its activity. The preferential binding mechanisms of the multi-catalytic sites of the Bcl2 enzyme have been a subject of debate in the literature. In addition to our previous report on the same protein, herein, we further investigate the preferential binding modes and the conformational implications of Venetoclax-JQ1 dual drug binding at both catalytic active sites of Bcl2. Structural analysis revealed asymmetric α3-α4 helices movement with the expansion of the distance between the α3 and α4 helix in Venetoclax-JQ1 dual inhibition by 15.2% and 26.3%, respectively when compared to JQ1 and Venetoclax individual drug inhibition-resulting in remarkable widening of hydrophobic groove. Moreso, a reciprocal enhanced binding effect was observed: Venetoclax increased the binding affinity of JQ1 by 11.5%, while the JQ1 fostered the binding affinity of Venetoclax by 16.3% compared with individual inhibition of each drug. This divergence has also resulted in higher protein stability, and prominent correlated motions were observed with the least fluctuations and multiple van der Waals interactions. Findings offer vital conformational dynamics and structural mechanisms of enzyme-single ligand and enzyme-dual ligand interactions, which could potentially shift the current therapeutic protocol of Waldenström macroglobulinemia.

Unveiling the Role of the Fatty Acid Binding Protein 4 in the Metabolic-Associated Fatty Liver Disease

Metabolic-associated fatty liver disease (MAFLD), the main cause of chronic liver disease worldwide, is a progressive disease ranging from fatty liver to steatohepatitis (metabolic-associated steatohepatitis; MASH). Nevertheless, it remains underdiagnosed due to the lack of effective non-invasive methods for its diagnosis and staging. Although MAFLD has been found in lean individuals, it is closely associated with obesity-related conditions. Adipose tissue is the main source of liver triglycerides and adipocytes act as endocrine organs releasing a large number of adipokines and pro-inflammatory mediators involved in MAFLD progression into bloodstream. Among the adipocyte-derived molecules, fatty acid binding protein 4 (FABP4) has been recently associated with fatty liver and additional features of advanced stages of MAFLD. Additionally, emerging data from preclinical studies propose FABP4 as a causal actor involved in the disease progression, rather than a mere biomarker for the disease. Therefore, the FABP4 regulation could be considered as a potential therapeutic strategy to MAFLD. Here, we review the current knowledge of FABP4 in MAFLD, as well as its potential role as a therapeutic target for this disease.

Impact of non-synonymous mutations on the structure and function of telomeric repeat binding factor 1

Telomeric repeat binding factor 1 (TRF1) is one of the major components of the shelterin complex. It directly binds to the telomere and controls its function by regulating the telomerase acting on it. Several variations are reported in the TRF1 gene; some are associated with variety of diseases. Here, we have studied the structural and functional significance of these variations in the TRFH domain of TRF1. We have used cutting-edge computational methods such as SIFT, PolyPhen-2, PROVEAN, Mutation Assessor, mCSM, SDM, STRUM, MAESTRO, and DUET to predict the effects of 124 mutations in the TRFH domain of TRF1. Out of 124 mutations, we have identified 12 deleterious mutations with high confidence based on their prediction. To see the impact of the finally selected mutations on the structure and stability of TRF1, all-atom molecular dynamics (MD) simulations on TRF1-Wild type (WT), L79R and P150R mutants for 200 ns were carried out. A significant conformational change in the structure of the P150R mutant was observed. Our integrated computational study provides a comprehensive understanding of structural changes in TRF1 incurred due to the mutations and subsequent function, leading to the progression of many diseases.Communicated by Ramaswamy H. Sarma.

Molecular mechanisms of inhibitor bindings to A-FABP deciphered by using molecular dynamics simulations and calculations of MM-GBSA

Adipocyte fatty-acid binding protein (A-FABP) plays a central role in many aspects of metabolic diseases. It is an important target in drug design for treatment of FABP-related diseases. In this study, molecular dynamics (MD) simulations followed by calculations of molecular mechanics generalized Born surface area (MM-GBSA) and principal components analysis (PCA) were implemented to decipher molecular mechanism correlating with binding of inhibitors 57Q, 57P and L96 to A-FABP. The results show that van der Waals interactions are the leading factors to control associations of 57Q, 57P, and L96 with A-FABP, which reveals an energetic basis for designing of clinically available inhibitors towards A-FABP. The information from PCA and cross-correlation analysis rationally unveils that inhibitor bindings affect conformational changes of A-FABP and change relative movements between residues. Decomposition of binding affinity into contributions of individual residues not only detects hot spots of inhibitor/A-FABP binding but also shows that polar interactions of the positively charged residue Arg126 with three inhibitors provide a significant contribution for stabilization of the inhibitor/A-FABP bindings. Furthermore, the binding strength of L96 to residues Ser55, Phe57 and Lys58 are stronger than that of inhibitors 57Q and 57P to these residues.

ABBV-744 as a potential inhibitor of SARS-CoV-2 main protease enzyme against COVID-19

A new pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide and become pandemic with thousands new deaths and infected cases globally. To address coronavirus disease (COVID-19), currently no effective drug or vaccine is available. This necessity motivated us to explore potential lead compounds by considering drug repurposing approach targeting main protease (Mpro) enzyme of SARS-CoV-2. This enzyme considered to be an attractive drug target as it contributes significantly in mediating viral replication and transcription. Herein, comprehensive computational investigations were performed to identify potential inhibitors of SARS-CoV-2 Mpro enzyme. The structure-based pharmacophore modeling was developed based on the co-crystallized structure of the enzyme with its biological active inhibitor. The generated hypotheses were applied for virtual screening based PhaseScore. Docking based virtual screening workflow was used to generate hit compounds using HTVS, SP and XP based Glide GScore. The pharmacological and physicochemical properties of the selected lead compounds were characterized using ADMET. Molecular dynamics simulations were performed to explore the binding affinities of the considered lead compounds. Binding energies revealed that compound ABBV-744 binds to the Mpro with strong affinity (ΔGbind -45.43 kcal/mol), and the complex is more stable in comparison with other protein-ligand complexes. Our study classified three best compounds which could be considered as promising inhibitors against main protease SARS-CoV-2 virus.

Molecular dynamics simulation of the heart type fatty acid binding protein in a crystal environment

Crystallographic data comes from a space-time average over all the unit cells within the crystal, so dynamic phenomena do not contribute significantly to the diffraction data. Many efforts have been made to reconstitute the movement of the macromolecules and explore the microstates that the confined proteins can adopt in the crystalline network. We explored different strategies to simulate a heart fatty acid binding protein (H-FABP) crystal by means of Molecular Dynamics (MD) simulations. We evaluate the effect of introducing restraints according to experimental isotropic B-factors and we analyzed the H-FABP motions in the crystal using Principal Component Analysis (PCA), isotropic and anisotropic B-factors. We compared the behavior of the protein simulated in the crystal confinement versus in solution, and we observed the effect of that confinement in the mobility of the protein residues. Restraining one-third of Cα atoms based on experimental B-factors produce lower B-factors than simulations without restraints, showing that the position restraint of the atoms with the lowest experimental B-factor is a good strategy to maintain the geometry of the crystal with an obvious decrease in the degrees of motion of the protein. PCA shows that, as position restraint reduces the conformational space explored by the system, the motion of the crystal is better recovered, for an essential subspace of the same size, in the simulations without restraints. Restraining only one Cα seems to be a good balance between giving flexibility to the system and preserving its structure. Communicated by Ramaswamy H. Sarma.

Molecular mechanism with regard to the binding selectivity of inhibitors toward FABP5 and FABP7 explored by multiple short molecular dynamics simulations and free energy analyses

Recently, fatty acid binding proteins 5 and 7 (FABP5 and FABP7) have been regarded as the prospective targets for clinically treating multiple diseases related to FABPs.

Evaluation of potential molecular interaction between quorum sensing receptor, LuxP and grouper fatty acids: in-silico screening and simulation

Pathologically relevant behaviors of Vibrio, such as the expression of virulence factors, biofilm production, and swarming motility, have been shown to be controlled by quorum sensing. The autoinducer-2 quorum sensing receptor protein LuxP is one of the target proteins for drug development to suppress the virulence of Vibrio. Here, we reported the potential molecular interaction of fatty acids identified in vibriosis-resistant grouper with LuxP. Fatty acid, 4-oxodocosahexaenoic acid (4R8) showed significant binding affinity toward LuxP (-6.0 kcal/mol) based on molecular docking analysis. The dynamic behavior of the protein-ligand complex was illustrated by molecular dynamic simulations. The fluctuation of the protein backbone, the stability of ligand binding, and hydrogen bond interactions were assessed, suggesting 4R8 possesses potential interaction with LuxP, which was supported by the low binding free energy (-29.144 kJ/mol) calculated using the molecular mechanics Poisson-Boltzmann surface area.

Exploring the flap dynamics of the South African HIV subtype C protease in presence of FDA-approved inhibitors: MD study

HIV-1 protease (HIV PR) is considered as one of the most attractive targets for the treatment of HIV and the impact of flap dynamics of HIV PR on the binding affinities of protease inhibitors (PIs) is a crucial ongoing research field. Recently, our research group evaluated the binding affinities of different FDA approved PIs against the South African HIV-1 subtype C (C-SA) protease (PR). The CSA-HIV PR displayed weaker binding affinity for most of the clinical PIs compared to HIV-1 B subtype for West and Central Europe, the Americas. In the current work, the flap dynamics of four different systems of HIV-1 C-SA PR complexed to FDA approved second generation PIs and its impact on binding was explored over the molecular dynamic trajectories. It was observed that the interactions of the selected drugs with the binding site residues of the protease may not be the major contributor for affinity towards PIs. Various post-MD analyses were performed, also entropic contributions, solvation free energies and hydrophobic core formation interactions were studied to assess how the flap dynamics of C-SA PR which is affected by such factors. From these contributions, large van der Waals interactions and low solvation free energies were found to be major factors for the higher activity of ATV against C-SA HIV PR. Furthermore, a comparatively stable hydrophobic core may be responsible for higher stability of the PR flaps of the ATV complex. The outcome of this study provides significant guidance to how the flap dynamics of C-SA PR is affected by various factors as a result of the binding affinity of various protease inhibitors. It will also assist with the design of potent inhibitors against C-SA HIV PR that apart from binding in the active site of PR can interacts with the flaps to prevent opening of the flaps resulting in inactivation of the protease.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein, retinoic acid ligand, and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.

Role of the fatty acid-binding protein 4 in heart failure and cardiovascular disease

Obesity and ectopic fat accumulation in non-adipose tissues are major contributors to heart failure (HF) and cardiovascular disease (CVD). Adipocytes act as endocrine organs by releasing a large number of bioactive molecules into the bloodstream, which participate in a communication network between white adipose tissue and other organs, including the heart. Among these molecules, fatty acid-binding protein 4 (FABP4) has recently been shown to increase cardiometabolic risk. Both clinical and experimental evidence have identified FABP4 as a relevant player in atherosclerosis and coronary artery disease, and it has been directly related to cardiac alterations such as left ventricular hypertrophy (LVH) and both systolic and diastolic cardiac dysfunction. The available interventional studies preclude the establishment of a direct causal role of this molecule in CVD and HF and propose FABP4 as a biomarker rather than as an aetiological factor. However, several experimental reports have suggested that FABP4 may act as a direct contributor to cardiac metabolism and physiopathology, and the pharmacological targeting of FABP4 may restore some of the metabolic alterations that are conducive to CVD and HF. Here, we review the current knowledge regarding FABP4 in the context of HF and CVD as well as the molecular basis by which this protein participates in the regulation of cardiac function.

Differential flap dynamics in l,d-transpeptidase2 from mycobacterium tuberculosis revealed by molecular dynamics

Despite the advances in tuberculosis treatment, TB is still one the most deadly infectious diseases and remains a major global health quandary. Mycobacterium tuberculosis (Mtb) is the only known mycobacterium with a high content of 3→3 crosslinks in the biosynthesis of peptidoglycan, which is negligible in most bacterial species. An Mtb lacking Ldt_Mt2 leads to alteration of the colony morphology and loss of virulence which makes this enzyme an attractive target. Regardless of the vital role of Ldt_Mt2 for cell wall survival, the impact of ligand binding on the dynamics of the β-hairpin flap is still unknown. Understanding the structural and dynamical behaviour of the flap regions provides clear insight into the design of the effective inhibitors against Ldt_Mt2. Carbapenems, an specific class of β-lactam family, have been shown to inactivate this enzyme. Herein a comprehensive investigation of the flap dynamics of Ldt_Mt2 complex with substrate and three carbapenems namely, ertapenem, imipenem and meropenem is discussed and analyzed for the first account using 140 ns molecular dynamics simulations. The structural features (RMSD, RMSF and R_g) derived by MD trajectories were analyzed. Distance analysis, particularly tip-tip SER135-ASN167 index, identified conformational changes in terms of flap opening and closure within binding process. Principal component analysis (PCA) was employed to qualitatively understand the divergent effects of different inhibitors on the dominant motion of each residue. To probe different internal dynamics induced by ligand binding, dynamic cross-correlation marix (DCCM) analysis was used. The binding free energies of the selected complexes were assessed using MM-GBSA method and per residue free energy decomposition analysis were performed to characterize the contribution of the key residues to the total binding free energies.

Large-scale molecular dynamics simulation: Effect of polarization on thrombin-ligand binding energy

Molecular dynamics (MD) simulations lasting 500 ns were performed in explicit water to investigate the effect of polarization on the binding of ligands to human α-thrombin based on the standard nonpolarizable AMBER force field and the quantum-derived polarized protein-specific charge (PPC). The PPC includes the electronic polarization effect of the thrombin-ligand complex, which is absent in the standard force field. A detailed analysis and comparison of the results of the MD simulation with experimental data provided strong evidence that intra-protein, protein-ligand hydrogen bonds and the root-mean-square deviation of backbone atoms were significantly stabilized through electronic polarization. Specifically, two critical hydrogen bonds between thrombin and the ligand were broken at approximately 190 ns when AMBER force field was used and the number of intra-protein backbone hydrogen bonds was higher under PPC than under AMBER. The thrombin-ligand binding energy was computed using the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, and the results were consistent with the experimental value obtained using PPC. Because hydrogen bonds were unstable, it was failed to predict the binding affinity under the AMBER force field. Furthermore, the results of the present study revealed that differences in the binding free energy between AMBER and PPC almost comes from the electrostatic interaction. Thus, this study provides evidence that protein polarization is critical to accurately describe protein-ligand binding.

The discovery of novel and selective fatty acid binding protein 4 inhibitors by virtual screening and biological evaluation

Adipocyte fatty acid binding protein (AFABP, FABP4) has been proven to be a potential therapeutic target for diabetes, atherosclerosis and inflammation-related diseases. In this study, a series of new scaffolds of small molecule inhibitors of FABP4 were identified by virtual screening and were validated by a bioassay. Fifty selected compounds were tested, which led to the discovery of seven hits. Structural similarity-based searches were then performed based on the hits and led to the identification of one high affinity compound 33b (Ki=0.29±0.07μM, ΔTm=8.5°C). This compound's effective blockade of inflammatory response was further validated by its ability to suppress pro-inflammatory cytokines induced by lipopolysaccharide (LPS) stimulation. Molecular dynamics simulation (MD) and mutagenesis studies validated key residues for its inhibitory potency and thus provide an important clue for the further development of drugs.

Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions

Molecular dynamics simulation and binding free energy calculations were used to probe drug resistance of HIV-1 protease mutations toward inhibitors.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Probing Difference in Binding Modes of Inhibitors to MDMX by Molecular Dynamics Simulations and Different Free Energy Methods

The p53-MDMX interaction has attracted extensive attention of anti-cancer drug development in recent years. This current work adopted molecular dynamics (MD) simulations and cross-correlation analysis to investigate conformation changes of MDMX caused by inhibitor bindings. The obtained information indicates that the binding cleft of MDMX undergoes a large conformational change and the dynamic behavior of residues obviously change by the presence of different structural inhibitors. Two different methods of binding free energy predictions were employed to carry out a comparable insight into binding mechanisms of four inhibitors PMI, pDI, WK23 and WW8 to MDMX. The data show that the main factor controlling the inhibitor bindings to MDMX arises from van der Waals interactions. The binding free energies were further divided into contribution of each residue and the derived information gives a conclusion that the hydrophobic interactions, such as CH-CH, CH-π and π-π interactions, are responsible for the inhibitor associations with MDMX.

Metabolic functions of FABPs--mechanisms and therapeutic implications

Intracellular and extracellular interactions with proteins enables the functional and mechanistic diversity of lipids. Fatty acid-binding proteins (FABPs) were originally described as intracellular proteins that can affect lipid fluxes, metabolism and signalling within cells. As the functions of this protein family have been further elucidated, it has become evident that they are critical mediators of metabolism and inflammatory processes, both locally and systemically, and therefore are potential therapeutic targets for immunometabolic diseases. In particular, genetic deficiency and small molecule-mediated inhibition of FABP4 (also known as aP2) and FABP5 can potently improve glucose homeostasis and reduce atherosclerosis in mouse models. Further research has shown that in addition to their intracellular roles, some FABPs are found outside the cells, and FABP4 undergoes regulated, vesicular secretion. The circulating form of FABP4 has crucial hormonal functions in systemic metabolism. In this Review we discuss the roles and regulation of both intracellular and extracellular FABP actions, highlighting new insights that might direct drug discovery efforts and opportunities for management of chronic metabolic diseases.

Revealing origin of decrease in potency of darunavir and amprenavir against HIV-2 relative to HIV-1 protease by molecular dynamics simulations

Clinical inhibitors Darunavir (DRV) and Amprenavir (APV) are less effective on HIV-2 protease (PR2) than on HIV-1 protease (PR1). To identify molecular basis associated with the lower inhibition, molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations were performed to investigate the effectiveness of the PR1 inhibitors DRV and APV against PR1/PR2. The rank of predicted binding free energies agrees with the experimental determined one. Moreover, our results show that two inhibitors bind less strongly to PR2 than to PR1, again in agreement with the experimental findings. The decrease in binding free energies for PR2 relative to PR1 is found to arise from the reduction of the van der Waals interactions induced by the structural adjustment of the triple mutant V32I, I47V and V82I. This result is further supported by the difference between the van der Waals interactions of inhibitors with each residue in PR2 and in PR1. The results from the principle component analysis suggest that inhibitor binding tends to make the flaps of PR2 close and the one of PR1 open. We expect that this study can theoretically provide significant guidance and dynamics information for the design of potent dual inhibitors targeting PR1/PR2.