Quantum and molecular dynamics study for binding of macrocyclic inhibitors to human alpha-thrombin

(1R,3S)-THCCA-Asn: To show the discovery of selective inhibitor of thrombin by successfully combining virtual screening and biological assay

Thrombin is the most potent platelet aggregator. To discover the selective inhibitor of thrombin that is important to curing platelet aggregation-related diseases, docking experiments were performed to dock (1R,3S)-2,3,4,9-tetrahydro-β-carboline-3- carboxylic acid, [(1R,3S)-THCCA], and (1S,3S)-2,3,4,9-tetrahydro-β-carboline-3- carboxylic acid, [(1S,3S)-THCCA], into the p pocket of bovine thrombin. The ideal match supported that (1R,3S)-THCCA could be used as a potential lead compound. In this case 20 natural amino acids were theoretically introduced into the 3-carboxyl of (1R,3S)-THCCA and 20 derivatives, (1R,3S)-THCCA-amino acids, were docked into p pocket of bovine thrombin to perform virtual screening. The screening revealed that comparing to (1R,3S)-THCCA itself the DockScores of 16 derivatives were higher, and (1R,3S)-THCCA-Asn (4j) got the highest DockScore. Thus, 16 derivatives were synthesized for experimental study. The in vitro anti-platelet aggregation assay showed that at 100 μM of concentration the 16 derivatives failed to inhibit the platelet aggregation induced by both adenosine diphosphate and arachidonic acid. On the other hand, however, the IC₅₀ value of the 16 derivatives inhibiting the platelet aggregation induced by platelet activating factor and thrombin ranged from 9.44 μM to 194.64 μM and from 0.07 μM to 9.56 μM, respectively. The in vitro anti-platelet aggregation assay suggested that the 16 derivatives selectively inhibited the platelet aggregation induced by thrombin. In particular, the IC₅₀ of (1R,3S)-THCCA-Asn (4j) had the lowest value. On rat model at 1 nmol/kg of dosage the 16 derivatives effectively prevented thrombus formation. It is worth pointing out that even at 0.01 nmol/kg of dosage, 4j still effectively prevented thrombus formation. 4j hardly has effects on the proliferation of mammalian cells and rat tail bleeding time. In conclusion, the combination of virtual screening and biological assays successfully lead to the discovery of 4j as a promising candidate of selective inhibitor of thrombin.

Fragmentation-Based Decomposition of a Metalloenzyme-Substrate Interaction: A Case Study for a Lytic Polysaccharide Monooxygenase

We present a novel decomposition scheme for electronic interaction energies based on the flexible formulation of fragmentation schemes through fragment combination ranges (FCRs; J. Chem. Phys., 2021, 155, 164105). We devise a clear additive decomposition with contribution of nondisjoint fragments and correction terms for overlapping fragments and apply this scheme to the metalloenzyme-substrate complex of a lytic polysaccharide monooxygenase (LPMO) with an oligosaccharide. By this, we further illustrate the straightforward adaptability of the FCR-based schemes to novel systems. Our calculations suggest that the description of the electronic structure is a larger error source than the fragmentation scheme. In particular, we find a large impact of the basis set size on the interaction energies. Still, the introduction of three-body interaction terms in the fragmentation setup improves the agreement to the supermolecular reference. Yet, the qualitative results for the decomposition scheme with two-body terms only largely agree within the investigated electronic-structure approaches and basis sets, which are B97-3c, DFT (TPSS and B3LYP), and MP2 methods. The overlap contributions are found to be small, allowing analysis of the interaction energy into individual amino acid residues: We find a particularly strong interaction between the substrate and the LPMO copper active site.

Exosite Binding in Thrombin: A Global Structural/Dynamic Overview of Complexes with Aptamers and Other Ligands

Thrombin is the key enzyme of the entire hemostatic process since it is able to exert both procoagulant and anticoagulant functions; therefore, it represents an attractive target for the developments of biomolecules with therapeutic potential. Thrombin can perform its many functional activities because of its ability to recognize a wide variety of substrates, inhibitors, and cofactors. These molecules frequently are bound to positively charged regions on the surface of protein called exosites. In this review, we carried out extensive analyses of the structural determinants of thrombin partnerships by surveying literature data as well as the structural content of the Protein Data Bank (PDB). In particular, we used the information collected on functional, natural, and synthetic molecular ligands to define the anatomy of the exosites and to quantify the interface area between thrombin and exosite ligands. In this framework, we reviewed in detail the specificity of thrombin binding to aptamers, a class of compounds with intriguing pharmaceutical properties. Although these compounds anchor to protein using conservative patterns on its surface, the present analysis highlights some interesting peculiarities. Moreover, the impact of thrombin binding aptamers in the elucidation of the cross-talk between the two distant exosites is illustrated. Collectively, the data and the work here reviewed may provide insights into the design of novel thrombin inhibitors.

Computational approach, scanning electron and fluorescence microscopies revealed insights into the action mechanisms of anticandidal peptide Mo-CBP3-PepIII

AIMS

The Candida genus is composed of opportunistic pathogens that threaten public health. Given the increase in resistance to current drugs, it is necessary to develop new drugs to treat infections by these pathogens. Antimicrobial peptides are promising alternative molecules with low cost, broad action spectrum and low resistance induction. This study aimed to clarify the action mechanisms of synthetic peptides against Candida albicans.

MAIN METHODS

The mode of action of the anticandidal peptides Mo-CBP₃-PepIII were analyzed through molecular dynamics and quantum biochemistry methods against Exo-β-1,3-glucanase (EXG), vital to cell wall metabolism. Furthermore, scanning electron (SEM) and fluorescence (FM) microscopies were employed to corroborate the in silico data.

KEY FINDINGS

Mo-CBP₃-PepIII strongly interacted with EXG (-122.2 kcal mol^-1) at the active site, higher than the commercial inhibitor pepstatin. Also, molecular dynamics revealed the insertion of Mo-CBP₃-PepIII into the yeast membrane. SEM analyses revealed that Mo-CBP₃-PepIII induced cracks and scars of the cell wall and FM analyses confirmed the pore formation on the Candida membrane.

SIGNIFICANCE

Mo-CBP₃-PepIII has strong potential as a new drug with a broad spectrum of action, given its different mode of action compared to conventional drugs.

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.

Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis

It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.

mTOR–mLST8 interaction: hot spot identification through quantum biochemistry calculations

Quantum calculation of mTOR–mLST8 interaction.

Effect of double mutations T790M/L858R on conformation and drug-resistant mechanism of epidermal growth factor receptor explored by molecular dynamics simulations

Epidermal growth factor receptor (EGFR) is one of the most promising targets for the treatment of cancers. Double mutations T790M/L858R lead to different degrees of drug resistance toward inhibitors. In this study, molecular dynamics (MD) simulations followed by principal component analysis are performed to study the conformational changes of EGFR induced by T790M/L858R. The results suggest that T790M/L858R cause obvious disturbance of the structural stability of EGFR. Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) and residue-based free energy decomposition methods are integrated to explore the drug-resistant mechanism of T790M/L858R toward inhibitors. The results show that the decrease in van der Waals interactions of inhibitors with the mutated EFGR relative to the wild-type (WT) one is the main force inducing drug resistance of T790M/L858R toward inhibitors TAK-285, while drug resistance toward W2P and HKI-272 is dominated by the decrease in van der Waals interactions and the increase in polar interactions. We expect that the information obtained from this study can aid rational design of effective drugs to relieve drug resistance of EGFR induced by T790M/L858R.

Quantum chemical studies on anion specificity of CαNN motif in functional proteins

Anion binding C^αNN motif is found in functionally important regions of protein structures. This motif based only on backbone atoms from three adjacent residues, recognizes free sulphate or phosphate ion as well as phosphate groups in nucleotides and in a variety of cofactors. The mode of anion recognition and microscopic picture of binding interaction remains unclear. Here we perform self-consistent quantum chemical calculations considering sulphate and phosphate bound C^αNN motif fragments from crystal structures of functional proteins in order to figure out microscopic basis of anion recognition. Our calculations indicate that stability and preference of the anion in the motif depends on the sequence of the motif. The stabilization energy is larger in case of polar residue containing motif fragment. Nitrogen atom of the polar residue of motif mainly participates in the coordination at the lowest energy levels. Anion replacement decreases stabilization energy along with coordination between motif atoms and oxygen atoms of anion shifted to higher energies, suggesting preference of the motif residues to specific anion. Our analysis may be helpful to understand microscopic basis of interaction between proteins and ionic species.

Molecular Dynamics Exploration of Selectivity of Dual Inhibitors 5M7, 65X, and 65Z toward Fatty Acid Binding Proteins 4 and 5

Designing highly selective inhibitors of fatty acid binding proteins 4 and 5 (FABP4 and FABP5) is of importance for treatment of some diseases related with inflammation, metabolism, and tumor growth. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were performed to probe binding selectivity of three inhibitors (5M7, 65X, and 65Z) to FABP4/FABP5 with K_i values of 0.022/0.50 μM, 0.011/0.086 μM, and 0.016/0.12 μM, respectively. The results not only suggest that all inhibitors associate more tightly with FABP4 than FABP5, but also prove that the main forces driving the selective bindings of inhibitors to FABP4 and FABP5 stem from the difference in the van der Waals interactions and polar interactions of inhibitors with two proteins. Meanwhile, a residue-based free energy decomposition method was applied to reveal molecular basis that inhibitors selectively interact with individual residues of two different proteins. The calculated results show that the binding difference of inhibitors to the residues (Phe16, Phe19), (Ala33, Gly36), (Phe57, Leu60), (Ala75, Ala78), (Arg126, Arg129), and (Tyr128, Tyr131) in (FABP4, FABP5) drive the selectivity of inhibitors toward FABP4 and FABP5. This study will provide great help for further design of effective drugs to protect against a series of metabolic diseases, arteriosclerosis, and inflammation.

Zinc ion-induced conformational changes in new Delphi metallo-β-lactamase 1 probed by molecular dynamics simulations and umbrella sampling

The hydrolysis of a β-lactam core ring caused by new Delphi metallo-β-lactamase 1 (NDM-1) with the help of two zinc cofactors induces significant resistance toward β-lactam antibiotics. Molecular dynamics (MD) simulations and the umbrella sampling method are integrated to study the conformational change mechanism of NDM-1 mediated by zinc ion binding. The statistical analyses of interaction contacts of the antibiotic ampicillin (AMP) with residues based on MD trajectories suggest that two Zn ions are essential for maintaining the binding of AMP with NDM-1. Umbrella sampling simulations further reveal that double-Zn coordination exerts strong restriction on the motions of loop L10 relative to loops L3 and L4. Principal component (PC) analysis also demonstrates that zinc ion binding totally inhibits the motion extent of NDM-1 and changes internal motion modes in NDM-1. We expect that the current study can provide significant dynamical information involving conformational changes of NDM-1 for the development of efficient inhibitors to decrease drug resistance of NDM-1 toward antibiotics.

Insight into binding mechanisms of inhibitors MKP56, MKP73, MKP86, and MKP97 to HIV-1 protease by using molecular dynamics simulation

HIV-1 protease (PR) has been a significant target for design of potent inhibitors curing acquired immunodeficiency syndrome. Molecular dynamics simulations coupled with molecular mechanics Poisson-Boltzmann surface area method were performed to study interaction modes of four inhibitors MKP56, MKP73, MKP86, and MKP97 with PR. The results suggest that the main force controlling interactions of inhibitors with PR should be contributed by van der Waals interactions between inhibitors and PR. The cross-correlation analyses based on MD trajectories show that inhibitor binding produces significant effect on the flap dynamics of PR. Hydrogen bond analyses indicate that inhibitors can form stable hydrogen bonding interactions with the residues from the catalytic strands of PR. The contributions of separate residues to inhibitor bindings are evaluated by using residue-based free energy decomposition method and the results demonstrate that the CH-π and CH-CH interactions between the hydrophobic groups of inhibitors with residues drive the associations of inhibitors with PR. We expect that this study can provide a significant theoretical aid for design of potent inhibitors targeting PR.

Energetic description of cilengitide bound to integrin

We employ quantum chemistry methods to investigate the binding energy features of the cyclic RGD pentapeptide cilengitide interacting with the integrin receptor α_Vβ₃.

Quantum chemical studies on nucleophilic sites in calcium ion bound zwitterionic calmodulin loops

Quantum chemical calculation on Ca²⁺ bound zwitterionic calmodulin-loops shows terminal capping contributions localized near HOMO and LUMO, which decay exponentially and presence of nucleophilic site at the phenyl-ring O_h of Y99.

Explaining RANKL inhibition by OPG through quantum biochemistry computations and insights into peptide-design for the treatment of osteoporosis

Quantum biochemistry computations are applied to precisely describe important protein–protein interactions, providing a basis for the design of inhibitory peptides against osteoporosis.

Quantum Chemical Studies on Stability and Chemical Activities in Calcium Ion Bound Calmodulin Loops

Quantum chemical (QC) calculations for macromolecules require truncation of the molecule, highlighting the portion of interest due to heavy computation cost. As a result, an estimation of the effects of truncation is important to interpret the energy spectrum of such calculations. We perform density functional theory based QC calculations on calcium ion bound EF-hand loops of Calmodulin isolated from the crystal structure in an implicit solvent. We find that the terminal contributions of neutral capping are negligible across the entire ground-state energy spectrum. The coordination energy range and the nature of hybridization of the coordination state molecular orbitals remain qualitatively similar across these loops. While the HOMO and LUMO of loops in the N-terminal domain are dominated by the acidic aspartates, and the polar/hydrophobic residues, respectively, these levels of the C-terminal domain loops show strong localized electron density on the phenyl rings of the tyrosines. The Fukui index calculation identifies the hydroxyl oxygen in the phenyl ring of Y99 as a potent nucleophile. Our analysis indicates a general way of interpreting the electronic energy spectra to understand stability and functions of large biomolecules where the truncation of the molecule and, hence, the terminal capping effects are inevitable.

A quantum biochemistry investigation of willardiine partial agonism in AMPA receptors

We employ quantum biochemistry methods based on the Density Functional Theory (DFT) approach to unveil detailed binding energy features of willardiines co-crystallized with the AMPA receptor.

Interaction of I50V mutant and I50L/A71V double mutant HIV-protease with inhibitor TMC114 (darunavir): molecular dynamics simulation and binding free energy studies

In the present work, the binding of inhibitor TMC114 (darunavir) to wild-type (WT), single (I50V) as well as double (I50L/A71V) mutant HIV-proteases (HIV-pr) was investigated with all-atom molecular dynamics (MD) simulations as well as molecular mechanic-Poisson-Boltzmann surface area (MM-PBSA) calculation. For both the apo and complexed HIV-pr, many intriguing effects due to double mutant, I50L/A71V, are observed. For example, the flap-flap distance and the distance from the active site to the flap residues in the apo I50L/A71V-HIV-pr are smaller than those of WT- and I50V-HIV-pr, probably making the active site smaller in volume and closer movement of flaps. For the complexed HIV-pr with TMC114, the double mutant I50L/A71V shows a less curling of the flap tips and less flexibility than WT and the single mutant I50V. As for the other previous studies, the present results also show that the single mutant I50V decreases the binding affinity of I50V-HIV-pr to TMC, resulting in a drug resistance; whereas the double mutant I50L/A71V increases the binding affinity, and as a result of the stronger binding, the I50L/A71V may be well adapted by the TMC114. The energy decomposition analysis suggests that the increase of the binding for the double mutant I50L/A71V-HIV-pr can be mainly attributed to the increase in electrostatic energy by -5.52 kacl/mol and van der Waals by -0.42 kcal/mol, which are canceled out in part by the increase of polar solvation energy of 1.99 kcal/mol. The I50L/A71V mutant directly increases the binding affinity by approximately -0.88 (Ile50 to Leu50) and -0.90 (Ile50' to Leu50') kcal/mol, accounting 45% for the total gain of the binding affinity. Besides the direct effects from the residues Leu50 and Leu50', the residue Gly49' increases the binding affinity of I50L/A71V-HIV-pr to the inhibitor by -0.74 kcal/mol, to which the electrostatic interaction of Leu50's backbone contributes by -1.23 kcal/mol. Another two residues Ile84 and Ile47' also increase the binding affinity by -0.22 and -0.29 kcal/mol, respectively, which can be mainly attributed to van der Waals terms (ΔT(vdw) = -0.21 and -0.39 kcal/mol).

Computational simulations of the interactions between acetyl-coenzyme-A carboxylase and clodinafop: resistance mechanism due to active and nonactive site mutations

Grass weed populations resistant to acetyl-CoA carboxylase-inhibiting (ACCase; EC 6.4.1.2) herbicides represent a major problem for the sustainable development of modern agriculture. In the present study, extensive computational simulations, including homology modeling, molecular dynamics (MD) simulations, and molecular mechanics-Poisson-Boltzmann surface area (MM/PBSA) calculations, have been carried out to uncover the detailed molecular mechanism of Alopecurus myosuroides resistance to clodinafop, a commercial herbicide targeting ACCase. All the computational model and energetic results indicated that W374C, I388N, D425G, and G443A mutations have great effects on the conformational change of the binding pocket and the hydrogen-bonding interactions. The pi-pi interaction between ligand and the residue of Phe377 and Tyr161', playing an important contribution to the binding affinity, were decreased after mutations. In addition, the hydrogen-bonding interactions between clodinafop and the residues (Ile158' and Ala54') disappeared or decreased significantly upon mutation. As a result, the mutant-type ACCase has a lower affinity for the inhibitor binding than the wild-type enzyme, which accounts for the molecular basis of herbicidal resistance. The structural and mechanistic insights obtained from the present study will provide a valuable clue for future designing of a promising inhibitor to reduce drug resistance associated with both active and nonactive site mutations.

Insights into drug resistance of mutations D30N and I50V to HIV-1 protease inhibitor TMC-114: free energy calculation and molecular dynamic simulation

The single mutations D30N and I50V are considered as the key residue mutations of the HIV-1 protease drug resistance to inhibitors in clinical use. In this work, molecular dynamics (MD) simulations combined with the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method have been performed to investigate the drug-resistant mechanisms of D30N and I50V to an inhibitor TMC-114. The analyses of absolute binding free energies using the separate trajectory approach suggests that the decrease in the van der Waals energy and electrostatic energy in the gas phase results in the drug resistance of D30N to TMC-114, while for I50V, the decrease in the electrostatic energy mainly drive its drug resistance to TMC-114. Detailed binding free energies between TMC-114 and individual protein residues are computed by using a per-residue basis decomposition method, which provides insights into the inhibitor-protein binding mechanism and also explains the drug-resistant mechanisms of mutations D30N and I50V to TMC-114. The study shows that the loss of the hydrogen bond between TMC-114 and the side chain of Asn30' is the main driving force of the resistance of D30N to TMC-114, and in the case of I50V, the increase in the polar solvation energies between TMC-114 and two residues Val50' and Asp30' definitively drives the resistance of I50V to TMC-114. We expect that this work can provide some helpful insights into the nature of mutational effect and aid the future design of better inhibitors.